cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A091602 Triangle: T(n,k) is the number of partitions of n such that some part is repeated k times and no part is repeated more than k times.

This page as a plain text file.
%I A091602 #29 Feb 13 2022 06:12:56
%S A091602 1,1,1,2,0,1,2,2,0,1,3,2,1,0,1,4,3,2,1,0,1,5,4,3,1,1,0,1,6,7,3,3,1,1,
%T A091602 0,1,8,8,6,3,2,1,1,0,1,10,12,7,5,3,2,1,1,0,1,12,15,11,6,5,2,2,1,1,0,1,
%U A091602 15,21,14,10,5,5,2,2,1,1,0,1,18,26,20,12,9,5,4,2,2,1,1,0,1,22,35,25,18,11,8,5,4,2,2,1,1,0,1
%N A091602 Triangle: T(n,k) is the number of partitions of n such that some part is repeated k times and no part is repeated more than k times.
%C A091602 From _Gary W. Adamson_, Mar 13 2010: (Start)
%C A091602 The triangle by rows = finite differences starting from the top, of an array in which row 1 = p(x)/p(x^2), row 2 = p(x)/p(x^3), ... row k = p(x)/p(x^k); such that p(x) = polcoeff A000041: (1 + x + 2x^2 + 3x^3 + 5x^4 + 7x^5 + ...)
%C A091602 Note that p(x)/p(x^2) = polcoeff A000009: (1 + x + x^2 + 2x^3 + 2x^4 + ...).
%C A091602 Refer to the example. (End)
%H A091602 Alois P. Heinz, <a href="/A091602/b091602.txt">Rows n = 1..141, flattened</a>
%F A091602 G.f.: G = G(t,x) = sum(k>=1, t^k*(prod(j>=1, (1-x^((k+1)*j))/(1-x^j) ) -prod(j>=1, (1-x^(k*j))/(1-x^j) ) ) ). - _Emeric Deutsch_, Mar 30 2006
%F A091602 Sum_{k=1..n} k * T(n,k) = A264397(n). - _Alois P. Heinz_, Nov 20 2015
%e A091602 Triangle starts:
%e A091602    1:  1;
%e A091602    2:  1,  1;
%e A091602    3:  2,  0,  1;
%e A091602    4:  2,  2,  0,  1;
%e A091602    5:  3,  2,  1,  0,  1;
%e A091602    6:  4,  3,  2,  1,  0,  1;
%e A091602    7:  5,  4,  3,  1,  1,  0,  1;
%e A091602    8:  6,  7,  3,  3,  1,  1,  0,  1;
%e A091602    9:  8,  8,  6,  3,  2,  1,  1,  0,  1;
%e A091602   10: 10, 12,  7,  5,  3,  2,  1,  1,  0,  1;
%e A091602   11: 12, 15, 11,  6,  5,  2,  2,  1,  1,  0,  1;
%e A091602   12: 15, 21, 14, 10,  5,  5,  2,  2,  1,  1,  0,  1;
%e A091602   13: 18, 26, 20, 12,  9,  5,  4,  2,  2,  1,  1,  0,  1;
%e A091602   14: 22, 35, 25, 18, 11,  8,  5,  4,  2,  2,  1,  1,  0,  1;
%e A091602   ...
%e A091602 In the partition 5+2+2+2+1+1, 2 is repeated 3 times, no part is repeated more than 3 times.
%e A091602 From _Gary W. Adamson_, Mar 13 2010: (Start)
%e A091602 First few rows of the array =
%e A091602   ...
%e A091602   1, 1, 1, 2, 2, 3,  4,  5,  6,  8, 10, ... = p(x)/p(x^2) = A000009
%e A091602   1, 1, 2, 2, 4, 5,  7,  9, 13, 16, 22, ... = p(x)/p(x^3)
%e A091602   1, 1, 2, 3, 4, 6,  9, 12, 16, 22, 29, ... = p(x)/p(x^4)
%e A091602   1, 1, 2, 3, 5, 6, 10, 13, 19, 25, 34, ... = p(x)/p(x^5)
%e A091602   1, 1, 2, 3, 5, 7, 10, 14, 20, 27, 37, ... = p(x)/p(x^6)
%e A091602   ...
%e A091602 Finally, taking finite differences from the top and deleting the first "1", we obtain triangle A091602 with row sums = A000041 starting with offset 1:
%e A091602   1;
%e A091602   1, 1;
%e A091602   2, 0, 1;
%e A091602   2, 2, 0, 1;
%e A091602   3, 2, 1, 0, 1;
%e A091602   4, 3, 2, 1, 0, 1;
%e A091602   ...
%e A091602 (End)
%p A091602 g:=sum(t^k*(product((1-x^((k+1)*j))/(1-x^j),j=1..50)-product((1-x^(k*j))/(1-x^j),j=1..50)),k=1..50): gser:=simplify(series(g,x=0,20)): for n from 1 to 13 do P[n]:=coeff(gser,x^n) od: for n from 1 to 13 do seq(coeff(P[n],t^j),j=1..n) od;
%p A091602 # yields sequence in triangular form - _Emeric Deutsch_, Mar 30 2006
%p A091602 b:= proc(n, i, k) option remember; `if`(n=0, 1,
%p A091602       `if`(i>n, 0, add(b(n-i*j, i+1, min(k,
%p A091602        iquo(n-i*j, i+1))), j=0..min(n/i, k))))
%p A091602     end:
%p A091602 T:= (n, k)-> b(n, 1, k) -`if`(k=0, 0, b(n, 1, k-1)):
%p A091602 seq(seq(T(n, k), k=1..n), n=1..20);
%p A091602 # _Alois P. Heinz_, Nov 27 2013
%t A091602 b[n_, i_, k_] := b[n, i, k] = If[n == 0, 1, If[i>n, 0, Sum[b[n-i*j, i+1, Min[k, Quotient[n-i*j, i+1]]], {j, 0, Min[n/i, k]}]]]; t[n_, k_] := b[n, 1, k] - If[k == 0, 0, b[n, 1, k-1]]; Table[t[n, k], {n, 1, 20}, {k, 1, n}] // Flatten (* _Jean-François Alcover_, Jan 17 2014, after _Alois P. Heinz_'s second Maple program *)
%Y A091602 Row sums: A000041. Inverse: A091603. Square: A091604.
%Y A091602 Columns 1-6: A000009, A091605-A091609. Convergent of columns: A002865.
%Y A091602 Cf. A000009. - _Gary W. Adamson_, Mar 13 2010
%Y A091602 T(2n,n) gives: A232697.
%Y A091602 Cf. A213177, A264397.
%K A091602 nonn,tabl
%O A091602 1,4
%A A091602 _Christian G. Bower_, Jan 23 2004