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

A052847 G.f.: 1 / Product_{k>=1} (1-x^k)^(k-1).

Original entry on oeis.org

1, 0, 1, 2, 4, 6, 12, 18, 33, 52, 88, 138, 229, 354, 568, 880, 1378, 2110, 3260, 4942, 7527, 11320, 17031, 25394, 37842, 55956, 82630, 121300, 177677, 258980, 376626, 545352, 787784, 1133764, 1627657, 2329020, 3324559, 4731396, 6717774, 9512060
Offset: 0

Views

Author

encyclopedia(AT)pommard.inria.fr, Jan 25 2000

Keywords

Comments

Euler transform of sequence [0,1,2,3,...]. - Michael Somos, Jul 02 2004
Number of partitions of n objects of 2 colors, where each part must contain at least one of each color. - Franklin T. Adams-Watters, Jan 23 2006
Number of partitions of n without 1s, one kind of 2s, two kinds of 3s, etc. - Joerg Arndt, Jul 31 2011
From Vaclav Kotesovec, Oct 17 2015: (Start)
In general, if v>=0 and g.f. = Product_{k>=1} 1/(1-x^(k+v))^k, then a(n) ~ d1(v) * n^(v^2/6 - 25/36) * exp(-Pi^4 * v^2 / (432*Zeta(3)) + 3*Zeta(3)^(1/3) * n^(2/3)/2^(2/3) - v * Pi^2 * n^(1/3) / (3 * 2^(4/3) * Zeta(3)^(1/3))) / (sqrt(3*Pi) * 2^(v^2/6 + 11/36) * Zeta(3)^(v^2/6 - 7/36)), where Zeta(3) = A002117.
d1(v) = exp(Integral_{x=0..infinity} (1/(x*exp((v-1)*x) * (exp(x)-1)^2) - (6*v^2-1) / (12*x*exp(x)) + v/x^2 - 1/x^3) dx).
d1(v) = (exp(Zeta'(-1) - v*Zeta'(0))) / Product_{j=0..v-1} j!, where Zeta'(0) = -A075700, Zeta'(-1) = A084448 and Product_{j=0..v-1} j! = A000178(v-1).
d1(v) = exp(1/12) * (2*Pi)^(v/2) / (A * G(v+1)), where A = A074962 is the Glaisher-Kinkelin constant and G is the Barnes G-function.
(End)

Examples

			1 + x^2 + 2*x^3 + 4*x^4 + 6*x^5 + 12*x^6 + 18*x^7 + 33*x^8 + 52*x^9 + ...
From _Gus Wiseman_, Jan 22 2019: (Start)
The partitions described in Franklin T. Adams-Watters's comment are (n = 2 through 6):
  {{12}}  {{112}}  {{1112}}    {{11112}}    {{111112}}
          {{122}}  {{1122}}    {{11122}}    {{111122}}
                   {{1222}}    {{11222}}    {{111222}}
                   {{12}{12}}  {{12222}}    {{112222}}
                               {{12}{112}}  {{122222}}
                               {{12}{122}}  {{112}{112}}
                                            {{112}{122}}
                                            {{12}{1112}}
                                            {{12}{1122}}
                                            {{12}{1222}}
                                            {{122}{122}}
                                            {{12}{12}{12}}
(End)
		

Crossrefs

Cf. A000219 (v=0), A052847 (v=1), A263358 (v=2), A263359 (v=3), A263360 (v=4), A263361 (v=5), A263362 (v=6), A263363 (v=7), A263364 (v=8).

