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 3899 results. Next

A218482 First differences of the binomial transform of the partition numbers (A000041).

Original entry on oeis.org

1, 1, 3, 8, 21, 54, 137, 344, 856, 2113, 5179, 12614, 30548, 73595, 176455, 421215, 1001388, 2371678, 5597245, 13166069, 30873728, 72185937, 168313391, 391428622, 908058205, 2101629502, 4853215947, 11183551059, 25718677187, 59030344851, 135237134812, 309274516740
Offset: 0

Views

Author

Paul D. Hanna, Oct 29 2012

Keywords

Comments

a(n) = A103446(n) for n>=1; here a(0) is set to 1 in accordance with the definition and other important generating functions.
From Gus Wiseman, Dec 12 2022: (Start)
Also the number of sequences of compositions (A133494) with weakly decreasing lengths and total sum n. For example, the a(0) = 1 through a(3) = 8 sequences are:
() ((1)) ((2)) ((3))
((11)) ((12))
((1)(1)) ((21))
((111))
((1)(2))
((2)(1))
((11)(1))
((1)(1)(1))
The case of constant lengths is A101509.
The case of strictly decreasing lengths is A129519.
The case of sequences of partitions is A141199.
The case of twice-partitions is A358831.
(End)

Examples

			G.f.: A(x) = 1 + x + 3*x^2 + 8*x^3 + 21*x^4 + 54*x^5 + 137*x^6 + 344*x^7 +...
The g.f. equals the product:
A(x) = (1-x)/((1-x)-x) * (1-x)^2/((1-x)^2-x^2) * (1-x)^3/((1-x)^3-x^3) * (1-x)^4/((1-x)^4-x^4) * (1-x)^5/((1-x)^5-x^5) * (1-x)^6/((1-x)^6-x^6) * (1-x)^7/((1-x)^7-x^7) *...
and also equals the series:
A(x) = 1  +  x*(1-x)/((1-x)-x)^2  +  x^4*(1-x)^2/(((1-x)-x)*((1-x)^2-x^2))^2  +  x^9*(1-x)^3/(((1-x)-x)*((1-x)^2-x^2)*((1-x)^3-x^3))^2  +  x^16*(1-x)^4/(((1-x)-x)*((1-x)^2-x^2)*((1-x)^3-x^3)*((1-x)^4-x^4))^2 +...
		

Crossrefs

Programs

  • Maple
    b:= proc(n) option remember;
          add(combinat[numbpart](k)*binomial(n,k), k=0..n)
        end:
    a:= n-> b(n)-b(n-1):
    seq(a(n), n=0..50);  # Alois P. Heinz, Aug 19 2014
  • Mathematica
    Flatten[{1, Table[Sum[Binomial[n-1,k]*PartitionsP[k+1],{k,0,n-1}],{n,1,30}]}] (* Vaclav Kotesovec, Jun 25 2015 *)
  • PARI
    {a(n)=sum(k=0,n,(binomial(n,k)-if(n>0,binomial(n-1,k)))*numbpart(k))}
    for(n=0,40,print1(a(n),", "))
    
  • PARI
    {a(n)=local(X=x+x*O(x^n));polcoeff(prod(k=1,n,(1-x)^k/((1-x)^k-X^k)),n)}
    
  • PARI
    {a(n)=local(X=x+x*O(x^n));polcoeff(sum(m=0,n,x^m*(1-x)^(m*(m-1)/2)/prod(k=1,m,((1-x)^k - X^k))),n)}
    
  • PARI
    {a(n)=local(X=x+x*O(x^n));polcoeff(sum(m=0,n,x^(m^2)*(1-X)^m/prod(k=1,m,((1-x)^k - x^k)^2)),n)}
    
  • PARI
    {a(n)=local(X=x+x*O(x^n));polcoeff(exp(sum(m=1,n+1,x^m/((1-x)^m-X^m)/m)),n)}
    
  • PARI
    {a(n)=local(X=x+x*O(x^n));polcoeff(exp(sum(m=1,n+1,sigma(m)*x^m/(1-X)^m/m)),n)}
    
  • PARI
    {a(n)=local(X=x+x*O(x^n));polcoeff(prod(k=1,n,(1 + x^k/(1-X)^k)^valuation(2*k,2)),n)}

Formula

G.f.: Product_{n>=1} (1-x)^n / ((1-x)^n - x^n).
G.f.: Sum_{n>=0} x^n * (1-x)^(n*(n-1)/2) / Product_{k=1..n} ((1-x)^k - x^k).
G.f.: Sum_{n>=0} x^(n^2) * (1-x)^n / Product_{k=1..n} ((1-x)^k - x^k)^2.
G.f.: exp( Sum_{n>=1} x^n/((1-x)^n - x^n) / n ).
G.f.: exp( Sum_{n>=1} sigma(n) * x^n/(1-x)^n / n ), where sigma(n) is the sum of divisors of n (A000203).
G.f.: Product_{n>=1} (1 + x^n/(1-x)^n)^A001511(n), where 2^A001511(n) is the highest power of 2 that divides 2*n.
a(n) ~ exp(Pi*sqrt(n/3) + Pi^2/24) * 2^(n-2) / (n*sqrt(3)). - Vaclav Kotesovec, Jun 25 2015

A056823 Number of compositions minus number of partitions: A011782(n) - A000041(n).

Original entry on oeis.org

0, 0, 0, 1, 3, 9, 21, 49, 106, 226, 470, 968, 1971, 3995, 8057, 16208, 32537, 65239, 130687, 261654, 523661, 1047784, 2096150, 4193049, 8387033, 16775258, 33551996, 67105854, 134214010, 268430891, 536865308, 1073734982, 2147475299, 4294957153, 8589922282
Offset: 0

Views

Author

Alford Arnold, Aug 29 2000

Keywords

Comments

