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.

Showing 1-5 of 5 results.

A115994 Triangle read by rows: T(n,k) is number of partitions of n with Durfee square of size k (n>=1; 1<=k<=floor(sqrt(n))).

Original entry on oeis.org

1, 2, 3, 4, 1, 5, 2, 6, 5, 7, 8, 8, 14, 9, 20, 1, 10, 30, 2, 11, 40, 5, 12, 55, 10, 13, 70, 18, 14, 91, 30, 15, 112, 49, 16, 140, 74, 1, 17, 168, 110, 2, 18, 204, 158, 5, 19, 240, 221, 10, 20, 285, 302, 20, 21, 330, 407, 34, 22, 385, 536, 59, 23, 440, 698, 94, 24, 506, 896, 149, 25
Offset: 1

Views

Author

Emeric Deutsch, Feb 11 2006

Keywords

Comments

Row n has floor(sqrt(n)) terms. Row sums yield A000041. Column 2 yields A006918. sum(k*T(n,k),k=1..floor(sqrt(n)))=A115995.
T(n,k) is number of partitions of n-k^2 into parts of 2 kinds with at most k of each kind.
The limit of the diagonals is A000712 (partitions into parts of two kinds). In particular, if 0<=m<=n, T(n(n+1)/2 + m, n) = A000712(m). These partitions in this range can be viewed as an equilateral right triangle of side n, with one partition appended on the top (at the left) and another appended on the right. - Franklin T. Adams-Watters, Jan 11 2006
Successive columns approach closer and closer to A000712. - N. J. A. Sloane, Mar 10 2007

Examples

			T(5,2) = 2 because the only partitions of 5 having Durfee square of size 2 are [3,2] and [2,2,1]; the other five partitions ([5], [4,1], [3,1,1], [2,1,1,1] and [1,1,1,1,1]) have Durfee square of size 1.
Triangle starts:
  1;
  2;
  3;
  4,  1;
  5,  2;
  6,  5;
  7,  8;
  8, 14;
  9, 20,  1;
  ...
		

References

  • G. E. Andrews, The Theory of Partitions, Addison-Wesley, 1976 (pp. 27-28).
  • G. E. Andrews and K. Eriksson, Integer Partitions, Cambridge Univ. Press, 2004 (pp. 75-78).

Crossrefs

For another version see A115720. Row lengths A000196.

