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

A003114 Number of partitions of n into parts 5k+1 or 5k+4.

Original entry on oeis.org

1, 1, 1, 1, 2, 2, 3, 3, 4, 5, 6, 7, 9, 10, 12, 14, 17, 19, 23, 26, 31, 35, 41, 46, 54, 61, 70, 79, 91, 102, 117, 131, 149, 167, 189, 211, 239, 266, 299, 333, 374, 415, 465, 515, 575, 637, 709, 783, 871, 961, 1065, 1174, 1299, 1429, 1579, 1735, 1913, 2100, 2311, 2533, 2785
Offset: 0

Views

Author

Keywords

Comments

Expansion of Rogers-Ramanujan function G(x) in powers of x.
Same as number of partitions into distinct parts where the difference between successive parts is >= 2.
As a formal power series, the limit of polynomials S(n,x): S(n,x)=sum(T(i,x),0<=i<=n); T(i,x)=S(i-2,x).x^i; T(0,x)=1,T(1,x)=x; S(n,1)=A000045(n+1), the Fibonacci sequence. - Claude Lenormand (claude.lenormand(AT)free.fr), Feb 04 2001
The Rogers-Ramanujan identity is 1 + Sum_{n >= 1} t^(n^2)/((1-t)*(1-t^2)*...*(1-t^n)) = Product_{n >= 1} 1/((1-t^(5*n-1))*(1-t^(5*n-4))).
Coefficients in expansion of permanent of infinite tridiagonal matrix:
1 1 0 0 0 0 0 0 ...
x 1 1 0 0 0 0 0 ...
0 x^2 1 1 0 0 0 ...
0 0 x^3 1 1 0 0 ...
0 0 0 x^4 1 1 0 ...
................... - Vladeta Jovovic, Jul 17 2004
Also number of partitions of n such that the smallest part is greater than or equal to number of parts. - Vladeta Jovovic, Jul 17 2004
Also number of partitions of n such that if k is the largest part, then each of {1, 2, ..., k-1} occur at least twice. Example: a(9)=5 because we have [3, 2, 2, 1, 1], [2, 2, 2, 1, 1, 1], [2, 2, 1, 1, 1, 1, 1], [2, 1, 1, 1, 1, 1, 1, 1] and [1, 1, 1, 1, 1, 1, 1, 1, 1]. - Emeric Deutsch, Feb 27 2006
Also number of partitions of n such that if k is the largest part, then k occurs at least k times. Example: a(9)=5 because we have [3, 3, 3], [2, 2, 2, 2, 1], [2, 2, 2, 1, 1, 1], [2, 2, 1, 1, 1, 1, 1] and [1, 1, 1, 1, 1, 1, 1, 1, 1]. - Emeric Deutsch, Apr 16 2006
a(n) = number of NW partitions of n, for n >= 1; see A237981.
For more about the generalized Rogers-Ramanujan series G[i](x) see the Andrews-Baxter and Lepowsky-Zhu papers. The present series is G[1](x). - N. J. A. Sloane, Nov 22 2015
Convolution of A109700 and A109697. - Vaclav Kotesovec, Jan 21 2017

Examples

			G.f. = 1 + x + x^2 + x^3 + 2*x^4 + 2*x^5 + 3*x^6 + 3*x^7 + 4*x^8 + 5*x^9 + ...
G.f. = 1/q + q^59 + q^119 + q^179 + 2*q^239 + 2*q^299 + 3*q^359 + 3*q^419 + ...
From _Joerg Arndt_, Dec 27 2012: (Start)
The a(16)=17 partitions of 16 where all parts are 1 or 4 (mod 5) are
  [ 1]  [ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ]
  [ 2]  [ 4 1 1 1 1 1 1 1 1 1 1 1 1 ]
  [ 3]  [ 4 4 1 1 1 1 1 1 1 1 ]
  [ 4]  [ 4 4 4 1 1 1 1 ]
  [ 5]  [ 4 4 4 4 ]
  [ 6]  [ 6 1 1 1 1 1 1 1 1 1 1 ]
  [ 7]  [ 6 4 1 1 1 1 1 1 ]
  [ 8]  [ 6 4 4 1 1 ]
  [ 9]  [ 6 6 1 1 1 1 ]
  [10]  [ 6 6 4 ]
  [11]  [ 9 1 1 1 1 1 1 1 ]
  [12]  [ 9 4 1 1 1 ]
  [13]  [ 9 6 1 ]
  [14]  [ 11 1 1 1 1 1 ]
  [15]  [ 11 4 1 ]
  [16]  [ 14 1 1 ]
  [17]  [ 16 ]
The a(16)=17 partitions of 16 where successive parts differ by at least 2 are
  [ 1]  [ 7 5 3 1 ]
  [ 2]  [ 8 5 3 ]
  [ 3]  [ 8 6 2 ]
  [ 4]  [ 9 5 2 ]
  [ 5]  [ 9 6 1 ]
  [ 6]  [ 9 7 ]
  [ 7]  [ 10 4 2 ]
  [ 8]  [ 10 5 1 ]
  [ 9]  [ 10 6 ]
  [10]  [ 11 4 1 ]
  [11]  [ 11 5 ]
  [12]  [ 12 3 1 ]
  [13]  [ 12 4 ]
  [14]  [ 13 3 ]
  [15]  [ 14 2 ]
  [16]  [ 15 1 ]
  [17]  [ 16 ]
(End)
		

References

  • G. E. Andrews, The Theory of Partitions, Addison-Wesley, 1976, p. 109, 238.
  • G. E. Andrews, R. Askey and R. Roy, Special Functions, Cambridge University Press, 1999; Exercise 6(e), p. 591.
  • Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 669.
  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 107.
  • G. H. Hardy, Ramanujan, AMS Chelsea Publ., Providence, RI, 2002, pp. 90-92.
  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, Fifth ed., Clarendon Press, Oxford, 2003, pp. 290-291.
  • H. P. Robinson, Letter to N. J. A. Sloane, Jan 04 1974.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A188216 (least part k occurs at least k times).
For the generalized Rogers-Ramanujan series G[1], G[2], G[3], G[4], G[5], G[6], G[7], G[8] see A003114, A003106, A006141, A264591, A264592, A264593, A264594, A264595. G[0] = G[1]+G[2] is given by A003113.
Row sums of A268187.

Programs

  • Haskell
    a003114 = p a047209_list where
       p _      0 = 1
       p ks'@(k:ks) m = if m < k then 0 else p ks' (m - k) + p ks m
    -- Reinhard Zumkeller, Jan 05 2011
    
  • Haskell
    a003114 = p 1 where
       p _ 0 = 1
       p k m = if k > m then 0 else p (k + 2) (m - k) + p (k + 1) m
    -- Reinhard Zumkeller, Feb 19 2013
  • Maple
    g:=sum(x^(k^2)/product(1-x^j,j=1..k),k=0..10): gser:=series(g,x=0,65): seq(coeff(gser,x,n),n=0..60); # Emeric Deutsch, Feb 27 2006
  • Mathematica
    CoefficientList[ Series[Sum[x^k^2/Product[1 - x^j, {j, 1, k}], {k, 0, 10}], {x, 0, 65}], x][[1 ;; 61]] (* Jean-François Alcover, Apr 08 2011, after Emeric Deutsch *)
    Table[Count[IntegerPartitions[n], p_ /; Min[p] >= Length[p]], {n, 0, 24}] (* Clark Kimberling, Feb 13 2014 *)
    a[ n_] := SeriesCoefficient[ 1 / (QPochhammer[ x^1, x^5] QPochhammer[ x^4, x^5]), {x, 0, n}]; (* Michael Somos, May 17 2015 *)
    a[ n_] := SeriesCoefficient[ Product[ (1 - x^k)^{-1, 0, 0, -1, 0}[[Mod[k, 5, 1]]], {k, n}], {x, 0, n}]; (* Michael Somos, May 17 2015 *)
    nmax = 60; kmax = nmax/5;
    s = Flatten[{Range[0, kmax]*5 + 1}~Join~{Range[0, kmax]*5 + 4}];
    Table[Count[IntegerPartitions@n, x_ /; SubsetQ[s, x]], {n, 0, nmax}] (* Robert Price, Aug 02 2020 *)
  • PARI
    {a(n) = my(t); if( n<0, 0, t = 1 + x * O(x^n); polcoeff( sum(k=1, sqrtint(n), t *= x^(2*k - 1) / (1 - x^k) * (1 + x * O(x^(n - k^2))), 1), n))}; /* Michael Somos, Oct 15 2008 */
    