Programs

  • Maple
    spec := [S,{B=Sequence(Z,1 <= card),C=Prod(B,B),S= Set(C)},unlabeled]: seq(combstruct[count](spec, size=n), n=0..20);
    with(numtheory): etr:= proc(p) local b; b:=proc(n) option remember; local d, j; if n=0 then 1 else add(add(d*p(d), d=divisors(j)) *b(n-j), j=1..n)/n fi end end: a:=etr(n-> n-1): seq(a(n), n=0..50); # Vaclav Kotesovec, Mar 04 2015 after Alois P. Heinz
  • Mathematica
    Clear[a]; a[n_]:= a[n] = 1/n*Sum[(DivisorSigma[2,k]-DivisorSigma[1,k])*a[n-k],{k,1,n}]; a[0]=1; Table[a[n],{n,0,100}] (* Vaclav Kotesovec, Mar 04 2015 *)
    nmax = 40; CoefficientList[Series[Product[1/(1-x^(k+1))^k, {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Oct 16 2015 *)
  • PARI
    {a(n) = if( n<0, 0, polcoeff( 1 / prod(k=1, n, (1 - x^k + x*O(x^n))^(k-1)), n))}

Formula

a(n) = 1/n*Sum_{k=1..n} (sigma[2](k)-sigma[1](k))*a(n-k).
G.f.: exp( Sum_{k>0} ( x^k / (1 - x^k) )^2 / k ).
G.f.: exp( sum(n>=0, (sigma[2](n)-sigma[1](n)) *x^n/n ) ). - Joerg Arndt, Jul 31 2011
a(n) ~ 2^(1/36) * Zeta(3)^(1/36) * exp(1/12 - Pi^4/(432*Zeta(3)) - Pi^2 * n^(1/3) / (3 * 2^(4/3) * Zeta(3)^(1/3)) + 3 * Zeta(3)^(1/3) * n^(2/3) / 2^(2/3)) / (A * 3^(1/2) * n^(19/36)), where A = A074962 = 1.2824271291... is the Glaisher-Kinkelin constant. - Vaclav Kotesovec, Mar 07 2015

Extensions

Edited by Vladeta Jovovic, Sep 10 2002

A321615 Triangle read by rows: T(n,k) is the number of k X k integer matrices with sum of elements n, with no zero rows or columns, up to row and column permutation.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 1, 2, 1, 0, 1, 6, 3, 1, 0, 1, 9, 13, 3, 1, 0, 1, 17, 38, 20, 3, 1, 0, 1, 23, 97, 82, 23, 3, 1, 0, 1, 36, 217, 311, 126, 24, 3, 1, 0, 1, 46, 453, 968, 624, 151, 24, 3, 1, 0, 1, 65, 868, 2825, 2637, 933, 162, 24, 3, 1, 0, 1, 80, 1585, 7394, 10098, 4942, 1132, 165, 24, 3, 1
Offset: 0

Views

Author

Andrew Howroyd, Nov 14 2018

Keywords

Comments

Also the number of non-isomorphic multiset partitions of weight n with k parts and k vertices, where the weight of a multiset partition is the sum of sizes of its parts. - Gus Wiseman, Nov 18 2018

Examples

			Triangle begins:
    1
    0  1
    0  1    1
    0  1    2    1
    0  1    6    3    1
    0  1    9   13    3    1
    0  1   17   38   20    3    1
    0  1   23   97   82   23    3    1
    0  1   36  217  311  126   24    3    1
    0  1   46  453  968  624  151   24    3    1
    0  1   65  868 2825 2637  933  162   24    3    1
		

Crossrefs

Programs

  • Mathematica
    (* See A318795 for M[m, n, k]. *)
    T[n_, k_] := M[k, k, n] - 2 M[k, k-1, n] + M[k-1, k-1, n];
    Table[T[n, k], {n, 1, 11}, {k, 1, n}] // Flatten (* Jean-François Alcover, Nov 24 2018, from PARI *)
  • PARI
    \\ See A318795 for M.
    T(n, k) = if(k==0, n==0, M(k, k, n) - 2*M(k, k-1, n) + M(k-1, k-1, n));
    
  • PARI
    \\ See A340652 for G.
    T(n)={[Vecrev(p) | p<-Vec(1 + sum(k=1, n, y^k*(polcoef(G(k, n, n, y), k, y) - polcoef(G(k-1, n, n, y), k, y))))]}
    { my(A=T(10)); for(i=1, #A, print(A[i])) } \\ Andrew Howroyd, Jan 16 2024

Extensions

Column k=0 inserted by Andrew Howroyd, Jan 17 2024

A323654 Number of non-isomorphic multiset partitions of weight n with no constant parts and only two distinct vertices.

Original entry on oeis.org

1, 0, 1, 1, 3, 3, 8, 9, 20, 26, 50, 69, 125, 177, 301, 440, 717, 1055, 1675, 2471, 3835, 5660, 8627, 12697, 19095, 27978, 41581, 60650, 89244, 129490, 188925, 272676, 394809, 566882, 815191, 1164510, 1664295, 2365698, 3361844, 4756030, 6723280, 9468138, 13319299
Offset: 0

Views

Author

Gus Wiseman, Jan 22 2019

Keywords

Comments

First differs from A304967 at a(10) = 50, A304967(10) = 49.
The weight of a multiset partition is the sum of sizes of its parts. Weight is generally not the same as number of vertices.
Also the number of positive integer matrices with only two columns and sum of entries equal to n, up to row and column permutations.

Examples

			Non-isomorphic representatives of the a(2) = 1 through a(7) = 9 multiset partitions:
  {{12}}  {{122}}  {{1122}}    {{11222}}    {{111222}}      {{1112222}}
                   {{1222}}    {{12222}}    {{112222}}      {{1122222}}
                   {{12}{12}}  {{12}{122}}  {{122222}}      {{1222222}}
                                            {{112}{122}}    {{112}{1222}}
                                            {{12}{1122}}    {{12}{11222}}
                                            {{12}{1222}}    {{12}{12222}}
                                            {{122}{122}}    {{122}{1122}}
                                            {{12}{12}{12}}  {{122}{1222}}
                                                            {{12}{12}{122}}
Inequivalent representatives of the a(8) = 20 matrices:
  [4 4] [3 5] [2 6] [1 7]
.
  [1 1] [1 1] [1 1] [2 1] [2 1] [1 2] [1 2] [3 1] [2 2] [2 2] [1 3]
  [3 3] [2 4] [1 5] [2 3] [1 4] [2 3] [1 4] [1 3] [2 2] [1 3] [1 3]
.
  [1 1] [1 1] [1 1] [1 1]
  [1 1] [1 1] [2 1] [1 2]
  [2 2] [1 3] [1 2] [1 2]
.
  [1 1]
  [1 1]
  [1 1]
  [1 1]
		

Crossrefs

Programs

  • PARI
    EulerT(v)={Vec(exp(x*Ser(dirmul(v,vector(#v,n,1/n))))-1, -#v)}
    seq(n)={concat(1,(EulerT(vector(n, k, k-1)) + EulerT(vector(n, k, if(k%2, 0, (k+2)\4))))/2)} \\ Andrew Howroyd, Aug 26 2019

Formula

a(2*n) = (A052847(2*n) + A003293(n))/2; a(2*n+1) = A052847(2*n+1)/2. - Andrew Howroyd, Aug 26 2019

Extensions

Terms a(11) and beyond from Andrew Howroyd, Aug 26 2019

A323655 Number of non-isomorphic multiset partitions of weight n with at most 2 distinct vertices, or with at most 2 (not necessarily distinct) edges.

Original entry on oeis.org

1, 1, 4, 7, 19, 35, 80, 149, 307, 566, 1092, 1974, 3643, 6447, 11498, 19947, 34636, 58974, 100182, 167713, 279659, 461056, 756562, 1230104, 1990255, 3195471, 5105540, 8103722, 12801925, 20107448, 31439978, 48907179, 75755094, 116797754, 179354540, 274253042
Offset: 0

Views

Author

Gus Wiseman, Jan 22 2019

Keywords

Comments

The weight of a multiset partition is the sum of sizes of its parts. Weight is generally not the same as number of vertices.
Also the number of nonnegative integer matrices with only one or two columns, no zero rows or columns, and sum of entries equal to n, up to row and column permutations.

Examples

			Non-isomorphic representatives of the a(1) = 1 through a(4) = 19 multiset partitions with at most 2 distinct vertices:
  {{1}}  {{11}}    {{111}}      {{1111}}
         {{12}}    {{122}}      {{1122}}
         {{1}{1}}  {{1}{11}}    {{1222}}
         {{1}{2}}  {{1}{22}}    {{1}{111}}
                   {{2}{12}}    {{11}{11}}
                   {{1}{1}{1}}  {{1}{122}}
                   {{1}{2}{2}}  {{11}{22}}
                                {{12}{12}}
                                {{1}{222}}
                                {{12}{22}}
                                {{2}{122}}
                                {{1}{1}{11}}
                                {{1}{1}{22}}
                                {{1}{2}{12}}
                                {{1}{2}{22}}
                                {{2}{2}{12}}
                                {{1}{1}{1}{1}}
                                {{1}{1}{2}{2}}
                                {{1}{2}{2}{2}}
Non-isomorphic representatives of the a(1) = 1 through a(4) = 19 multiset partitions with at most 2 edges:
  {{1}}  {{11}}    {{111}}    {{1111}}
         {{12}}    {{122}}    {{1122}}
         {{1}{1}}  {{123}}    {{1222}}
         {{1}{2}}  {{1}{11}}  {{1233}}
                   {{1}{22}}  {{1234}}
                   {{1}{23}}  {{1}{111}}
                   {{2}{12}}  {{11}{11}}
                              {{1}{122}}
                              {{11}{22}}
                              {{12}{12}}
                              {{1}{222}}
                              {{12}{22}}
                              {{1}{233}}
                              {{12}{33}}
                              {{1}{234}}
                              {{12}{34}}
                              {{13}{23}}
                              {{2}{122}}
                              {{3}{123}}
Inequivalent representatives of the a(4) = 19 matrices:
  [4] [2 2] [1 3]
.
  [1] [1 0] [1 0] [0 1] [2] [2 0] [1 1] [1 1]
  [3] [1 2] [0 3] [1 2] [2] [0 2] [1 1] [0 2]
.
  [1] [1 0] [1 0] [1 0] [0 1]
  [1] [1 0] [0 1] [0 1] [0 1]
  [2] [0 2] [1 1] [0 2] [1 1]
.
  [1] [1 0] [1 0]
  [1] [1 0] [0 1]
  [1] [0 1] [0 1]
  [1] [0 1] [0 1]
		

Crossrefs

Programs

  • PARI
    EulerT(v)={Vec(exp(x*Ser(dirmul(v,vector(#v,n,1/n))))-1, -#v)}
    seq(n)={concat(1, (EulerT(vector(n, k, k+1)) + EulerT(vector(n, k, if(k%2, 0, (k+6)\4))))/2)} \\ Andrew Howroyd, Aug 26 2019

Formula

a(2*n) = (A005380(2*n) + A005986(n))/2; a(2*n+1) = A005380(2*n+1)/2. - Andrew Howroyd, Aug 26 2019

Extensions

Terms a(11) and beyond from Andrew Howroyd, Aug 26 2019

A323656 Number of non-isomorphic multiset partitions of weight n with exactly 2 distinct vertices, or with exactly 2 (not necessarily distinct) edges.

Original entry on oeis.org

0, 0, 2, 4, 14, 28, 69, 134, 285, 536, 1050, 1918, 3566, 6346, 11363, 19771, 34405, 58677, 99797, 167223, 279032, 460264, 755560, 1228849, 1988680, 3193513, 5103104, 8100712, 12798207, 20102883, 31434374, 48900337, 75746745, 116787611, 179342230, 274238159
Offset: 0

Views

Author

Gus Wiseman, Jan 22 2019

Keywords

Comments

The weight of a multiset partition is the sum of sizes of its parts. Weight is generally not the same as number of vertices.
Also the number of nonnegative integer matrices with only two columns, no zero rows or columns, and sum of entries equal to n, up to row and column permutations.

Examples

			Non-isomorphic representatives of the a(2) = 2 through a(4) = 14 multiset partitions with exactly 2 distinct vertices:
  {{12}}    {{122}}      {{1122}}
  {{1}{2}}  {{1}{22}}    {{1222}}
            {{2}{12}}    {{1}{122}}
            {{1}{2}{2}}  {{11}{22}}
                         {{12}{12}}
                         {{1}{222}}
                         {{12}{22}}
                         {{2}{122}}
                         {{1}{1}{22}}
                         {{1}{2}{12}}
                         {{1}{2}{22}}
                         {{2}{2}{12}}
                         {{1}{1}{2}{2}}
                         {{1}{2}{2}{2}}
Non-isomorphic representatives of the a(2) = 2 through a(4) = 14 multiset partitions with exactly 2 edges:
  {{1}{1}}  {{1}{11}}  {{1}{111}}
  {{1}{2}}  {{1}{22}}  {{11}{11}}
            {{1}{23}}  {{1}{122}}
            {{2}{12}}  {{11}{22}}
                       {{12}{12}}
                       {{1}{222}}
                       {{12}{22}}
                       {{1}{233}}
                       {{12}{33}}
                       {{1}{234}}
                       {{12}{34}}
                       {{13}{23}}
                       {{2}{122}}
                       {{3}{123}}
Inequivalent representatives of the a(4) = 14 matrices:
  [2 2] [1 3]
.
  [1 0] [1 0] [0 1] [2 0] [1 1] [1 1]
  [1 2] [0 3] [1 2] [0 2] [1 1] [0 2]
.
  [1 0] [1 0] [1 0] [0 1]
  [1 0] [0 1] [0 1] [0 1]
  [0 2] [1 1] [0 2] [1 1]
.
  [1 0] [1 0]
  [1 0] [0 1]
  [0 1] [0 1]
  [0 1] [0 1]
		

Crossrefs

Programs

  • PARI
    EulerT(v)={Vec(exp(x*Ser(dirmul(v,vector(#v,n,1/n))))-1, -#v)}
    seq(n)={concat(0, (EulerT(vector(n, k, k+1)) + EulerT(vector(n, k, if(k%2, 0, (k+6)\4))))/2 - EulerT(vector(n,k,1)))} \\ Andrew Howroyd, Aug 26 2019

Formula

a(n) = A323655(n) - A000041(n). - Andrew Howroyd, Aug 26 2019

Extensions

Terms a(11) and beyond from Andrew Howroyd, Aug 26 2019

A055007 Number of nonnegative integer 4 X 4 matrices with no zero rows or columns and with sum of elements equal to n.

Original entry on oeis.org

1, 0, 0, 0, 24, 528, 4648, 26224, 112666, 401424, 1246000, 3476368, 8905432, 21266208, 47875272, 102482048, 210000931, 414160240, 789572072, 1460372624, 2628456428, 4615495808, 7924479264, 13328517504, 21997272036, 35674700896, 56926058920, 89477437120
Offset: 0

Views

Author

Vladeta Jovovic, May 30 2000

Keywords

Crossrefs

Formula

Number of nonnegative integer p X q matrices with no zero rows or columns and with sum of elements equal to n is Sum_{k=0...q} (-1)^k*C(q, k)*m(p, q-k, n) where m(p, q, n)=Sum_{k=0..p} (-1)^k*C(p, k)*C((p-k)*q+n-1, n).
For p = q = 4 we get a(n) = (1/15!)*(n^15 + 120*n^14 + 6580*n^13 + 218400*n^12 + 4637542*n^11 + 61261200*n^10 + 423591740*n^9 + 164392800*n^8 - 17247717487*n^7 - 47940252360*n^6 + 346941238280*n^5 + 557885764800*n^4 - 4897231459056*n^3 + 8643549191040*n^2 - 5894285241600*n + 1307674368000).
G.f.: -(16*x^15 -192*x^14 +1040*x^13 -3356*x^12 +7200*x^11 -10952*x^10 +12544*x^9 -11712*x^8 +9664*x^7 -7088*x^6 +4224*x^5 -1844*x^4 +560*x^3 -120*x^2 +16*x -1) / (x -1)^16. - Colin Barker, Jul 11 2013

Extensions

More terms from James Sellers, May 31 2000

A055005 Number of nonnegative integer 3 X 3 matrices with no zero rows or columns and with sum of elements equal to n.

Original entry on oeis.org

1, 0, 0, 6, 63, 306, 1038, 2844, 6750, 14437, 28521, 52911, 93258, 157509, 256581, 405171, 622719, 934542, 1373158, 1979820, 2806281, 3916812, 5390496, 7323822, 9833604, 13060251, 17171415, 22366045, 28878876, 36985383, 47007231
Offset: 0

Views

Author

Vladeta Jovovic, May 30 2000

Keywords

Crossrefs

Programs

  • PARI
    a(n)=([0,1,0,0,0,0,0,0,0; 0,0,1,0,0,0,0,0,0; 0,0,0,1,0,0,0,0,0; 0,0,0,0,1,0,0,0,0; 0,0,0,0,0,1,0,0,0; 0,0,0,0,0,0,1,0,0; 0,0,0,0,0,0,0,1,0; 0,0,0,0,0,0,0,0,1; 1,-9,36,-84,126,-126,84,-36,9]^n*[1;0;0;6;63;306;1038;2844;6750])[1,1] \\ Charles R Greathouse IV, Aug 14 2023

Formula

Number of nonnegative integer p X q matrices with no zero rows or columns and with sum of elements equal to n is Sum_{k=0...q} (-1)^k*C(q, k)*m(p, q-k, n) where m(p, q, n)=Sum_{k=0..p} (-1)^k*C(p, k)*C((p-k)*q+n-1, n).
For p=q=3 we get a(n)=C(n + 8, 8) - 6*C(n + 5, 5) + 9*C(n + 3, 3) + 6*C(n + 2, 2) - 18*C(n + 1, 1) + 9=(1/8!)*(n^8 + 36*n^7 + 546*n^6 + 2520*n^5 - 7791*n^4 - 43596*n^3 + 148364*n^2 - 140400*n + 40320).
G.f.: -(9*x^8-54*x^7+132*x^6-171*x^5+135*x^4-78*x^3+36*x^2-9*x+1) / (x-1)^9. - Colin Barker, Jul 13 2013
Showing 1-7 of 7 results.