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-3 of 3 results.

A054225 Triangle read by rows: row n (n>=0) gives the number of partitions of (n,0), (n-1,1), (n-2,2), ..., (0,n) respectively into sums of pairs.

Original entry on oeis.org

1, 1, 1, 2, 2, 2, 3, 4, 4, 3, 5, 7, 9, 7, 5, 7, 12, 16, 16, 12, 7, 11, 19, 29, 31, 29, 19, 11, 15, 30, 47, 57, 57, 47, 30, 15, 22, 45, 77, 97, 109, 97, 77, 45, 22, 30, 67, 118, 162, 189, 189, 162, 118, 67, 30, 42, 97, 181, 257, 323, 339, 323, 257, 181, 97, 42, 56, 139, 267, 401, 522, 589, 589, 522, 401, 267, 139, 56
Offset: 0

Views

Author

Marc LeBrun, Feb 04 2000

Keywords

Comments

By analogy with ordinary partitions (A000041). The empty partition gives T(0,0)=1 by definition. A054225 and A201377 give partitions of pairs into sums of distinct pairs. Parts (i,j) are "positive" in the sense that min {i,j} >= 0 and max {i,j} >0. The empty partition of (0,0) is counted as 1.
Or, triangle T(n,k) of bipartite partitions of n objects, k of which are black.
Or, number of ways to factor p^(n-k)*q^k where p and q are distinct primes.
In the paper by F. C. Auluck: "On partitions of bipartite numbers", p.74, in the formula for fixed m there should be factor 1/m!. The correct asymptotic formula is p(m, n) ~ (sqrt(6*n)/Pi)^m * exp(Pi*sqrt(2*n/3)) / (4*sqrt(3)*m!*n). - Vaclav Kotesovec, Feb 01 2016
T(n,k)=T(n,k-n) is the number of multiset partitions of the multiset {1^k, 2^(n-k)}, see example link. - Joerg Arndt, Jan 01 2024
Let R be the ring of power series in two countably infinite sets of variables x_1,y_1,x_2,y_2,... that are invariant under the diagonal action (i.e, the group S of permutations of positive integers acts by w(x_i)=x_{w(i)} and w(y_i)=y_{w(i)}). Then T(n,k) is the dimension of the (n,k)-bigraded piece of R, i.e., the bihomogeneous power series of degree n in the x-variables and k in the y-variables that are S-invariant. - Jeremy L. Martin, Nov 27 2024

Examples

			The second row (n=1) is 1,1 since (1,0) and (0,1) each have a single partition.
The third row (n=2) is 2, 2, 2 from (2,0) = (1,0)+(1,0), (1,1) = (1,0)+(0,1), (0,2) = (0,1)+(0,1).
In the fourth row (n=3), T(2,1)=4 from (2,1) = (2,0)+(0,1) = (1,0)+(1,1) = (1,0)+(1,0)+(0,1).
The triangle begins:
   1;
   1,  1;
   2,  2,  2;
   3,  4,  4,  3;
   5,  7,  9,  7,   5;
   7, 12, 16, 16,  12,  7;
  11, 19, 29, 31,  29, 19, 11;
  15, 30, 47, 57,  57, 47, 30, 15;
  22, 45, 77, 97, 109, 97, 77, 45, 22;
  ...
A further example: T(2,2) = 9:
[(2,2)],
[(2,1),(0,1)],
[(2,0),(0,2)],
[(2,0),(0,1),(0,1)],
[(1,2),(1,0)],
[(1,1),(1,1)],
[(1,1),(1,0),(0,1)],
[(1,0),(1,0),(0,2)],
[(1,0),(1,0),(0,1),(0,1)].
		

References

  • M. S. Cheema, Tables of partitions of Gaussian integers, National Institute of Sciences of India, New Delhi, 1956.
  • D. E. Knuth, The Art of Computer Programming, Vol. 4A, Table A-1, page 778. - N. J. A. Sloane, Dec 30 2018

Crossrefs