Formula

G.f.: Sum_{k>=0} x^(k^2)/(Product_{i=1..k} 1-x^i).
The g.f. above is the special case D=2 of sum(n>=0, x^(D*n*(n+1)/2 - (D-1)*n) / prod(k=1..n, 1-x^k) ), the g.f. for partitions into distinct part where the difference between successive parts is >= D. - Joerg Arndt, Mar 31 2014
G.f.: 1 + sum(i=1, oo, x^(5i+1)/prod(j=1 or 4 mod 5 and j<=5i+1, 1-x^j) + x^(5i+4)/prod(j=1 or 4 mod 5 and j<=5i+4, 1-x^j)). - Jon Perry, Jul 06 2004
G.f.: (Product_{k>0} 1 + x^(2*k)) * (Sum_{k>=0} x^(k^2) / (Product_{i=1..k} 1 - x^(4*i))). - Michael Somos, Oct 19 2006
Euler transform of period 5 sequence [ 1, 0, 0, 1, 0, ...]. - Michael Somos, Oct 15 2008
Expansion of f(-x^5) / f(-x^1, -x^4) in powers of x where f(,) is the Ramanujan general theta function. - Michael Somos, May 17 2015
Expansion of f(-x^2, -x^3) / f(-x) in powers of x where f(,) is the Ramanujan general theta function. - Michael Somos, Jun 13 2015
a(n) ~ phi^(1/2) * exp(2*Pi*sqrt(n/15)) / (2 * 3^(1/4) * 5^(1/2) * n^(3/4)) * (1 - (3*sqrt(15)/(16*Pi) + Pi/(60*sqrt(15))) / sqrt(n)), where phi = A001622 = (1+sqrt(5))/2 is the golden ratio. - Vaclav Kotesovec, Aug 23 2015, extended Jan 24 2017
a(n) = (1/n)*Sum_{k=1..n} A284150(k)*a(n-k), a(0) = 1. - Seiichi Manyama, Mar 21 2017

A047993 Number of balanced partitions of n: the largest part equals the number of parts.

Original entry on oeis.org

1, 0, 1, 1, 1, 1, 3, 2, 4, 4, 6, 7, 11, 11, 16, 19, 25, 29, 40, 45, 60, 70, 89, 105, 134, 156, 196, 232, 285, 336, 414, 485, 591, 696, 839, 987, 1187, 1389, 1661, 1946, 2311, 2702, 3201, 3731, 4400, 5126, 6018, 6997, 8195, 9502, 11093, 12849, 14949, 17281, 20062
Offset: 1

Views

Author

Keywords

Comments

Useful in the creation of plane partitions with C3 or C3v symmetry.
The function T[m,a,b] used here gives the partitions of m whose Ferrers plot fits within an a X b box.
Central terms of triangle in A063995: a(n) = A063995(n,0). - Reinhard Zumkeller, Jul 24 2013
Sequence enumerates the collection of partitions of size n that are in the monoid of Dyson rank=0, or balanced partitions, under the binary operation A*B = (a1,a2,...,a[k-1],k)*(b1,...,b[n-1,n) = (a1*b1,...,a1*n,a2*b1,...,a2*n,...,k*b1,...,k*n), where A is a partition with k parts and B is a partition with n parts, and A*B is a partition with k*n parts. Note that the rank of A*B is 0, as required. For example, the product of the rank 0 partitions (1,2,3) of 6 and (1,1,3) of 5 is the rank 0 partition (1,1,2,2,3,3,3,6,9) of 30. There is no rank zero partition of 2, as shown in the sequence. It can be seen that any element of the monoid that partitions an odd prime p or a composite number of form 2p cannot be a product of smaller nontrivial partitions, whether in this monoid or not. - Richard Locke Peterson, Jul 15 2018
The "multiplication" given above was noted earlier by Franklin T. Adams-Watters in A122697. - Richard Peterson, Jul 19 2023
The Heinz numbers of these integer partitions are given by A106529. - Gus Wiseman, Mar 09 2019

Examples

			From _Joerg Arndt_, Oct 08 2012: (Start)
a(12) = 7 because the partitions of 12 where the largest part equals the number of parts are
   2 + 3 + 3 + 4,
   2 + 2 + 4 + 4,
   1 + 3 + 4 + 4,
   1 + 2 + 2 + 2 + 5,
   1 + 1 + 2 + 3 + 5,
   1 + 1 + 1 + 4 + 5, and
   1 + 1 + 1 + 1 + 2 + 6.
(End)
From _Gus Wiseman_, Mar 09 2019: (Start)
The a(1) = 1 through a(13) = 11 integer partitions:
  1  21  22  311  321  322   332   333    4222   4322    4332    4333
                       331   4211  4221   4321   4331    4422    4432
                       4111        4311   4411   4421    4431    4441
                                   51111  52111  52211   52221   52222
                                                 53111   53211   53221
                                                 611111  54111   53311
                                                         621111  54211
                                                                 55111
                                                                 622111
                                                                 631111
                                                                 7111111
(End)
		

Crossrefs

Programs

  • Haskell
    a047993 = flip a063995 0  -- Reinhard Zumkeller, Jul 24 2013
  • Maple
    A047993 := proc(n)
         a := 0 ;
         for p in combinat[partition](n) do
            r := max(op(p))-nops(p) ;
            if r = 0 then
                 a := a+1 ;
            end if;
         end do:
         a ;
     end proc:
    seq(A047993(n),n=1..20) ; # Emeric Deutsch, Dec 11 2004
  • Mathematica
    Table[ Count[Partitions[n], par_List/; First[par]===Length[par]], {n, 12}] or recur: Sum[T[n-(2m-1), m-1, m-1], {m, Ceiling[Sqrt[n]], Floor[(n+1)/2]}] with T[m_, a_, b_]/; b < a := T[m, b, a]; T[m_, a_, b_]/; m > a*b := 0; T[m_, a_, b_]/; (2m > a*b) := T[a*b-m, a, b]; T[m_, 1, b_] := If[b < m, 0, 1]; T[0, , ] := 1; T[m_, a_, b_] := T[m, a, b]=Sum[T[m-a*i, a-1, b-i], {i, 0, Floor[m/a]}];
    Table[Sum[ -(-1)^k*(p[n-(3*k^2-k)/2] - p[n-(3*k^2+k)/2]), {k, 1, Floor[(1+Sqrt[1+24*n])/6]}] /. p -> PartitionsP, {n, 1, 64}] (* Wouter Meeussen *)
    (* also *)
    Table[Count[IntegerPartitions[n], q_ /; Max[q] == Length[q]], {n, 24}]
    (* Clark Kimberling, Feb 13 2014 *)
    nmax = 100; p = 1; s = 1; Do[p = Normal[Series[p*x^2*(1 - x^(2*k - 1))*(1 + x^k)/(1 - x^k), {x, 0, nmax}]]; s += p;, {k, 1, nmax + 1}]; Take[CoefficientList[s, x], nmax] (* Vaclav Kotesovec, Oct 16 2024 *)
  • PARI
    N=66;  q='q + O('q^N );
    S=2+2*ceil(sqrt(N));
    gf= sum(k=1, S,  (-1)^k * ( q^((3*k^2+k)/2) - q^((3*k^2-k)/2) ) ) / prod(k=1,N, 1-q^k );
    /* Joerg Arndt, Oct 08 2012 */
    
  • PARI
    my(N=66, x='x+O('x^N)); Vec(sum(k=1, N, x^(2*k-1)*prod(j=1, k-1, (1-x^(k+j-1))/(1-x^j)))) \\ Seiichi Manyama, Jan 24 2022
    

