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 21-30 of 87 results. Next

A114921 Number of unimodal compositions of n+2 where the maximal part appears exactly twice.

Original entry on oeis.org

1, 0, 1, 2, 4, 6, 11, 16, 27, 40, 63, 92, 141, 202, 299, 426, 614, 862, 1222, 1694, 2362, 3242, 4456, 6054, 8229, 11072, 14891, 19872, 26477, 35050, 46320, 60866, 79827, 104194, 135703, 176008, 227791, 293702, 377874, 484554, 620011, 790952, 1006924
Offset: 0

Views

Author

Michael Somos, Jan 07 2006

Keywords

Comments

Old name was: Expansion of a q-series.
a(n) is also the number of 2-colored partitions of n with the same number of parts in each color. - Shishuo Fu, May 30 2017
From Gus Wiseman, Mar 25 2021: (Start)
Also the number of even-length compositions of n with alternating parts weakly decreasing. Allowing odd lengths also gives A342528. The version with alternating parts strictly decreasing appears to be A064428. The a(2) = 1 through a(7) = 16 compositions are:
(1,1) (1,2) (1,3) (1,4) (1,5) (1,6)
(2,1) (2,2) (2,3) (2,4) (2,5)
(3,1) (3,2) (3,3) (3,4)
(1,1,1,1) (4,1) (4,2) (4,3)
(1,2,1,1) (5,1) (5,2)
(2,1,1,1) (1,2,1,2) (6,1)
(1,3,1,1) (1,3,1,2)
(2,1,2,1) (1,4,1,1)
(2,2,1,1) (2,2,1,2)
(3,1,1,1) (2,2,2,1)
(1,1,1,1,1,1) (2,3,1,1)
(3,1,2,1)
(3,2,1,1)
(4,1,1,1)
(1,2,1,1,1,1)
(2,1,1,1,1,1)
(End)

Examples

			From _Joerg Arndt_, Jun 10 2013: (Start)
There are a(7)=16 such compositions of 7+2=9 where the maximal part appears twice:
  01:  [ 1 1 1 1 1 2 2 ]
  02:  [ 1 1 1 1 2 2 1 ]
  03:  [ 1 1 1 2 2 1 1 ]
  04:  [ 1 1 1 3 3 ]
  05:  [ 1 1 2 2 1 1 1 ]
  06:  [ 1 1 3 3 1 ]
  07:  [ 1 2 2 1 1 1 1 ]
  08:  [ 1 2 3 3 ]
  09:  [ 1 3 3 1 1 ]
  10:  [ 1 3 3 2 ]
  11:  [ 1 4 4 ]
  12:  [ 2 2 1 1 1 1 1 ]
  13:  [ 2 3 3 1 ]
  14:  [ 3 3 1 1 1 ]
  15:  [ 3 3 2 1 ]
  16:  [ 4 4 1 ]
(End)
		

Crossrefs

Cf. A226541 (max part appears three times), A188674 (max part m appears m times), A001523 (max part appears any number of times).
Column k=2 of A247255.
A000041 counts weakly increasing (or weakly decreasing) compositions.
A000203 adds up divisors.
A002843 counts compositions with all adjacent parts x <= 2y.
A003242 counts anti-run compositions.
A034008 counts even-length compositions.
A065608 counts even-length compositions with alternating parts equal.
A342528 counts compositions with alternating parts weakly decreasing.
A342532 counts even-length compositions with alternating parts unequal.

