OFFSET
1,2
COMMENTS
The asymptotic density of this sequence is d = zeta(6) * Product_{p prime} (1 - 1/p^2 + 1/p^5 - 2/p^6 + 1/p^7) = 0.62636565734301010494... .
The asymptotic density of the squarefree numbers within this sequence is 1/(zeta(2) * d) = 0.97056263338702466237... .
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
FORMULA
Sum_{n>=1} 1/a(n)^s = zeta(6*s) * Product_{p prime} (1 + 1/p^s + 1/p^(5*s) - 1/p^(6*s)).
MATHEMATICA
q[n_] := AllTrue[FactorInteger[n][[;; , 2]], CoprimeQ[6, #] &]; Select[Range[105], q]
PROG
(PARI) isok(k) = vecsum(apply(x -> if(gcd(x, 6) == 1, 0, 1), factor(k)[, 2])) == 0;
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Amiram Eldar, Nov 05 2025
STATUS
approved