Formula

a(n) = p(n-1) - p(n-2) - p(n-5) + p(n-7) + ... + (-1)^k*(p(n-(3*k^2-k)/2) - p(n-(3*k^2+k)/2)) + ..., where p() is A000041(). E.g., a(20) = p(19) - p(18) - p(15) + p(13) + p(8) - p(5) = 490 - 385 - 176 + 101 + 22 - 7 = 45. - Vladeta Jovovic, Aug 04 2004
G.f.: ( Sum_{k>=1} (-1)^k * ( x^((3*k^2+k)/2) - x^((3*k^2-k)/2) ) ) / Product_{k>=1} (1-x^k). - Vladeta Jovovic, Aug 05 2004
a(n) ~ exp(Pi*sqrt(2*n/3))*Pi / (48*sqrt(2)*n^(3/2)) ~ p(n) * Pi / (4*sqrt(6*n)), where p(n) is the partition function A000041. - Vaclav Kotesovec, Oct 06 2016
G.f.: Sum_{k>=1} x^(2*k-1) * Product_{j=1..k-1} (1-x^(k+j-1))/(1-x^j). - Seiichi Manyama, Jan 24 2022

A001156 Number of partitions of n into squares.

Original entry on oeis.org

1, 1, 1, 1, 2, 2, 2, 2, 3, 4, 4, 4, 5, 6, 6, 6, 8, 9, 10, 10, 12, 13, 14, 14, 16, 19, 20, 21, 23, 26, 27, 28, 31, 34, 37, 38, 43, 46, 49, 50, 55, 60, 63, 66, 71, 78, 81, 84, 90, 98, 104, 107, 116, 124, 132, 135, 144, 154, 163, 169, 178, 192, 201, 209, 220, 235, 247, 256
Offset: 0

Views

Author

Keywords

Comments

Number of partitions of n such that number of parts equal to k is multiple of k for all k. - Vladeta Jovovic, Aug 01 2004
Of course p_{4*square}(n)>0. In fact p_{4*square}(32n+28)=3 times p_{4*square}(8n+7) and p_{4*square}(72n+69) is even. These seem to be the only arithmetic properties the function p_{4*square(n)} possesses. Similar results hold for partitions into positive squares, distinct squares and distinct positive squares. - Michael David Hirschhorn, May 05 2005
The Heinz numbers of these partitions are given by A324588. - Gus Wiseman, Mar 09 2019

Examples

			p_{4*square}(23)=1 because 23 = 3^2 + 3^2 + 2^2 + 1^2 and there is no other partition of 23 into squares.
G.f.: A(x) = 1 + x + x^2 + x^3 + 2*x^4 + 2*x^5 + 2*x^6 + 2*x^7 +...
such that the g.f. A(x) satisfies the identity [_Paul D. Hanna_]:
A(x) = 1/((1-x)*(1-x^4)*(1-x^9)*(1-x^16)*(1-x^25)*...)
A(x) = 1 + x/(1-x) + x^4/((1-x)*(1-x^4)) + x^9/((1-x)*(1-x^4)*(1-x^9)) + x^16/((1-x)*(1-x^4)*(1-x^9)*(1-x^16)) + ...
From _Gus Wiseman_, Mar 09 2019: (Start)
The a(14) = 6 integer partitions into squares are:
  (941)
  (911111)
  (44411)
  (44111111)
  (41111111111)
  (11111111111111)
while the a(14) = 6 integer partitions in which the multiplicity of k is a multiple of k for all k are:
  (333221)
  (33311111)
  (22222211)
  (2222111111)
  (221111111111)
  (11111111111111)
(End)
		

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

Cf. A000041, A000161 (partitions into 2 squares), A000290, A033461, A131799, A218494, A285218, A304046.
Cf. A078134 (first differences).
Row sums of A243148.
Euler trans. of A010052 (see also A308297).

