%I #34 Sep 13 2020 14:43:16
%S 1,5,5,15,5,25,5,35,15,25,5,75,5,25,25,70,5,75,5,75,25,25,5,175,15,25,
%T 35,75,5,125,5,126,25,25,25,225,5,25,25,175,5,125,5,75,75,25,5,350,15,
%U 75,25,75,5,175,25,175,25,25,5,375,5,25,75,210,25,125,5,75,25,125,5
%N tau_5(n) = number of ordered 5-factorizations of n.
%H Seiichi Manyama, <a href="/A061200/b061200.txt">Table of n, a(n) for n = 1..10000</a> (terms 1..1000 from Enrique PĂ©rez Herrero)
%F tau_k(n) = |{(x_1,x_2,...,x_k): x_1*x_2*...*x_k=n}|, number of ordered k-factorizations of n.
%F tau_k(p^m) = (-1)^(k-1)*binomial(-m-1,k-1), p prime.
%F limit(tau_k(n)/n^epsilon, n=infinity) = 0, for any epsilon>0.
%F tau_k(n) = Sum_{d|n} tau_(k-1)(d), tau_1(n)=1.
%F Dirichlet g.f.: (zeta(s))^k.
%F For explicit formula, see A007425.
%F G.f.: Sum_{k>=1} tau_4(k)*x^k/(1 - x^k). - _Ilya Gutkovskiy_, Oct 30 2018
%t tau[n_, 1] = 1; tau[n_, k_] := tau[n, k] = Plus @@ (tau[ #, k - 1] & /@ Divisors[n]); Table[ tau[n, 5], {n, 77}] (* _Robert G. Wilson v_ *)
%t tau[1, k_] := 1; tau[n_, k_] := Times @@ (Binomial[Last[#]+k-1, k-1]& /@ FactorInteger[n]); Table[tau[n, 5], {n, 1, 100}] (* _Amiram Eldar_, Sep 13 2020 *)
%o (PARI) for(n=1,100,print1(sumdiv(n,k,sumdiv(k,x,sumdiv(x,y,numdiv(y)))),","))
%o (PARI) a(n)=sumdivmult(n,k,sumdivmult(k,x,sumdivmult(x,y,numdiv(y)))) \\ _Charles R Greathouse IV_, Sep 09 2014
%o (PARI) a(n, f=factor(n))=f=f[, 2]; prod(i=1, #f, binomial(f[i]+4, 4)) \\ _Charles R Greathouse IV_, Oct 28 2017
%Y Cf. tau_2(n): A000005, tau_3(n): A007425, tau_4(n): A007426, tau_6(n): A034695, (unordered) 2-factorization of n: A038548, (unordered) 3-factorization of n: A034836, A001055, A006218, A061201, A061202, A061203, A061204.
%Y Column k=5 of A077592.
%K easy,nonn,mult
%O 1,2
%A _Vladeta Jovovic_, Apr 21 2001