rMod[m,n_x]:=rMod[m,n]=Mod[分子[m]*PowerMod[分母[m],-1,n],n,-n/2]
p[n_u]:=p[n]=素数[n];
Do[r=0;Do[Do[If[Mod[p[g]^(Part[Divisors[p[n]-1],i])-1,p[n]]==0,Goto[aa]],{i,1,长度[除数[p[n]-1]]-1}];
Do[If[Mod[p[g]^(Part[Divisors[p[p[g]]-1],j])-1,p[p[g]]]==0,Goto[aa]],{j,1,长度[除数[p[p[g]]-1]]-1};
r=r+1;标签[aa],{g,1,n-1}];打印[n,“”,r],{n,1,80}]
|