Programs

  • Maple
    g:=sum(t^k*q^(k^2)/product((1-q^j)^2,j=1..k),k=1..40): gser:=series(g,q=0,32): for n from 1 to 27 do P[n]:=coeff(gser,q^n) od: for n from 1 to 27 do seq(coeff(P[n],t^j),j=1..floor(sqrt(n))) od; # yields sequence in triangular form
    # second Maple program:
    b:= proc(n, i) option remember;
          `if`(n=0, 1, `if`(i<1, 0, b(n, i-1)+`if`(i>n, 0, b(n-i, i))))
        end:
    T:= (n, k)-> add(b(m, k)*b(n-k^2-m, k), m=0..n-k^2):
    seq(seq(T(n, k), k=1..floor(sqrt(n))), n=1..30); # Alois P. Heinz, Apr 09 2012
  • Mathematica
    Map[Select[#,#>0&]&,Drop[Transpose[Table[CoefficientList[ Series[x^(n^2)/Product[1-x^i,{i,1,n}]^2,{x,0,nn}],x],{n,1,10}]],1]] //Grid (* Geoffrey Critzer, Sep 27 2013 *)
    b[n_, i_] := b[n, i] = If[n==0, 1, If[i<1, 0, b[n, i-1] + If[i>n, 0, b[n-i, i]]]]; T[n_, k_] := Sum[b[m, k]*b[n-k^2-m, k], {m, 0, n-k^2}]; Table[T[n, k], {n, 1, 30}, {k, 1, Sqrt[n]}] // Flatten (* Jean-François Alcover, Dec 25 2015, after Alois P. Heinz *)

Formula

G.f.: sum(k>=1, t^k*q^(k^2)/product(j=1..k, (1-q^j)^2 ) ).
T(n,k) = Sum_{i=0}^{n-k^2} P*(i,k)*P*(n-k^2-i), where P*(n,k) = P(n+k,k) is the number of partitions of n objects into at most k parts.

Extensions

Edited and verified by Franklin T. Adams-Watters, Mar 11 2006

A115720 Triangle T(n,k) is the number of partitions of n with Durfee square k.

Original entry on oeis.org

1, 0, 1, 0, 2, 0, 3, 0, 4, 1, 0, 5, 2, 0, 6, 5, 0, 7, 8, 0, 8, 14, 0, 9, 20, 1, 0, 10, 30, 2, 0, 11, 40, 5, 0, 12, 55, 10, 0, 13, 70, 18, 0, 14, 91, 30, 0, 15, 112, 49, 0, 16, 140, 74, 1, 0, 17, 168, 110, 2, 0, 18, 204, 158, 5, 0, 19, 240, 221, 10, 0, 20, 285, 302, 20, 0, 21, 330, 407
Offset: 0

Views

Author

Keywords

Comments

T(n,k) is number of partitions of n-k^2 into parts of 2 kinds with at most k of each kind.

Examples

			Triangle starts:
  1;
  0,  1;
  0,  2;
  0,  3;
  0,  4,  1;
  0,  5,  2;
  0,  6,  5;
  0,  7,  8;
  0,  8, 14;
  0,  9, 20,  1;
  0, 10, 30,  2;
From _Gus Wiseman_, Apr 12 2019: (Start)
Row n = 9 counts the following partitions:
  (9)          (54)       (333)
  (81)         (63)
  (711)        (72)
  (6111)       (432)
  (51111)      (441)
  (411111)     (522)
  (3111111)    (531)
  (21111111)   (621)
  (111111111)  (3222)
               (3321)
               (4221)
               (4311)
               (5211)
               (22221)
               (32211)
               (33111)
               (42111)
               (222111)
               (321111)
               (2211111)
(End)
		

Crossrefs

For a version without zeros see A115994. Row lengths are A003059. Row sums are A000041. Column k = 2 is A006918. Column k = 3 is A117485.
Related triangles are A096771, A325188, A325189, A325192, with Heinz-encoded versions A263297, A325169, A065770, A325178.

Programs

  • Maple
    b:= proc(n, i) option remember;
          `if`(n=0, 1, `if`(i<1, 0, b(n, i-1)+`if`(i>n, 0, b(n-i, i))))
        end:
    T:= (n, k)-> add(b(m, k)*b(n-k^2-m, k), m=0..n-k^2):
    seq(seq(T(n, k), k=0..floor(sqrt(n))), n=0..30); # Alois P. Heinz, Apr 09 2012
  • Mathematica
    b[n_, i_] := b[n, i] = If[n==0, 1, If[i<1, 0, b[n, i-1] + If[i>n, 0, b[n-i, i]]]]; T[n_, k_] := Sum[b[m, k]*b[n-k^2-m, k], {m, 0, n-k^2}]; Table[ T[n, k], {n, 0, 30}, {k, 0, Sqrt[n]}] // Flatten (* Jean-François Alcover, Dec 03 2015, after Alois P. Heinz *)
    durf[ptn_]:=Length[Select[Range[Length[ptn]],ptn[[#]]>=#&]];
    Table[Length[Select[IntegerPartitions[n],durf[#]==k&]],{n,0,10},{k,0,Sqrt[n]}] (* Gus Wiseman, Apr 12 2019 *)

Formula

T(n,k) = Sum_{i=0..n-k^2} P*(i,k)*P*(n-k^2-i), where P*(n,k) = P(n+k,k) is the number of partitions of n objects into at most k parts.

A115995 Sum of the sizes of the Durfee squares of all partitions of n.

Original entry on oeis.org

0, 1, 2, 3, 6, 9, 16, 23, 36, 52, 76, 106, 152, 207, 286, 386, 522, 691, 920, 1202, 1576, 2038, 2636, 3373, 4320, 5478, 6944, 8738, 10984, 13717, 17116, 21232, 26308, 32441, 39944, 48977, 59970, 73147, 89090, 108151, 131090, 158417, 191166, 230049, 276444
Offset: 0

Views

Author

Emeric Deutsch, Feb 11 2006

Keywords

Comments

Also sum of positive cranks of all partitions of n, n>1; see A064391. - Vladeta Jovovic, Oct 20 2006
This sequence, its author and the author of the above comment were mentioned in the Andrews-Chan-Kim paper, where it is called C_1 (see the remark on page 6). - Omar E. Pol, Apr 06 2012

Examples

			a(4) = 6 because the partitions [4], [3,1], [2,2], [2,1,1] and [1,1,1,1] of 4 have Durfee squares of sizes 1,1,2,1 and 1, respectively.
		

References

  • G. E. Andrews, The Theory of Partitions, Addison-Wesley, 1976 (pp. 27-28).
  • G. E. Andrews and K. Eriksson, Integer Partitions, Cambridge Univ. Press, 2004 (pp. 75-78).

Crossrefs

Programs

  • Maple
    g:= add(k*z^(k^2)/mul((1-z^j)^2,j=1..k),k=1..10): gser:=series(g,z=0,56): seq(coeff(gser,z,n), n=0..52);
    # second Maple program:
    b:= proc(n, i) option remember;
          `if`(n=0, 1, `if`(i<1, 0, b(n, i-1)+`if`(i>n, 0, b(n-i, i))))
        end:
    a:= n-> add(add(b(k, d)*b(n-d^2-k, d), k=0..n-d^2)*d, d=1..isqrt(n)):
    seq(a(n), n=0..70);  # Alois P. Heinz, Apr 09 2012
    # Third Maple program, based on Theorem 1 of Andrews-Chan-Kim:
    M:=101;
    qinf:=mul(1-q^i,i=1..M);
    qinf:=series(qinf,q,M);
    C1:=add((-1)^(n+1)*q^(n*(n+1)/2)/(1-q^n),n=1..M);
    C1:=series(C1/qinf,q,M);
    seriestolist(%); # N. J. A. Sloane, Sep 04 2012
  • Mathematica
    b[n_, i_] := b[n, i] = If[n == 0, 1, If[i < 1, 0, b[n, i - 1] + If[i > n, 0, b[n - i, i]]]] ; a[n_] := Sum[ Sum[b[k, d]*b[n - d^2 - k, d], {k, 0, n - d^2}]*d, {d, 1, Sqrt[n]}]; Table [a[n], {n, 0, 70}] (* Jean-François Alcover, Jan 16 2015, after Alois P. Heinz *)
  • PARI
    N=66; x='x+O('x^N); concat([0], Vec( sum(n=0,N, n*x^(n^2) / prod(k=1,n, 1-x^k)^2))) \\ Joerg Arndt, Mar 26 2014
    
  • Sage
    [sum(p.frobenius_rank() for p in Partitions(n)) for n in range(45)] # Peter Luschny, Sep 15 2014

Formula

G.f.: Sum_{k>=1} (k*z^(k^2) / Product_{j=1..k} (1 - z^j)^2 ).
a(n) = Sum_{k=1..floor(sqrt(n))} k*A115994(n,k).
Convolution of A067742 and A000041. - Vladeta Jovovic, Oct 20 2006
a(n) = A195012(n) + A209616(n), n >= 1. - Omar E. Pol, Apr 06 2012
a(n) ~ log(2) * exp(Pi*sqrt(2*n/3)) / (2^(3/2)*Pi*sqrt(n)). - Vaclav Kotesovec, Jan 02 2019

Extensions

Edited and verified by Franklin T. Adams-Watters, Mar 11 2006

A115722 Table of Durfee square of partitions in Mathematica order.

Original entry on oeis.org

0, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 1, 2, 1, 2, 2, 1, 2, 2, 1, 1, 1, 1, 2, 1, 2, 2, 1, 2, 2, 2, 1, 2, 2, 1, 1, 1, 1, 2, 1, 2, 2, 1, 2, 2, 2, 2, 1, 2, 2, 2, 2, 1, 2, 2, 2, 1, 1, 1, 1, 2, 1, 2, 2, 1, 2, 2, 2, 2, 1, 2, 2, 2, 2, 2, 1, 3, 2, 2, 2, 2, 2, 1, 2, 2, 2, 1, 1, 1, 1, 2, 1, 2, 2, 1, 2
Offset: 0

Views

Author

Keywords

Examples

			First few rows:
0;
1,1;
1,1,1;
1,1,2,1,1;
1,1,2,1,2,1,1;
		

Crossrefs

Row lengths A000041, totals A115995.

Programs

  • Mathematica
    {0}~Join~Array[Map[Block[{k = Length@ #}, While[Nand[k > 0, AllTrue[Take[#, k], # >= k &]], k--]; k] &, IntegerPartitions@ #] &, 10] // Flatten (* Michael De Vlieger, Jan 17 2020 *)

Formula

If partition is laid out in descending order p(1),p(2),...,p(k) without repetition factors (e.g. [3,2,2,1,1,1]), a(P) = max_k min(k,p(k)).

A330376 Irregular triangle read by rows: T(n,k) is the total number of parts in all partitions of n with Durfee square of size k (n>=1; 1<=k<=floor(sqrt(n))).

Original entry on oeis.org

1, 3, 6, 10, 2, 15, 5, 21, 14, 28, 26, 36, 50, 45, 80, 3, 55, 130, 7, 66, 190, 19, 78, 280, 41, 91, 385, 80, 105, 532, 143, 120, 700, 248, 136, 924, 399, 4, 153, 1176, 627, 9, 171, 1500, 949, 24, 190, 1860, 1397, 51, 210, 2310, 2003, 107, 231, 2805, 2823, 193
Offset: 1

Views

Author

Omar E. Pol, Dec 22 2019

Keywords

Examples

			Triangle begins:
   1;
   3;
   6;
  10,  2;
  15,  5;
  21, 14;
  28, 26;
  36, 50;
  45, 80, 3;
		

Crossrefs

Row sums give A006128, n >= 1.
Column 1 gives A000217, n >= 1.
Cf. A330369.

Programs

  • PARI
    \\ by enumeration over partitions.
    ds(p)={for(i=2, #p, if(p[#p+1-i]Andrew Howroyd, Feb 02 2022
    
  • PARI
    \\ by generating function.
    P(n,k,y)={1/prod(j=1, k, 1 - y*x^j + O(x*x^n))}
    T(n,k)={my(r=n-k^2); if(r<0, 0, subst(deriv(polcoef(y^k*P(r,k,1)*P(r,k,y), r)), y, 1))}
    { for(n=1, 10, print(vector(sqrtint(n), k, T(n,k)))) } \\ Andrew Howroyd, Feb 02 2022

Extensions

Terms a(10) and beyond from Andrew Howroyd, Feb 02 2022
Showing 1-5 of 5 results.