Programs

  • Mathematica
    max = 50; s = (1+Sum[2*(-1)^k*q^(k(k+1)/2), {k, 1, max}])/QPochhammer[q]^2+ O[q]^max; CoefficientList[s, q] (* Jean-François Alcover, Nov 30 2015, from 1st g.f. *)
    wdw[q_]:=And@@Table[q[[i]]>=q[[i+2]],{i,Length[q]-2}];
    Table[Length[Select[Join@@Permutations/@Select[IntegerPartitions[n],EvenQ[Length[#]]&],wdw]],{n,0,15}] (* Gus Wiseman, Mar 25 2021 *)
  • PARI
    {a(n) = if( n<0, 0, polcoeff( sum(k=0, n\2, x^(2*k) / prod(i=1, k, 1 - x^i, 1 + x * O(x^n))^2), n))};
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( sum(k=1, sqrtint(8*n + 1)\2, 2*(-1)^k * x^((k^2+k)/2), 1 + A) / eta(x + A)^2, n))};

Formula

G.f.: 1 + Sum_{k>0} (x^k / ((1-x)(1-x^2)...(1-x^k)))^2 = (1 + Sum_{k>0} 2 (-1)^k x^((k^2+k)/2) ) / (Product_{k>0} (1 - x^k))^2.
G.f.: 1 + x*(1 - G(0))/(1-x) where G(k) = 1 - x/(1-x^(k+1))^2/(1-x/(x-1/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Jan 23 2013
a(n) = A006330(n) - A001523(n). - Vaclav Kotesovec, Jun 22 2015
a(n) ~ Pi * exp(2*Pi*sqrt(n/3)) / (16 * 3^(5/4) * n^(7/4)). - Vaclav Kotesovec, Oct 24 2018

Extensions

New name from Joerg Arndt, Jun 10 2013

A174713 Triangle read by rows, A173305 (A000009 shifted down twice) * A174712 (diagonalized variant of A000041).

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 2, 1, 2, 3, 2, 2, 4, 2, 2, 3, 5, 3, 4, 3, 6, 4, 4, 3, 5, 8, 5, 6, 6, 5, 10, 6, 8, 6, 5, 7, 12, 8, 10, 9, 10, 7, 15, 10, 12, 12, 10, 7, 11, 18, 12, 16, 15, 15, 14, 11, 22, 15, 20, 18, 20, 14, 11, 15
Offset: 0

Views

Author

Gary W. Adamson, Mar 27 2010

Keywords

Comments

Row sums = A000041, the partition numbers.
The current triangle is the 2nd in an infinite set, followed by A174714 (k=3), and A174715, (k=4); in which row sums of each triangle = A000041.
k-th triangle in the infinite set can be defined as having the sequence:
"Euler transform of ones: (1,1,1,...) interleaved with (k-1) zeros"; shifted down k times (except column 0) in successive columns, then multiplied * triangle A174712, the diagonalized variant of A000041, A174713 begins with A000009 shifted down twice (triangle A173305); where A000009 = the Euler transform of period 2 sequence: [1,0,1,0,...].
Similarly, triangle A174714 begins with A000716 shifted down thrice; where A000716 = the Euler transform of period 3 series: [1,1,0,1,1,0,...]. Then multiply the latter as an infinite lower triangular matrix * A174712, the diagonalized variant of A000041, obtaining triangle A174714 with row sums = A000041.
Case k=4 = triangle A174715 which begins with the Euler transform of period 4 series: [1,1,1,0,1,1,1,0,...], shifted down 4 times in successive columns then multiplied * A174712, the diagonalized variant of A000041.
All triangles in the infinite set have row sums = A000041.
The sequences: "Euler transform of ones interleaved with (k-1) zeros" have the following properties, beginning with k=2:
...
k=2, A000009: = Euler transform of [1,0,1,0,1,0,...] and satisfies
.....A000009. = p(x)/p(x^2), where p(x) = polcoeff A000041; and A000041 =
.....A000009(x) = r(x), then p(x) = r(x) * r(x^2) * r(x^4) * r(x^8) * ...
...
k=3, A000726: = Euler transform of [1,1,0,1,1,0,...] and satisfies
.....A000726(x): = p(x)/p(x^3), and given s(x) = polcoeff A000726, we get
.....A000041(x) = p(x) = s(x) * s(x^3) * s(x^9) * s(x^27) * ...
...
k=4, A001935: = Euler transform of [1,1,1,0,1,1,1,0,...] and satisfies
.....A001935(x) = p(x)/p(x^4) and given t(x) = polcoeff A001935, we get
.....A000041(x) = p(x) = t(x) * t(x^4) * t(x^16) * t(x^64) * ...
...
Also the number of integer partitions of n whose even parts sum to k, for k an even number from zero to n. The version including odd k is A113686. - Gus Wiseman, Oct 23 2023

Examples

			First few rows of the triangle =
1;
1;
1, 1;
2, 1;
2, 1, 2;
3, 2, 2;
4, 2, 2, 3;
5, 3, 4, 3;
6, 4, 4, 3, 5;
8, 5, 6, 6, 5;
10, 6, 8, 6, 5, 7;
12, 8, 10, 9, 10, 7;
15, 10, 12, 12, 10, 7, 11;
18, 12, 16, 15, 15, 14, 11;
22, 15, 20, 18, 20, 14, 11, 15;
...
From _Gus Wiseman_, Oct 23 2023: (Start)
Row n = 9 counts the following partitions:
  (9)          (72)        (54)       (63)      (81)
  (711)        (5211)      (522)      (6111)    (621)
  (531)        (3321)      (4311)     (432)     (441)
  (51111)      (321111)    (411111)   (42111)   (4221)
  (333)        (21111111)  (32211)    (3222)    (22221)
  (33111)                  (2211111)  (222111)
  (3111111)
  (111111111)
(End)
		

Crossrefs

Row sums are A000041.
The odd version is A365067.
The corresponding rank statistic is A366531, odd version A366528.
A000009 counts partitions into odd parts, ranks A066208.
A113685 counts partitions by sum of odd parts, even version A113686.
A239261 counts partitions with (sum of odd parts) = (sum of even parts).

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],Total[Select[#,EvenQ]]==k&]],{n,0,15},{k,0,n,2}] (* Gus Wiseman, Oct 23 2023 *)

Formula

As infinite lower triangular matrices, A173305 * A174712.
T(n,k) = A000009(n-2k) * A000041(k). - Gus Wiseman, Oct 23 2023

A261775 Expansion of Product_{k>=1} (1 - x^(8*k))/(1 - x^k).

Original entry on oeis.org

1, 1, 2, 3, 5, 7, 11, 15, 21, 29, 40, 53, 72, 94, 124, 161, 208, 266, 341, 431, 545, 684, 856, 1064, 1322, 1631, 2009, 2464, 3014, 3672, 4467, 5411, 6543, 7888, 9489, 11383, 13632, 16280, 19409, 23088, 27415, 32483, 38430, 45371, 53485, 62939, 73950, 86742
Offset: 0

Views

Author

Vaclav Kotesovec, Aug 31 2015

Keywords

Comments

Number of partitions in which no part occurs more than 7 times. - Ilya Gutkovskiy, May 31 2017

Crossrefs

Number of r-regular partitions for r = 2 through 12: A000009, A000726, A001935, A035959, A219601, A035985, A261775, A104502, A261776, A328545, A328546.

Programs

  • Maple
    a:= proc(n) option remember; `if`(n=0, 1, add(a(n-j)*add(d*
         signum(irem(d, 8)), d=numtheory[divisors](j)), j=1..n)/n)
        end:
    seq(a(n), n=0..50);  # Alois P. Heinz, Aug 07 2022
  • Mathematica
    nmax = 50; CoefficientList[Series[Product[(1 - x^(8*k))/(1 - x^k), {k, 1, nmax}], {x, 0, nmax}], x]
    Table[Count[IntegerPartitions@n, x_ /; ! MemberQ [Mod[x, 8], 0, 2] ], {n, 0, 47}] (* Robert Price, Jul 28 2020 *)
  • PARI
    Vec(prod(k=1, 51, (1 - x^(8*k))/(1 - x^k)) + O(x^51)) \\ Indranil Ghosh, Mar 25 2017

Formula

a(n) ~ Pi*sqrt(7) * BesselI(1, sqrt(7*(24*n + 7)/8) * Pi/6) / (4*sqrt(24*n + 7)) ~ exp(Pi*sqrt(7*n/3)/2) * 7^(1/4) / (2^(7/2) * 3^(1/4) * n^(3/4)) * (1 + (7^(3/2)*Pi/(96*sqrt(3)) - 3*sqrt(3)/(4*Pi*sqrt(7))) / sqrt(n) + (343*Pi^2/55296 - 45/(224*Pi^2) - 35/128) / n). - Vaclav Kotesovec, Aug 31 2015, extended Jan 14 2017
a(n) = (1/n)*Sum_{k=1..n} A284341(k)*a(n-k), a(0) = 1. - Seiichi Manyama, Mar 25 2017
G.f.: A(x)*A(x^2)*A(x^4) where A(x) is the o.g.f. for A000009. (see Flajolet, Sedgewick link) - Geoffrey Critzer, Aug 07 2022

A232432 Number of compositions of n avoiding the pattern 111.

Original entry on oeis.org

1, 1, 2, 3, 7, 11, 21, 34, 59, 114, 178, 284, 522, 823, 1352, 2133, 3739, 5807, 9063, 14074, 23639, 36006, 56914, 87296, 131142, 214933, 324644, 487659, 739291, 1108457, 1724673, 2558386, 3879335, 5772348, 8471344, 12413666, 19109304, 27886339, 40816496
Offset: 0

Views

Author

Alois P. Heinz, Nov 23 2013

Keywords

Comments

Number of compositions of n into parts with multiplicity <= 2.

Examples

			a(4) = 7: [4], [3,1], [2,2], [1,3], [2,1,1], [1,2,1], [1,1,2].
a(5) = 11: [5], [4,1], [3,2], [2,3], [1,4], [3,1,1], [2,2,1], [1,3,1], [2,1,2], [1,2,2], [1,1,3].
a(6) = 21: [6], [4,2], [3,3], [5,1], [2,4], [1,5], [2,1,3], [1,2,3], [1,1,4], [4,1,1], [3,2,1], [2,3,1], [1,4,1], [3,1,2], [1,3,2], [1,2,2,1], [2,1,1,2], [1,2,1,2], [1,1,2,2], [2,2,1,1], [2,1,2,1].
		

Crossrefs

Cf. A000726 (partitions avoiding 111), A032020 (pattern 11), A128695 (adjacent pattern 111).
Column k=2 of A243081.
The case of partitions is ranked by A004709.
The version for patterns is A080599.
(1,1,1,1)-avoiding partitions are counted by A232464.
The (1,1,1)-matching version is A335455.
Patterns matched by compositions are counted by A335456.
The version for prime indices is A335511.
(1,1,1)-avoiding compositions are ranked by A335513.

Programs

  • Maple
    b:= proc(n, i, p) option remember; `if`(n=0, p!, `if`(i<1, 0,
          add(b(n-i*j, i-1, p+j)/j!, j=0..min(n/i, 2))))
        end:
    a:= n-> b(n$2, 0):
    seq(a(n), n=0..50);
  • Mathematica
    f[list_]:=Apply[And,Table[Count[list,i]<3,{i,1,Max[list]}]];
    g[list_]:=Length[list]!/Apply[Times,Table[Count[list,i]!,{i,1,Max[list]}]];
    a[n_] := If[n == 0, 1, Total[Map[g, Select[IntegerPartitions[n], f]]]];
    Table[a[n], {n, 0, 35}] (* Geoffrey Critzer, Nov 25 2013, updated by Jean-François Alcover, Nov 20 2023 *)

A137569 Expansion of f(-x) / f(-x^3) in powers of x where f() is a Ramanujan theta function.

Original entry on oeis.org

1, -1, -1, 1, -1, 0, 2, -1, -1, 3, -2, -1, 4, -3, -2, 5, -4, -2, 8, -6, -4, 10, -7, -4, 14, -10, -6, 18, -13, -7, 24, -17, -10, 30, -21, -12, 40, -28, -17, 49, -35, -19, 63, -44, -26, 78, -55, -31, 98, -69, -40, 120, -84, -47, 150, -105, -61, 182, -127, -71
Offset: 0

Views

Author

Michael Somos, Jan 26 2008

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).

Examples

			G.f. = 1 - x - x^2 + x^3 - x^4 + 2*x^6 - x^7 - x^8 + 3*x^9 - 2*x^10 - x^11 + ...
G.f. = 1/q - q^11 - q^23 + q^35 - q^47 + 2*q^71 - q^83 - q^95 + 3*q^107 + ...
		

Crossrefs

Programs

  • Mathematica
    a[ n_] := SeriesCoefficient[ QPochhammer[ x] / QPochhammer[ x^3], {x, 0, n}]; (* Michael Somos, Oct 08 2015 *)
    a[ n_] := SeriesCoefficient[ QPochhammer[ x, x^3] QPochhammer[ x^2, x^3], {x, 0, n}]; (* Michael Somos, Oct 08 2015 *)
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x + A) / eta(x^3 + A), n))};

Formula

Expansion of q^(1/12) eta(q) / eta(q^3) in powers of q.
Euler transform of period 3 sequence [ -1, -1, 0, ...].
Given g.f. A(x) then B(q) = A(q^6)^2 / q satisfies 0 = f(B(q), B(q^2), B(q^4)) where f(u, v, w) = 4*v^2 + (u^2 - v) * (w^2 + v).
G.f. is a period 1 Fourier series which satisfies f(-1 / (432 t)) = 3^(1/2) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A000726.
G.f.: Product_{k>0} (1 - x^(3*k-1)) * (1 - x^(3*k-2)).
a(3*n) = A035943(n). a(3*n + 1) = - A035941(n). a(3*n + 2) = - A035940(n).
Convolution inverse of A000726.
Convolution square is A112157. Convolution 4th power is A058095. - Michael Somos, Oct 08 2015
a(2*n) = A263050(n). a(2*n + 1) = - A263051(n). - Michael Somos, Oct 08 2015
G.f.: (Product_{k>0} (1 + x^k + x^(2*k)))^-1. - Michael Somos, Oct 08 2015
a(n) = -(1/n)*Sum_{k=1..n} A046913(k)*a(n-k), a(0) = 1. - Seiichi Manyama, Mar 25 2017

A224958 Number of compositions [p(1), p(2), ..., p(k)] of n such that p(j) != p(j-2).

Original entry on oeis.org

1, 1, 2, 3, 6, 9, 18, 29, 53, 91, 162, 277, 495, 855, 1508, 2625, 4618, 8049, 14130, 24675, 43255, 75621, 132475, 231697, 405751, 709887, 1242824, 2174763, 3806989, 6662291, 11661737, 20409409, 35723307, 62521919, 109431810, 191527623, 335225350, 586717615
Offset: 0

Views

Author

Joerg Arndt, Apr 21 2013

Keywords

Examples

			The a(6) = 18 such compositions of 6 are
01:  [ 1 1 2 2 ]
02:  [ 1 1 4 ]
03:  [ 1 2 2 1 ]
04:  [ 1 2 3 ]
05:  [ 1 3 2 ]
06:  [ 1 5 ]
07:  [ 2 1 1 2 ]
08:  [ 2 1 3 ]
09:  [ 2 2 1 1 ]
10:  [ 2 3 1 ]
11:  [ 2 4 ]
12:  [ 3 1 2 ]
13:  [ 3 2 1 ]
14:  [ 3 3 ]
15:  [ 4 1 1 ]
16:  [ 4 2 ]
17:  [ 5 1 ]
18:  [ 6 ]
		

Crossrefs

Cf. A000726 (partitions such that p(j) != p(j-2)), A003242, A241902.

Programs

  • Maple
    b:= proc(n, i, j) option remember; `if`(n=0, 1, add(`if`(k=j, 0,
          b(n-k, `if`(n-k b(n, 0, 0):
    seq(a(n), n=0..50);  # Alois P. Heinz, May 02 2013
  • Mathematica
    b[n_, i_, j_] := b[n, i, j] = If[n==0, 1, Sum[If[k==j, 0, b[n-k, If[n-k < k, 0, k], If[n-k < i, 0, i]]], {k, 1, n}]]; a[n_] := b[n, 0, 0]; Table[a[n], {n, 0, 50}] (* Jean-François Alcover, Apr 08 2015, after Alois P. Heinz *)

Formula

a(n) ~ c * d^n, where d = 1.7502412917183090312497386246... (see A241902) and c = 0.5940298439978189763822100914... - Vaclav Kotesovec, May 01 2014

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

Original entry on oeis.org

1, 1, 0, 2, 0, 1, 3, 1, 1, 2, 6, 1, 4, 2, 5, 10, 5, 4, 9, 7, 8, 21, 9, 13, 13, 19, 13, 27, 32, 23, 29, 33, 27, 45, 37, 45, 79, 49, 57, 68, 82, 67, 101, 83, 109, 155, 124, 113, 174, 148, 171, 196, 215, 198, 262, 310, 269, 330, 314, 342, 414, 430, 393, 536, 493
Offset: 0

Views

Author

Vaclav Kotesovec, Jan 03 2016

Keywords

Crossrefs

Programs

  • Mathematica
    nmax = 80; CoefficientList[Series[Product[1+x^k-x^(2*k), {k, 1, nmax}], {x, 0, nmax}], x]
    nmax = 100; p = ConstantArray[0, nmax + 1]; p[[1]] = 1; p[[2]] = 1; p[[3]] = -1; Do[Do[p[[j+1]] = p[[j+1]] + p[[j - k + 1]] - If[j < 2*k, 0, p[[j - 2*k + 1]]], {j, nmax, k, -1}];, {k, 2, nmax}]; p (* Vaclav Kotesovec, May 10 2018 *)

Formula

a(n) ~ sqrt(log(phi)) * phi^sqrt(8*n) / (2^(3/4)*sqrt(Pi)*n^(3/4)), where phi = A001622 = (1+sqrt(5))/2 is the golden ratio. - Vaclav Kotesovec, Jan 03 2016

A328545 Number of 11-regular partitions of n (no part is a multiple of 11).

Original entry on oeis.org

1, 1, 2, 3, 5, 7, 11, 15, 22, 30, 42, 55, 76, 99, 132, 171, 224, 286, 370, 468, 597, 750, 945, 1177, 1472, 1820, 2255, 2772, 3410, 4165, 5092, 6185, 7515, 9085, 10978, 13207, 15884, 19025, 22774, 27170, 32388, 38489, 45705, 54120, 64030, 75569, 89100
Offset: 0

Views

Author

N. J. A. Sloane, Oct 19 2019

Keywords

References

  • Kathiravan, T., and S. N. Fathima. "On L-regular bipartitions modulo L." The Ramanujan Journal 44.3 (2017): 549-558.

Crossrefs

Number of r-regular partitions for r = 2 through 12: A000009, A000726, A001935, A035959, A219601, A035985, A261775, A104502, A261776, A328545, A328546.

Programs

  • Maple
    f:=(k,M) -> mul(1-q^(k*j),j=1..M);
    LRP := (L,M) -> f(L,M)/f(1,M);
    s := L -> seriestolist(series(LRP(L,80),q,60));
    s(11);
  • Mathematica
    Table[Count[IntegerPartitions@n, x_ /; ! MemberQ [Mod[x, 11], 0, 2] ], {n, 0, 46}]

Formula

a(n) ~ exp(Pi*sqrt(2*n*(s-1)/(3*s))) * (s-1)^(1/4) / (2 * 6^(1/4) * s^(3/4) * n^(3/4)) * (1 + ((s-1)^(3/2)*Pi/(24*sqrt(6*s)) - 3*sqrt(6*s) / (16*Pi * sqrt(s-1))) / sqrt(n) + ((s-1)^3*Pi^2/(6912*s) - 45*s/(256*(s-1)*Pi^2) - 5*(s-1)/128) / n), set s=11. - Vaclav Kotesovec, Aug 01 2022

A328546 Number of 12-regular partitions of n (no part is a multiple of 12).

Original entry on oeis.org

1, 1, 2, 3, 5, 7, 11, 15, 22, 30, 42, 56, 76, 100, 133, 173, 226, 290, 374, 475, 605, 762, 960, 1199, 1497, 1856, 2299, 2831, 3482, 4261, 5208, 6337, 7700, 9321, 11266, 13572, 16325, 19578, 23444, 27999, 33389, 39721, 47185, 55929, 66199, 78199, 92246
Offset: 0

Views

Author

N. J. A. Sloane, Oct 19 2019

Keywords

References

  • Kathiravan, T., and S. N. Fathima. "On L-regular bipartitions modulo L." The Ramanujan Journal 44.3 (2017): 549-558.

Crossrefs

Number of r-regular partitions for r = 2 through 12: A000009, A000726, A001935, A035959, A219601, A035985, A261775, A104502, A261776, A328545, A328546.

Programs

  • Maple
    f:=(k,M) -> mul(1-q^(k*j),j=1..M);
    LRP := (L,M) -> f(L,M)/f(1,M);
    s := L -> seriestolist(series(LRP(L,80),q,60));
    s(12);
  • Mathematica
    Table[Count[IntegerPartitions@n, x_ /; ! MemberQ [Mod[x, 12], 0, 2] ], {n, 0, 46}] (* Robert Price, Jul 28 2020 *)

Formula

a(n) ~ exp(Pi*sqrt(2*n*(s-1)/(3*s))) * (s-1)^(1/4) / (2 * 6^(1/4) * s^(3/4) * n^(3/4)) * (1 + ((s-1)^(3/2)*Pi/(24*sqrt(6*s)) - 3*sqrt(6*s) / (16*Pi * sqrt(s-1))) / sqrt(n) + ((s-1)^3*Pi^2/(6912*s) - 45*s/(256*(s-1)*Pi^2) - 5*(s-1)/128) / n), set s=12. - Vaclav Kotesovec, Aug 01 2022

A209318 Number T(n,k) of partitions of n with k parts in which no part occurs more than twice; triangle T(n,k), n>=0, 0<=k<=A055086(n), read by rows.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 1, 1, 0, 1, 2, 1, 0, 1, 2, 2, 0, 1, 3, 2, 1, 0, 1, 3, 4, 1, 0, 1, 4, 5, 3, 0, 1, 4, 6, 4, 1, 0, 1, 5, 8, 6, 2, 0, 1, 5, 10, 8, 3, 0, 1, 6, 11, 12, 5, 1, 0, 1, 6, 14, 14, 8, 1, 0, 1, 7, 16, 19, 11, 3, 0, 1, 7, 18, 23, 16, 5
Offset: 0

Views

Author

Alois P. Heinz, Jan 19 2013

Keywords

Examples

			T(8,3) = 5: [6,1,1], [5,2,1], [4,3,1], [4,2,2], [3,3,2].
T(8,4) = 3: [4,2,1,1], [3,3,1,1], [3,2,2,1].
T(9,3) = 6: [7,1,1], [6,2,1], [5,3,1], [4,4,1], [5,2,2], [4,3,2].
T(9,4) = 4: [5,2,1,1], [4,3,1,1], [4,2,2,1], [3,3,2,1].
T(9,5) = 1: [3,2,2,1,1].
Triangle begins:
  1;
  0, 1;
  0, 1, 1;
  0, 1, 1;
  0, 1, 2, 1;
  0, 1, 2, 2;
  0, 1, 3, 2, 1;
  0, 1, 3, 4, 1;
  0, 1, 4, 5, 3;
  0, 1, 4, 6, 4, 1;
		

Crossrefs

Columns k=0-10 give: A000007, A057427, A004526, A230059 (conjectured), A320592, A320593, A320594, A320595, A320596, A320597, A320598.
Row sums give: A000726.
Row lengths give: A000267.
Cf. A002620, A008289 (no part more than once), A055086, A117147 (no part more than 3 times).

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
          add(expand(b(n-i*j, i-1)*x^j), j=0..min(2, n/i))))
        end:
    T:= n-> (p-> seq(coeff(p, x, i), i=0..degree(p)))(b(n$2)):
    seq(T(n), n=0..20);
  • Mathematica
    max = 15; g = -1+Product[1+t*x^j+t^2*x^(2j), {j, 1, max}]; t[n_, k_] := SeriesCoefficient[g, {x, 0, n}, {t, 0, k}]; t[0, 0] = 1; Table[Table[t[n, k], {k, 0, n}] /. {a__, 0 ..} -> {a}, {n, 0, max}] // Flatten (* Jean-François Alcover, Jan 08 2014 *)
Previous Showing 21-30 of 87 results. Next