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-10 of 82 results. Next

A117488 The number 1 followed by 2*k+1 terms from column k of table A115994.

Original entry on oeis.org

1, 1, 2, 3, 1, 2, 5, 8, 14, 1, 2, 5, 10, 18, 30, 49, 1, 2, 5, 10, 20, 34, 59, 94, 149, 1, 2, 5, 10, 20, 36, 63, 104, 169, 264, 405, 1, 2, 5, 10, 20, 36, 65, 108, 179, 284, 445, 676, 1017, 1, 2, 5, 10, 20, 36, 65, 110, 183, 294, 465, 716, 1089, 1622, 2387, 1, 2, 5, 10, 20, 36, 65
Offset: 1

Views

Author

Alford Arnold, Mar 22 2006

Keywords

Examples

			Column two of table A115994 begins 1 2 5 8 14 20 30 ... A006918
so row three of A117488 is 1 2 5 8 14.
Triangle starts
1
1 2 3
1 2 5 8 14
1 2 5 10 18 30 49
1 2 5 10 20 34 59 94 149
1 2 5 10 20 36 63 104 169 264 405
1 2 5 10 20 36 65 108 179 284 445 676 1017
1 2 5 10 20 36 65 110 183 294 465 716 1089 1622 2387
		

Crossrefs

Programs

  • Maple
    A026820 := proc(n,k) if k > n then combinat[numbpart](n,n) ; else combinat[numbpart](n,k) ; fi ; end: A115994 := proc(n,k) local i ; add(A026820(i,k)*A026820(n-k^2-i,k),i=0..n-k^2) ; end: A117488 := proc(n,k) if k >= 2*n then 0 ; else if n = 1 then 1; else A115994(k+n^2-2*n,n-1) ; fi ; fi ; end: for n from 1 to 10 do for k from 1 to 2*n-1 do printf("%d ",A117488(n,k)) ; od ; od ; # R. J. Mathar, Feb 22 2007

Extensions

More terms from R. J. Mathar, Feb 22 2007

A128626 A triangular array distributing the values of sequence A072213 (cf. A115994).

Original entry on oeis.org

1, 4, 1, 9, 20, 1, 16, 140, 74, 1, 25, 572, 1136, 224, 1, 36, 1785, 8866, 6685, 604, 1, 49, 4600, 47152, 88380, 31851, 1492, 1, 64, 10416, 194282, 737059, 665542, 130808, 3458, 1, 81, 21320, 665769, 4512584, 8211274, 4105870, 479826, 7602, 1, 100, 40425
Offset: 0

Views

Author

Alford Arnold, Mar 15 2007

Keywords

Examples

			A115994 distributes the numeric partition sequence A000041 and A072213 records partition values for sequence A000290 (the squares).
Therefore the table begins:
   1;
   4,    1;
   9,   20,     1;
  16,  140,    74,     1;
  25,  572,  1136,   224,     1;
  36, 1785,  8866,  6685,   604,    1;
  49, 4600, 47152, 88380, 31851, 1492, 1;
  ...
		

Crossrefs

Programs

  • Maple
    nn:=8: g:=sum(t^k*q^(k^2)/product((1-q^j)^2, j=1..k), k=1..nn): gser:=series(g, q=0, nn^2+1): for n from 1 to nn do P[n]:=coeff(gser, q^(n^2)) od: for n from 1 to nn do seq(coeff(P[n], t^j), j=1..n); od; # Nathaniel Johnston, Apr 30 2011

A006918 a(n) = binomial(n+3, 3)/4 for odd n, n*(n+2)*(n+4)/24 for even n.

Original entry on oeis.org

0, 1, 2, 5, 8, 14, 20, 30, 40, 55, 70, 91, 112, 140, 168, 204, 240, 285, 330, 385, 440, 506, 572, 650, 728, 819, 910, 1015, 1120, 1240, 1360, 1496, 1632, 1785, 1938, 2109, 2280, 2470, 2660, 2870, 3080, 3311, 3542, 3795, 4048, 4324, 4600, 4900, 5200, 5525, 5850, 6201, 6552, 6930
Offset: 0

Views

Author

Keywords

Comments

Maximal number of inconsistent triples in a tournament on n+2 nodes [Kac]. - corrected by Leen Droogendijk, Nov 10 2014
a(n-4) is the number of aperiodic necklaces (Lyndon words) with 4 black beads and n-4 white beads.
a(n-3) is the maximum number of squares that can be formed from n lines, for n>=3. - Erich Friedman; corrected by Leen Droogendijk, Nov 10 2014
Number of trees with diameter 4 where at most 2 vertices 1 away from the graph center have degree > 2. - Jon Perry, Jul 11 2003
a(n+1) is the number of partitions of n into parts of two kinds, with at most two parts of each kind. Also a(n-3) is the number of partitions of n with Durfee square of size 2. - Franklin T. Adams-Watters, Jan 27 2006
Factoring the g.f. as x/(1-x)^2 times 1/(1-x^2)^2 we find that the sequence equals (1, 2, 3, 4, ...) convolved with (1, 0, 2, 0, 3, 0, 4, ...), A000027 convolved with its aerated variant. - Gary W. Adamson, May 01 2009
Starting with "1" = triangle A171238 * [1,2,3,...]. - Gary W. Adamson, Dec 05 2009
The Kn21, Kn22, Kn23, Fi2 and Ze2 triangle sums, see A180662 for their definitions, of the Connell-Pol triangle A159797 are linear sums of shifted versions of this sequence, e.g., Kn22(n) = a(n+1) + a(n) + 2*a(n-1) + a(n-2) and Fi2(n) = a(n) + 4*a(n-1) + a(n-2). - Johannes W. Meijer, May 20 2011
For n>3, a(n-4) is the number of (w,x,y,z) having all terms in {1,...,n} and w+x+y+z=|x-y|+|y-z|. - Clark Kimberling, May 23 2012
a(n) is the number of (w,x,y) having all terms in {0,...,n} and w+x+y < |w-x|+|x-y|. - Clark Kimberling, Jun 13 2012
For n>0 number of inequivalent (n-1) X 2 binary matrices, where equivalence means permutations of rows or columns or the symbol set. - Alois P. Heinz, Aug 17 2014
Number of partitions p of n+5 such that p[3] = 2. Examples: a(1)=1 because we have (2,2,2); a(2)=2 because we have (2,2,2,1) and (3,2,2); a(3)=5 because we have (2,2,2,1,1), (2,2,2,2), (3,2,2,1), (3,3,2), and (4,2,2). See the R. P. Stanley reference. - Emeric Deutsch, Oct 28 2014
Sum over each antidiagonal of A243866. - Christopher Hunt Gribble, Apr 02 2015
Number of nonisomorphic outer planar graphs of order n>=3, size n+2, and maximum degree 3. - Christian Barrientos and Sarah Minion, Feb 27 2018
a(n) is the number of 2413-avoiding odd Grassmannian permutations of size n+1. - Juan B. Gil, Mar 09 2023

