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[5000], CoprimeQ[DivisorSigma[0, #], DivisorSigma[3, #]] &] (* Amiram Eldar, Aug 08 2019 *)
PROG
(PARI) isok(n) = gcd(numdiv(n), sigma(n, 3)) == 1; \\ Michel Marcus, Sep 24 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved