OFFSET
1,1
COMMENTS
LINKS
Max Alekseyev, PARI/GP Scripts for Miscellaneous Math Problems (invphi.gp).
David M. Bressoud, A Course in Computational Number Theory (web page), CNT.m, Computational Number Theory Mathematica package.
EXAMPLE
MATHEMATICA
q[n_] := Count[Join @@ PhiInverse[PhiInverse[2*n]], _?(IntegerQ[PrimitiveRoot[#]] &)] == 0; Select[Range[120], q] (* Amiram Eldar, Jan 28 2025, using David M. Bressoud's CNT.m *)
PROG
(PARI) isA033948(n) = {my(f = factor(n)); lcm(znstar(f)[2]) == eulerphi(f); }
isok(k) = {my(v = invphi(2*k), w, c = 0); for(i = 1, #v, c += vecsum(apply(x -> isA033948(x), invphi(v[i])))); c == 0; } \\ Amiram Eldar, Jan 28 2025, using Max Alekseyev's invphi.gp
CROSSREFS
KEYWORD
nonn
AUTHOR
David James Sycamore, Jan 28 2025
STATUS
approved