OFFSET
0,19
COMMENTS
Let b(n,k,m) be the number of partitions of n into k distinct parts containing the part m. Define G_{k,m}(q) = Sum_{n>=0} b(n,k,m) * q^n.
Then G_{k,m}(q) = q^(k*(k+1)/2) * Sum_{i=1..min(k,m)} q^((k-i+1)*(m-i)) * q_binomial(m-1,i-1) / Product_{j=1..k-i} (1-q^j).
Also, G_{k,m}(q) = Sum_{i=1..k} (-1)^(i-1) * q^(m*i) * Product_{j=1..k-i} q^j/(1-q^j).
LINKS
Index entries for linear recurrences with constant coefficients, signature (1,1,0,0,-2,0,0,1,1,-1).
FORMULA
G.f.: G_{5,5}(q) = q^15 * (1 + q^2*(1+q+q^2+q^3)/(1-q) + q^6*(1+q+2*q^2+q^3+q^4)/((1-q)*(1-q^2)) + q^12*(1+q+q^2+q^3)/((1-q)*(1-q^2)*(1-q^3)) + q^20/((1-q)*(1-q^2)*(1-q^3)*(1-q^4))).
G.f.: Sum_{j=1..5} (-1)^(j-1) * q^(5*j) * Product_{k=1..5-j} q^k/(1-q^k).
a(n) = a(n-1) + a(n-2) - 2*a(n-5) + a(n-8) + a(n-9) - a(n-10) for n > 35.
a(n) = [n>14]*floor((n^3 - 42*n^2 + 672*n - 4176)/144 + (n mod 2)*(n-12)/16 + [n<16] + [n<21] + [n=25]), where "[.]" denote the Iverson brackets. - Hoang Xuan Thanh, May 16 2026
PROG
(PARI) my(N=70, q='q+O('q^N)); concat([0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], Vec(sum(j=1, 5, (-1)^(j-1)*q^(5*j)*prod(k=1, 5-j, q^k/(1-q^k)))))
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, May 15 2026
STATUS
approved
