OFFSET
1,2
COMMENTS
(tau<=)_k(n) = |{(x_1,x_2,...,x_k): x_1*x_2*...*x_k<=n}|, i.e. (tau<=)_k(n) is number of solutions to x_1*x_2*...*x_k<=n, x_i>0.
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 1..10000
Vaclav Kotesovec, Graph - The asymptotic ratio (100000 terms)
FORMULA
(tau<=)_k(n)=Sum_{i=1..n} tau_k(i). a(n)=partial sums of A034695.
a(n) = Sum_{k=1..n} tau_{5}(k) * floor(n/k), where tau_{5} is A061200. - Enrique Pérez Herrero, Jan 23 2013
a(n) ~ n*(log(n)^5/120 + (g/4 - 1/24)*log(n)^4 + (5*g^2/2 - g - g1 + 1/6)*log(n)^3 + (10*g^3 - 15*g^2/2 + (3 - 15*g1)*g + 3*g1 + 3*g2/2 - 1/2)*log(n)^2 + (15*g^4 - 20*g^3 + (15 - 60*g1)*g^2 + (30*g1 + 15*g2 - 6)*g + 15*g1^2 - 6*g1 - 3*g2 - g3 + 1)*log(n) + 6*g^5 - 15*g^4 + (20 - 60*g1)*g^3 + (60*g1 + 30*g2 - 15)*g^2 + (60*g1^2 - 30*g1 - 15*g2 - 5*g3 + 6)*g - 15*g1^2 + g1*(6 - 15*g2) + 3*g2 + g3 + g4/4 - 1), where g is the Euler-Mascheroni constant A001620 and g1, g2, g3, g4 are the Stieltjes constants, see A082633, A086279, A086280 and A086281. - Vaclav Kotesovec, Sep 10 2018
MATHEMATICA
nmax = 50; tau4 = Table[DivisorSum[n, DivisorSigma[0, n/#]*DivisorSigma[0, #] &], {n, 1, nmax}]; tau5 = Table[Sum[tau4[[d]], {d, Divisors[n]}], {n, nmax}]; Accumulate[Table[Sum[tau5[[d]], {d, Divisors[n]}], {n, nmax}]] (* Vaclav Kotesovec, Sep 10 2018 *)
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Vladeta Jovovic, Apr 21 2001
STATUS
approved