OFFSET
1,5
COMMENTS
It seems that a(10^n) = (6, 168, 2754, 38561, 495569, ...) ~ 1.1*(n-0.5)*10^n; otherwise said, a(n) ~ 1.1*(log_10(n)-0.5)*n, asymptotically.
The exact value of the coefficient above is 3*log(10)/(2*Pi) = 1.09940339... . - Amiram Eldar, Dec 01 2023
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
FORMULA
a(n) = A147807(n) - n.
a(n) ~ c * n * log(n), where c = 3/(2*Pi) = 0.477464... (A093582). - Amiram Eldar, Dec 01 2023
MATHEMATICA
Accumulate[DivisorSigma[0, Range[72]^2 + 1]/2 - 1] (* Amiram Eldar, Oct 25 2019 *)
PROG
(PARI) s=0; a147806=vector(99, n, s+=numdiv(n^2+1))/2
A147806(n)=sum(p=1, n, numdiv(n^2+1))/2-n
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
M. F. Hasler, Dec 13 2008
STATUS
approved