20/09/07 00:26:33.91 CMvmi5X2.net
>>498
Ｃのソースありがとうございます。
アルゴリズムが理解できないままＲに移植してみました。
> lo=1e6
> up=1e7
> for(a in lo:up){
+ b=a; c=3; d=1
+ while(b>0){
+ if((b%%2)==1) d=(d*c)%%a
+ c=(c*c)%%a
+ b=b%/%2
+ }
+ if((d+1)==a) cat(a,'\n')
+ a=a+1
+ }
1476250
2125250
2550250
3156250
3751250
3906250
5964050
7381250
メモリ不足にならずに計算できたのに驚き