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.

Previous Showing 11-20 of 46 results. Next

A284593 Square array read by antidiagonals: T(n,k) = the number of pairs of partitions of n and k respectively, such that each partition is composed of distinct parts and the pair of partitions have no part in common.

Original entry on oeis.org

1, 1, 1, 1, 0, 1, 2, 1, 1, 2, 2, 1, 0, 1, 2, 3, 1, 1, 1, 1, 3, 4, 2, 2, 2, 2, 2, 4, 5, 2, 2, 2, 2, 2, 2, 5, 6, 3, 2, 3, 2, 3, 2, 3, 6, 8, 3, 3, 4, 3, 3, 4, 3, 3, 8, 10, 5, 4, 6, 5, 6, 5, 6, 4, 5, 10, 12, 5, 5, 6, 5, 6, 6, 5, 6, 5, 5, 12, 15, 7, 6, 8, 7, 8, 8, 8, 7, 8, 6, 7, 15
Offset: 0

Views

Author

Peter Bala, Mar 30 2017

Keywords

Comments

Compare with A284592.

Examples

			Square array begins
  n\k| 0  1  2  3  4  5  6   7   8   9  10  11  12  13
- - - - - - - - - - - - - - - - - - - - - - - - - - - -
  0  | 1  1  1  2  2  3  4   5   6   8  10  12  15  18: A000009
  1  | 1  0  1  1  1  2  2   3   3   5   5   7   8  10: A096765
  2  | 1  1  0  1  2  2  2   3   4   5   6   7   9  11: A015744
  3  | 2  1  1  2  2  3  4   6   6   8   9  12  15  18
  4  | 2  1  2  2  2  3  5   5   7   9  10  14  15  19
  5  | 3  2  2  3  3  6  6   8   9  12  16  19  22  28
  6  | 4  2  2  4  5  6  8   9  11  16  18  22  27  33
  7  | 5  3  3  6  5  8  9  14  16  20  23  29  34  41
  ...
T(3,7) = 6: the six pairs of partitions of 3 and 7 into distinct parts and with no parts in common are (3, 7), (3, 6 + 1), (3, 5 + 2), (3, 4 + 2 + 1), (2 + 1, 7) and (2 + 1, 4 + 3).
		

Crossrefs

Rows n=0..2 give A000009, A096765, A015744.
Main diagonal gives A365662.
Antidiagonal sums give A032302.

