A340952 - OEIS

%I #8 Jan 31 2021 20:22:03

%S 1,0,7,0,21,7,35,42,35,112,42,182,112,210,260,217,462,252,651,399,728,

%T 777,672,1232,749,1533,1127,1659,1617,1792,2289,1890,2926,2212,3339,

%U 2990,3584,3654,4046,4613,4263,5754,4487,6636,5733,6825,7014,7203,8617,7560,10087,8302

%N Number of ways to write n as an ordered sum of 7 nonzero triangular numbers.

%H Alois P. Heinz, <a href="/A340952/b340952.txt">Table of n, a(n) for n = 7..10000</a>

%F G.f.: (theta_2(sqrt(x)) / (2 * x^(1/8)) - 1)^7, where theta_2() is the Jacobi theta function.

%p b:= proc(n, k) option remember; local r, t, d; r, t, d:= $0..2;

%p       if n=0 then `if`(k=0, 1, 0) else

%p       while t<=n do r:= r+b(n-t, k-1); t, d:= t+d, d+1 od; r fi

%p     end:

%p a:= n-> b(n, 7):

%p seq(a(n), n=7..58);  # _Alois P. Heinz_, Jan 31 2021

%t nmax = 58; CoefficientList[Series[(EllipticTheta[2, 0, Sqrt[x]]/(2 x^(1/8)) - 1)^7, {x, 0, nmax}], x] // Drop[#, 7] &

%Y Cf. A000217, A010054, A053603, A053604, A226252, A319817, A340949, A340950, A340951, A340953, A340954, A340955.

%K nonn

%O 7,3

%A _Ilya Gutkovskiy_, Jan 31 2021