OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
Max Alekseyev, PARI/GP Scripts for Miscellaneous Math Problems (invphi.gp).
PROG
(PARI) isok(k) = if(k == 1, 1, my(v = invsigma(sigma(k))); #v > 1 && gcd(v) == 1); \\ Amiram Eldar, May 28 2025, using Max Alekseyev's invphi.gp (see links).
CROSSREFS
KEYWORD
nonn
AUTHOR
Max Alekseyev, Jan 19 2025
STATUS
approved
