OFFSET
1,2
LINKS
Seiichi Manyama, Table of n, a(n) for n = 1..10000 (terms 1..1000 from Enrique PĂ©rez Herrero)
FORMULA
G.f.: Sum_{k>=1} tau_8(k)*x^k/(1 - x^k). - Ilya Gutkovskiy, Oct 30 2018
Multiplicative with a(p^e) = binomial(e+8,8). - Amiram Eldar, Sep 13 2020
MATHEMATICA
tau[n_, 1] = 1; tau[n_, k_] := tau[n, k] = Plus @@ (tau[ #, k - 1] & /@ Divisors[n]); Table[ tau[n, 9], {n, 65}] (* Robert G. Wilson v, Nov 02 2005 *)
tau[1, k_] := 1; tau[n_, k_] := Times @@ (Binomial[Last[#]+k-1, k-1]& /@ FactorInteger[n]); Table[tau[n, 9], {n, 1, 100}] (* Amiram Eldar, Sep 13 2020 *)
PROG
(PARI) for(n=1, 100, print1(sumdiv(n, i, sumdiv(i, j, sumdiv(j, k, sumdiv(k, l, sumdiv(l, m, sumdiv(m, o, sumdiv(o, x, numdiv(x)))))))), ", "))
(PARI) a(n, f=factor(n))=f=f[, 2]; prod(i=1, #f, binomial(f[i]+8, 8)) \\ Charles R Greathouse IV, Oct 28 2017
CROSSREFS
KEYWORD
mult,nonn
AUTHOR
Gerald McGarvey, Oct 25 2005
STATUS
approved