See A201376 for the same triangle formatted in a different way.
Row sums: A005380. a(2n, n): A002774. a(n, [n/2]): A091437. Cf. A060244.
The outer edges are T(n,0) = T(0,n) = A000041(n).
A054242 gives partitions into sums of distinct pairs.

Programs

  • Haskell
    see Zumkeller link.
  • Maple
    read transforms; t1 := mul( mul( 1/(1-x^(i-j)*y^j), j=0..i), i=1..11): SERIES2(t1,x,y,6);
  • Mathematica
    rows = 11; se = Series[ Product[ 1/(1-x^(n-k)*y^k), {n, 1, rows}, {k, 0, n}], {x, 0, rows}, {y, 0, rows}]; coes = CoefficientList[ se, {x, y}]; Flatten[ Table[ coes[[n-k+1, k]], {n, 1, rows+1}, {k, 1, n}]] (* Jean-François Alcover, Nov 21 2011, after g.f. *)
    p = 2; q = 3; b[n_, k_] := b[n, k] = If[n>k, 0, 1] + If[PrimeQ[n], 0, Sum[If[d>k, 0, b[n/d, d]], {d, DeleteCases[Divisors[n], 1|n]}]]; t[n_, k_] := b[p^(n-k)*q^k, p^(n-k)*q^k]; Table[t[n, k], {n, 0, 11}, {k, 0, n}] // Flatten (* Jean-François Alcover, Mar 13 2014, after Alois P. Heinz *)
  • PARI
    {T(n, k) = if( n<0 || k<0, 0, polcoeff( polcoeff( prod( i=1, n, prod( j=0, i, 1 / (1 - x^i * y^j), 1 + x * O(x^n))),n),k))} /* Michael Somos, Apr 19 2005 */
    

Formula

G.f.: Product_{i>=1, j=0..i} 1/(1-x^(i-j)*y^j).
Series ends ... + 7*x^5 + 12*x^4*y + 16*x^3*y^2 + 16*x^2*y^3 + 12*x*y^4 + 7*y^5 + 5*x^4 + 7*x^3*y + 9*x^2*y^2 + 7*x*y^3 + 5*y^4 + 3*x^3 + 4*x^2*y + 4*x*y^2 + 3*y^3 + 2*x^2 + 2*x*y + 2*y^2 + x + y + 1.

Extensions

Entry revised by N. J. A. Sloane, Nov 30 2011, to incorporate corrections provided by Reinhard Zumkeller, who also contributed the alternative version A201376. Once the errors were corrected, this sequence coincided with A060243, due to N. J. A. Sloane, Mar 22 2001, which included edits by Vladeta Jovovic, Mar 23 2001, and Christian G. Bower, Jan 08 2004. The two entries have now been merged.

A002774 Number of bipartite partitions of n white objects and n black ones.

Original entry on oeis.org

1, 2, 9, 31, 109, 339, 1043, 2998, 8406, 22652, 59521, 151958, 379693, 927622, 2224235, 5236586, 12130780, 27669593, 62229990, 138095696, 302673029, 655627975, 1404599867, 2977831389, 6251060785, 12999299705, 26791990052, 54750235190, 110977389012
Offset: 0

Views

Author

Keywords

Comments

Number of ways to factor p^n * q^n where p and q are distinct primes.

References

  • M. S. Cheema and H. Gupta, Tables of Partitions of Gaussian Integers. National Institute of Sciences of India, Mathematical Tables, Vol. 1, New Delhi, 1956, p. 1.
  • D. E. Knuth, The Art of Computer Programming, Vol. 4A, Table A-1, see p(n,n), page 778. - N. J. A. Sloane, Dec 30 2018
  • A. Murthy, Generalization of partition function, introducing Smarandache factor partitions. Smarandache Notions Journal, Vol. 11, No. 1-2-3, Spring 2000.
  • A. Murthy, Program for finding out the number of Smarandache factor partitions. (To be published in Smarandache Notions Journal).
  • Amarnath Murthy and Charles Ashbacher, Generalized Partitions and Some New Ideas on Number Theory and Smarandache Sequences, Hexis, Phoenix; USA 2005. See Section 1.4, 1.14.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A005380.
