OFFSET
1,1
COMMENTS
Apart from 1, all terms of A211656 belong here since the solutions to sigma(x)=sigma(n) form a singleton and thus their GCD is n itself.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
Max Alekseyev, PARI/GP Scripts for Miscellaneous Math Problems (invphi.gp).
FORMULA
Numbers k such that A241479(k) is not equal to 1.
PROG
(PARI) sigv(n) = select(i->sigma(i) == n, vector(n, i, i));
isok(n) = (gcd(sigv(sigma(n))) != 1);
(PARI) is(k) = gcd(invsigma(sigma(k))) > 1; \\ Amiram Eldar, Dec 19 2024, using Max Alekseyev's invphi.gp
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Marcus, Apr 23 2014
STATUS
approved