24/04/16 22:25:34.50 7gGe0Okf.net
>>238
その前提で組んだけど、(n-1)乗固定ではないのなら、ちょっとだけ改変
count=0;
Do[
For[flag=1;k=1,flag==1 && k<n,k++,If[Mod[k^m,n]==1,Null,flag=0]];
If[flag==1,count++;Print[{count,n,m,Prime[count]}],Null];
,{n,2,18000},{m,1,n-1}]
与えられた、mとnに対し、k=1,2,3,...,n-1と変化しても、常に、Mod[k^m,n]==1なら、出力