Programs

  • Maple
    # A284593 as a square array
    ser := taylor(taylor(mul(1 + x^j + y^j, j = 1..10), x, 11), y, 11):
    convert(ser, polynom):
    s := convert(%, polynom):
    with(PolynomialTools):
    for n from 0 to 10 do CoefficientList(coeff(s, y, n), x) end do;
    # second Maple program:
    b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
          b(n, i-1)+expand((x^i+1)*b(n-i, min(n-i, i-1)))))
        end:
    T:= (n, k)-> coeff(b(n+k$2), x, k):
    seq(seq(T(n, d-n), n=0..d), d=0..14);  # Alois P. Heinz, Aug 24 2019
  • Mathematica
    nmax = 12; M = CoefficientList[#, y][[;; nmax+1]]& /@ (Product[1 + x^j + y^j, {j, 1, nmax}] + O[x]^(nmax+1) // CoefficientList[#, x]& // Expand);
    T[n_, k_] := M[[n+1, k+1]];
    Table[T[n-k, k], {n, 0, nmax}, {k, 0, n}] // Flatten (* Jean-François Alcover, Dec 07 2019 *)

Formula

O.g.f. Product_{j >= 1} (1 + x^j + y^j) = Sum_{n,k >= 0} T(n,k)*x^n*y^k (see Wilf, Example 7).
Antidiagonal sums are A032302.

A365662 Number of ordered pairs of disjoint strict integer partitions of n.

Original entry on oeis.org

1, 0, 0, 2, 2, 6, 8, 14, 18, 32, 42, 66, 92, 136, 190, 280, 374, 532, 744, 1014, 1366, 1896, 2512, 3384, 4526, 6006, 7910, 10496, 13648, 17842, 23338, 30116, 38826, 50256, 64298, 82258, 105156, 133480, 169392, 214778, 270620, 340554, 428772, 536302, 670522
Offset: 0

Views

Author

Gus Wiseman, Sep 19 2023

Keywords

Comments

Also the number of ways to first choose a strict partition of 2n, then a subset of it summing to n.

Examples

			The a(0) = 1 through a(7) = 14 pairs:
  ()()  .  .  (21)(3)  (31)(4)  (32)(5)   (42)(6)   (43)(7)
              (3)(21)  (4)(31)  (41)(5)   (51)(6)   (52)(7)
                                (5)(32)   (6)(42)   (61)(7)
                                (5)(41)   (6)(51)   (7)(43)
                                (32)(41)  (321)(6)  (7)(52)
                                (41)(32)  (42)(51)  (7)(61)
                                          (51)(42)  (421)(7)
                                          (6)(321)  (43)(52)
                                                    (43)(61)
                                                    (52)(43)
                                                    (52)(61)
                                                    (61)(43)
                                                    (61)(52)
                                                    (7)(421)
		

Crossrefs

For subsets instead of partitions we have A000244, non-disjoint A000302.
If the partitions can have different sums we get A032302.
The non-strict version is A054440, non-disjoint A001255.
The unordered version is A108796, non-strict A260669.
A000041 counts integer partitions, strict A000009.
A000124 counts distinct possible sums of subsets of {1..n}.
A000712 counts distinct submultisets of partitions.
A002219 and A237258 count partitions of 2n including a partition of n.
A304792 counts subset-sums of partitions, positive A276024, strict A284640.
A364272 counts sum-full strict partitions, sum-free A364349.

Programs

  • Mathematica
    Table[Length[Select[Tuples[Select[IntegerPartitions[n], UnsameQ@@#&],2], Intersection@@#=={}&]], {n,0,15}]
    Table[SeriesCoefficient[Product[(1 + x^k + y^k), {k, 1, n}], {x, 0, n}, {y, 0, n}], {n, 0, 50}] (* Vaclav Kotesovec, Apr 24 2025 *)

Formula

a(n) = 2*A108796(n) for n > 1.
a(n) = [(x*y)^n] Product_{k>=1} (1 + x^k + y^k). - Ilya Gutkovskiy, Apr 24 2025

A303346 Expansion of Product_{n>=1} ((1 + 2*x^n)/(1 - 2*x^n))^(1/2).

Original entry on oeis.org

1, 2, 4, 10, 18, 38, 72, 142, 260, 510, 940, 1814, 3362, 6490, 12112, 23466, 44114, 85766, 162516, 317190, 604806, 1184682, 2271248, 4461514, 8591784, 16916490, 32696708, 64496130, 125037142, 247007142, 480077432, 949510526, 1849375796, 3661330398, 7144215452
Offset: 0

Views

Author

Seiichi Manyama, Apr 22 2018

Keywords

Crossrefs

Programs

  • Mathematica
    nmax = 30; CoefficientList[Series[Product[((1 + 2*x^k)/(1 - 2*x^k))^(1/2), {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Apr 22 2018 *)
    nmax = 30; CoefficientList[Series[Sqrt[-QPochhammer[-2, x] / (3*QPochhammer[2, x])], {x, 0, nmax}], x] (* Vaclav Kotesovec, Apr 22 2018 *)
  • PARI
    N=66; x='x+O('x^N); Vec(prod(k=1, N, ((1+2*x^k)/(1-2*x^k))^(1/2)))

Formula

a(n) ~ 2^n / sqrt(c*Pi*n), where c = A048651 * A083864 = 1/2 * Product_{j>=1} (2^j-1)/(2^j+1) = 0.06056210400129025123042464659093375290492912341... - Vaclav Kotesovec, Apr 22 2018

A070877 Expansion of Product_{k>=1} (1 - 2x^k).

Original entry on oeis.org

1, -2, -2, 2, 2, 6, -2, 2, -6, -10, -2, -6, -6, 6, 22, -6, 26, 14, 22, -6, -14, -2, -10, -46, -46, -50, -18, 18, -78, 22, 14, 82, 42, 166, 14, 42, 170, 118, 6, 106, -150, -66, -122, -118, -62, -370, -282, -350, -126, -354, -2, -94, 226, -250, 30, 450, 730, 342, 894, 474, 890, 358, 758, 58, 1210, 782, -778, 26, -270, -1250
Offset: 0

Views

Author

Sharon Sela (sharonsela(AT)hotmail.com), May 24 2002

Keywords

Examples

			G.f. = 1 - 2*x - 2*x^2 + 2*x^3 + 2*x^4 + 6*x^5 - 2*x^6 + 2*x^7 - 6*x^8 - 10*x^9 + ...
		

Crossrefs

Programs

  • Mathematica
    CoefficientList[ Series[ Product[(1 - 2t^k), {k, 1, 80}], {t, 0, 80}], t]
    a[ n_] := SeriesCoefficient[ -QPochhammer[2, x], {x, 0, n}]; (* Michael Somos, Mar 11 2014 *)
  • PARI
    N=66; q='q+O('q^N); Vec(sum(n=0, N, (-2)^n*q^(n*(n+1)/2) / prod(k=1, n, 1-q^k ) )) \\ Joerg Arndt, Mar 09 2014
    
  • PARI
    N=66; q='q+O('q^N); t2=Vec( prod(k=1, N, 1-2*q^k) ) \\ Joerg Arndt, Mar 11 2014

Extensions

Edited by Robert G. Wilson v, May 26 2002
Corrected by Vincenzo Librandi, Mar 11 2014

A261569 Expansion of Product_{k>=1} (1 + 5*x^k).

Original entry on oeis.org

1, 5, 5, 30, 30, 55, 180, 205, 330, 480, 1230, 1380, 2255, 3030, 4530, 8555, 10680, 15330, 21330, 29730, 39480, 67380, 81505, 116280, 153030, 210930, 270805, 370080, 534330, 675480, 900480, 1180380, 1544130, 1997280, 2597280, 3304805, 4581180, 5653080
Offset: 0

Views

Author

Vaclav Kotesovec, Aug 24 2015

Keywords

Comments

In general, for a fixed integer m > 0, if g.f. = Product_{k>=1} (1 + m*x^k) then a(n) ~ c^(1/4) * exp(2*sqrt(c*n)) / (2*sqrt((m+1)*Pi)*n^(3/4)), where c = Pi^2/6 + log(m)^2/2 + polylog(2, -1/m) = -polylog(2, -m). - Vaclav Kotesovec, Jan 04 2016

Crossrefs

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(i*(i+1)/2n, 0, 5*b(n-i, i-1))))
        end:
    a:= n-> b(n$2):
    seq(a(n), n=0..60);  # Alois P. Heinz, Aug 24 2015
  • Mathematica
    nmax = 40; CoefficientList[Series[Product[1 + 5*x^k, {k, 1, nmax}], {x, 0, nmax}], x]
    nmax = 40; CoefficientList[Series[Exp[Sum[(-1)^(k+1)*5^k/k*x^k/(1-x^k), {k, 1, nmax}]], {x, 0, nmax}], x] (* Vaclav Kotesovec, Aug 25 2015 *)
    (QPochhammer[-5, x]/6 + O[x]^58)[[3]] (* Vladimir Reshetnikov, Nov 20 2015 *)

Formula

a(n) ~ c^(1/4) * exp(2*sqrt(c*n)) / (2*sqrt(6*Pi)*n^(3/4)), where c = Pi^2/6 + log(5)^2/2 + polylog(2, -1/5) = 2.74927912606080829002558751537626864449... . - Vaclav Kotesovec, Jan 04 2016
G.f.: Sum_{i>=0} 5^i*x^(i*(i+1)/2)/Product_{j=1..i} (1 - x^j). - Ilya Gutkovskiy, Apr 12 2018

A268498 Expansion of Product_{k>=1} ((1 + 2*x^k) / (1 + x^k)).

Original entry on oeis.org

1, 1, 0, 3, -1, 3, 3, 3, 0, 4, 12, 0, 9, -3, 21, 12, 17, -3, 33, 0, 33, 36, 36, 27, 21, 52, 24, 90, 72, 99, 24, 138, 21, 207, 0, 261, 149, 267, 45, 333, 174, 339, 174, 345, 411, 654, 330, 456, 657, 535, 684, 483, 1233, 489, 1353, 882, 1803, 720, 1902, 756
Offset: 0

Views

Author

Vaclav Kotesovec, Feb 06 2016

Keywords

Comments

It appears that this sequence contains only finitely many nonpositive terms, namely at indices {2, 4, 8, 11, 13, 17, 19, 34}. - Gus Wiseman, Jan 23 2019

Crossrefs

Programs

  • Mathematica
    nmax = 100; CoefficientList[Series[Product[(1+2*x^k)/(1+x^k), {k, 1, nmax}], {x, 0, nmax}], x]

Formula

a(n) ~ c^(1/4) * exp(sqrt(c*n)) / (2*sqrt(3*Pi)*n^(3/4)), where c = Pi^2/3 + 2*log(2)^2 + 4*polylog(2, -1/2) = 2.4571173338382709125... .
a(n) = Sum_{k = 0...n} (-1)^k * A133121(n,k). - Gus Wiseman, Jan 23 2019
G.f.: Product_{k>=1} (1 - Sum_{j>=1} (-1)^j * x^(k*j)). - Ilya Gutkovskiy, Nov 06 2019

A264686 Expansion of Product_{k>=1} (1 + 2*x^k)/(1 - x^k).

Original entry on oeis.org

1, 3, 6, 15, 27, 51, 93, 159, 264, 432, 696, 1086, 1683, 2553, 3837, 5700, 8367, 12147, 17505, 24972, 35361, 49728, 69402, 96243, 132657, 181782, 247692, 335838, 453042, 608289, 813102, 1082256, 1434519, 1894215, 2491644, 3265869, 4265973, 5553771, 7207167
Offset: 0

Views

Author

Vaclav Kotesovec, Nov 21 2015

Keywords

Comments

Convolution of A000041 and A032302.

Crossrefs

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(i*(i+1)/2n, 0, 2*b(n-i, i-1))))
        end:
    a:= n-> add(b(i$2)*combinat[numbpart](n-i), i=0..n):
    seq(a(n), n=0..60);  # Alois P. Heinz, Dec 22 2017
  • Mathematica
    nmax = 40; CoefficientList[Series[Product[(1 + 2*x^k)/(1 - x^k), {k, 1, nmax}], {x, 0, nmax}], x]
  • PARI
    { my(n=40); Vec(prod(k=1, n, 3/(1-x^k) - 2 + O(x*x^n))) } \\ Andrew Howroyd, Dec 22 2017

Formula

a(n) ~ sqrt(c) * exp(sqrt(2*c*n)) / (4*Pi*sqrt(3)*n), where c = 2*Pi^2/3 + log(2)^2 + 2*polylog(2, -1/2) = 6.163360867463814765670634381079217086937812673723341... . - Vaclav Kotesovec, Jan 04 2016

A286957 Square array A(n,k), n>=0, k>=0, read by antidiagonals, where column k is the expansion of Product_{j>=1} (1 + k*x^j).

Original entry on oeis.org

1, 1, 0, 1, 1, 0, 1, 2, 1, 0, 1, 3, 2, 2, 0, 1, 4, 3, 6, 2, 0, 1, 5, 4, 12, 6, 3, 0, 1, 6, 5, 20, 12, 10, 4, 0, 1, 7, 6, 30, 20, 21, 18, 5, 0, 1, 8, 7, 42, 30, 36, 48, 22, 6, 0, 1, 9, 8, 56, 42, 55, 100, 57, 30, 8, 0, 1, 10, 9, 72, 56, 78, 180, 116, 84, 42, 10, 0, 1, 11, 10, 90, 72, 105, 294, 205, 180, 120, 66, 12, 0
Offset: 0

Views

Author

Ilya Gutkovskiy, May 17 2017

Keywords

Comments

A(n,k) is the number of partitions of n into distinct parts of k sorts: the parts are unordered, but not the sorts.

Examples

			Square array begins:
1,  1,   1,   1,   1,   1,  ...
0,  1,   2,   3,   4,   5,  ...
0,  1,   2,   3,   4,   5,  ...
0,  2,   6,  12,  20,  30,  ...
0,  2,   6,  12,  20,  30,  ...
0,  3,  10,  21,  36,  55,  ...
		

Crossrefs

Columns k=0-5 give: A000007, A000009, A032302, A032308, A261568, A261569.
Main diagonal gives A291698.
Cf. A246935.

Programs

  • Mathematica
    Table[Function[k, SeriesCoefficient[Product[(1 + k x^i), {i, 1, Infinity}], {x, 0, n}]][j - n], {j, 0, 12}, {n, 0, j}] // Flatten
    Table[Function[k, SeriesCoefficient[QPochhammer[-k, x]/(1 + k), {x, 0, n}]][j - n], {j, 0, 12}, {n, 0, j}] // Flatten

Formula

G.f. of column k: Product_{j>=1} (1 + k*x^j).

A370709 a(n) = 2^n * [x^n] Product_{k>=1} (1 + 2*x^k)^(1/2).

Original entry on oeis.org

1, 2, 2, 20, 6, 108, 148, 776, -186, 5964, -4, 51032, -89700, 512120, -1259416, 6406032, -19733434, 78363148, -268823572, 1047941688, -3800035916, 14327505832, -52766730600, 199492430192, -746479735524, 2811936761016, -10588174502568, 40092283176560, -151796846803592
Offset: 0

Views

Author

Vaclav Kotesovec, Feb 27 2024

Keywords

Crossrefs

Programs

  • Mathematica
    nmax = 30; CoefficientList[Series[Product[(1 + 2*x^k), {k, 1, nmax}]^(1/2), {x, 0, nmax}], x] * 2^Range[0, nmax]
    nmax = 30; CoefficientList[Series[Product[(1 + 2*(2*x)^k), {k, 1, nmax}]^(1/2), {x, 0, nmax}], x]
    nmax = 30; CoefficientList[Series[Sqrt[QPochhammer[-2, x]/3], {x, 0, nmax}], x] * 2^Range[0, nmax]

Formula

G.f.: Product_{k>=1} (1 + 2*(2*x)^k)^(1/2).
a(n) ~ (-1)^(n+1) * c * 4^n / n^(3/2), where c = QPochhammer(-1/2)^(1/2) / (2*sqrt(Pi)) = 0.31039710860287467176143051675437...

A370716 a(n) = 3^(2*n) * [x^n] Product_{k>=1} (1 + 2*x^k)^(1/3).

Original entry on oeis.org

1, 6, 18, 1170, -1890, 133326, 101250, 20498994, -164656314, 3778220862, -28085954094, 771567716970, -10691904063114, 183594050113518, -2711145260068326, 49416883617381354, -789899109743435994, 13176840267952166070, -216403389726994588086, 3681309971143060236810
Offset: 0

Views

Author

Vaclav Kotesovec, Feb 27 2024

Keywords

Crossrefs

Programs

  • Mathematica
    nmax = 20; CoefficientList[Series[Product[(1 + 2*x^k), {k, 1, nmax}]^(1/3), {x, 0, nmax}], x] * 3^(2*Range[0, nmax])
    nmax = 20; CoefficientList[Series[Product[(1 + 2*(9*x)^k), {k, 1, nmax}]^(1/3), {x, 0, nmax}], x]
    nmax = 20; CoefficientList[Series[(QPochhammer[-2, x]/3)^(1/3), {x, 0, nmax}], x] * 3^(2*Range[0, nmax])

Formula

G.f.: Product_{k>=1} (1 + 2*(9*x)^k)^(1/3).
a(n) ~ (-1)^(n+1) * c * 18^n / n^(4/3), where c = QPochhammer(-1/2)^(1/3) / (3*Gamma(2/3)) = 0.2623638446186535909018671540030519...
Previous Showing 11-20 of 46 results. Next