OFFSET
1,1
LINKS
Ivan Neretin, Table of n, a(n) for n = 1..10000
Max Alekseyev, PARI/GP Scripts for Miscellaneous Math Problems (invphi.gp).
MATHEMATICA
mx = 100; Complement[Range[mx], Table[Nest[DivisorSigma[1, #] &, n, 2], {n, mx}]] (* Ivan Neretin, Feb 12 2018 *)
PROG
(PARI) is(k) = {my(v = invsigma(k)); for(i = 1, #v, if(invsigmaNum(v[i]) > 0, return(0))); 1; } \\ Amiram Eldar, Dec 26 2024, using Max Alekseyev's invphi.gp
CROSSREFS
KEYWORD
nonn
AUTHOR
Jaroslav Krizek, Jul 17 2009
STATUS
approved