24/04/16 15:02:12.74 7gGe0Okf.net
>>220
UnionもLengthもMinもいらない。
O(n^2)のメモリもいらない。
count=0;
Do[
For[flag=1;k=1,flag==1 && k<n,k++,If[Mod[k^(n-1),n]==1,Null,flag=0]];
If[flag==1,count++;Print[{count,n,Prime[count]}],Null];
,{n,2,18000}]
k=1,2,3,...,n-1に対し、Mod[k^(n-1),n]==1 となるような n を見つけたら、
countを1アップして、countとnとPrime[count]を表示する。
これでもかなり遅いが、正常には動く。