ｌｃさんより，追加メールを頂いた．

『Wagstaff の原論文のようなものが見つかりました。

「https://www.ams.org/journals/mcom/1983-40-161/S0025-5718-1983-0679454-X/home.html」

この論文の式は、

f(x) = exp(Euler)/log_e(2) * log_e( log_e(x) )

= exp(Euler) * log_2( log_e(x) )

に間違いないです。

47番に当てはめると、あまり合ってないですが。

f(2^43112609) = 44.229...

参考までに、下記に48番候補までの計算結果があり、計算自体は、おおむね正確のようです。

「http://www.cscanada.net/index.php/sms/article/download

/3976/6166」』

Lenstra-Pomerance-Wagstaff の式は大ざっぱな近似なので、log_eとlog_2の違いは、大した問題ではないのかもしれませんが、その違いに本質的な意味があるのなら、コラムに書いていただければとも思います。

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