Cf. A219554. Column k=2 of A219727. - Alois P. Heinz, Nov 26 2012
Main diagonal of A054225 if that entry is drawn as a square array. - N. J. A. Sloane, Dec 30 2018

Programs

  • Maple
    with(numtheory):
    b:= proc(n, k) option remember; `if`(n>k, 0, 1) +`if`(isprime(n), 0,
          add(`if`(d>k, 0, b(n/d, d)), d=divisors(n) minus {1, n}))
        end:
    a:= n-> b(6^n$2):
    seq(a(n), n=0..30);  # Alois P. Heinz, Jun 27 2013
  • Mathematica
    max = 26; se = Series[ Sum[ Log[1 - x^(n-k)*y^k], {n, 1, 2max}, {k, 0, n}], {x, 0, 2max}, {y, 0, 2max}]; coes = CoefficientList[ Series[ Exp[-se], {x, 0, 2max}, {y, 0, 2max}], {x, y}]; a[n_] := coes[[n+1, n+1]]; Table[a[n], {n, 0, max} ](* Jean-François Alcover, Dec 06 2011 *)

Formula

a(n) = A054225(2n, n) = A091437(2n).
a(n) ~ Zeta(3)^(19/36) * exp(3*Zeta(3)^(1/3) * n^(2/3) + Pi^2 * n^(1/3) / (6*Zeta(3)^(1/3)) + Zeta'(-1) - Pi^4/(432*Zeta(3))) / (sqrt(3) * (2*Pi)^(3/2) * n^(55/36)). - Vaclav Kotesovec, Jan 30 2016
Formula (25) in the article by Auluck is incorrect. The correct formula is: p(n,n) ~ c^(19/12) * exp(3*c*n^(2/3) + 3*d*n^(1/3) + Zeta'(-1) - 3*d^2/(4*c)) / (sqrt(3) * (2*Pi)^(3/2) * n^(55/36)), where c = Zeta(3)^(1/3), d = Zeta(2)/(3*c). Also formula (24) is incorrect. - Vaclav Kotesovec, Jan 30 2016
From Vaclav Kotesovec, Feb 04 2016: (Start)
The correct formula (24) is p(m,n) ~ c^(7/4)/(2*Pi*sqrt(3)) * exp(3*c*(m*n)^(1/3) + 3*d*(m+n)/(2*(m*n)^(1/3)) - 19*log(m*n)/24 - ((m/n - 2*n/m)*log(m) + (n/m - 2*m/n)*log(n))/36 - (m/n + n/m)*(log(c)/12 + Zeta'(-1) - 1/12 + 3*d^2/(4*c)) + 3*d^2/(4*c) - 3*log(2*Pi)/4 + fi((n/m)^(1/2))),
where m and n are of the same order, c = Zeta(3)^(1/3), d = Zeta(2)/(3*c) and fi(alfa) = Integral_{t=0..infinity} (1/t)*(1/(exp(alfa*t)-1)/(exp(t/alfa)-1) - (alfa/t)/(exp(alfa*t)-1) - ((1/alfa)/t)/(exp(t/alfa)-1) + 1/t^2 + (1/4)/(exp(alfa*t)-1) + (1/4)/(exp(t/alfa)-1) - (alfa/4)/t - ((1/4)/alfa)/t).
If m = n then alfa = 1 and fi(1) = 3*Zeta'(-1) + log(2*Pi)/4 - 1/6.
For the asymptotic formula for fixed m see A054225.
(End)

Extensions

Corrected using A000491.
Edited by Christian G. Bower, Jan 08 2004

A201376 Triangle read by rows: T(n,k) (0 <= k <= n) is the number of partitions of (n,k) into a sum of pairs.

Original entry on oeis.org

1, 1, 2, 2, 4, 9, 3, 7, 16, 31, 5, 12, 29, 57, 109, 7, 19, 47, 97, 189, 339, 11, 30, 77, 162, 323, 589, 1043, 15, 45, 118, 257, 522, 975, 1752, 2998, 22, 67, 181, 401, 831, 1576, 2876, 4987, 8406, 30, 97, 267, 608, 1279, 2472, 4571, 8043, 13715, 22652, 42, 139, 392, 907, 1941, 3804, 7128, 12693, 21893, 36535, 59521
Offset: 0

Views

Author

Reinhard Zumkeller, Nov 30 2011

Keywords

Comments

By analogy with ordinary partitions (A000041). The empty partition gives T(0,0)=1 by definition. A201377 and A054225 give partitions of pairs into sums of distinct pairs.
Parts (i,j) are "positive" in the sense that min {i,j} >= 0 and max {i,j} >0. The empty partition of (0,0) is counted as 1.

Examples

			Partitions of (3,1) into positive pairs, T(3,1) = 7:
(3,1),
(3,0) + (0,1),
(2,1) + (1,0),
(2,0) + (1,1),
(2,0) + (1,0) + (0,1),
(1,1) + (1,0) + (1,0),
(1,0) + (1,0) + (1,0) + (0,1).
First ten rows of triangle:
0:                      1
1:                    1  2
2:                  2  4  9
3:                3  7  16  31
4:              5  12  29  57  109
5:            7  19  47  97  189  339
6:          11  30  77  162  323  589  1043
7:        15  45  118  257  522  975  1752  2998
8:      22  67  181  401  831  1576  2876  4987  8406
9:    30  97  267  608  1279  2472  4571  8043  13715  22652
X:  42  139  392  907  1941  3804  7128  12693  21893  36535  59521
		

Crossrefs

T(n,0) = A000041(n);
T(1,k) = A000070(k), k <= 1; T(n,1) = A000070(n), n > 1;
T(2,k) = A000291(k), k <= 2; T(n,2) = A000291(n), n > 2;
T(3,k) = A000412(k), k <= 3; T(n,3) = A000412(n), n > 3;
T(4,k) = A000465(k), k <= 4; T(n,4) = A000465(n), n > 4;
T(5,k) = A000491(k), k <= 5; T(n,5) = A000491(n), n > 5;
T(6,k) = A002755(k), k <= 6; T(n,6) = A002755(n), n > 6;
T(7,k) = A002756(k), k <= 7; T(n,7) = A002756(n), n > 7;
T(8,k) = A002757(k), k <= 8; T(n,8) = A002757(n), n > 8;
T(9,k) = A002758(k), k <= 9; T(n,9) = A002758(n), n > 9;
T(10,k) = A002759(n), k <= 10; T(n,10) = A002759(n), n > 10;
T(n,n) = A002774(n).
See A054225 for another version.

Programs

  • Haskell
    -- see link.
  • Mathematica
    max = 10; se = Series[ Sum[ Log[1 - x^(n-k)*y^k], {n, 1, 2max }, {k, 0, n}], {x, 0, 2max }, {y, 0, 2max }]; coes = CoefficientList[ Series[ Exp[-se], {x, 0, 2max }, {y, 0, 2max }], {x, y}]; t[n_, k_] := coes[[n+1, k+1]]; Flatten[ Table[ t[n, k], {n, 0, max}, {k, 0, n}]] (* Jean-François Alcover, Dec 05 2011 *)
    p = 2; q = 3; b[n_, k_] := b[n, k] = If[n > k, 0, 1] + If[PrimeQ[n], 0, Sum[If[d > k, 0, b[n/d, d]], {d, DeleteCases[Divisors[n] , 1|n]}]]; t[n_, k_] := b[p^n*q^k, p^n*q^k]; Table[t[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Jean-François Alcover, Mar 13 2014, after Alois P. Heinz *)

Formula

For references, programs and g.f. see A054225.

Extensions

Entry revised by N. J. A. Sloane, Nov 30 2011
Showing 1-3 of 3 results.