Previous name was: Counts members of A056808 by number of factors.
A056808 relates to least prime signatures (cf. A025487)
a(n) is also the number of compositions of n that are not partitions of n. - Omar E. Pol, Jan 31 2009, Oct 14 2013
a(n) is the number of compositions of n into positive parts containing pattern [1,2]. - Bob Selcoe, Jul 08 2014

Examples

			A011782 begins     1 1 2 4 8 16 32 64 128 256 ...;
A000041 begins     1 1 2 3 5  7 11 15  22  30 ...;
so sequence begins 0 0 0 1 3  9 21 49 106 226 ... .
For n = 3 the factorizations are 8=2*2*2, 12=2*2*3, 18=2*3*3 and 30=2*3*5.
a(5) = 9: {[1,1,1,2], [1,1,2,1], [1,1,3], [1,2,1,1], [1,2,2], [1,3,1], [1,4], [2,1,2], [2,3]}. - _Bob Selcoe_, Jul 08 2014
		

Crossrefs

The version for patterns is A002051.
(1,2)-avoiding compositions are just partitions A000041.
The (1,1)-matching version is A261982.
The version for prime indices is A335447.
(1,2)-matching compositions are ranked by A335485.
Patterns matched by compositions are counted by A335456.

Programs

  • Maple
    a:= n-> ceil(2^(n-1))-combinat[numbpart](n):
    seq(a(n), n=0..37);  # Alois P. Heinz, Jan 30 2020
  • Mathematica
    Table[Length[Select[Join@@Permutations/@IntegerPartitions[n],!GreaterEqual@@#&]],{n,0,10}] (* Gus Wiseman, Jun 24 2020 *)
    a[n_] := If[n == 0, 0, 2^(n-1) - PartitionsP[n]];
    a /@ Range[0, 37] (* Jean-François Alcover, May 23 2021 *)

Formula

a(n) = A011782(n) - A000041(n).
a(n) = 2*a(n-1) + A117989(n-1). - Bob Selcoe, Apr 11 2014
G.f.: (1 - x) / (1 - 2*x) - Product_{k>=1} 1 / (1 - x^k). - Ilya Gutkovskiy, Jan 30 2020

Extensions

More terms from James Sellers, Aug 31 2000
New name from Joerg Arndt, Sep 02 2013

A340035 Irregular triangle read by rows T(n,k) in which row n lists n blocks, where the m-th block consists of A000041(n-m) copies of the divisors of m, with 1 <= m <= n.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 1, 2, 1, 3, 1, 1, 1, 1, 2, 1, 2, 1, 3, 1, 2, 4, 1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 2, 1, 3, 1, 3, 1, 2, 4, 1, 5, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 3, 1, 3, 1, 3, 1, 2, 4, 1, 2, 4, 1, 5, 1, 2, 3, 6, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 2, 1, 2
Offset: 1

Views

Author

Omar E. Pol, Dec 26 2020

Keywords

Comments

For further information about the correspondence divisor/part see A338156.

Examples

			Triangle begins:
  1;
  1, 1, 2;
  1, 1, 1, 2, 1, 3;
  1, 1, 1, 1, 2, 1, 2, 1, 3, 1, 2, 4;
  1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 2, 1, 3, 1, 3, 1, 2, 4, 1, 5;
  ...
Written as an irregular tetrahedron the first five slices are:
  1;
  --
  1,
  1, 2;
  -----
  1,
  1,
  1, 2
  1, 3;
  -----
  1,
  1,
  1,
  1, 2,
  1, 2,
  1, 3,
  1, 2, 4;
  --------
  1,
  1,
  1,
  1,
  1,
  1, 2,
  1, 2,
  1, 2,
  1, 3,
  1, 3,
  1, 2, 4,
  1, 5;
--------
The slices of the tetrahedron appear in the upper zone of the following table (formed by three zones) which shows the correspondence between divisors and parts (n = 1..5):
.
|---|---------|-----|-------|---------|-----------|-------------|
| n |         |  1  |   2   |    3    |     4     |      5      |
|---|---------|-----|-------|---------|-----------|-------------|
|   | A027750 |     |       |         |           |  1          |
|   | A027750 |     |       |         |           |  1          |
|   | A027750 |     |       |         |           |  1          |
|   | A027750 |     |       |         |           |  1          |
| D | A027750 |     |       |         |           |  1          |
| I |---------|-----|-------|---------|-----------|-------------|
| V | A027750 |     |       |         |  1        |  1 2        |
| I | A027750 |     |       |         |  1        |  1 2        |
| S | A027750 |     |       |         |  1        |  1 2        |
| O |---------|-----|-------|---------|-----------|-------------|
| R | A027750 |     |       |  1      |  1 2      |  1   3      |
| S | A027750 |     |       |  1      |  1 2      |  1   3      |
|   |---------|-----|-------|---------|-----------|-------------|
|   | A027750 |     |  1    |  1 2    |  1   3    |  1 2   4    |
|   |---------|-----|-------|---------|-----------|-------------|
|   | A027750 |  1  |  1 2  |  1   3  |  1 2   4  |  1       5  |
|---|---------|-----|-------|---------|-----------|-------------|
.
|---|---------|-----|-------|---------|-----------|-------------|
|   | A138785 |  1  |  2 2  |  4 2 3  |  7 6 3 4  | 12 8 6 4 5  |
|   |         |  =  |  = =  |  = = =  |  = = = =  |  = = = = =  |
| L | A002260 |  1  |  1 2  |  1 2 3  |  1 2 3 4  |  1 2 3 4 5  |
| I |         |  *  |  * *  |  * * *  |  * * * *  |  * * * * *  |
| N | A066633 |  1  |  2 1  |  4 1 1  |  7 3 1 1  | 12 4 2 1 1  |
| K |         |  |  |  |\|  |  |\|\|  |  |\|\|\|  |  |\|\|\|\|  |
|   | A181187 |  1  |  3 1  |  6 2 1  | 12 5 2 1  | 20 8 4 2 1  |
|---|---------|-----|-------|---------|-----------|-------------|
.
|---|---------|-----|-------|---------|-----------|-------------|
| P |         |  1  |  1 1  |  1 1 1  |  1 1 1 1  |  1 1 1 1 1  |
| A |         |     |  2    |  2 1    |  2 1 1    |  2 1 1 1    |
| R |         |     |       |  3      |  3 1      |  3 1 1      |
| T |         |     |       |         |  2 2      |  2 2 1      |
| I |         |     |       |         |  4        |  4 1        |
| T |         |     |       |         |           |  3 2        |
| I |         |     |       |         |           |  5          |
| O |         |     |       |         |           |             |
| N |         |     |       |         |           |             |
| S |         |     |       |         |           |             |
|---|---------|-----|-------|---------|-----------|-------------|
.
The table is essentially the same table of A340032 but here, in the upper zone, every row is A027750 instead of A127093.
Also the above table is the table of A338156 upside down.
The connection with the tower described in A221529 is as follows (n = 7):
|--------|------------------------|
| Level  |                        |
| in the | 7th slice of divisors  |
| tower  |                        |
|--------|------------------------|
|  11    |   1,                   |
|  10    |   1,                   |
|   9    |   1,                   |
|   8    |   1,                   |
|   7    |   1,                   |
|   6    |   1,                   |
|   5    |   1,                   |
|   4    |   1,                   |
|   3    |   1,                   |
|   2    |   1,                   |
|   1    |   1,                   |
|--------|------------------------|
|   7    |   1, 2,                |
|   6    |   1, 2,                |
|   5    |   1, 2,                |
|   4    |   1, 2,                |
|   3    |   1, 2,                |
|   2    |   1, 2,                |
|   1    |   1, 2,                |
|--------|------------------------|
|   5    |   1,    3,             |
|   4    |   1,    3,             |
|   3    |   1,    3,             |
|   2    |   1,    3,             |      Level
|   1    |   1,    3,             |             _
|--------|------------------------|       11   | |
|   3    |   1, 2,    4,          |       10   | |
|   2    |   1, 2,    4,          |        9   | |
|   1    |   1, 2,    4,          |        8   |_|_
|--------|------------------------|        7   |   |
|   2    |   1,          5,       |        6   |_ _|_
|   1    |   1,          5,       |        5   |   | |
|--------|------------------------|        4   |_ _|_|_
|   1    |   1, 2, 3,       6,    |        3   |_ _ _| |_
|--------|------------------------|        2   |_ _ _|_ _|_ _
|   1    |   1,                7; |        1   |_ _ _ _|_|_ _|
|--------|------------------------|
             Figure 1.                            Figure 2.
                                                Lateral view
                                                of the tower.
.
                                                _ _ _ _ _ _ _
                                               |_| | | | |   |
                                               |_ _|_| | |   |
                                               |_ _|  _|_|   |
                                               |_ _ _|    _ _|
                                               |_ _ _|  _|
                                               |       |
                                               |_ _ _ _|
.
                                                  Figure 3.
                                                  Top view
                                                of the tower.
.
Figure 1 shows the terms of the 7th row of the triangle arranged as the 7th slice of the tetrahedron. The left hand column (see figure 1) gives the level of the sum of the divisors in the tower (see figures 2 and 3).
		

Crossrefs

Programs

  • Mathematica
    A340035row[n_]:=Flatten[Array[ConstantArray[Divisors[#],PartitionsP[n-#]]&,n]];
    nrows=7;Array[A340035row,nrows] (* Paolo Xausa, Jun 20 2022 *)

A218481 Binomial transform of the partition numbers (A000041).

Original entry on oeis.org

1, 2, 5, 13, 34, 88, 225, 569, 1425, 3538, 8717, 21331, 51879, 125474, 301929, 723144, 1724532, 4096210, 9693455, 22859524, 53733252, 125919189, 294232580, 685661202, 1593719407, 3695348909, 8548564856, 19732115915, 45450793102, 104481137953, 239718272765
Offset: 0

Views

Author

Paul D. Hanna, Oct 29 2012

Keywords

Comments

Partial sums of A218482.
From Vaclav Kotesovec, Nov 02 2023: (Start)
Let 0 < p < 1, r > 0, v > 0, f(n) = v*exp(r*n^p)/n^b, then
Sum_{k=0..n} binomial(n,k) * f(k) ~ f(n/2) * 2^n * exp(g(n)), where
g(n) = p^2 * r^2 * n^p / (2^(1+2*p)*n^(1-p) + p*r*(1-p)*2^(1+p)).
Special cases:
p < 1/2, g(n) = 0
p = 1/2, g(n) = r^2/16
p = 2/3, g(n) = r^2 * n^(1/3) / (9 * 2^(1/3)) - r^3/81
p = 3/4, g(n) = 9*r^2*sqrt(n)/(64*sqrt(2)) - 27*r^3*n^(1/4)/(2048*2^(1/4)) + 81*r^4/65536
p = 3/5, g(n) = 9*r^2*n^(1/5)/(100*2^(1/5))
p = 4/5, g(n) = 2^(7/5)*r^2*n^(3/5)/25 - 4*2^(3/5)*r^3*n^(2/5)/625 + 8*2^(4/5)*r^4*n^(1/5)/15625 - 32*r^5/390625
(End)

Examples

			G.f.: A(x) = 1 + 2*x + 5*x^2 + 13*x^3 + 34*x^4 + 88*x^5 + 225*x^6 + 569*x^7 +...
The g.f. equals the product:
A(x) = 1/((1-x)-x) * (1-x)^2/((1-x)^2-x^2) * (1-x)^3/((1-x)^3-x^3) * (1-x)^4/((1-x)^4-x^4) * (1-x)^5/((1-x)^5-x^5) * (1-x)^6/((1-x)^6-x^6) * (1-x)^7/((1-x)^7-x^7) *...
and also equals the series:
A(x) = 1/(1-x) * (1  +  x*(1-x)/((1-x)-x)^2  +  x^4*(1-x)^2/(((1-x)-x)*((1-x)^2-x^2))^2  +  x^9*(1-x)^3/(((1-x)-x)*((1-x)^2-x^2)*((1-x)^3-x^3))^2  +  x^16*(1-x)^4/(((1-x)-x)*((1-x)^2-x^2)*((1-x)^3-x^3)*((1-x)^4-x^4))^2 +...).
The terms begin:
a(0) = 1*1,
a(1) = 1*1 + 1*1 = 2;
a(2) = 1*1 + 2*1 + 1*2 = 5;
a(3) = 1*1 + 3*1 + 3*2 + 1*3 = 13;
a(4) = 1*1 + 4*1 + 6*2 + 4*3 + 1*5 = 34; ...
		

Crossrefs

Programs

  • Mathematica
    Table[Sum[Binomial[n,k]*PartitionsP[k],{k,0,n}],{n,0,30}] (* Vaclav Kotesovec, Jun 25 2015 *)
    nmax = 30; CoefficientList[Series[Sum[PartitionsP[k] * x^k / (1-x)^(k+1), {k, 0, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Jul 31 2022 *)
  • PARI
    {a(n)=sum(k=0,n,binomial(n,k)*numbpart(k))}
    for(n=0,40,print1(a(n),", "))
    
  • PARI
    {a(n)=local(X=x+x*O(x^n));polcoeff(1/(1-X)*prod(k=1,n,(1-x)^k/((1-x)^k-X^k)),n)}
    
  • PARI
    {a(n)=local(X=x+x*O(x^n));polcoeff(1/(1-X)*sum(m=0,n,x^m*(1-x)^(m*(m-1)/2)/prod(k=1,m,((1-x)^k - X^k))),n)}
    
  • PARI
    {a(n)=local(X=x+x*O(x^n));polcoeff(1/(1-X)*sum(m=0,n,x^(m^2)*(1-X)^m/prod(k=1,m,((1-x)^k - x^k)^2)),n)}
    
  • PARI
    {a(n)=local(X=x+x*O(x^n));polcoeff(1/(1-X)*exp(sum(m=1,n+1,x^m/((1-x)^m-X^m)/m)),n)}
    
  • PARI
    {a(n)=local(X=x+x*O(x^n));polcoeff(1/(1-X)*exp(sum(m=1,n+1,sigma(m)*x^m/(1-X)^m/m)),n)}
    
  • PARI
    {a(n)=local(X=x+x*O(x^n));polcoeff(1/(1-X)*prod(k=1,n,(1 + x^k/(1-X)^k)^valuation(2*k,2)),n)}

Formula

G.f.: 1/(1-x)*Product_{n>=1} (1-x)^n / ((1-x)^n - x^n).
G.f.: 1/(1-x)*Sum_{n>=0} x^n * (1-x)^(n*(n-1)/2) / Product_{k=1..n} ((1-x)^k - x^k).
G.f.: 1/(1-x)*Sum_{n>=0} x^(n^2) * (1-x)^n / Product_{k=1..n} ((1-x)^k - x^k)^2.
G.f.: 1/(1-x)*exp( Sum_{n>=1} x^n/((1-x)^n - x^n) / n ).
G.f.: 1/(1-x)*exp( Sum_{n>=1} sigma(n) * x^n/(1-x)^n / n ), where sigma(n) is the sum of divisors of n (A000203).
G.f.: 1/(1-x)*Product_{n>=1} (1 + x^n/(1-x)^n)^A001511(n), where 2^A001511(n) is the highest power of 2 that divides 2*n.
Logarithmic derivative yields A222115.
a(n) ~ exp(Pi*sqrt(n/3) + Pi^2/24) * 2^(n-1) / (n*sqrt(3)). - Vaclav Kotesovec, Jun 25 2015

A109085 G.f. A(x) satisfies: A(x) = P(x*A(x)) where P(x) = A(x/P(x)) is the g.f. of the partition numbers A000041.

Original entry on oeis.org

1, 1, 3, 10, 38, 153, 646, 2816, 12585, 57343, 265401, 1244256, 5896512, 28200365, 135935424, 659754072, 3221354296, 15812501100, 77985955410, 386254209762, 1920391362054, 9580985321554, 47951223856445, 240680464689600
Offset: 0

Views

Author

Paul D. Hanna, Jun 18 2005

Keywords

Comments

a(n) = Sum[Product(1 + n/h(v)^2)]/(n+1), where the product is over all boxes v in the Ferrers diagram of a partition L of n, h(v) is the hook length of v and the summation is over all partitions L of n. Example: a(3)=10 because for the partitions L=(3), (2,1), (1,1,1) of n=3 the hook length multi-sets are {3,2,1}, {3,1,1},{3,2,1}, respectively, the products are (1+3/9)(1+3/4)(1+3/1)=28/3, (1+3/9)(1+3/1)(1+3/1)=64/3, (1+3/9)(1+3/4)(1+3/1)=28/3 and now a(3)=(1/4)(28+64+28)/3=10. - Emeric Deutsch, May 15 2008

Examples

			G.f.: A(x) = 1 + x + 3*x^2 + 10*x^3 + 38*x^4 + 153*x^5 + 646*x^6 + ...
G.f. satisfies: P(x*A(x)) = A(x) where P(x) is the partition function:
P(x) = 1 + x + 2*x^2 + 3*x^3 + 5*x^4 + 7*x^5 + 11*x^6 + 15*x^7 + 22*x^8 + ...
The g.f. A = A(x) also satisfies the identities:
(1) A(x) = 1/((1-x*A) * (1-x^2*A^2) * (1-x^3*A^3) * (1-x^4*A^4) * ...).
(2) A(x) = 1 + x*A/(1-x*A) + x^2*A^2/((1-x*A)*(1-x^2*A^2)) + x^3*A^3/((1-x*A)*(1-x^2*A^2)*(1-x^3*A^3)) + ...
(3) A(x) = 1 + x*A/(1-x*A)^2 + x^4*A^4/((1-x*A)*(1-x^2*A^2))^2 + x^9*A^9/((1-x*A)*(1-x^2*A^2)*(1-x^3*A^3))^2 + ...
The logarithm of the g.f. is given by:
log(A(x)) = x*A(x)/(1-x*A(x)) + x^2*A(x)^2/(2*(1-x^2*A(x)^2)) + x^3*A(x)^3/(3*(1-x^3*A(x)^3)) + x^4*A(x)^4/(4*(1-x^4*A(x)^4)) + x^5*A(x)^5/(5*(1-x^5*A(x)^5)) + ...
Explicitly,
log(A(x)) = x + 5*x^2/2 + 22*x^3/3 + 105*x^4/4 + 506*x^5/5 + 2492*x^6/6 + 12405*x^7/7 + 62337*x^8/8 + ... + A008485(n)*x^n/n + ...
		

Crossrefs

Programs

  • Mathematica
    A109085list[n_] := Module[{m = 1, A = 1 + x}, For[i = 1, i <= n, i++, A = 1/Product[(1 - x^k*(A + x*O[x]^n)^k), {k, 1, n}]]; CoefficientList[A, x][[1 ;; n]]]; A109085list[24] (* Jean-François Alcover, Apr 21 2016, adapted from PARI *)
    InverseSeries[x QPochhammer[x] + O[x]^25][[3]] (* Vladimir Reshetnikov, Nov 17 2016 *)
    Table[SeriesCoefficient[1/QPochhammer[x, x]^(n+1), {x, 0, n}]/(n+1), {n, 0, 24}] (* Vladimir Reshetnikov, Nov 20 2016 *)
  • PARI
    {a(n)=polcoeff(1/x*serreverse(x*eta(x+x*O(x^n))),n)}
    for(n=0,30,print1(a(n),", "))
    
  • PARI
    {a(n)=local(A=1+x); for(i=1, n, A=1/prod(k=1, n, (1-x^k*(A+x*O(x^n))^k))); polcoeff(A, n)}
    
  • PARI
    {a(n)=local(A=1+x); for(i=1, n, A=sum(m=0,n,x^m*A^m/prod(k=1, m, (1-x^k*(A+x*O(x^n))^k)))); polcoeff(A, n)} \\ Paul D. Hanna, Nov 24 2012
    for(n=0,30,print1(a(n),", "))
    
  • PARI
    {a(n)=local(A=1+x); for(i=1, n, A=sum(m=0,sqrtint(n+1),(x*A)^(m^2)/prod(k=1, m, (1-x^k*(A+x*O(x^n))^k)^2))); polcoeff(A, n)} \\ Paul D. Hanna, Nov 24 2012
    for(n=0,30,print1(a(n),", "))
    
  • PARI
    {a(n)=local(A=1+x);for(i=1,n,A=exp(sum(m=1,n,(x^m*A^m/m)/(1-x^m*A^m+x*O(x^n)) )));polcoeff(A,n)} \\ Paul D. Hanna, Jun 01 2011
    
  • PARI
    {A008485(n)=polcoeff(prod(k=1,n,1/(1-x^k +x*O(x^n))^n),n)}
    {a(n)=polcoeff(exp(sum(m=1,n,A008485(m)*x^m/m)+x*O(x^n)),n)} \\ Paul D. Hanna, Feb 06 2012

Formula

G.f. A(x) satisfies:
(1) A(x) = (1/x)*Series_Reversion(x*eta(x)), where eta(x) is Dedekind's eta(q) function without the q^(1/24) factor.
(2) A(x) = 1/G(x) where G(x) is g.f. of A109084.
(3) A(x) = 1 / Product_{n>=1} (1 - x^n*A(x)^n).
(4) A(x) = Sum_{n>=0} x^n*A(x)^n / Product_{k=1..n} (1-x^k*A(x)^k).
(5) A(x) = Sum_{n>=0} (x*A(x))^(n^2) / Product_{k=1..n} (1-x^k*A(x)^k)^2.
(6) A(x) = exp( Sum_{n>=1} (x^n/n) * A(x)^n/(1 - x^n*A(x)^n) ). - Paul D. Hanna, Jun 01 2011
Logarithmic derivative yields A008485, where A008485(n) is the number of partitions of n into parts of n kinds. - Paul D. Hanna, Feb 06 2012
a(n) ~ c * d^n / n^(3/2), where d = A270915 = 5.3527013334866426877724... and c = A366022 = 0.489635226684303373081541660578468619322416625... . - Vaclav Kotesovec, Nov 21 2016

A058698 a(n) = p(P(n)), P = primes (A000040), p = partition numbers (A000041).

Original entry on oeis.org

2, 3, 7, 15, 56, 101, 297, 490, 1255, 4565, 6842, 21637, 44583, 63261, 124754, 329931, 831820, 1121505, 2679689, 4697205, 6185689, 13848650, 23338469, 49995925, 133230930, 214481126, 271248950, 431149389, 541946240, 851376628, 3913864295, 5964539504, 11097645016
Offset: 1

Views

Author

N. J. A. Sloane, Dec 31 2000

Keywords

Comments

Number of partitions of n-th prime. - Omar E. Pol, Aug 05 2011

Examples

			a(2) = 3 because the second prime is 3 and there are three partitions of 3: {1, 1, 1}, {1, 2}, {3}.
		

Crossrefs

Programs

  • Haskell
    import Data.MemoCombinators (memo2, integral)
    a058698 n = a058698_list !! (n-1)
    a058698_list = map (pMemo 1) a000040_list where
       pMemo = memo2 integral integral p
       p _ 0 = 1
       p k m | m < k     = 0
             | otherwise = pMemo k (m - k) + pMemo (k + 1) m
    -- Reinhard Zumkeller, Aug 09 2015
  • Mathematica
    Table[PartitionsP[Prime[n]], {n, 30}] (* Vladimir Joseph Stephan Orlovsky, Dec 05 2008 *)

Formula

a(n) = A000041(A000040(n)). - Omar E. Pol, Aug 05 2011

A339278 Irregular triangle read by rows T(n,k), (n >= 1, k >= 1), in which the partition number A000041(n-1) is the length of row n and every column k is A000203, the sum of divisors function.

Original entry on oeis.org

1, 3, 4, 1, 7, 3, 1, 6, 4, 3, 1, 1, 12, 7, 4, 3, 3, 1, 1, 8, 6, 7, 4, 4, 3, 3, 1, 1, 1, 1, 15, 12, 6, 7, 7, 4, 4, 3, 3, 3, 3, 1, 1, 1, 1, 13, 8, 12, 6, 6, 7, 7, 4, 4, 4, 4, 3, 3, 3, 3, 1, 1, 1, 1, 1, 1, 1, 18, 15, 8, 12, 12, 6, 6, 7, 7, 7, 7, 4, 4, 4, 4, 3, 3, 3, 3, 3, 3, 3, 1, 1, 1, 1, 1, 1, 1, 1
Offset: 1

Views

Author

Omar E. Pol, Nov 29 2020

Keywords

Comments

The sum of row n equals A138879(n), the sum of all parts in the last section of the set of partitions of n.
T(n,k) is also the number of cubic cells (or cubes) added at the n-th stage in the k-th level starting from the base in the tower described in A221529, assuming that the tower is an object under construction (see the example). - Omar E. Pol, Jan 20 2022

Examples

			Triangle begins:
   1;
   3;
   4,  1;
   7,  3,  1;
   6,  4,  3, 1, 1;
  12,  7,  4, 3, 3, 1, 1;
   8,  6,  7, 4, 4, 3, 3, 1, 1, 1, 1;
  15, 12,  6, 7, 7, 4, 4, 3, 3, 3, 3, 1, 1, 1, 1;
  13,  8, 12, 6, 6, 7, 7, 4, 4, 4, 4, 3, 3, 3, 3, 1, 1, 1, 1, 1, 1, 1;
...
From _Omar E. Pol_, Jan 13 2022: (Start)
Illustration of the first six rows of triangle showing the growth of the symmetric tower described in A221529:
    Level k: 1              2         3        4       5      6     7
Stage
  n   _ _ _ _ _ _ _ _
     |            _  |
  1  |           |_| |
     |_ _ _ _ _ _ _ _|
     |          _    |
     |         | |_  |
  2  |         |_ _| |
     |_ _ _ _ _ _ _ _|_ _ _ _ _ _
     |        _      |        _  |
     |       | |     |       |_| |
  3  |       |_|_ _  |           |
     |         |_ _| |           |
     |_ _ _ _ _ _ _ _|_ _ _ _ _ _|_ _ _ _ _
     |      _        |      _    |      _  |
     |     | |       |     | |_  |     |_| |
  4  |     | |_      |     |_ _| |         |
     |     |_  |_ _  |           |         |
     |       |_ _ _| |           |         |
     |_ _ _ _ _ _ _ _|_ _ _ _ _ _|_ _ _ _ _|_ _ _ _ _ _ _ _
     |    _          |    _      |    _    |    _  |    _  |
     |   | |         |   | |     |   | |_  |   |_| |   |_| |
     |   | |         |   |_|_ _  |   |_ _| |       |       |
  5  |   |_|_        |     |_ _| |         |       |       |
     |       |_ _ _  |           |         |       |       |
     |       |_ _ _| |           |         |       |       |
     |_ _ _ _ _ _ _ _|_ _ _ _ _ _|_ _ _ _ _|_ _ _ _|_ _ _ _|_ _ _ _ _ _
     |  _            |  _        |  _      |  _    |  _    |  _  |  _  |
     | | |           | | |       | | |     | | |_  | | |_  | |_| | |_| |
     | | |           | | |_      | |_|_ _  | |_ _| | |_ _| |     |     |
     | | |_ _        | |_  |_ _  |   |_ _| |       |       |     |     |
  6  | |_    |       |   |_ _ _| |         |       |       |     |     |
     |   |_  |_ _ _  |           |         |       |       |     |     |
     |     |_ _ _ _| |           |         |       |       |     |     |
     |_ _ _ _ _ _ _ _|_ _ _ _ _ _|_ _ _ _ _|_ _ _ _|_ _ _ _|_ _ _|_ _ _|
.
Every cell in the diagram of the symmetric representation of sigma represents a cubic cell or cube.
For n = 6 and k = 3 we add four cubes at 6th stage in the third level of the structure of the tower starting from the base so T(6,3) = 4.
For n = 9 another connection with the tower is as follows:
First we take the columns from the above triangle and build a new triangle in which all columns start at row 1 as shown below:
.
   1,  1,  1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1;
   3,  3,  3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3;
   4,  4,  4, 4, 4, 4, 4, 4, 4, 4, 4;
   7,  7,  7, 7, 7, 7, 7;
   6,  6,  6, 6, 6;
  12, 12, 12;
   8,  8;
  15;
  13;
.
Then we rotate the triangle by 90 degrees as shown below:
                                       _
  1;                                  | |
  1;                                  | |
  1;                                  | |
  1;                                  | |
  1;                                  | |
  1;                                  | |
  1;                                  |_|_
  1, 3;                               |   |
  1, 3;                               |   |
  1, 3;                               |   |
  1, 3;                               |_ _|_
  1, 3, 4;                            |   | |
  1, 3, 4;                            |   | |
  1, 3, 4;                            |   | |
  1, 3, 4;                            |_ _|_|_
  1, 3, 4, 7;                         |     | |
  1, 3, 4, 7;                         |_ _ _| |_
  1, 3, 4, 7, 6;                      |     |   |
  1, 3, 4, 7, 6;                      |_ _ _|_ _|_
  1, 3, 4, 7, 6, 12;                  |_ _ _ _| | |_
  1, 3, 4, 7, 6, 12, 8;               |_ _ _ _|_|_ _|_ _
  1, 3, 4, 7, 6, 12, 8, 15; 13;       |_ _ _ _ _|_ _|_ _|
.
                                         Lateral view
                                         of the tower
.                                      _ _ _ _ _ _ _ _ _
                                      |_| | | | | | |   |
                                      |_ _|_| | | | |   |
                                      |_ _|  _|_| | |   |
                                      |_ _ _|    _|_|   |
                                      |_ _ _|  _|    _ _|
                                      |_ _ _ _|     |
                                      |_ _ _ _|  _ _|
                                      |         |
                                      |_ _ _ _ _|
.
                                           Top view
                                         of the tower
.
The sum of the m-th row of the new triangle equals A024916(j) where j is the length of the m-th row, equaling the number of cubic cells in the m-th level of the tower. For example: the last row of triangle has 9 terms and the sum of the last row is 1 + 3 + 4 + 7 + 6 + 12 + 8 + 15 + 13 = A024916(9) = 69, equaling the number of cubes in the base of the tower. (End)
		

Crossrefs

Sum of divisors of A336811.
Row n has length A000041(n-1).
Every column gives A000203.
The length of the m-th block in row n is A187219(m), m >= 1.
Row sums give A138879.
Cf. A337209 (another version).
Cf. A272172 (analog for the stepped pyramid described in A245092).

Programs

  • Mathematica
    A339278[rowmax_]:=Table[Flatten[Table[ConstantArray[DivisorSigma[1,n-m],PartitionsP[m]-PartitionsP[m-1]],{m,0,n-1}]],{n,rowmax}];
    A339278[15] (* Generates 15 rows *) (* Paolo Xausa, Feb 17 2023 *)
  • PARI
    f(n) = numbpart(n-1);
    T(n, k) = {if (k > f(n), error("invalid k")); if (k==1, return (sigma(n))); my(s=0); while (k <= f(n-1), s++; n--;); sigma(1+s);}
    tabf(nn) = {for (n=1, nn, for (k=1, f(n), print1(T(n,k), ", ");); print;);} \\ Michel Marcus, Jan 13 2021
    
  • PARI
    A339278(rowmax)=vector(rowmax,n,concat(vector(n,m,vector(numbpart(m-1)-numbpart(m-2),i,sigma(n-m+1)))));
    A339278(15) \\ Generates 15 rows \\ Paolo Xausa, Feb 17 2023

Formula

a(m) = A000203(A336811(m)).
T(n,k) = A000203(A336811(n,k)).

A211971 Column 0 of square array A211970 (in which column 1 is A000041).

Original entry on oeis.org

1, 1, 2, 4, 6, 10, 16, 24, 36, 54, 78, 112, 160, 224, 312, 432, 590, 802, 1084, 1452, 1936, 2568, 3384, 4440, 5800, 7538, 9758, 12584, 16160, 20680, 26376, 33520, 42468, 53644, 67552, 84832, 106246, 132706, 165344, 205512, 254824, 315256, 389168, 479368
Offset: 0

Views

Author

Omar E. Pol, Jun 10 2012

Keywords

Comments

Partial sums give A015128. - Omar E. Pol, Jan 09 2014

Crossrefs

Programs

  • Mathematica
    Flatten[{1, Differences[Table[Sum[PartitionsP[n-k]*PartitionsQ[k], {k, 0, n}], {n, 0, 60}]]}] (* Vaclav Kotesovec, Oct 25 2016 *)
    CoefficientList[Series[(1 - x)/EllipticTheta[4, 0, x], {x, 0, 43}], x] (* Robert G. Wilson v, Mar 06 2018 *)

Formula

a(n) ~ exp(Pi*sqrt(n))*Pi / (16*n^(3/2)) * (1 - (3/Pi + Pi/4)/sqrt(n) + (3/2 + 3/Pi^2+ Pi^2/24)/n). - Vaclav Kotesovec, Oct 25 2016, extended Nov 04 2016
G.f.: (1 - x)/theta_4(x), where theta_4() is the Jacobi theta function. - Ilya Gutkovskiy, Mar 05 2018

A061260 G.f.: Product_{k>=1} (1-y*x^k)^(-numbpart(k)), where numbpart(k) = number of partitions of k, cf. A000041.

Original entry on oeis.org

1, 2, 1, 3, 2, 1, 5, 6, 2, 1, 7, 11, 6, 2, 1, 11, 23, 15, 6, 2, 1, 15, 40, 32, 15, 6, 2, 1, 22, 73, 67, 37, 15, 6, 2, 1, 30, 120, 134, 79, 37, 15, 6, 2, 1, 42, 202, 255, 172, 85, 37, 15, 6, 2, 1, 56, 320, 470, 348, 187, 85, 37, 15, 6, 2, 1, 77, 511, 848, 697, 397, 194, 85, 37, 15, 6, 2, 1
Offset: 1

Views

Author

Vladeta Jovovic, Apr 23 2001

Keywords

Comments

Multiset transformation of A000041. - R. J. Mathar, Apr 30 2017
Number of orderless twice-partitions of n of length k. A twice-partition of n is a choice of a partition of each part in a partition of n. The T(5,3) = 6 orderless twice-partitions: (3)(1)(1), (21)(1)(1), (111)(1)(1), (2)(2)(1), (2)(11)(1), (11)(11)(1). - Gus Wiseman, Mar 23 2018

Examples

			:  1;
:  2,   1;
:  3,   2,   1;
:  5,   6,   2,   1;
:  7,  11,   6,   2,  1;
: 11,  23,  15,   6,  2,  1;
: 15,  40,  32,  15,  6,  2,  1;
: 22,  73,  67,  37, 15,  6,  2, 1;
: 30, 120, 134,  79, 37, 15,  6, 2, 1;
: 42, 202, 255, 172, 85, 37, 15, 6, 2, 1;
		

Crossrefs

Row sums: A001970, first column: A000041.
T(2,n) gives A061261,

Programs

  • Maple
    b:= proc(n, i, p) option remember; `if`(p>n, 0, `if`(n=0, 1,
          `if`(min(i, p)<1, 0, add(b(n-i*j, i-1, p-j)*binomial(
           combinat[numbpart](i)+j-1, j), j=0..min(n/i, p)))))
        end:
    T:= (n, k)-> b(n$2, k):
    seq(seq(T(n, k), k=1..n), n=1..14);  # Alois P. Heinz, Apr 13 2017
  • Mathematica
    b[n_, i_, p_] := b[n, i, p] = If[p > n, 0, If[n == 0, 1, If[Min[i, p] < 1, 0, Sum[b[n - i*j, i - 1, p - j]*Binomial[PartitionsP[i] + j - 1, j], {j, 0, Min[n/i, p]}]]]];
    T[n_, k_] := b[n, n, k];
    Table[T[n, k], {n, 1, 14}, {k, 1, n}] // Flatten (* Jean-François Alcover, May 17 2018, after Alois P. Heinz *)

A167137 E.g.f.: P(exp(x)-1) where P(x) is the g.f. of the partition numbers (A000041).

Original entry on oeis.org

1, 1, 5, 31, 257, 2551, 30065, 407191, 6214577, 105530071, 1972879025, 40213910551, 886979957297, 21044674731991, 534313527291185, 14448883517785111, 414475305054698417, 12568507978358276311, 401658204472560090545, 13490011548122407566871, 474964861088609044357937, 17491333169997896126211031
Offset: 0

Views

Author

Paul D. Hanna, Nov 03 2009

Keywords

Comments

CONJECTURE: Sum_{n>=0} a(n)^m * log(1+x)^n/n! is an integer series in x for all integer m>0; see A167138 and A167139 for examples.
From Peter Bala, Jul 07 2022: (Start)
Conjecture: Let k be a positive integer. The sequence obtained by reducing a(n) modulo k is eventually periodic with the period dividing phi(k) = A000010(k). For example, modulo 16 we obtain the sequence [1, 1, 5, 15, 1, 7, 1, 7, 1, 7, ...], with an apparent period of 2 beginning at a(4).
More generally, we conjecture that the same property holds for integer sequences having an e.g.f. of the form G(exp(x) - 1), where G(x) is an integral power series. (End)

Examples

			E.g.f.: A(x) = 1 + x + 5*x^2/2! + 31*x^3/3! + 257*x^4/4! +...
A(log(1+x)) = P(x) = 1 + x + 2*x^2 + 3*x^3 + 5*x^4 + 7*x^5 +...
		

Crossrefs

Programs

  • Mathematica
    Table[Sum[PartitionsP[k]*StirlingS2[n, k]*k!, {k, 0, n}], {n, 0, 20}] (* Vaclav Kotesovec, May 21 2018 *)
    nmax = 20; CoefficientList[Series[1/QPochhammer[E^x - 1], {x, 0, nmax}], x] * Range[0, nmax]! (* Vaclav Kotesovec, Nov 22 2021 *)
  • PARI
    {a(n)=if(n==0,1,n!*polcoeff(exp(sum(m=1,n,sigma(m)*(exp(x+x*O(x^n))-1)^m/m) ),n))}
    
  • PARI
    {Stirling2(n, k)=if(k<0||k>n, 0, sum(i=0, k, (-1)^i*binomial(k, i)/k!*(k-i)^n))}
    {a(n)=sum(k=0,n,numbpart(k)*Stirling2(n, k)*k!)}
    
  • PARI
    x='x+O('x^66); Vec( serlaplace( 1/eta(exp(x)-1) ) ) \\ Joerg Arndt, Sep 18 2013

Formula

a(n) = Sum_{k=0..n} A000041(k)*Stirling2(n,k)*k! where A000041 is the partition numbers.
E.g.f.: exp( Sum_{n>=1} sigma(n)*[exp(x)-1]^n/n ).
Sum_{n>=0} a(n) * log(1+x)^n/n! = g.f. of the partition numbers (A000041).
Sum_{n>=0} a(n)^2*log(1+x)^n/n! = g.f. of A167138.
From Peter Bala, Sep 18 2013: (Start)
Sum {n >= 0} (-1)^n*a(n)*(log(1 - x))^n/n! = 1 + x + 3*x^2 + 8*x^3 + 21*x^4 + ... is the o.g.f. of A218482.
a(n) is always odd. Congruences for n >= 1: a(2*n) = 2 (mod 3); a(4*n) = 2 (mod 5); a(6*n) = 0 (mod 7); a(10*n) = 7 (mod 11); a(12*n) = 5 (mod 13); a(16*n) = 0 (mod 17). (End)
From Vaclav Kotesovec, Jun 17 2018: (Start)
a(n) ~ n! * exp((1/log(2) - 1) * Pi^2 / 24 + Pi*sqrt(n/(3*log(2)))) / (4 * sqrt(3) * n * (log(2))^n).
a(n) ~ sqrt(Pi) * exp((1/log(2) - 1) * Pi^2 / 24 + Pi*sqrt(n/(3*log(2))) - n) * n^(n + 1/2) / (2^(3/2) * sqrt(3) * n * (log(2))^n). (End)
Previous Showing 21-30 of 3899 results. Next