高校数学の質問スレ Part434

235 ：１３２人目の素数さん：2024/04/16(火) 15:02:12.74 ID: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]を表示する。
これでもかなり遅いが、正常には動く。