- プログラミングのお題スレ Part6©2ch.net
230 ：デフォルトの名無しさん[sage]：2015/01/04(日) 15:42:07.21 ID:ulO+nk3L - >>220 Io
powerMod := method(a, b, m,
r := 1
for(i, 0, b log2 floor, if(b at(i) == 1, r = r * a % m); a = a * a % m)
r
)
isPrime := method(n,
if(n isEven or n < 2, return(n == 2))
for(i, 3, n sqrt floor, 2, if(n % i == 0, return(false)))
true
)
isPrime2 := method(n,
if(n == 2, (2 ** n - 2) % n == 0, powerMod(2, n, n) == 2)
)
f := method(z,
for(i, 2, z, if(isPrime(i) != isPrime2(i), i println))
)
Io> f(1000)
341
561
645
|