OFFSET
2,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 2..70
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.
FORMULA
a(n) >= A007374(n).
MATHEMATICA
seq[len_] := Module[{v = Table[0, {len+1}], c = 0, k = 2, s, ns}, While[c < len, s = PhiInverse[k]; ns = Length[s]; If[0 < ns <= len + 1 && AllTrue[s, PhiMultiplicity[#] == 0 &] && v[[ns]] == 0, v[[ns]] = k; c++]; k += 2]; Rest[v]]; seq[10] (* using David M. Bressoud's CNT.m *)
PROG
(PARI) lista(len) = {my(v = vector(len+1), c = 0, k = 2, s, ns, ans); while(c < len, s = invphi(k); ns = #s; ans = 1; for(i = 1, ns, if(istotient(s[i]), ans = 0; break)); if(ans && ns > 0 && ns <= len + 1 && v[ns] == 0, c++; v[ns] = k); k += 2); vecextract(v, "^1"); } \\ using Max Alekseyev's invphi.gp
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Nov 29 2024
STATUS
approved
