fix bug in 153.pi: i*j<=N