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

A100883 Number of partitions of n in which the sequence of frequencies of the summands is nondecreasing.

Original entry on oeis.org

1, 1, 2, 3, 5, 6, 11, 13, 19, 26, 36, 43, 64, 77, 102, 129, 169, 205, 268, 323, 413, 504, 629, 751, 947, 1131, 1384, 1661, 2024, 2393, 2919, 3442, 4136, 4884, 5834, 6836, 8162, 9531, 11262, 13155, 15493, 17981, 21138, 24472, 28571, 33066, 38475, 44305
Offset: 0

Views

Author

David S. Newman, Nov 21 2004

Keywords

Comments

From Gus Wiseman, Jan 21 2019: (Start)
Also the number of semistandard Young tableaux where the rows are constant and the entries sum to n. For example, the a(8) = 19 tableaux are:
8 44 2222 11111111
.
1 2 11 3 111 22 1111 11 11111 1111 111111
7 6 6 5 5 4 4 33 3 22 2
.
1 1 11 111
2 3 2 2
5 4 4 3
(End)

Examples

			a(5) = 6 because, of the 7 unrestricted partitions of 5, only one, 2 + 2 + 1, has a decreasing sequence of frequencies. Two is used twice, but 1 is used only once.
		

Crossrefs

Programs

  • Maple
    b:= proc(n, i, t) option remember; `if`(n<0, 0, `if`(n=0, 1,
          `if`(i=1, `if`(n>=t, 1, 0), `if`(i=0, 0, b(n, i-1, t)+
           add(b(n-i*j, i-1, j), j=t..floor(n/i))))))
        end:
    a:= n-> b(n$2, 1):
    seq(a(n), n=0..60);  # Alois P. Heinz, Jul 03 2014
  • Mathematica
    b[n_, i_, t_] := b[n, i, t] = If[n<0, 0, If[n == 0, 1, If[i == 1, If[n >= t, 1, 0], If[i == 0, 0, b[n, i-1, t] + Sum[b[n-i*j, i-1, j], {j, t, Floor[n/i]}]]]]]; a[n_] := b[n, n, 1]; Table[a[n], {n, 0, 60}] (* Jean-François Alcover, Mar 16 2015, after Alois P. Heinz *)
    Table[Length[Select[IntegerPartitions[n],OrderedQ[Length/@Split[#]]&]],{n,20}] (* Gus Wiseman, Jan 21 2019 *)

Extensions

More terms from Vladeta Jovovic, Nov 23 2004

A006951 Number of conjugacy classes in GL(n,2).

Original entry on oeis.org

1, 1, 3, 6, 14, 27, 60, 117, 246, 490, 1002, 1998, 4053, 8088, 16284, 32559, 65330, 130626, 261726, 523374, 1047690, 2095314, 4192479, 8384808, 16773552, 33546736, 67101273, 134202258, 268420086, 536839446, 1073710914, 2147420250, 4294904430, 8589807438
Offset: 0

Views

Author

Keywords

Comments

Unlabeled permutations of sets. - Christian G. Bower, Jan 29 2004
From Joerg Arndt, Jan 02 2013: (Start)
Set q=2 and f(m)=q^(m-1)*(q-1), then a(n) is the sum over all partitions P of n over all products Product_{k=1..L} f(m_k) where L is the number of different parts in the partition P=[p_1^m_1, p_2^m_2, ..., p_L^m_L], see the Macdonald reference.
Setting q to a prime power gives the sequence "Number of conjugacy classes in GL(n,q)":
q=3: A006952, q=4: A049314, q=5: A049315, q=7: A049316, q=8: A182603,
q=9: A182604, q=11: A182605, q=13: A182606, q=16: A182607, q=17: A182608,
q=19: A182609, q=23: A182610, q=25: A182611, q=27: A182612.
Sequences where q is not a prime power are:
q=6: A221578, q=10: A221579, q=12: A221580,
q=14: A221581, q=15: A221582, q=18: A221583, q=20: A221584.
(End)
From Gus Wiseman, Jan 21 2019: (Start)
Also the number of ways to split an integer partition of n into consecutive constant subsequences. For example, the a(5) = 27 ways (subsequences shown as rows) are:
5 11111
.
4 3 3 22 2 1111 1 111 11
1 2 11 1 111 1 1111 11 111
.
3 2 2 2 111 1 1 11 11 1
1 2 11 1 1 111 1 11 1 11
1 1 1 11 1 1 111 1 11 11
.
2 11 1 1 1
1 1 11 1 1
1 1 1 11 1
1 1 1 1 11
.
1
1
1
1
1
(End)

Examples

			For the 5 partitions of 4 (namely [1^4]; [2,1^2]; [2^2]; [3,1]; [4]) we have
(f(m) = 2^(m-1)*(2-1) = 2^(m-1) and)
f([1^4]) = 2^3 = 8,
f([2,1^2]) = 1*2^1 = 2,
f([2^2]) = 2^1 = 2,
f([3,1]) = 1*1 = 1,
f([4]) = 1,
the sum is 8+2+2+1+1 = 14 = a(4).
- _Joerg Arndt_, Jan 02 2013
		

References

  • W. D. Smith, personal communication.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Magma
    /* The program does not work for n>19: */
    [1] cat [NumberOfClasses(GL(n,2)): n in [1..19]]; // Sergei Haller (sergei(AT)sergei-haller.de), Dec 21 2006; edited by Vincenzo Librandi Jan 24 2013
    
  • Maple
    with(numtheory):
    b:= n-> add(phi(d)*2^(n/d), d=divisors(n))/n-1:
    a:= proc(n) option remember; `if`(n=0, 1,
           add(add(d*b(d), d=divisors(j)) *a(n-j), j=1..n)/n)
        end:
    seq(a(n), n=0..40);  # Alois P. Heinz, Oct 20 2012
  • Mathematica
    b[n_] := Sum[EulerPhi[d]*2^(n/d), {d, Divisors[n]}]/n-1; a[n_] := a[n] = If[n == 0, 1, Sum[Sum[d*b[d], {d, Divisors[j]}]*a[n-j], {j, 1, n}]/n]; Table[a[n], {n, 0, 40}] (* Jean-François Alcover, Feb 17 2014, after Alois P. Heinz *)
    Table[Sum[2^(Length[ptn]-Length[Split[ptn]]),{ptn,IntegerPartitions[n]}],{n,30}] (* Gus Wiseman, Jan 21 2019 *)
  • PARI
    N=66; x='x+O('x^N);
    gf=prod(n=1,N, (1-x^n)/(1-2*x^n)  );
    v=Vec(gf)
    /* Joerg Arndt, Jan 02 2013 */

Formula

G.f.: Product_{n>=1} (1-x^n)/(1-2*x^n). - Joerg Arndt, Jan 02 2013
The number a(n) of conjugacy classes in the group GL(n, q) is the coefficient of t^n in Product_{k>=1} (1-t^k)/(1-q*t^k). - Noam Katz (noamkj(AT)hotmail.com), Mar 30 2001
Euler transform of A008965. - Christian G. Bower, Jan 29 2004
a(n) ~ 2^n - (1+sqrt(2) + (-1)^n*(1-sqrt(2))) * 2^(n/2-1). - Vaclav Kotesovec, Nov 21 2015
G.f.: exp(Sum_{k>=1} ( Sum_{d|k} d*(2^(k/d) - 1) ) * x^k/k). - Ilya Gutkovskiy, Sep 27 2018

Extensions

More terms from Christian G. Bower, Jan 29 2004

A075900 Expansion of g.f.: Product_{n>0} 1/(1 - 2^(n-1)*x^n).

Original entry on oeis.org

1, 1, 3, 7, 19, 43, 115, 259, 659, 1523, 3731, 8531, 20883, 47379, 113043, 259219, 609683, 1385363, 3245459, 7344531, 17028499, 38579603, 88585619, 199845267, 457864595, 1028904339, 2339763603, 5256820115, 11896157587, 26626389395
Offset: 0

Views

Author

N. J. A. Sloane, Oct 15 2002

Keywords

Comments

Number of compositions of partitions of n. a(3) = 7: 3, 21, 12, 111, 2|1, 11|1, 1|1|1. - Alois P. Heinz, Sep 16 2019
Also the number of ways to split an integer composition of n into consecutive subsequences with weakly decreasing (or increasing) sums. - Gus Wiseman, Jul 13 2020
This sequence is obtained from the generalized Euler transform in A266964 by taking f(n) = 1, g(n) = 2^(n-1). - Seiichi Manyama, Aug 22 2020

Examples

			From _Gus Wiseman_, Jul 13 2020: (Start)
The a(0) = 1 through a(4) = 19 splittings:
  ()  (1)  (2)      (3)          (4)
           (1,1)    (1,2)        (1,3)
           (1),(1)  (2,1)        (2,2)
                    (1,1,1)      (3,1)
                    (2),(1)      (1,1,2)
                    (1,1),(1)    (1,2,1)
                    (1),(1),(1)  (2,1,1)
                                 (2),(2)
                                 (3),(1)
                                 (1,1,1,1)
                                 (1,1),(2)
                                 (1,2),(1)
                                 (2),(1,1)
                                 (2,1),(1)
                                 (1,1),(1,1)
                                 (1,1,1),(1)
                                 (2),(1),(1)
                                 (1,1),(1),(1)
                                 (1),(1),(1),(1)
(End)
		

Crossrefs

Row sums of A327549.
The strict case is A304961.
Partitions of partitions are A001970.
Splittings with equal sums are A074854.
Splittings of compositions are A133494.
Splittings of partitions are A323583.
Splittings with distinct sums are A336127.
Starting with a reversed partition gives A316245.
Starting with a partition instead of composition gives A336136.

Programs

  • Magma
    m:=80;
    R:=PowerSeriesRing(Integers(), m);
    Coefficients(R!( 1/(&*[1-2^(j-1)*x^j: j in [1..m+2]]) )); // G. C. Greubel, Jan 25 2024
    
  • Maple
    oo := 101; t1 := mul(1/(1-x^n/2),n=1..oo): t2 := series(t1,x,oo-1): t3 := seriestolist(t2): A075900 := n->2^n*t3[n+1];
    with(combinat); A075900 := proc(n) local i,t1,t2,t3; t1 := partition(n); t2 := 0; for i from 1 to nops(t1) do t3 := t1[i]; t2 := t2+2^(n-nops(t3)); od: t2; end;
  • Mathematica
    b[n_]:= b[n]= Sum[d*2^(n - n/d), {d, Divisors[n]}];
    a[0]= 1; a[n_]:= a[n]= 1/n*Sum[b[k]*a[n-k], {k,n}];
    Table[a[n], {n,0,30}] (* Jean-François Alcover, Mar 20 2014, after Vladeta Jovovic, fixed by Vaclav Kotesovec, Mar 08 2018 *)
  • Maxima
    s(m,n):=if nVladimir Kruchinin, Sep 06 2014 */
    
  • PARI
    {a(n)=polcoeff(prod(k=1,n,1/(1-2^(k-1)*x^k+x*O(x^n))),n)} \\ Paul D. Hanna, Jan 13 2013
    
  • PARI
    {a(n)=polcoeff(exp(sum(k=1,n+1,x^k/(k*(1-2^k*x^k)+x*O(x^n)))),n)} \\ Paul D. Hanna, Jan 13 2013
    
  • SageMath
    m=80;
    def A075900_list(prec):
        P. = PowerSeriesRing(QQ, prec)
        return P( 1/product(1-2^(j-1)*x^j for j in range(1,m+1)) ).list()
    A075900_list(m) # G. C. Greubel, Jan 25 2024

Formula

a(n) = Sum_{ partitions n = c_1 + ... + c_k } 2^(n-k). If p(n, m) = number of partitions of n into m parts, a(n) = Sum_{m=1..n} p(n, m)*2^(n-m).
G.f.: Sum_{n>=0} (a(n)/2^n)*x^n = Product_{n>0} 1/(1-x^n/2). - Vladeta Jovovic, Feb 11 2003
a(n) = 1/n*Sum_{k=1..n} A080267(k)*a(n-k). - Vladeta Jovovic, Feb 11 2003
G.f.: exp( Sum_{n>=1} x^n / (n*(1 - 2^n*x^n)) ). - Paul D. Hanna, Jan 13 2013
a(n) = s(1,n), a(0)=1, where s(m,n) = Sum_{k=m..n/2} 2^(k-1)*s(k, n-k) + 2^(n-1), s(n,n) = 2^(n-1), s(m,n)=0, m>. - Vladimir Kruchinin, Sep 06 2014
a(n) ~ 2^(n-2) * (Pi^2 - 6*log(2)^2)^(1/4) * exp(sqrt((Pi^2 - 6*log(2)^2)*n/3)) / (3^(1/4) * sqrt(Pi) * n^(3/4)). - Vaclav Kotesovec, Mar 09 2018

Extensions

More terms from Vladeta Jovovic, Feb 11 2003

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

Original entry on oeis.org

1, 1, 2, 6, 12, 32, 72, 176, 384, 960, 2112, 4992, 11264, 26112, 58368, 136192, 301056, 688128, 1548288, 3489792, 7766016, 17596416, 38993920, 87293952, 194248704, 432537600, 957349888, 2132803584, 4699717632, 10406068224, 23001563136, 50683969536, 111434268672, 245819768832
Offset: 0

Views

Author

Ilya Gutkovskiy, May 22 2018

Keywords

Comments

Number of compositions of partitions of n into distinct parts. a(3) = 6: 3, 21, 12, 111, 2|1, 11|1. - Alois P. Heinz, Sep 16 2019
Also the number of ways to split a composition of n into contiguous subsequences with strictly decreasing sums. - Gus Wiseman, Jul 13 2020
This sequence is obtained from the generalized Euler transform in A266964 by taking f(n) = -1, g(n) = (-1) * 2^(n-1). - Seiichi Manyama, Aug 22 2020

Examples

			From _Gus Wiseman_, Jul 13 2020: (Start)
The a(0) = 1 through a(4) = 12 splittings:
  ()  (1)  (2)    (3)        (4)
           (1,1)  (1,2)      (1,3)
                  (2,1)      (2,2)
                  (1,1,1)    (3,1)
                  (2),(1)    (1,1,2)
                  (1,1),(1)  (1,2,1)
                             (2,1,1)
                             (3),(1)
                             (1,1,1,1)
                             (1,2),(1)
                             (2,1),(1)
                             (1,1,1),(1)
(End)
		

Crossrefs

The non-strict version is A075900.
Starting with a reversed partition gives A323583.
Starting with a partition gives A336134.
Partitions of partitions are A001970.
Splittings with equal sums are A074854.
Splittings of compositions are A133494.
Splittings with distinct sums are A336127.

Programs

  • Mathematica
    nmax = 33; CoefficientList[Series[Product[(1 + 2^(k - 1) x^k), {k, 1, nmax}], {x, 0, nmax}], x]
  • PARI
    N=40; x='x+O('x^N); Vec(prod(k=1, N, 1+2^(k-1)*x^k)) \\ Seiichi Manyama, Aug 22 2020

Formula

G.f.: Product_{k>=1} (1 + A011782(k)*x^k).
a(n) ~ 2^n * exp(2*sqrt(-polylog(2, -1/2)*n)) * (-polylog(2, -1/2))^(1/4) / (sqrt(6*Pi) * n^(3/4)). - Vaclav Kotesovec, Sep 19 2019

A323583 Number of ways to split an integer partition of n into consecutive subsequences.

Original entry on oeis.org

1, 1, 3, 7, 17, 37, 83, 175, 373, 773, 1603, 3275, 6693, 13557, 27447, 55315, 111397, 223769, 449287, 900795, 1805465, 3615929, 7240327, 14491623, 29001625, 58027017, 116093259, 232237583, 464558201, 929224589, 1858623819, 3717475031, 7435314013, 14871103069
Offset: 0

Views

Author

Gus Wiseman, Jan 19 2019

Keywords

Examples

			The a(3) = 7 ways to split an integer partition of 3 into consecutive subsequences are (3), (21), (2)(1), (111), (11)(1), (1)(11), (1)(1)(1).
		

Crossrefs

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0, 1/2, `if`(i<1, 0,
          b(n, i-1) +`if`(i>n, 0, 2*b(n-i, i))))
        end:
    a:= n-> ceil(b(n$2)):
    seq(a(n), n=0..33);  # Alois P. Heinz, Jan 01 2023
  • Mathematica
    Table[Sum[2^(Length[ptn]-1),{ptn,IntegerPartitions[n]}],{n,40}]
    (* Second program: *)
    (1/2) CoefficientList[1 - 1/QPochhammer[2, x] + O[x]^100 , x] (* Jean-François Alcover, Jan 02 2022, after Vladimir Reshetnikov in A070933 *)

Formula

a(n) = A070933(n)/2.
O.g.f.: (1/2)*Product_{n >= 1} 1/(1 - 2*x^n).
G.f.: 1 + Sum_{k>=1} 2^(k - 1) * x^k / Product_{j=1..k} (1 - x^j). - Ilya Gutkovskiy, Jan 28 2020

A229915 Number of espalier polycubes of a given volume in dimension 3.

Original entry on oeis.org

1, 1, 3, 5, 10, 14, 26, 34, 57, 76, 116, 150, 227, 284, 408, 520, 718, 895, 1226, 1508, 2018, 2487, 3248, 3968, 5160, 6235, 7970, 9653, 12179, 14630, 18367, 21924, 27241, 32506, 39985, 47492, 58203, 68752, 83613, 98730, 119269, 140224, 168799, 197758, 236753, 277052, 329867, 384852, 457006, 531500, 628338
Offset: 0

Views

Author

Matthieu Deneufchâtel, Oct 03 2013

Keywords

Comments

A pyramid polycube is obtained by gluing together horizontal plateaux (parallelepipeds of height 1) in such a way that (0,0,0) belongs to the first plateau and each cell with coordinates (0,b,c) belongs to the first plateau such that b,c >= 0. If the cell with coordinates (a,b,c) belongs to the (a+1)-st plateau (a>0), then the cell with coordinates (a-1, b, c) belongs to the a-th plateau.
An espalier polycube is a special pyramid such that each plateau contains the cell with coordinates (a,0,0).

Crossrefs

Formula

The generating function for the numbers of espaliers of height h and volumes v_1 , ... v_h is x_1^{n_1} * ... x_h^{n_h} / ((1-x_1^{n_1}) *(1-x_1^{n_1}*x_2^{n_2}) *... *(1-x_1^{n_1}*x_2^{n_2}*...x_h^{n_h})).
This sequence is obtained with x_1 = ... = x_h = p by summing over n_1>= ... >= n_h>=1 and then over h.

Extensions

a(0)=1 prepended by Seiichi Manyama, Aug 20 2020
Showing 1-6 of 6 results.