OFFSET
1,1
COMMENTS
The numbers of terms that do not exceed 10^k, for k = 3, 4, ..., are 5, 47, 459, 4655, 46733, 460693, 4612685, 46177602, ... . Apparently, the asymptotic density of this sequence exists and equals 0.00461... .
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
Select[Range[10000], DivisorSigma[-1, #^2] > 2 && DivisorSigma[-1, (#+1)^2] > 2 &]
PROG
(PARI) is1(k) = {my(f = factor(k)); prod(i = 1, #f~, f[i, 2] *= 2); sigma(f, -1) > 2; }
list(lim) = {my(q1 = is1(1), q2); for(k = 2, lim, q2 = is1(k); if(q1 && q2, print1(k-1, ", ")); q1 = q2); }
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Apr 25 2025
STATUS
approved
