18/08/17 01:41:51.87 Smlos65L.net
>>1の論文はともかくf(pr)が2m+1の倍数になるのか検証。
Prelude> mapM_ print $ [[(flip mod (2*m+1)) $div ((2*pr-1)^(4*m+2) - 1) (pr-1) | m<-[1,2..10]]|pr<-[3,7,19,31,37,79,97,139,157,199,211,229,271,307]]
[0,2,5,0,1,12,12,12,12,0]
[0,3,0,0,6,2,3,11,9,0]
[0,1,6,0,10,11,6,8,0,0]
[0,0,5,0,3,7,0,5,10,0]
[0,3,1,0,5,5,3,12,15,0]
[0,1,1,0,8,0,6,10,12,0]
[0,3,3,0,3,11,3,14,8,0]
[0,1,3,0,6,10,6,12,5,0]
[0,3,5,0,1,0,3,16,1,0]
[0,1,5,0,0,3,6,14,17,0]
[0,0,0,0,8,12,0,11,8,0]
[0,1,6,0,3,6,6,15,0,0]
[0,0,6,0,6,5,0,13,1,0]
[0,3,3,0,7,6,3,0,12,0]
ダメやね。
pr, 2pr-1がともに素数になる数でやってみたけどこれだけの条件では2m+1は((2pr-1)^(4m+2)-1)/(pr-1)の約数とは言えない。
ので2m+1の因子に2yの因子でないものが現れる可能性を否定できないね。