160326(2)

Ruby


Wieferich prime

「数論〈未解決問題〉の事典」の「A3」によれば、

1.25 × 10^15 までの素数で

2^p - 2 がp^2 で割り切れるものはただ２つ

らしい。

10^5 までなら、一瞬で答えが出る。

require 'prime'

Prime.each(10 ** 5).each{|i|
  m = 2 ** i - 2
  p i if m % (i * i) == 0
}

出力結果
1093
3511