OFFSET
0,4
COMMENTS
Also number of ordered pairs of positive integers (i, j) such that i+j <= n and i*j <= n. - Seiichi Manyama, Sep 04 2021
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..10000
Eric Weisstein's World of Mathematics, Quadratic Equation
FORMULA
a(n) = A091626(n) - n - 1. a(n) = a(n-1) + ceiling(tau(n)/2), n>1. Partial sums of A038548. - Vladeta Jovovic, Jun 12 2004
G.f.: (1/(1 - x)) * (-x + Sum_{k>=1} x^(k^2)/(1 - x^k)). - Seiichi Manyama, Sep 04 2021
MATHEMATICA
Accumulate[ Join[{0, 0}, Table[ Ceiling[ DivisorSigma[0, n]/2], {n, 2, 64}]]] (* Jean-François Alcover, Oct 23 2012, after Vladeta Jovovic *)
PROG
(PARI) a(n) = sum(i=1, n, sum(j=i, n-i, i*j<=n)); \\ Seiichi Manyama, Sep 04 2021
(PARI) N=66; x='x+O('x^N); concat([0, 0], Vec((-x+sum(k=1, sqrtint(N), x^k^2/(1-x^k)))/(1-x))) \\ Seiichi Manyama, Sep 04 2021
(Python)
from math import isqrt
def A091627(n):
m = isqrt(n)
return 0 if n == 0 else sum(n//k for k in range(1, m+1))-m*(m-1)//2-1 # Chai Wah Wu, Oct 07 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Eric W. Weisstein, Jan 24 2004
STATUS
approved