Programs

  • Haskell
    a001156 = p (tail a000290_list) where
       p _          0 = 1
       p ks'@(k:ks) m = if m < k then 0 else p ks' (m - k) + p ks m
    -- Reinhard Zumkeller, Oct 31 2012, Aug 14 2011
    
  • Magma
    m:=70; R:=PowerSeriesRing(Integers(), m); Coefficients(R!( (&*[1/(1-x^(k^2)): k in [1..(m+2)]]) )); // G. C. Greubel, Nov 11 2018
  • Maple
    b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
          b(n, i-1)+ `if`(i^2>n, 0, b(n-i^2, i))))
        end:
    a:= n-> b(n, isqrt(n)):
    seq(a(n), n=0..120);  # Alois P. Heinz, May 30 2014
  • Mathematica
    CoefficientList[ Series[Product[1/(1 - x^(m^2)), {m, 70}], {x, 0, 68}], x] (* Or *)
    Join[{1}, Table[Length@PowersRepresentations[n, n, 2], {n, 68}]] (* Robert G. Wilson v, Apr 12 2005, revised Sep 27 2011 *)
    f[n_] := Length@ IntegerPartitions[n, All, Range@ Sqrt@ n^2]; Array[f, 67] (* Robert G. Wilson v, Apr 14 2013 *)
    b[n_, i_] := b[n, i] = If[n==0, 1, If[i<1, 0, b[n, i-1] + If[i^2>n, 0, b[n-i^2, i]]]]; a[n_] := b[n, Sqrt[n]//Floor]; Table[a[n], {n, 0, 120}] (* Jean-François Alcover, Nov 02 2015, after Alois P. Heinz *)
  • PARI
    {a(n)=polcoeff(1/prod(k=1, sqrtint(n+1), 1-x^(k^2)+x*O(x^n)), n)} \\ Paul D. Hanna, Mar 09 2012
    
  • PARI
    {a(n)=polcoeff(1+sum(m=1, sqrtint(n+1), x^(m^2)/prod(k=1, m, 1-x^(k^2)+x*O(x^n))), n)} \\ Paul D. Hanna, Mar 09 2012
    

Formula

G.f.: Product_{m>=1} 1/(1-x^(m^2)).
G.f.: Sum_{n>=0} x^(n^2) / Product_{k=1..n} (1 - x^(k^2)). - Paul D. Hanna, Mar 09 2012
a(n) = (1/n)*Sum_{k=1..n} A035316(k)*a(n-k). - Vladeta Jovovic, Nov 20 2002
a(n) = f(n,1,3) with f(x,y,z) = if xReinhard Zumkeller, Nov 08 2009
Conjecture (Jan Bohman, Carl-Erik Fröberg, Hans Riesel, 1979): a(n) ~ c * n^(-alfa) * exp(beta*n^(1/3)), where c = 1/18.79656, beta = 3.30716, alfa = 1.16022. - Vaclav Kotesovec, Aug 19 2015
From Vaclav Kotesovec, Dec 29 2016: (Start)
Correct values of these constants are:
1/c = sqrt(3) * (4*Pi)^(7/6) / Zeta(3/2)^(2/3) = 17.49638865935104978665...
alfa = 7/6 = 1.16666666666666666...
beta = 3/2 * (Pi/2)^(1/3) * Zeta(3/2)^(2/3) = 3.307411783596651987...
a(n) ~ 3^(-1/2) * (4*Pi*n)^(-7/6) * Zeta(3/2)^(2/3) * exp(2^(-4/3) * 3 * Pi^(1/3) * Zeta(3/2)^(2/3) * n^(1/3)). [Hardy & Ramanujan, 1917]
(End)

Extensions

More terms from Eric W. Weisstein
More terms from Gh. Niculescu (ghniculescu(AT)yahoo.com), Oct 08 2006

A006141 Number of integer partitions of n whose smallest part is equal to the number of parts.

Original entry on oeis.org

1, 0, 0, 1, 1, 1, 1, 1, 2, 2, 3, 3, 4, 4, 5, 6, 7, 8, 10, 11, 13, 15, 17, 19, 23, 25, 29, 33, 38, 42, 49, 54, 62, 69, 78, 87, 99, 109, 123, 137, 154, 170, 191, 211, 236, 261, 290, 320, 357, 392, 435, 479, 530, 582, 644, 706, 779, 854, 940, 1029, 1133, 1237, 1358, 1485
Offset: 1

Views

Author

Keywords

Comments

Or, number of partitions of n in which number of largest parts is equal to the largest part.
a(n) is the number of partitions of n-1 without parts that differ by less than 2 and which have no parts less than three. [MacMahon]
There are two conflicting choices for the offset in this sequence. For the definition given here the offset is 1, and that is what we shall adopt. On the other hand, if one arrives at this sequence via the Rogers-Ramanujan identities (see the next comment), the natural offset is 0.
Related to Rogers-Ramanujan identities: Let G[1](q) and G[2](q) be the generating functions for the two Rogers-Ramanujan identities of A003114 and A003106, starting with the constant term 1. The g.f. for the present sequence is G[3](q) = (G[1](q) - G[2](q))/q = 1+q^3+q^4+q^5+q^6+q^7+2*q^8+2*q^9+3*q^10+.... - Joerg Arndt, Oct 08 2012; N. J. A. Sloane, Nov 18 2015
For more about the generalized Rogers-Ramanujan series G[i](x) see the Andrews-Baxter and Lepowsky-Zhu papers. The present series is G[3](x). - N. J. A. Sloane, Nov 22 2015
From Wolfdieter Lang, Oct 31 2016: (Start)
From Hardy (H) p. 94, eq. (6.12.1) and Hardy-Wright (H-W), p. 293, eq. (19.14.3) for H_2(a,x) - H_1(a,x) = a*H_1(a*x,x) one finds from the result for H_1(a,x) (in (H) on top on p. 95), after putting a=x, the o.g.f. of a(n) = A003114(n) - A003106(n), n >= 0, with a(0) = 0 as Sum_{m>=0} x^((m+1)^2) / Product_{j=1..m} (1 - x^j). The m=0 term is 1*x^1. See the formula given by Joerg Arndt, Jan 29 2011.
This formula has a combinatorial interpretation (found similar to the one given in (H) section 6.0, pp. 91-92 or (H-W) pp. 290-291): a(n) is the number of partitions of n with parts differing by at least 2 and part 1 present. See the example for a(15) below. (End)
The Heinz numbers of these integer partitions are given by A324522. - Gus Wiseman, Mar 09 2019

Examples

			G.f. = x + x^4 + x^5 + x^6 + x^7 + x^8 + 2*x^9 + 2*x^10 + 3*x^11 + 3*x^12 + ...
a(15) = 5 because the partitions of 15 where the smallest part equals the number of parts are
3 + 6 + 6,
3 + 5 + 7,
3 + 4 + 8,
3 + 3 + 9, and
2 + 13.
- _Joerg Arndt_, Oct 08 2012
a(15) = 5 because the partitions of 15 with parts differing by at least 2 and part 1 present are: [14,1] obtained from the partition of 11 with one part, [11], added to the first part of the special partition [3,1] of 4 and  [11,3,1], [10,4,1], [9,5,1], [8,6,1] from adding all partition of 15 - 9 = 6 with one part, [6], and those with two parts, [5,1], [4,1], [3,3], to the special partition [5,3,1] of 9. - _Wolfdieter Lang_, Oct 31 2016
a(15) = 5 because the partitions of 14 with parts >= 3 and parts differing by at least 2 are [14], [11,3], [10,4], [9,5] and [8,6]. See the second [MacMahon] comment. This follows from the g.f. G[3](q) given in Andrews - Baxter, eq. (5.1) for i=3, (using summation index  m) and  m*(m+2) = 3 + 5 + ... + (2*m+1). - _Wolfdieter Lang_, Nov 02 2016
From _Gus Wiseman_, Mar 09 2019: (Start)
The a(8) = 1 through a(15) = 5 integer partitions:
  (6,2)  (7,2)    (8,2)    (9,2)    (10,2)   (11,2)   (12,2)   (13,2)
         (3,3,3)  (4,3,3)  (4,4,3)  (5,4,3)  (5,5,3)  (6,5,3)  (6,6,3)
                           (5,3,3)  (6,3,3)  (6,4,3)  (7,4,3)  (7,5,3)
                                             (7,3,3)  (8,3,3)  (8,4,3)
                                                               (9,3,3)
(End)
		

References

  • G. H. Hardy, Ramanujan, AMS Chelsea Publ., Providence, RI, 2002, pp. 92-95.
  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, Fifth ed., Clarendon Press, Oxford, 2003, pp. 292-294.
  • P. A. MacMahon, Combinatory Analysis, Cambridge Univ. Press, London and New York, Vol. 1, 1915 and Vol. 2, 1916; see vol. 2, p 45, Section 293.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

For the generalized Rogers-Ramanujan series G[1], G[2], G[3], G[4], G[5], G[6], G[7], G[8] see A003114, A003106, A006141, A264591, A264592, A264593, A264594, A264595. G[0] = G[1]+G[2] is given by A003113.
A003106 counts partitions with minimum > length.
A003114 counts partitions with minimum >= length.
A026794 counts partitions by minimum.
A039899 counts partitions with minimum < length.
A039900 counts partitions with minimum <= length.
A239950 counts partitions with minimum equal to number of distinct parts.
Sequences related to balance:
- A010054 counts balanced strict partitions.
- A047993 counts balanced partitions.
- A098124 counts balanced compositions.
- A106529 ranks balanced partitions.
- A340596 counts co-balanced factorizations.
- A340598 counts balanced set partitions.
- A340599 counts alt-balanced factorizations.
- A340600 counts unlabeled balanced multiset partitions.
- A340653 counts balanced factorizations.

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n<0, 0, `if`(n=0, 1,
          `if`(i<1, 0, b(n, i-1)+`if`(i>n, 0, b(n-i,i)))))
        end:
    a:= n-> add(b(n-j^2, j-1), j=0..isqrt(n)):
    seq(a(n), n=1..80);  # Alois P. Heinz, Oct 08 2012
  • Mathematica
    b[n_, i_] := b[n, i] = If[n<0, 0, If[n == 0, 1, If[i<1, 0, b[n, i-1] + If[i>n, 0, b[n-i, i]]]]]; a[n_] := Sum[b[n-j^2, j-1], {j, 0, Sqrt[n]}]; Table[a[n], {n, 1, 80}] (* Jean-François Alcover, Mar 17 2014, after Alois P. Heinz *)
    Table[Length[Select[IntegerPartitions[n],Min[#]==Length[#]&]],{n,30}] (* Gus Wiseman, Mar 09 2019 *)
  • PARI
    {a(n) = if( n<1, 0, polcoeff( sum(k=1, sqrtint(n), x^k^2 / prod(j=1, k-1, 1 - x^j, 1 + O(x ^ (n - k^2 + 1) ))), n))} /* Michael Somos, Jan 22 2008 */

Formula

G.f.: Sum_{m>=1} (x^(m^2)-x^(m*(m+1))) / Product_{i=1..m} (1-x^i) .
G.f.: Sum_{n>=1} x^(n^2)/Product_{k=1..n-1} (1-x^k). - Joerg Arndt, Jan 29 2011
a(n) = A003114(n) - A003106(n) = A039900(n) - A039899(n), (offset 1). - Vladeta Jovovic, Jul 17 2004
Plouffe in his 1992 dissertation conjectured that this has g.f. = (1+z+z^4+2*z^5-z^3-z^8+3*z^10-z^7+z^9)/(1+z-z^4-2*z^3-z^8+z^10), but Michael Somos pointed out on Jan 22 2008 that this is false.
Expansion of ( f(-x^2, -x^3) - f(-x, -x^4) ) / f(-x) in powers of x where f(, ) is Ramanujan's general theta function. - Michael Somos, Jan 22 2007
a(n) ~ sqrt(1/sqrt(5) - 2/5) * exp(2*Pi*sqrt(n/15)) / (2*3^(1/4)*n^(3/4)). - Vaclav Kotesovec, Nov 01 2016

Extensions

More terms from Kok Seng Chua (chuaks(AT)ihpc.nus.edu.sg), Jun 20 2000
Better description from Naohiro Nomoto, Feb 06 2002
Name shortened by Gus Wiseman, Apr 07 2021 (balanced partitions are A047993).

A064174 Number of partitions of n with nonnegative rank.

Original entry on oeis.org

1, 1, 2, 3, 4, 6, 9, 12, 17, 23, 31, 42, 56, 73, 96, 125, 161, 207, 265, 336, 426, 536, 672, 840, 1046, 1296, 1603, 1975, 2425, 2970, 3628, 4417, 5367, 6503, 7861, 9482, 11412, 13702, 16423, 19642, 23447, 27938, 33231, 39453, 46767, 55342, 65386, 77135
Offset: 1

Views

Author

Vladeta Jovovic, Sep 20 2001

Keywords

Comments

The rank of a partition is the largest summand minus the number of summands.
This sequence (up to proof) equals "partitions of 2n with even number of parts, ending in 1, with max descent of 1, where the number of odd parts in odd places equals the number of odd parts in even places. (See link and 2nd Mathematica line.) - Wouter Meeussen, Mar 29 2013
Number of partitions p of n such that max(max(p), number of parts of p) is a part of p. - Clark Kimberling, Feb 28 2014
From Gus Wiseman, Mar 09 2019: (Start)
Also the number of integer partitions of n with maximum part greater than or equal to the number of parts. The Heinz numbers of these integer partitions are given by A324521. For example, the a(1) = 1 through a(8) = 12 partitions are:
(1) (2) (3) (4) (5) (6) (7) (8)
(21) (22) (32) (33) (43) (44)
(31) (41) (42) (52) (53)
(311) (51) (61) (62)
(321) (322) (71)
(411) (331) (332)
(421) (422)
(511) (431)
(4111) (521)
(611)
(4211)
(5111)
Also the number of integer partitions of n with maximum part less than or equal to the number of parts. The Heinz numbers of these integer partitions are given by A324562. For example, the a(1) = 1 through a(8) = 12 partitions are:
(1) (11) (21) (22) (221) (222) (322) (332)
(111) (211) (311) (321) (331) (2222)
(1111) (2111) (2211) (2221) (3221)
(11111) (3111) (3211) (3311)
(21111) (4111) (4211)
(111111) (22111) (22211)
(31111) (32111)
(211111) (41111)
(1111111) (221111)
(311111)
(2111111)
(11111111)
(End)

Examples

			a(20) = p(19) - p(15) + p(8) = 490 - 176 + 22 = 336.
		

Crossrefs

Programs

  • Maple
    f:= n -> add((-1)^(k+1)*combinat:-numbpart(n-(3*k^2-k)/2),k=1..floor((1+sqrt(24*n+1))/6)):
    map(f, [$1..100]); # Robert Israel, Aug 03 2015
  • Mathematica
    Table[Count[IntegerPartitions[n], q_ /; First[q] >= Length[q]], {n, 16}]
    (* also *)
    Table[Count[IntegerPartitions[2n],q_/;Last[q]===1 && Max[q-PadRight[Rest[q],Length[q]]]<=1 && Count[First/@Partition[q,2],?OddQ]==Count[Last/@Partition[q,2],?OddQ]],{n,16}]
    (* also *)
    Table[Count[IntegerPartitions[n], p_ /; MemberQ[p, Max[Max[p], Length[p]]]], {n, 50}] (* Clark Kimberling, Feb 28 2014 *)
  • PARI
    {a(n) = my(A=1); A = sum(m=0,n,x^m*prod(k=1,m,(1-x^(m+k-1))/(1-x^k +x*O(x^n)))); polcoeff(A,n)}
    for(n=1,60,print1(a(n),", ")) \\ Paul D. Hanna, Aug 03 2015
    
  • PARI
    my(N=50, x='x+O('x^N)); Vec(1/prod(k=1, N, 1-x^k)*sum(k=1, N, (-1)^(k-1)*x^(k*(3*k-1)/2))) \\ Seiichi Manyama, May 21 2023

Formula

a(n) = (A000041(n) + A047993(n))/2.
a(n) = p(n-1) - p(n-5) + p(n-12) - ... -(-1)^k*p(n-(3*k^2-k)/2) + ..., where p() is A000041(). - Vladeta Jovovic, Aug 04 2004
G.f.: Sum_{n>=1} x^n * Product_{k=1..n} (1 - x^(n+k-1))/(1 - x^k). - Paul D. Hanna, Aug 03 2015
A064173(n) + a(n) = A000041(n). - R. J. Mathar, Feb 22 2023
G.f.: (1/Product_{k>=1} (1-x^k)) * Sum_{k>=1} (-1)^(k-1) * x^(k*(3*k-1)/2). - Seiichi Manyama, May 21 2023

Extensions

Mathematica programs modified by Clark Kimberling, Feb 12 2014

A087897 Number of partitions of n into odd parts greater than 1.

Original entry on oeis.org

1, 0, 0, 1, 0, 1, 1, 1, 1, 2, 2, 2, 3, 3, 4, 5, 5, 6, 8, 8, 10, 12, 13, 15, 18, 20, 23, 27, 30, 34, 40, 44, 50, 58, 64, 73, 83, 92, 104, 118, 131, 147, 166, 184, 206, 232, 256, 286, 320, 354, 394, 439, 485, 538, 598, 660, 730, 809, 891, 984, 1088, 1196, 1318, 1454, 1596, 1756
Offset: 0

Views

Author

N. J. A. Sloane, Dec 04 2003

Keywords

Comments

Also number of partitions of n into distinct parts which are not powers of 2.
Also number of partitions of n into distinct parts such that the two largest parts differ by 1.
Also number of partitions of n such that the largest part occurs an odd number of times that is at least 3 and every other part occurs an even number of times. Example: a(10) = 2 because we have [2,2,2,1,1,1,1] and [2,2,2,2,2]. - Emeric Deutsch, Mar 30 2006
Also difference between number of partitions of 1+n into distinct parts and number of partitions of n into distinct parts. - Philippe LALLOUET, May 08 2007
In the Berndt reference replace {a -> -x, q -> x} in equation (3.1) to get f(x). G.f. is 1 - x * (1 - f(x)).
Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
Also number of symmetric unimodal compositions of n+3 where the maximal part appears three times. - Joerg Arndt, Jun 11 2013
Let c(n) = number of palindromic partitions of n whose greatest part has multiplicity 3; then c(n) = a(n-3) for n>=3. - Clark Kimberling, Mar 05 2014
From Gus Wiseman, Aug 22 2021: (Start)
Also the number of integer partitions of n - 1 whose parts cover an interval of positive integers starting with 2. These partitions are ranked by A339886. For example, the a(6) = 1 through a(16) = 5 partitions are:
32 222 322 332 432 3322 3332 4332 4432 5432 43332
2222 3222 22222 4322 33222 33322 33332 44322
32222 222222 43222 43322 333222
322222 332222 432222
2222222 3222222
(End)

Examples

			1 + x^3 + x^5 + x^6 + x^7 + x^8 + 2*x^9 + 2*x^10 + 2*x^11 + 3*x^12 + 3*x^13 + ...
q + q^73 + q^121 + q^145 + q^169 + q^193 + 2*q^217 + 2*q^241 + 2*q^265 + ...
a(10)=2 because we have [7,3] and [5,5].
From _Joerg Arndt_, Jun 11 2013: (Start)
There are a(22)=13 symmetric unimodal compositions of 22+3=25 where the maximal part appears three times:
01:  [ 1 1 1 1 1 1 1 1 3 3 3 1 1 1 1 1 1 1 1 ]
02:  [ 1 1 1 1 1 1 2 3 3 3 2 1 1 1 1 1 1 ]
03:  [ 1 1 1 1 1 5 5 5 1 1 1 1 1 ]
04:  [ 1 1 1 1 2 2 3 3 3 2 2 1 1 1 1 ]
05:  [ 1 1 1 2 5 5 5 2 1 1 1 ]
06:  [ 1 1 2 2 2 3 3 3 2 2 2 1 1 ]
07:  [ 1 1 3 5 5 5 3 1 1 ]
08:  [ 1 1 7 7 7 1 1 ]
09:  [ 1 2 2 5 5 5 2 2 1 ]
10:  [ 1 4 5 5 5 4 1 ]
11:  [ 2 2 2 2 3 3 3 2 2 2 2 ]
12:  [ 2 3 5 5 5 3 2 ]
13:  [ 2 7 7 7 2 ]
(End)
From _Gus Wiseman_, Feb 16 2021: (Start)
The a(7) = 1 through a(19) = 8 partitions are the following (A..J = 10..19). The Heinz numbers of these partitions are given by A341449.
  7  53  9    55  B    75    D    77    F      97    H      99      J
         333  73  533  93    553  95    555    B5    755    B7      775
                       3333  733  B3    753    D3    773    D5      955
                                  5333  933    5533  953    F3      973
                                        33333  7333  B33    5553    B53
                                                     53333  7533    D33
                                                            9333    55333
                                                            333333  73333
(End)
		

References

  • J. W. L. Glaisher, Identities, Messenger of Mathematics, 5 (1876), pp. 111-112. see Eq. I

Crossrefs

The ordered version is A000931.
Partitions with no ones are counted by A002865, ranked by A005408.
The even version is A035363, ranked by A066207.
The version for factorizations is A340101.
Partitions whose only even part is the smallest are counted by A341447.
The Heinz numbers of these partitions are given by A341449.
A000009 counts partitions into odd parts, ranked by A066208.
A025147 counts strict partitions with no 1's.
A025148 counts strict partitions with no 1's or 2's.
A026804 counts partitions whose smallest part is odd, ranked by A340932.
A027187 counts partitions with even length/maximum, ranks A028260/A244990.
A027193 counts partitions with odd length/maximum, ranks A026424/A244991.
A058695 counts partitions of odd numbers, ranked by A300063.
A058696 counts partitions of even numbers, ranked by A300061.
A340385 counts partitions with odd length and maximum, ranked by A340386.

Programs

  • Haskell
    a087897 = p [3,5..] where
       p [] _ = 0
       p _  0 = 1
       p ks'@(k:ks) m | m < k     = 0
                      | otherwise = p ks' (m - k) + p ks m
    -- Reinhard Zumkeller, Aug 12 2011
    
  • Maple
    To get 128 terms: t4 := mul((1+x^(2^n)),n=0..7); t5 := mul((1+x^k),k=1..128): t6 := series(t5/t4,x,100); t7 := seriestolist(t6);
    # second Maple program:
    b:= proc(n, i) option remember; `if`(n=0, 1,
          `if`(i<3, 0, b(n, i-2)+`if`(i>n, 0, b(n-i, i))))
        end:
    a:= n-> b(n, n-1+irem(n, 2)):
    seq(a(n), n=0..80);  # Alois P. Heinz, Jun 11 2013
  • Mathematica
    max = 65; f[x_] := Product[ 1/(1 - x^(2k+1)), {k, 1, max}]; CoefficientList[ Series[f[x], {x, 0, max}], x] (* Jean-François Alcover, Dec 16 2011, after Emeric Deutsch *)
    b[n_, i_] := b[n, i] = If[n==0, 1, If[i<3, 0, b[n, i-2]+If[i>n, 0, b[n-i, i]]] ]; a[n_] := b[n, n-1+Mod[n, 2]]; Table[a[n], {n, 0, 80}] (* Jean-François Alcover, Apr 01 2015, after Alois P. Heinz *)
    Flatten[{1, Table[PartitionsQ[n+1] - PartitionsQ[n], {n, 0, 80}]}] (* Vaclav Kotesovec, Dec 01 2015 *)
    Table[Length[Select[IntegerPartitions[n],FreeQ[#,1]&&OddQ[Times@@#]&]],{n,0,30}] (* Gus Wiseman, Feb 16 2021 *)
  • PARI
    {a(n) = local(A); if( n<0, 0, A = x * O(x^n); polcoeff( (1 - x) * eta(x^2 + A) / eta(x + A), n))} /* Michael Somos, Nov 13 2011 */
    
  • Python
    from functools import lru_cache
    @lru_cache(maxsize=None)
    def A087897_T(n,k):
        if n==0: return 1
        if k<3 or n<0: return 0
        return A087897_T(n,k-2)+A087897_T(n-k,k)
    def A087897(n): return A087897_T(n,n-(n&1^1)) # Chai Wah Wu, Sep 23 2023, after Alois P. Heinz

Formula

Expansion of q^(-1/24) * (1 - q) * eta(q^2) / eta(q) in powers of q.
Expansion of (1 - x) / chi(-x) in powers of x where chi() is a Ramanujan theta function.
G.f.: 1 + x^3 + x^5*(1 + x) + x^7*(1 + x)*(1 + x^2) + x^9*(1 + x)*(1 + x^2)*(1 + x^3) + ... [Glaisher 1876]. - Michael Somos, Jun 20 2012
G.f.: Product_{k >= 1} 1/(1-x^(2*k+1)).
G.f.: Product_{k >= 1, k not a power of 2} (1+x^k).
G.f.: Sum_{k >= 1} x^(3*k)/Product_{j = 1..k} (1 - x^(2*j)). - Emeric Deutsch, Mar 30 2006
a(n) ~ exp(Pi*sqrt(n/3)) * Pi / (8 * 3^(3/4) * n^(5/4)) * (1 - (15*sqrt(3)/(8*Pi) + 11*Pi/(48*sqrt(3)))/sqrt(n) + (169*Pi^2/13824 + 385/384 + 315/(128*Pi^2))/n). - Vaclav Kotesovec, Aug 30 2015, extended Nov 04 2016
G.f.: 1/(1 - x^3) * Sum_{n >= 0} x^(5*n)/Product_{k = 1..n} (1 - x^(2*k)) = 1/((1 - x^3)*(1 - x^5)) * Sum_{n >= 0} x^(7*n)/Product_{k = 1..n} (1 - x^(2*k)) = ..., extending Deutsch's result dated Mar 30 2006. - Peter Bala, Jan 15 2021
G.f.: Sum_{n >= 0} x^(n*(2*n+1))/Product_{k = 2..2*n+1} (1 - x^k). (Set z = x^3 and q = x^2 in Mc Laughlin et al., Section 1.3, Entry 7.) - Peter Bala, Feb 02 2021
a(2*n+1) = Sum{j>=1} A008284(n+1-j,2*j - 1) and a(2*n) = Sum{j>=1} A008284(n-j, 2*j). - Gregory L. Simay, Sep 22 2023

A168659 Number of partitions of n such that the number of parts is divisible by the greatest part. Also number of partitions of n such that the greatest part is divisible by the number of parts.

Original entry on oeis.org

1, 1, 2, 2, 3, 3, 6, 6, 8, 9, 14, 16, 22, 25, 33, 39, 51, 60, 79, 92, 116, 137, 174, 204, 254, 300, 368, 435, 530, 625, 760, 896, 1076, 1267, 1518, 1780, 2121, 2484, 2946, 3444, 4070, 4749, 5594, 6514, 7637, 8879, 10384, 12043, 14040, 16255
Offset: 1

Views

Author

Vladeta Jovovic, Dec 02 2009

Keywords

Examples

			a(5)=3 because in the partitions [1,1,1,1,1], [1,1,1,2], [1,1,3] the number of parts is divisible by the greatest part; not true for the partitions [1,2,2],[2,3], [1,4], and [5]. - _Emeric Deutsch_, Dec 04 2009
From _Gus Wiseman_, Feb 08 2021: (Start)
The a(1) = 1 through a(10) = 9 partitions of the first type:
  1  11  21   22    311    321     322      332       333        4222
         111  1111  2111   2211    331      2222      4221       4321
                    11111  111111  2221     4211      4311       4411
                                   4111     221111    51111      52111
                                   211111   311111    222111     222211
                                   1111111  11111111  321111     322111
                                                      21111111   331111
                                                      111111111  22111111
                                                                 1111111111
The a(1) = 1 through a(11) = 14 partitions of the second type (A=10, B=11):
  1   2   3    4    5     6     7      8      9       A       B
          21   22   41    42    43     44     63      64      65
                    311   321   61     62     81      82      83
                                322    332    333     622     A1
                                331    611    621     631     632
                                4111   4211   4221    4222    641
                                              4311    4321    911
                                              51111   4411    4322
                                                      52111   4331
                                                              4421
                                                              8111
                                                              52211
                                                              53111
                                                              611111
(End)
		

Crossrefs

Note: A-numbers of Heinz-number sequences are in parentheses below.
The case of equality is A047993 (A106529).
The Heinz numbers of these partitions are A340609/A340610.
If all parts (not just the greatest) are divisors we get A340693 (A340606).
The strict case in the second interpretation is A340828 (A340856).
A006141 = partitions whose length equals their minimum (A324522).
A067538 = partitions whose length/max divides their sum (A316413/A326836).
A200750 = partitions with length coprime to maximum (A340608).
Row sums of A350879.

Programs

  • Maple
    a := proc (n) local pn, ct, j: with(combinat): pn := partition(n): ct := 0: for j to numbpart(n) do if `mod`(nops(pn[j]), max(seq(pn[j][i], i = 1 .. nops(pn[j])))) = 0 then ct := ct+1 else end if end do: ct end proc: seq(a(n), n = 1 .. 50); # Emeric Deutsch, Dec 04 2009
  • Mathematica
    Table[Length[Select[IntegerPartitions[n],Divisible[Length[#],Max[#]]&]],{n,30}] (* Gus Wiseman, Feb 08 2021 *)
    nmax = 100; s = 0; Do[s += Normal[Series[Sum[x^((m+1)*k - 1) * Product[(1 - x^(m*k + j - 1))/(1 - x^j), {j, 1, k-1}], {k, 1, (1 + nmax)/(1 + m) + 1}], {x, 0, nmax}]], {m, 1, nmax}]; Rest[CoefficientList[s, x]] (* Vaclav Kotesovec, Oct 18 2024 *)

Formula

G.f.: Sum_{i>=1} Sum_{j>=1} x^((i+1)*j-1) * Product_{k=1..j-1} (1-x^(i*j+k-1))/(1-x^k). - Seiichi Manyama, Jan 24 2022
a(n) ~ c * exp(Pi*sqrt(2*n/3)) / n^(3/2), where c = 0.04628003... - Vaclav Kotesovec, Nov 16 2024

Extensions

Extended by Emeric Deutsch, Dec 04 2009

A064173 Number of partitions of n with positive rank.

Original entry on oeis.org

0, 1, 1, 2, 3, 5, 6, 10, 13, 19, 25, 35, 45, 62, 80, 106, 136, 178, 225, 291, 366, 466, 583, 735, 912, 1140, 1407, 1743, 2140, 2634, 3214, 3932, 4776, 5807, 7022, 8495, 10225, 12313, 14762, 17696, 21136, 25236, 30030, 35722, 42367, 50216, 59368, 70138, 82665
Offset: 1

Views

Author

Vladeta Jovovic, Sep 19 2001

Keywords

Comments

The rank of a partition is the largest summand minus the number of summands.
Also number of partitions of n with negative rank. - Omar E. Pol, Mar 05 2012
Column 1 of A208478. - Omar E. Pol, Mar 11 2012
Number of partitions p of n such that max(max(p), number of parts of p) is not a part of p. - Clark Kimberling, Feb 28 2014
The sequence enumerates the semigroup of partitions of positive rank for each number n. The semigroup is a subsemigroup of the monoid of partitions of nonnegative rank under the binary operation "*": Let A be the positive rank partition (a1,...,ak) where ak > k, and let B=(b1,...bj) with bj > j. Then let A*B be the partition (a1b1,...,a1bj,...,akb1,...,akbj), which has akbj > kj, thus having positive rank. For example, the partition (2,3,4) of 9 has rank 1, and its product with itself is (4,6,6,8,8,9,12,12,16) of 81, which has rank 7. A similar situation holds for partitions of negative rank--they are a subsemigroup of the monoid of nonpositive rank partitions. - Richard Locke Peterson, Jul 15 2018

Examples

			a(20) = p(18) - p(13) + p(5) = 385 - 101 + 7 = 291.
From _Gus Wiseman_, Feb 09 2021: (Start)
The a(2) = 1 through a(9) = 13 partitions of positive rank:
  (2)  (3)  (4)   (5)   (6)    (7)    (8)     (9)
            (31)  (32)  (33)   (43)   (44)    (54)
                  (41)  (42)   (52)   (53)    (63)
                        (51)   (61)   (62)    (72)
                        (411)  (421)  (71)    (81)
                               (511)  (422)   (432)
                                      (431)   (441)
                                      (521)   (522)
                                      (611)   (531)
                                      (5111)  (621)
                                              (711)
                                              (5211)
                                              (6111)
(End)
		

Crossrefs

Note: A-numbers of ranking sequences are in parentheses below.
The negative-rank version is also A064173 (A340788).
The case of odd positive rank is A101707 (A340604).
The case of even positive rank is A101708 (A340605).
These partitions are ranked by (A340787).
A063995/A105806 count partitions by rank.
A072233 counts partitions by sum and length.
A168659 counts partitions whose length is a multiple of the greatest part.
A200750 counts partitions whose length and greatest part are coprime.
- Rank -
A064174 counts partitions of nonnegative/nonpositive rank (A324562/A324521).
A101198 counts partitions of rank 1 (A325233).
A257541 gives the rank of the partition with Heinz number n.
A340601 counts partitions of even rank (A340602).
A340692 counts partitions of odd rank (A340603).
- Balance -
A047993 counts balanced partitions (A106529).
A340599 counts alt-balanced factorizations.
A340653 counts balanced factorizations.

Programs

  • Maple
    A064173 := proc(n)
        a := 0 ;
        for p in combinat[partition](n) do
            r := max(op(p))-nops(p) ;
            if r > 0 then
                a := a+1 ;
            end if;
        end do:
        a ;
    end proc:
    seq(A064173(n),n=0..40) ;# Emeric Deutsch, Dec 11 2004
  • Mathematica
    Table[Count[IntegerPartitions[n], q_ /; First[q] > Length[q]], {n, 24}] (* Clark Kimberling, Feb 12 2014 *)
    Table[Count[IntegerPartitions[n], p_ /; ! MemberQ[p, Max[Max[p], Length[p]]]], {n, 20}] (* Clark Kimberling, Feb 28 2014 *)
    P = PartitionsP;
    a[n_] := (P[n] - Sum[-(-1)^k (P[n - (3k^2 - k)/2] - P[n - (3k^2 + k)/2]), {k, 1, Floor[(1 + Sqrt[1 + 24n])/6]}])/2;
    a /@ Range[48] (* Jean-François Alcover, Jan 11 2020, after Wouter Meeussen in A047993 *)
  • PARI
    my(N=66, x='x+O('x^N)); concat(0, Vec(sum(k=1, N, x^k*prod(j=1, k, (1-x^(k+j-2))/(1-x^j))))) \\ Seiichi Manyama, Jan 25 2022

Formula

a(n) = (A000041(n) - A047993(n))/2.
a(n) = p(n-2) - p(n-7) + p(n-15) - ... - (-1)^k*p(n-(3*k^2+k)/2) + ..., where p() is A000041(). - Vladeta Jovovic, Aug 04 2004
G.f.: Product_{k>=1} (1/(1-q^k)) * Sum_{k>=1} ( (-1)^k * (-q^(3*k^2/2+k/2))) (conjectured). - Thomas Baruchel, May 12 2018
G.f.: Sum_{k>=1} x^k * Product_{j=1..k} (1-x^(k+j-2))/(1-x^j). - Seiichi Manyama, Jan 25 2022
a(n)+A064174(n) = A000041(n). - R. J. Mathar, Feb 22 2023

A340601 Number of integer partitions of n of even rank.

Original entry on oeis.org

1, 1, 0, 3, 1, 5, 3, 11, 8, 18, 16, 34, 33, 57, 59, 98, 105, 159, 179, 262, 297, 414, 478, 653, 761, 1008, 1184, 1544, 1818, 2327, 2750, 3480, 4113, 5137, 6078, 7527, 8899, 10917, 12897, 15715, 18538, 22431, 26430, 31805, 37403, 44766, 52556, 62620, 73379
Offset: 0

Views

Author

Gus Wiseman, Jan 21 2021

Keywords

Comments

The Dyson rank of a nonempty partition is its maximum part minus its number of parts. For this sequence, the rank of an empty partition is 0.

Examples

			The a(1) = 1 through a(9) = 18 partitions (empty column indicated by dot):
  (1)  .  (3)    (22)  (5)      (42)    (7)        (44)      (9)
          (21)         (41)     (321)   (43)       (62)      (63)
          (111)        (311)    (2211)  (61)       (332)     (81)
                       (2111)           (322)      (521)     (333)
                       (11111)          (331)      (2222)    (522)
                                        (511)      (4211)    (531)
                                        (2221)     (32111)   (711)
                                        (4111)     (221111)  (4221)
                                        (31111)              (4311)
                                        (211111)             (6111)
                                        (1111111)            (32211)
                                                             (33111)
                                                             (51111)
                                                             (222111)
                                                             (411111)
                                                             (3111111)
                                                             (21111111)
                                                             (111111111)
		

Crossrefs

Note: Heinz numbers are given in parentheses below.
The positive case is A101708 (A340605).
The Heinz numbers of these partitions are A340602.
The odd version is A340692 (A340603).
- Rank -
A047993 counts partitions of rank 0 (A106529).
A072233 counts partitions by sum and length.
A101198 counts partitions of rank 1 (A325233).
A101707 counts partitions of odd positive rank (A340604).
A101708 counts partitions of even positive rank (A340605).
A257541 gives the rank of the partition with Heinz number n.
A340653 counts factorizations of rank 0.
- Even -
A024430 counts set partitions of even length.
A027187 counts partitions of even length (A028260).
A027187 (also) counts partitions of even maximum (A244990).
A034008 counts compositions of even length.
A035363 counts partitions into even parts (A066207).
A052841 counts ordered set partitions of even length.
A058696 counts partitions of even numbers (A300061).
A067661 counts strict partitions of even length (A030229).
A236913 counts even-length partitions of even numbers (A340784).
A339846 counts factorizations of even length.

Programs

  • Maple
    b:= proc(n, i, r) option remember; `if`(n=0, 1-max(0, r),
          `if`(i<1, 0, b(n, i-1, r) +b(n-i, min(n-i, i), 1-
          `if`(r<0, irem(i, 2), r))))
        end:
    a:= n-> b(n$2, -1):
    seq(a(n), n=0..55);  # Alois P. Heinz, Jan 22 2021
  • Mathematica
    Table[If[n==0,1,Length[Select[IntegerPartitions[n],EvenQ[Max[#]-Length[#]]&]]],{n,0,30}]
    (* Second program: *)
    b[n_, i_, r_] := b[n, i, r] = If[n == 0, 1 - Max[0, r], If[i < 1, 0, b[n, i - 1, r] + b[n - i, Min[n - i, i], 1 - If[r < 0, Mod[i, 2], r]]]];
    a[n_] := b[n, n, -1];
    a /@ Range[0, 55] (* Jean-François Alcover, May 10 2021, after Alois P. Heinz *)
  • PARI
    p_q(k) = {prod(j=1, k, 1-q^j); }
    GB_q(N, M)= {if(N>=0 && M>=0,  p_q(N+M)/(p_q(M)*p_q(N)), 0 ); }
    A_q(N) = {my(q='q+O('q^N), g=1+sum(i=1,N, sum(j=1,N/i, q^(i*j) * ( ((1/2)*(1+(-1)^(i+j))) + sum(k=1,N-(i*j), ((q^k)*GB_q(k,i-2)) * ((1/2)*(1+(-1)^(i+j+k)))))))); Vec(g)}
    A_q(50) \\ John Tyler Rascoe, Apr 15 2024

Formula

G.f.: 1 + Sum_{i, j>0} q^(i*j) * ( (1+(-1)^(i+j))/2 + Sum_{k>0} q^k * q_binomial(k,i-2) * (1+(-1)^(i+j+k))/2 ). - John Tyler Rascoe, Apr 15 2024
a(n) ~ exp(Pi*sqrt(2*n/3)) / (8*n*sqrt(3)). - Vaclav Kotesovec, Apr 17 2024

A117144 Partitions of n in which each part k occurs at least k times.

Original entry on oeis.org

1, 1, 1, 1, 2, 2, 3, 3, 4, 5, 6, 6, 8, 9, 10, 12, 15, 16, 19, 21, 25, 28, 32, 34, 41, 46, 51, 55, 64, 70, 79, 86, 97, 106, 119, 129, 146, 159, 175, 190, 214, 232, 256, 277, 306, 334, 367, 394, 434, 472, 515, 556, 607, 654, 714, 770, 836, 901, 978, 1048, 1140, 1226, 1322
Offset: 0

Views

Author

Emeric Deutsch, Mar 06 2006

Keywords

Comments

The Heinz numbers of these integer partitions are given by A324525. - Gus Wiseman, Mar 09 2019

Examples

			a(9)=5 because we have [3,3,3], [2,2,2,2,1], [2,2,2,1,1,1], [2,2,1,1,1,1,1] and [1,1,1,1,1,1,1,1,1].
From _Gus Wiseman_, Mar 09 2019: (Start)
The a(1) = 1 through a(9) = 5 integer partitions:
  1  11  111  22    221    222     2221     2222      333
              1111  11111  2211    22111    22211     22221
                           111111  1111111  221111    222111
                                            11111111  2211111
                                                      111111111
(End)
		

Crossrefs

Programs

  • Maple
    g:=product((1-x^k+x^(k^2))/(1-x^k),k=1..100): gser:=series(g,x=0,70): seq(coeff(gser,x,n),n=0..66);
    # second Maple program:
    b:= proc(n,i) option remember; `if`(n=0, 1, `if`(i<1, 0,
          b(n, i-1) +add(b(n-i*j, i-1), j=i..n/i)))
        end:
    a:= n-> b(n$2):
    seq(a(n), n=0..80);  # Alois P. Heinz, Dec 28 2016
  • Mathematica
    b[n_, i_] := b[n, i] = If[n==0, 1, If[i<1, 0, b[n, i-1] + Sum[b[n-i*j, i-1], {j, i, n/i}]]]; a[n_] := b[n, n]; Table[a[n], {n, 0, 80}] (* Jean-François Alcover, Feb 03 2017, after Alois P. Heinz *)
    Table[Length[Select[IntegerPartitions[n],And@@Table[Count[#,i]>=i,{i,Union[#]}]&]],{n,0,30}] (* Gus Wiseman, Mar 09 2019 *)
    nmax = 100; CoefficientList[Series[Product[(1-x^k+x^(k^2))/(1-x^k), {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Jan 28 2024 *)

Formula

G.f.: Product_{k>=1} (1-x^k+x^(k^2))/(1-x^k).
Showing 1-10 of 28 results. Next