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..5000
MATHEMATICA
Select[Range[300000], CoprimeQ[DivisorSigma[0, #], DivisorSigma[4, #]]&] (* Harvey P. Dale, May 30 2012 *)
PROG
(PARI) isok(n) = gcd(sigma(n, 4), numdiv(n)) == 1; \\ Michel Marcus, Sep 24 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved