OFFSET
1,2
COMMENTS
It can be shown that this is a subsequence of A028982.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
Select[Range[91000], CoprimeQ[DivisorSigma[0, #], DivisorSigma[2, #]]&] (* Harvey P. Dale, May 11 2011 *)
PROG
(PARI) isok(n) = gcd(sigma(n, 2), numdiv(n)) == 1; \\ Michel Marcus, Sep 24 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
a(1)=1 added by Amiram Eldar, Sep 24 2019
STATUS
approved