OFFSET
1,1
FORMULA
a(n) = 10^n - A071541(n). - Amiram Eldar, Jul 27 2024
MATHEMATICA
f[n_] := Plus @@ Select[ Divisors[n], Sqrt[n/2] <= # < Sqrt[n*2] &]; s = 0; k = 0; Do[ While[k < 10^n, k++; If[ f[k] == 0, s++ ]]; Print[s], {n, 1, 7}]
PROG
(PARI) lista(nmax) = {my(c = 0, pow = 10); for(k = 1, 10^nmax, if(sumdiv(k, d, if(d^2 >= k/2 && d^2 < 2*k, d, 0)) == 0, c++); if(k == pow, print1(c, ", "); pow *= 10)); } \\ Amiram Eldar, Jul 27 2024
CROSSREFS
KEYWORD
nonn,hard,more
AUTHOR
Robert G. Wilson v, May 31 2002
EXTENSIONS
a(8) from Sean A. Irvine, Jul 24 2024
a(9)-a(10) from Amiram Eldar, Jul 27 2024
STATUS
approved