Examples

			G.f. = x + 2*x^2 + 5*x^3 + 8*x^4 + 14*x^5 + 20*x^6 + 30*x^7 + 40*x^8 + 55*x^9 + ...
From _Gus Wiseman_, Apr 06 2019: (Start)
The a(4 - 3) = 1 through a(8 - 3) = 14 integer partitions with Durfee square of length 2 are the following (see Franklin T. Adams-Watters's second comment). The Heinz numbers of these partitions are given by A325164.
  (22)  (32)   (33)    (43)     (44)
        (221)  (42)    (52)     (53)
               (222)   (322)    (62)
               (321)   (331)    (332)
               (2211)  (421)    (422)
                       (2221)   (431)
                       (3211)   (521)
                       (22111)  (2222)
                                (3221)
                                (3311)
                                (4211)
                                (22211)
                                (32111)
                                (221111)
The a(0 + 1) = 1 through a(4 + 1) = 14 integer partitions of n into parts of two kinds with at most two parts of each kind are the following (see Franklin T. Adams-Watters's first comment).
  ()()  ()(1)  ()(2)   ()(3)    ()(4)
        (1)()  (2)()   (3)()    (4)()
               ()(11)  (1)(2)   (1)(3)
               (1)(1)  ()(21)   ()(22)
               (11)()  (2)(1)   (2)(2)
                       (21)()   (22)()
                       (1)(11)  ()(31)
                       (11)(1)  (3)(1)
                                (31)()
                                (11)(2)
                                (1)(21)
                                (2)(11)
                                (21)(1)
                                (11)(11)
The a(6 - 5) = 1 through a(10 - 5) = 14 integer partitions whose third part is 2 are the following (see Emeric Deutsch's comment). The Heinz numbers of these partitions are given by A307373.
  (222)  (322)   (332)    (432)     (442)
         (2221)  (422)    (522)     (532)
                 (2222)   (3222)    (622)
                 (3221)   (3321)    (3322)
                 (22211)  (4221)    (4222)
                          (22221)   (4321)
                          (32211)   (5221)
                          (222111)  (22222)
                                    (32221)
                                    (33211)
                                    (42211)
                                    (222211)
                                    (322111)
                                    (2221111)
(End)
		

References

  • J. M. Borwein, D. H. Bailey and R. Girgensohn, Experimentation in Mathematics, A K Peters, Ltd., Natick, MA, 2004. x+357 pp. See p. 147.
  • M. Kac, An example of "counting without counting", Philips Res. Reports, 30 (1975), 20*-22* [Special issue in honour of C. J. Bouwkamp].
  • E. V. McLaughlin, Numbers of factorizations in non-unique factorial domains, Senior Thesis, Allegeny College, Meadville, PA, 2004.
  • K. B. Reid and L. W. Beineke "Tournaments", pp. 169-204 in L. W. Beineke and R. J. Wilson, editors, Selected Topics in Graph Theory, Academic Press, NY, 1978, p. 186, Theorem 6.11.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 1, 2nd ed., 2012, Exercise 4.16, pp. 530, 552.
  • W. A. Whitworth, DCC Exercises in Choice and Chance, Stechert, NY, 1945, p. 33.

Crossrefs

Cf. A000031, A001037, A028723, A051168. a(n) = T(n,4), array T as in A051168.
Cf. A000094.
Cf. A171238. - Gary W. Adamson, Dec 05 2009
Row sums of A173997. - Gary W. Adamson, Mar 05 2010
Column k=2 of A242093. Column k=2 of A115720 and A115994.

Programs

  • Haskell
    a006918 n = a006918_list !! n
    a006918_list = scanl (+) 0 a008805_list
    -- Reinhard Zumkeller, Feb 01 2013
    
  • Magma
    [Floor(Binomial(n+4, 4)/(n+4))-Floor((n+2)/8)*(1+(-1)^n)/2: n in [0..60]]; // Vincenzo Librandi, Nov 10 2014
  • Maple
    with(combstruct):ZL:=[st,{st=Prod(left,right),left=Set(U,card=r),right=Set(U,card=r),U=Sequence(Z,card>=3)}, unlabeled]: subs(r=1,stack): seq(count(subs(r=2,ZL),size=m),m=11..58) ; # Zerinvary Lajos, Mar 09 2007
    A006918 := proc(n)
        if type(n,'even') then
            n*(n+2)*(n+4)/24 ;
        else
            binomial(n+3,3)/4 ;
        fi ;
    end proc: # R. J. Mathar, May 17 2016
  • Mathematica
    f[n_]:=If[EvenQ[n],(n(n+2)(n+4))/24,Binomial[n+3,3]/4]; Join[{0},Array[f,60]]  (* Harvey P. Dale, Apr 20 2011 *)
    durf[ptn_]:=Length[Select[Range[Length[ptn]],ptn[[#]]>=#&]];
    Table[Length[Select[IntegerPartitions[n],durf[#]==2&]],{n,0,30}] (* Gus Wiseman, Apr 06 2019 *)
  • PARI
    { parttrees(n)=local(pt,k,nk); if (n%2==0, pt=(n/2+1)^2, pt=ceil(n/2)*(ceil(n/2)+1)); pt+=floor(n/2); for (x=1,floor(n/2),pt+=floor(x/2)+floor((n-x)/2)); if (n%2==0 && n>2, pt-=floor(n/4)); k=1; while (3*k<=n, for (x=k,floor((n-k)/2), pt+=floor(k/2); if (x!=k, pt+=floor(x/2)); if ((n-x-k)!=k && (n-x-k)!=x, pt+=floor((n-x-k)/2))); k++); pt }
    
  • PARI
    {a(n) = n += 2; (n^3 - n * (2-n%2)^2) / 24}; /* Michael Somos, Aug 15 2009 */
    

Formula

G.f.: x/((1-x)^2*(1-x^2)^2) = x/((1+x)^2*(1-x)^4).
0, 0, 0, 1, 2, 5, 8, 14, ... has a(n) = (Sum_{k=0..n} floor(k(n-k)/2))/2. - Paul Barry, Sep 14 2003
0, 0, 0, 0, 0, 1, 2, 5, 8, 14, 20, 30, 40, 55, ... has a(n) = binomial(floor(1/2 n), 3) + binomial(floor(1/2 n + 1/2), 3) [Eke]. - N. J. A. Sloane, May 12 2012
a(0)=0, a(1)=1, a(n) = (2/(n-1))*a(n-1) + ((n+3)/(n-1))*a(n-2). - Benoit Cloitre, Jun 28 2004
a(n) = floor(binomial(n+4, 4)/(n+4)) - floor((n+2)/8)(1+(-1)^n)/2. - Paul Barry, Jan 01 2005
a(n+1) = a(n) + binomial(floor(n/2)+2,2), i.e., first differences are A008805. Convolution of A008619 with itself, then shifted right (or A004526 with itself, shifted left by 3). - Franklin T. Adams-Watters, Jan 27 2006
a(n+1) = (A027656(n) + A003451(n+5))/2 with a(1)=0. - Yosu Yurramendi, Sep 12 2008
Linear recurrence: a(n) = 2a(n-1) + a(n-2) - 4a(n-3) + a(n-4) + 2a(n-5) - a(n-6). - Jaume Oliver Lafont, Dec 05 2008
Euler transform of length 2 sequence [2, 2]. - Michael Somos, Aug 15 2009
a(n) = -a(-4-n) for all n in Z.
a(n+1) + a(n) = A002623(n). - Johannes W. Meijer, May 20 2011
a(n) = (n+2)*(2*n*(n+4)-3*(-1)^n+3)/48. - Bruno Berselli, May 21 2011
a(2n) = A007290(n+2). - Jon Perry, Nov 10 2014
G.f.: (1/(1-x)^4-1/(1-x^2)^2)/4. - Herbert Kociemba, Oct 23 2016
E.g.f.: (x*(18 + 9*x + x^2)*cosh(x) + (6 + 15*x + 9*x^2 + x^3)*sinh(x))/24. - Stefano Spezia, Dec 07 2021
From Amiram Eldar, Mar 20 2022: (Start)
Sum_{n>=1} 1/a(n) = 75/4 - 24*log(2).
Sum_{n>=1} (-1)^(n+1)/a(n) = 69/4 - 24*log(2). (End)

A001522 Number of n-stacks with strictly receding walls, or the number of Type A partitions of n in the sense of Auluck (1951).

Original entry on oeis.org

1, 1, 1, 1, 2, 3, 5, 7, 10, 14, 19, 26, 35, 47, 62, 82, 107, 139, 179, 230, 293, 372, 470, 591, 740, 924, 1148, 1422, 1756, 2161, 2651, 3244, 3957, 4815, 5844, 7075, 8545, 10299, 12383, 14859, 17794, 21267, 25368, 30207, 35902, 42600, 50462, 59678, 70465, 83079, 97800, 114967, 134956, 158205, 185209, 216546, 252859
Offset: 0

Views

Author

Keywords

Comments

Also number of partitions of n with positive crank (n>=2), cf. A064391. - Vladeta Jovovic, Sep 30 2001
Number of smooth weakly unimodal compositions of n into positive parts such that the first and last part are 1 (smooth means that successive parts differ by at most one), see example. Dropping the requirement for unimodality gives A186085. - Joerg Arndt, Dec 09 2012
Number of weakly unimodal compositions of n where the maximal part m appears at least m times, see example. - Joerg Arndt, Jun 11 2013
Also weakly unimodal compositions of n with first part 1, maximal up-step 1, and no consecutive up-steps; see example. The smooth weakly unimodal compositions are recovered by shifting all rows above the bottom row to the left by one position with respect to the next lower row. - Joerg Arndt, Mar 30 2014
It would seem from Stanley that he regards a(0)=0 for this sequence and A001523. - Michael Somos, Feb 22 2015
From Gus Wiseman, Mar 30 2021: (Start)
Also the number of odd-length compositions of n with alternating parts strictly decreasing. These are finite odd-length sequences q of positive integers summing to n such that q(i) > q(i+2) for all possible i. The even-length version is A064428. For example, the a(1) = 1 through a(9) = 14 compositions are:
(1) (2) (3) (4) (5) (6) (7) (8) (9)
(211) (221) (231) (241) (251) (261)
(311) (312) (322) (332) (342)
(321) (331) (341) (351)
(411) (412) (413) (423)
(421) (422) (432)
(511) (431) (441)
(512) (513)
(521) (522)
(611) (531)
(612)
(621)
(711)
(32211)
(End)
In the Ferrers diagram of a partition x of n, count the dots in each diagonal parallel to the main diagonal (starting at the top-right, say). The result diag(x) is a smooth weakly unimodal composition of n into positive parts such that the first and last part are 1. For example, diag(5541) = 11233221. The function diag is many-to-one; the size of its codomain as a set is a(n). If diag(x) = diag(y), each hook of x can be slid by the same amount past the main diagonal to get y. For example, diag(5541) = diag(44331). - George Beck, Sep 26 2021
From Gus Wiseman, May 23 2022: (Start)
Conjecture: Also the number of integer partitions y of n with a fixed point y(i) = i. These partitions are ranked by A352827. The conjecture is stated at A238395, but Resta tells me he may not have had a proof. The a(1) = 1 through a(8) = 10 partitions are:
(1) (11) (111) (22) (32) (42) (52) (62)
(1111) (221) (222) (322) (422)
(11111) (321) (421) (521)
(2211) (2221) (2222)
(111111) (3211) (3221)
(22111) (4211)
(1111111) (22211)
(32111)
(221111)
(11111111)
Note that these are not the same partitions (compare A352827 to A352874), only the same count (apparently).
(End)
The above conjecture is true. See Section 4 of the Blecher-Knopfmacher paper in the Links section. - Jeremy Lovejoy, Sep 26 2022

Examples

			For a(6)=5 we have the following stacks:
.x... ..x.. ...x. .xx.
xxxxx xxxxx xxxxx xxxx xxxxxx
.
From _Joerg Arndt_, Dec 09 2012: (Start)
There are a(9) = 14 smooth weakly unimodal compositions of 9:
01:   [ 1 1 1 1 1 1 1 1 1 ]
02:   [ 1 1 1 1 1 1 2 1 ]
03:   [ 1 1 1 1 1 2 1 1 ]
04:   [ 1 1 1 1 2 1 1 1 ]
05:   [ 1 1 1 1 2 2 1 ]
06:   [ 1 1 1 2 1 1 1 1 ]
07:   [ 1 1 1 2 2 1 1 ]
08:   [ 1 1 2 1 1 1 1 1 ]
09:   [ 1 1 2 2 1 1 1 ]
10:   [ 1 1 2 2 2 1 ]
11:   [ 1 2 1 1 1 1 1 1 ]
12:   [ 1 2 2 1 1 1 1 ]
13:   [ 1 2 2 2 1 1 ]
14:   [ 1 2 3 2 1 ]
(End)
From _Joerg Arndt_, Jun 11 2013: (Start)
There are a(9) = 14 weakly unimodal compositions of 9 where the maximal part m appears at least m times:
01:  [ 1 1 1 1 1 1 1 1 1 ]
02:  [ 1 1 1 1 1 2 2 ]
03:  [ 1 1 1 1 2 2 1 ]
04:  [ 1 1 1 2 2 1 1 ]
05:  [ 1 1 1 2 2 2 ]
06:  [ 1 1 2 2 1 1 1 ]
07:  [ 1 1 2 2 2 1 ]
08:  [ 1 2 2 1 1 1 1 ]
09:  [ 1 2 2 2 1 1 ]
10:  [ 1 2 2 2 2 ]
11:  [ 2 2 1 1 1 1 1 ]
12:  [ 2 2 2 1 1 1 ]
13:  [ 2 2 2 2 1 ]
14:  [ 3 3 3 ]
(End)
From _Joerg Arndt_, Mar 30 2014: (Start)
There are a(9) = 14 compositions of 9 with first part 1, maximal up-step 1, and no consecutive up-steps:
01:  [ 1 1 1 1 1 1 1 1 1 ]
02:  [ 1 1 1 1 1 1 1 2 ]
03:  [ 1 1 1 1 1 1 2 1 ]
04:  [ 1 1 1 1 1 2 1 1 ]
05:  [ 1 1 1 1 1 2 2 ]
06:  [ 1 1 1 1 2 1 1 1 ]
07:  [ 1 1 1 1 2 2 1 ]
08:  [ 1 1 1 2 1 1 1 1 ]
09:  [ 1 1 1 2 2 1 1 ]
10:  [ 1 1 1 2 2 2 ]
11:  [ 1 1 2 1 1 1 1 1 ]
12:  [ 1 1 2 2 1 1 1 ]
13:  [ 1 1 2 2 2 1 ]
14:  [ 1 1 2 2 3 ]
(End)
G.f. = 1 + x + x^2 + x^3 + 2*x^4 + 3*x^5 + 5*x^6 + 7*x^7 + 10*x^8 + 14*x^9 + ...
		

References

  • G. E. Andrews, The reasonable and unreasonable effectiveness of number theory in statistical mechanics, pp. 21-34 of S. A. Burr, ed., The Unreasonable Effectiveness of Number Theory, Proc. Sympos. Appl. Math., 46 (1992). Amer. Math. Soc.
  • G. E. Andrews, Three-quadrant Ferrers graphs, Indian J. Math., 42 (No. 1, 2000), 1-7.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 1, 1999; see section 2.5 on page 76.

Crossrefs

A version for permutations is A002467, complement A000166.
The case of zero crank is A064410, ranked by A342192.
The case of nonnegative crank is A064428, ranked by A352873.
A strict version is A352829, complement A352828.
Conjectured to be column k = 1 of A352833.
These partitions (positive crank) are ranked by A352874.
A000700 counts self-conjugate partitions, ranked by A088902.
A064391 counts partitions by crank.
A115720 and A115994 count partitions by their Durfee square.
A257989 gives the crank of the partition with Heinz number n.
Counting compositions: A003242, A114921, A238351, A342527, A342528, A342532.
Fixed points of reversed partitions: A238352, A238394, A238395, A352822, A352830, A352872.

Programs

  • Maple
    b:= proc(n, i, t) option remember; `if`(n<=0, `if`(i=1, 1, 0),
          `if`(n<0 or i<1, 0, b(n-i, i, t)+b(n-(i-1), i-1, false)+
          `if`(t, b(n-(i+1), i+1, t), 0)))
        end:
    a:= n-> b(n-1, 1, true):
    seq(a(n), n=0..70);  # Alois P. Heinz, Feb 26 2014
    # second Maple program:
    A001522 := proc(n)
        local r,a;
        a := 0 ;
        if n = 0 then
            return 1 ;
        end if;
        for r from 1 do
            if r*(r+1) > 2*n then
                return a;
            else
                a := a-(-1)^r*combinat[numbpart](n-r*(r+1)/2) ;
            end if;
        end do:
    end proc: # R. J. Mathar, Mar 07 2015
  • Mathematica
    max = 50; f[x_] := 1 + Sum[-(-1)^k*x^(k*(k+1)/2), {k, 1, max}] / Product[(1-x^k), {k, 1, max}]; CoefficientList[ Series[ f[x], {x, 0, max}], x] (* Jean-François Alcover, Dec 27 2011, after g.f. *)
    b[n_, i_, t_] := b[n, i, t] = If[n <= 0, If[i == 1, 1, 0], If[n<0 || i<1, 0, b[n-i, i, t] + b[n - (i-1), i-1, False] + If[t, b[n - (i+1), i+1, t], 0]]]; a[n_] := b[n-1, 1, True]; Table[a[n], {n, 0, 70}] (* Jean-François Alcover, Dec 01 2015, after Alois P. Heinz *)
    Flatten[{1, Table[Sum[(-1)^(j-1)*PartitionsP[n-j*((j+1)/2)], {j, 1, Floor[(Sqrt[8*n + 1] - 1)/2]}], {n, 1, 60}]}] (* Vaclav Kotesovec, Sep 26 2016 *)
    ici[q_]:=And@@Table[q[[i]]>q[[i+2]],{i,Length[q]-2}];
    Table[If[n==0,1,Length[Select[Join@@Permutations/@Select[IntegerPartitions[n],OddQ@*Length],ici]]],{n,0,15}] (* Gus Wiseman, Mar 30 2021 *)
  • PARI
    {a(n) = if( n<1, n==0, polcoeff( sum(k=1, (sqrt(1+8*n) - 1)\2, -(-1)^k * x^((k + k^2)/2)) / eta(x + x * O(x^n)), n))}; /* Michael Somos, Jul 22 2003 */
    
  • PARI
    N=66; q='q+O('q^N);
    Vec( 1 + sum(n=1, N, q^(n^2)/(prod(k=1,n-1,1-q^k)^2*(1-q^n)) ) ) \\ Joerg Arndt, Dec 09 2012
    
  • Sage
    def A001522(n):
        if n < 4: return 1
        return (number_of_partitions(n) - [p.crank() for p in Partitions(n)].count(0))/2
    [A001522(n) for n in range(30)]  # Peter Luschny, Sep 15 2014

Formula

a(n) = (A000041(n) - A064410(n)) / 2 for n>=2.
G.f.: 1 + ( Sum_{k>=1} -(-1)^k * x^(k*(k+1)/2) ) / ( Product_{k>=1} 1-x^k ).
G.f.: 1 + ( Sum_{n>=1} q^(n^2) / ( ( Product_{k=1..n-1} 1-q^k )^2 * (1-q^n) ) ). - Joerg Arndt, Dec 09 2012
a(n) ~ exp(Pi*sqrt(2*n/3)) / (8*sqrt(3)*n) [Auluck, 1951]. - Vaclav Kotesovec, Sep 26 2016
a(n) = A000041(n) - A064428(n). - Gus Wiseman, Mar 30 2021
a(n) = A064428(n) - A064410(n). - Gus Wiseman, May 23 2022

Extensions

a(0) changed from 0 to 1 by Joerg Arndt, Mar 30 2014
Edited definition. - N. J. A. Sloane, Mar 31 2021

A115720 Triangle T(n,k) is the number of partitions of n with Durfee square k.

Original entry on oeis.org

1, 0, 1, 0, 2, 0, 3, 0, 4, 1, 0, 5, 2, 0, 6, 5, 0, 7, 8, 0, 8, 14, 0, 9, 20, 1, 0, 10, 30, 2, 0, 11, 40, 5, 0, 12, 55, 10, 0, 13, 70, 18, 0, 14, 91, 30, 0, 15, 112, 49, 0, 16, 140, 74, 1, 0, 17, 168, 110, 2, 0, 18, 204, 158, 5, 0, 19, 240, 221, 10, 0, 20, 285, 302, 20, 0, 21, 330, 407
Offset: 0

Views

Author

Keywords

Comments

T(n,k) is number of partitions of n-k^2 into parts of 2 kinds with at most k of each kind.

Examples

			Triangle starts:
  1;
  0,  1;
  0,  2;
  0,  3;
  0,  4,  1;
  0,  5,  2;
  0,  6,  5;
  0,  7,  8;
  0,  8, 14;
  0,  9, 20,  1;
  0, 10, 30,  2;
From _Gus Wiseman_, Apr 12 2019: (Start)
Row n = 9 counts the following partitions:
  (9)          (54)       (333)
  (81)         (63)
  (711)        (72)
  (6111)       (432)
  (51111)      (441)
  (411111)     (522)
  (3111111)    (531)
  (21111111)   (621)
  (111111111)  (3222)
               (3321)
               (4221)
               (4311)
               (5211)
               (22221)
               (32211)
               (33111)
               (42111)
               (222111)
               (321111)
               (2211111)
(End)
		

Crossrefs

For a version without zeros see A115994. Row lengths are A003059. Row sums are A000041. Column k = 2 is A006918. Column k = 3 is A117485.
Related triangles are A096771, A325188, A325189, A325192, with Heinz-encoded versions A263297, A325169, A065770, A325178.

Programs

  • Maple
    b:= proc(n, i) option remember;
          `if`(n=0, 1, `if`(i<1, 0, b(n, i-1)+`if`(i>n, 0, b(n-i, i))))
        end:
    T:= (n, k)-> add(b(m, k)*b(n-k^2-m, k), m=0..n-k^2):
    seq(seq(T(n, k), k=0..floor(sqrt(n))), n=0..30); # Alois P. Heinz, Apr 09 2012
  • Mathematica
    b[n_, i_] := b[n, i] = If[n==0, 1, If[i<1, 0, b[n, i-1] + If[i>n, 0, b[n-i, i]]]]; T[n_, k_] := Sum[b[m, k]*b[n-k^2-m, k], {m, 0, n-k^2}]; Table[ T[n, k], {n, 0, 30}, {k, 0, Sqrt[n]}] // Flatten (* Jean-François Alcover, Dec 03 2015, after Alois P. Heinz *)
    durf[ptn_]:=Length[Select[Range[Length[ptn]],ptn[[#]]>=#&]];
    Table[Length[Select[IntegerPartitions[n],durf[#]==k&]],{n,0,10},{k,0,Sqrt[n]}] (* Gus Wiseman, Apr 12 2019 *)

Formula

T(n,k) = Sum_{i=0..n-k^2} P*(i,k)*P*(n-k^2-i), where P*(n,k) = P(n+k,k) is the number of partitions of n objects into at most k parts.

A000701 One half of number of non-self-conjugate partitions; also half of number of asymmetric Ferrers graphs with n nodes.

Original entry on oeis.org

0, 0, 1, 1, 2, 3, 5, 7, 10, 14, 20, 27, 37, 49, 66, 86, 113, 146, 190, 242, 310, 392, 497, 623, 782, 973, 1212, 1498, 1851, 2274, 2793, 3411, 4163, 5059, 6142, 7427, 8972, 10801, 12989, 15572, 18646, 22267, 26561, 31602, 37556, 44533, 52743, 62338, 73593
Offset: 0

Views

Author

Keywords

Comments

Also number of cycle types of odd permutations.
Also number of partitions of n with an odd number of even parts. There is no restriction on the odd parts. - N. Sato, Jul 20 2005. E.g., a(6)=5 because we have [6],[4,1,1],[3,2,1],[2,2,2] and [2,1,1,1,1]. - Emeric Deutsch, Mar 02 2006
Also number of partitions of n with largest part not congruent to n modulo 2: a(2*n)=A027193(2*n), a(2*n+1)=A027187(2*n+1); a(n)=A000041(n)-A046682(n). - Reinhard Zumkeller, Apr 22 2006
From Gus Wiseman, Mar 31 2022: (Start)
Also the number of integer partitions of n with Heinz number greater than that of their conjugate, where the Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). These partitions are ranked by A352490. The complement is counted by A046682. For example, the a(n) partitions for n = 2...8 are:
(11) (111) (211) (221) (222) (331) (2222)
(1111) (2111) (2211) (2221) (3221)
(11111) (3111) (3211) (3311)
(21111) (22111) (22211)
(111111) (31111) (32111)
(211111) (41111)
(1111111) (221111)
(311111)
(2111111)
(11111111)
Also the number of integer partitions of n with Heinz number less than that of their conjugate, ranked by A352487. For example, the a(n) partitions for n = 2...8 are:
(2) (3) (4) (5) (6) (7) (8)
(31) (32) (33) (43) (44)
(41) (42) (52) (53)
(51) (61) (62)
(411) (322) (71)
(421) (422)
(511) (431)
(521)
(611)
(5111)
(End)

Examples

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

References

  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

A000700 counts self-conjugate partitions, ranked by A088902.
A330644 counts non-self-conjugate partitions, ranked by A352486.
Heinz number (rank) and partition:
- A122111 = rank of conjugate.
- A296150 = parts of partition, conjugate A321649.
- A352487 = rank less than conjugate, counted by A000701.
- A352488 = rank greater than or equal to conjugate, counted by A046682.
- A352489 = rank less than or equal to conjugate, counted by A046682.
- A352490 = rank greater than conjugate, counted by A000701.
- A352491 = rank minus conjugate.

Programs

  • Maple
    with(combinat); A000701 := n->(numbpart(n)-A000700(n))/2;
  • Mathematica
    a41 = PartitionsP; a700[n_] := SeriesCoefficient[ Product[1 + x^k, {k, 1, n, 2}], {x, 0, n}]; a[0] = 0; a[n_] := (a41[n] - a700[n])/2; Table[a[n], {n, 0, 48}] (* Jean-François Alcover, Feb 21 2012, after first formula *)
    a[ n_] := SeriesCoefficient[ (1 / QPochhammer[ x] - 1 / QPochhammer[ x, -x]) / 2, {x, 0, n}]; (* Michael Somos, Aug 25 2015 *)
    a[ n_] := SeriesCoefficient[ (1 - EllipticTheta[ 4, 0, x^2]) / (2 QPochhammer[ x]), {x, 0, n}]; (* Michael Somos, Aug 25 2015 *)
    a[ n_] := SeriesCoefficient[ QPochhammer[ -x, x] Sum[ x^(2 k) / QPochhammer[ x^2, x^2, k], {k, 1, n/2, 2}], {x, 0, n}] (* Michael Somos, Aug 25 2015 *)
    a[ n_] := If[ n < 0, 0, SeriesCoefficient[ Sum[ (1 / QPochhammer[ x, x, k]^2 - 1 / QPochhammer[ x^2, x^2, k]) x^k^2, {k, Sqrt@n}] / 2, {x, 0, n}]]; (* Michael Somos, Aug 25 2015 *)
    conj[y_]:=If[Length[y]==0,y,Table[Length[Select[y,#>=k&]],{k,1,Max[y]}]];
    Table[Length[Select[IntegerPartitions[n],Times@@Prime/@#>Times@@Prime/@conj[#]&]],{n,0,15}] (* Gus Wiseman, Mar 31 2022 *)
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( (1 - eta(x^2 + A)^2 / eta(x^4 + A) ) / (2 * eta(x + A)), n))}; /* Michael Somos, Aug 25 2015 */
    
  • PARI
    q='q+O('q^60); concat([0, 0], Vec((1-eta(q^2)^2/eta(q^4))/(2*eta(q)))) \\ Altug Alkan, Sep 26 2018

Formula

a(n) = (A000041(n) - A000700(n))/2.
From Bill Gosper, Aug 08 2005: (Start)
Sum a(n) q^n = q^2 + q^3 + 2 q^4 + 3 q^5 + 5 q^6 + 7 q^7 + ...
= -( Sum_{n>=1} (-q^2)^(n^2) ) / ( Sum_{ n = -oo..oo } (-1)^n q^(n(3n-1)/2) )
= (- q; q){oo} Sum{n>=1} q^(2(2n-1))/(q^2;q^2)_{2n-1}
= (1/(q;q)_oo - 1/(q;-q)_oo)/2
= (1/(q;q)_oo - (-q;q^2)_oo)/2
= Sum{k>=0} ( 1/((q;q)_k)^2 - 1/(q^2;q^2)_k ) q^(k^2)/2
using the "q-Pochhammer" notation (a;q)n := Product{k=0..n-1} (1 - a*q^k).
(End)
a(n) = p(n-2) - p(n-8) + p(n-18) - p(n-32) + ... + (-1)^(k+1)*p(n-2*k^2) + ..., where p() is A000041(). E.g., a(20) = p(18) - p(12) + p(2) = 385 - 77 + 2 = 310. - Vladeta Jovovic, Aug 08 2004
G.f.: (1/2)*(1 - Product_{j>=1} (1-x^(2j))/(1+x^(2j)))/Product_{j>=1} (1 - x^j). - Emeric Deutsch, Mar 02 2006
a(2*n) = A236559(n). a(2*n + 1) = A236914(n). - Michael Somos, Aug 25 2015
a(n) = A330644(n)/2. - Omar E. Pol, Jan 10 2020
a(n) = A000041(n) - A046682(n) = A046682(n) - A000700(n). - Gus Wiseman, Mar 31 2022

Extensions

Better description and more terms from Christian G. Bower, Apr 27 2000

A064428 Number of partitions of n with nonnegative crank.

Original entry on oeis.org

1, 0, 1, 2, 3, 4, 6, 8, 12, 16, 23, 30, 42, 54, 73, 94, 124, 158, 206, 260, 334, 420, 532, 664, 835, 1034, 1288, 1588, 1962, 2404, 2953, 3598, 4392, 5328, 6466, 7808, 9432, 11338, 13632, 16326, 19544, 23316, 27806, 33054, 39273, 46534, 55096, 65076, 76808
Offset: 0

Views

Author

Vladeta Jovovic, Sep 30 2001

Keywords

Comments

For a partition p, let l(p) = largest part of p, w(p) = number of 1's in p, m(p) = number of parts of p larger than w(p). The crank of p is given by l(p) if w(p) = 0, otherwise m(p)-w(p).
From Gus Wiseman, Mar 30 2021 and May 21 2022: (Start)
Also the number of even-length compositions of n with alternating parts strictly decreasing, or properly 2-colored partitions (proper = no equal parts of the same color) with the same number of parts of each color, or ordered pairs of strict partitions of the same length with total n. The odd-length case is A001522, and there are a total of A000041 compositions with alternating parts strictly decreasing (see A342528 for a bijective proof). The a(2) = 1 through a(7) = 8 ordered pairs of strict partitions of the same length 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)
(4)(1) (4)(2) (4)(3)
(5)(1) (5)(2)
(21)(21) (6)(1)
(21)(31)
(31)(21)
Conjecture: Also the number of integer partitions y of n without a fixed point y(i) = i, ranked by A352826. This is stated at A238394, but Resta tells me he may not have had a proof. The a(2) = 1 through a(7) = 8 partitions without a fixed point are:
(2) (3) (4) (5) (6) (7)
(21) (31) (41) (33) (43)
(211) (311) (51) (61)
(2111) (411) (331)
(3111) (511)
(21111) (4111)
(31111)
(211111)
The version for permutations is A000166, complement A002467.
The version for compositions is A238351.
This is column k = 0 of A352833.
A238352 counts reversed partitions by fixed points, rank statistic A352822.
A238394 counts reversed partitions without a fixed point, ranked by A352830.
A238395 counts reversed partitions with a fixed point, ranked by A352872. (End)
The above conjecture is true. See Section 4 of the Blecher-Knopfmacher paper in the Links section. - Jeremy Lovejoy, Sep 26 2022

Examples

			G.f. = 1 + x^2 + 2*x^3 + 3*x^4 + 4*x^5 + 6*x^6 + 8*x^7 + 12*x^8 + 16*x^9 + 23*x^10 + ... - _Michael Somos_, Jan 15 2018
From _Gus Wiseman_, May 21 2022: (Start)
The a(0) = 1 through a(8) = 12 partitions with nonnegative crank:
  ()  .  (2)  (3)   (4)   (5)    (6)    (7)     (8)
              (21)  (22)  (32)   (33)   (43)    (44)
                    (31)  (41)   (42)   (52)    (53)
                          (221)  (51)   (61)    (62)
                                 (222)  (322)   (71)
                                 (321)  (331)   (332)
                                        (421)   (422)
                                        (2221)  (431)
                                                (521)
                                                (2222)
                                                (3221)
                                                (3311)
(End)
		

References

  • B. C. Berndt, Ramanujan's Notebooks Part III, Springer-Verlag, see p. 18 Entry 9 Corollary (i).
  • G. E. Andrews, B. C. Berndt, Ramanujan's Lost Notebook Part I, Springer, see p. 169 Entry 6.7.1.

Crossrefs

These are the row-sums of the right (or left) half of A064391, inclusive.
The case of crank 0 is A064410, ranked by A342192.
The strict case is A352828.
These partitions are ranked by A352873.
A000700 = self-conjugate partitions, ranked by A088902, complement A330644.
A001522 counts partitions with positive crank, ranked by A352874.
A034008 counts even-length compositions.
A115720 and A115994 count partitions by their Durfee square.
A224958 counts compositions w/ alternating parts unequal (even: A342532).
A257989 gives the crank of the partition with Heinz number n.
A342527 counts compositions w/ alternating parts equal (even: A065608).
A342528 = compositions w/ alternating parts weakly decr. (even: A114921).

Programs

  • Mathematica
    a[ n_] := If[ n < 0, 0, SeriesCoefficient[ Sum[ (-1)^k x^(k (k + 1)/2) , {k, 0, (Sqrt[1 + 8 n] - 1)/2}] / QPochhammer[ x], {x, 0, n}]]; (* Michael Somos, Jan 15 2018 *)
    a[ n_] := If[ n < 0, 0, SeriesCoefficient[ Sum[  x^(k (k + 1)) / QPochhammer[ x, x, k]^2 , {k, 0, (Sqrt[1 + 4 n] - 1)/2}], {x, 0, n}]]; (* Michael Somos, Jan 15 2018 *)
    ck[y_]:=With[{w=Count[y,1]},If[w==0,If[y=={},0,Max@@y],Count[y,?(#>w&)]-w]];Table[Length[Select[IntegerPartitions[n],ck[#]>=0&]],{n,0,30}] (* _Gus Wiseman, Mar 30 2021 *)
    ici[q_]:=And@@Table[q[[i]]>q[[i+2]],{i,Length[q]-2}];
    Table[Length[Select[Join@@Permutations/@Select[IntegerPartitions[n], EvenQ@*Length],ici]],{n,0,15}] (* Gus Wiseman, Mar 30 2021 *)
  • PARI
    {a(n) = if( n<0, 0, polcoeff( sum(k=0, (sqrtint(1 + 8*n) -1)\2, (-1)^k * x^((k+k^2)/2)) / eta( x + x * O(x^n)), n))}; /* Michael Somos, Jul 28 2003 */

Formula

a(n) = (A000041(n) + A064410(n)) / 2, n>1. - Michael Somos, Jul 28 2003
G.f.: (Sum_{k>=0} (-1)^k * x^(k(k+1)/2)) / (Product_{k>0} 1-x^k). - Michael Somos, Jul 28 2003
G.f.: Sum_{i>=0} x^(i*(i+1)) / (Product_{j=1..i} 1-x^j )^2. - Jon Perry, Jul 18 2004
a(n) ~ exp(Pi*sqrt(2*n/3)) / (8*n*sqrt(3)). - Vaclav Kotesovec, Sep 26 2016
G.f.: (Sum_{i>=0} x^i / (Product_{j=1..i} 1-x^j)^2 ) * (Product_{k>0} 1-x^k). - Li Han, May 23 2020
a(n) = A000041(n) - A001522(n). - Gus Wiseman, Mar 30 2021
a(n) = A064410(n) + A001522(n). - Gus Wiseman, May 21 2022

A114088 Triangle read by rows: T(n,k) is number of partitions of n whose tail below its Durfee square has k parts (n >= 1; 0 <= k <= n-1).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 1, 3, 3, 2, 1, 1, 1, 3, 4, 3, 2, 1, 1, 1, 4, 5, 5, 3, 2, 1, 1, 1, 5, 6, 6, 5, 3, 2, 1, 1, 1, 6, 8, 8, 7, 5, 3, 2, 1, 1, 1, 7, 10, 10, 9, 7, 5, 3, 2, 1, 1, 1, 9, 13, 13, 12, 10, 7, 5, 3, 2, 1, 1, 1, 10, 16, 17, 15, 13, 10, 7, 5, 3, 2, 1, 1, 1, 12, 20, 22, 20, 17
Offset: 1

Views

Author

Emeric Deutsch, Feb 12 2006

Keywords

Comments

From Gus Wiseman, May 21 2022: (Start)
Also the number of integer partitions of n with k parts below the diagonal. For example, the partition (3,2,2,1) has two parts (at positions 3 and 4) below the diagonal (1,2,3,4). Row n = 8 counts the following partitions:
8 71 611 5111 41111 311111 2111111 11111111
44 332 2222 22211 221111
53 422 3221 32111
62 431 3311
521 4211
Indices of parts below the diagonal are also called strong nonexcedances.
(End)

Examples

			T(7,2)=3 because we have [5,1,1], [3,2,1,1] and [2,2,2,1] (the bottom tails are [1,1], [1,1] and [2,1], respectively).
Triangle starts:
  1;
  1, 1;
  1, 1, 1;
  2, 1, 1, 1;
  2, 2, 1, 1, 1;
  3, 3, 2, 1, 1, 1;
  3, 4, 3, 2, 1, 1, 1;
		

References

  • G. E. Andrews, The Theory of Partitions, Addison-Wesley, 1976 (pp. 27-28).
  • G. E. Andrews and K. Eriksson, Integer Partitions, Cambridge Univ. Press, 2004 (pp. 75-78).

Crossrefs

Row sums: A000041.
Column k = 0: A003114.
Weak opposite: A115994.
Permutations: A173018, weak A123125.
Ordered: A352521, rank stat A352514, weak A352522.
Opposite ordered: A352524, first col A008930, rank stat A352516.
Weak opposite ordered: A352525, first col A177510, rank stat A352517.
Weak: A353315.
Opposite: A353318.
A000700 counts self-conjugate partitions, ranked by A088902.
A115720 counts partitions by Durfee square, rank stat A257990.
A352490 gives the (strong) nonexcedance set of A122111, counted by A000701.

Programs

  • Maple
    g:=sum(z^(k^2)/product((1-z^j)*(1-t*z^j),j=1..k),k=1..20): gserz:=simplify(series(g,z=0,30)): for n from 1 to 14 do P[n]:=coeff(gserz,z^n) od: for n from 1 to 14 do seq(coeff(t*P[n],t^j),j=1..n) od; # yields sequence in triangular form
  • Mathematica
    subdiags[y_]:=Length[Select[Range[Length[y]],#>y[[#]]&]];
    Table[Length[Select[IntegerPartitions[n],subdiags[#]==k&]],{n,1,15},{k,0,n-1}] (* Gus Wiseman, May 21 2022 *)
  • PARI
    T_qt(max_row) = {my(N=max_row+1, q='q+O('q^N), h = sum(k=1,N, q^(k^2)/prod(j=1,k, (1-q^j)*(1-t*q^j))) ); for(i=1, N-1, print(Vecrev(polcoef(h, i))))}
    T_qt(10) \\ John Tyler Rascoe, Oct 24 2024

Formula

G.f. = Sum_{k>=1} q^(k^2) / Product_{j=1..k} (1 - q^j)*(1 - t*q^j).
Sum_{k=0..n-1} k*T(n,k) = A114089(n).

A046682 Number of cycle types of conjugacy classes of all even permutations of n elements.

Original entry on oeis.org

1, 1, 1, 2, 3, 4, 6, 8, 12, 16, 22, 29, 40, 52, 69, 90, 118, 151, 195, 248, 317, 400, 505, 632, 793, 985, 1224, 1512, 1867, 2291, 2811, 3431, 4186, 5084, 6168, 7456, 9005, 10836, 13026, 15613, 18692, 22316, 26613, 31659, 37619, 44601, 52815, 62416, 73680, 86809, 102162
Offset: 0

Views

Author

Keywords

Comments

Also number of partitions of n with even number of even parts. There is no restriction on the odd parts.
a(n) = u(n) + v(n), n >= 2, of the Osima reference, p. 383.
Also number of partitions of n with largest part congruent to n modulo 2: a(2*n) = A027187(2*n), a(2*n-1) = A027193(2*n-1); a(n) = A000041(n) - A000701(n). - Reinhard Zumkeller, Apr 22 2006
Equivalently, number of partitions of n with number of parts having the same parity as n. - Olivier Gérard, Apr 04 2012
Also number of distinct free Young diagrams (Ferrers graphs with n nodes). Free Young diagrams are distinct when none is a rigid transformation (translation, rotation, reflection or glide reflection) of another. - Jani Melik, May 08 2016
Let the cycle type of an even permutation be represented by the partition A=(O1,O2,...,Oi,E1,E2,...,E2j), where the Os are parts with odd length and the Es are parts with even lengths, and where j may be zero, using Reinhard Zumkeller's observation that the partition associated with a cycle type of an even permutation has an even number of even parts. The set of even cycle types enumerated here can be considered a monoid under a binary operation *: Let A be as above and B=(o1,o2,...,ok,e1,e2,...,e2m). A*B is the partition (O1o1,O1o2,...,O1ok,O1e1,...,O1e2m,O2o1,...,O2e2m,...,Oio1,...,Oie2m,E1o1,...,E1e2m,...,E2je2m). This product has 2im+2jk+4jm even parts, so it represents the cycle type of an even permutation. - Richard Locke Peterson, Aug 20 2018
From Gus Wiseman, Mar 31 2022: (Start)
Also the number of integer partitions of n with Heinz number greater than or equal to that of their conjugate, where the Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). These partitions are ranked by A352488. The complement is counted by A000701. For example, the a(n) partitions for n = 1...7 are:
(1) (11) (21) (22) (221) (222) (331)
(111) (211) (311) (321) (2221)
(1111) (2111) (2211) (3211)
(11111) (3111) (4111)
(21111) (22111)
(111111) (31111)
(211111)
(1111111)
Also the number of integer partitions of n with Heinz number less than or equal to their conjugate, ranked by A352489. For example, the a(n) partitions for n = 1...7 are:
(1) (2) (3) (4) (5) (6) (7)
(21) (22) (32) (33) (43)
(31) (41) (42) (52)
(311) (51) (61)
(321) (322)
(411) (421)
(511)
(4111)
(End)

Examples

			1 + x + x^2 + 2*x^3 + 3*x^4 + 4*x^5 + 6*x^6 + 8*x^7 + 12*x^8 + 16*x^9 + ...
a(3)=2 since cycle types of even permutations of 3 elements is (.)(.)(.), (...).
a(4)=3 since cycle types of even permutations of 4 elements is (.)(.)(.)(.), (...)(.), (..)(..).
a(5)=4 (free Young diagrams):
  XXXXX XXXX. XXX.. XXX..
  ..... X.... XX... X....
  ..... ..... ..... X....
  ..... ..... ..... .....
  ..... ..... ..... .....
		

Crossrefs

For the number of conjugacy classes of the alternating group A_n, n>=2, see A000702.
Cf. A118301.
A000041 counts integer partitions.
A000700 counts self-conjugate partitions, ranked by A088902.
A330644 counts non-self-conjugate partitions, ranked by A352486.
Heinz number (rank) and partition:
- A122111 = rank of conjugate.
- A296150 = parts of partition, conjugate A321649.
- A352487 = rank less than conjugate, counted by A000701.
- A352488 = rank greater than or equal to conjugate, counted by A046682.
- A352489 = rank less than or equal to conjugate, counted by A046682.
- A352490 = rank greater than conjugate, counted by A000701.
- A352491 = rank minus conjugate.

Programs

  • Maple
    seq(add((-1)^(n-k)*combinat:-numbpart(n,k),k=0..n),n=0..48); # Peter Luschny, Aug 03 2015
  • Mathematica
    max = 48; f[q_] := Sum[(-q^2)^n^2, {n, 0, max}]/Product[1-q^n, {n, 1, max}]; CoefficientList[ Series[f[q], {q, 0, max}], q] (* Jean-François Alcover, Oct 18 2011, after g.f. *)
    conj[y_]:=If[Length[y]==0,y,Table[Length[Select[y,#>=k&]],{k,1,Max[y]}]];
    Table[Length[Select[IntegerPartitions[n],Times@@Prime/@#>=Times@@Prime/@conj[#]&]],{n,0,15}] (* Gus Wiseman, Mar 31 2022 *)
  • PARI
    list(lim)=my(q='q);Vec(sum(n=0,sqrt(lim),(-q^2)^(n^2))/prod(n=1,lim,1-q^n)+O(q^(lim\1+1))) \\ Charles R Greathouse IV, Oct 18 2011
    
  • PARI
    {a(n) = if( n<0, 0, (numbpart(n) + polcoeff( 1 / prod( k=1, n, 1 + (-x)^k, 1 + x * O(x^n)), n)) / 2)} /* Michael Somos, Jul 24 2012 */

Formula

G.f.: Sum_{n>=0} (-q^2)^(n^2) / Product_{m>=1} (1-q^m ) = ( 1/Product_{m>=1} (1-q^m) + Product_{m>=1} (1+q^(2*m-1) ) ) / 2. - Mamuka Jibladze, Sep 07 2003
a(n) = (A000041(n) + A000700(n)) / 2.
a(n) = A000041(n) - A000701(n). - Gus Wiseman, Mar 31 2022

A325169 Origin-to-boundary graph-distance of the Young diagram of the integer partition with Heinz number n.

Original entry on oeis.org

0, 1, 1, 1, 1, 2, 1, 1, 2, 2, 1, 2, 1, 2, 2, 1, 1, 2, 1, 2, 2, 2, 1, 2, 2, 2, 2, 2, 1, 3, 1, 1, 2, 2, 2, 2, 1, 2, 2, 2, 1, 3, 1, 2, 3, 2, 1, 2, 2, 3, 2, 2, 1, 2, 2, 2, 2, 2, 1, 3, 1, 2, 3, 1, 2, 3, 1, 2, 2, 3, 1, 2, 1, 2, 3, 2, 2, 3, 1, 2, 2, 2, 1, 3, 2, 2, 2
Offset: 1

Views

Author

Gus Wiseman, Apr 05 2019

Keywords

Comments

The origin-to-boundary graph-distance of a Young diagram is the minimum number of unit steps left or down from the upper-left square to a nonsquare in the lower-right quadrant. It is also the side-length of the minimum triangular partition contained inside the diagram.
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    otb[ptn_]:=Min@@MapIndexed[#1+#2[[1]]-1&,Append[ptn,0]];
    Table[otb[Reverse[primeMS[n]]],{n,100}]

Formula

A257990(n) <= a(n) <= 2 * A257990(n).
Showing 1-10 of 82 results. Next