|
MATHEMATICA
|
G[n_, k_] := G[n, k] = Mod[Sum[PowerMod[i, k, n], {i, 1, n}], n];
Select[Range[2000], G[#, EulerPhi[#]] == 2 &]
fa=FactorInteger;
se[n_, k_] := Select[Transpose[fa[n]][[1]], IntegerQ[k/(# - 1)] &];
sumlis[li_] := Sum[li[[i]], {i, 1, Length[li]}]
Table[If[Mod[-n/se[n, EulerPhi[n]] // sumlis, n] == 2, n], {n, 1, 1000000}] // Union
|