cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Previous Showing 11-20 of 49 results. Next

A303810 Mirror image of the triangle A026794.

Original entry on oeis.org

1, 1, 1, 1, 0, 2, 1, 0, 1, 3, 1, 0, 0, 1, 5, 1, 0, 0, 1, 2, 7, 1, 0, 0, 0, 1, 2, 11, 1, 0, 0, 0, 1, 1, 4, 15, 1, 0, 0, 0, 0, 1, 2, 4, 22, 1, 0, 0, 0, 0, 1, 1, 2, 7, 30, 1, 0, 0, 0, 0, 0, 1, 1, 3, 8, 42, 1, 0, 0, 0, 0, 0, 1, 1, 2, 4, 12, 56, 1, 0, 0, 0, 0, 0, 0, 1, 1, 2, 5, 14, 77, 1, 0, 0, 0, 0, 0, 0, 1, 1, 1, 3
Offset: 1

Views

Author

Omar E. Pol, May 05 2018

Keywords

Examples

			Triangle begins:
1;
1,  1;
1,  0,  2;
1,  0,  1,  3;
1,  0,  0,  1,  5;
1,  0,  0,  1,  2,  7;
1,  0,  0,  0,  1,  2, 11;
1,  0,  0,  0,  1,  1,  4, 15;
1,  0,  0,  0,  0,  1,  2,  4, 22;
1,  0,  0,  0,  0,  1,  1,  2,  7, 30;
1,  0,  0,  0,  0,  0,  1,  1,  3,  8, 42;
1,  0,  0,  0,  0,  0,  1,  1,  2,  4, 12, 56;
1,  0,  0,  0,  0,  0,  0,  1,  1,  2,  5, 14, 77;
1,  0,  0,  0,  0,  0,  0,  1,  1,  1,  3,  6, 21, 101;
1,  0,  0,  0,  0,  0,  0,  0,  1,  1,  2,  3,  9,  24, 135;
...
		

Crossrefs

Leading diagonal gives A000041.
Second diagonal gives A002865.
Row sums give A000041, n >= 1.
Columns 1-4: A000012, A000007, A033322, A278105.
Cf. A026794.

A000070 a(n) = Sum_{k=0..n} p(k) where p(k) = number of partitions of k (A000041).

Original entry on oeis.org

1, 2, 4, 7, 12, 19, 30, 45, 67, 97, 139, 195, 272, 373, 508, 684, 915, 1212, 1597, 2087, 2714, 3506, 4508, 5763, 7338, 9296, 11732, 14742, 18460, 23025, 28629, 35471, 43820, 53963, 66273, 81156, 99133, 120770, 146785, 177970, 215308, 259891, 313065, 376326, 451501
Offset: 0

Views

Author

Keywords

Comments

Also the total number of all different integers in all partitions of n + 1. E.g., a(3) = 7 because the partitions of 4 comprise the sets {1},{1, 2},{2},{1, 3},{4} of different integers and their total number is 7. - Thomas Wieder, Apr 10 2004
With offset 1, also the number of 1's in all partitions of n. For example, 3 = 2+1 = 1+1+1, a(3) = (zero 1's) + (one 1's) + (three 1's), so a(3) = 4. - Naohiro Nomoto, Jan 09 2002. See the Riordan reference p. 184, last formula, first term, for a proof based on Fine's identity given in Riordan, p. 182 (20).
Also, number of partitions of n into parts when there are two kinds of parts of size one.
Also number of graphical forest partitions of 2n+2.
a(n) = count 2 for each partition of n and 1 for each decrement. E.g., the partitions of 4 are 4 (2), 31 (3), 22 (2), 211 (3) and 1111 (2). 2 + 3 + 2 + 3 + 2 = 12. This is related to the Ferrers representation. We can see that taking the Ferrers diagram for each partition of n and adding a new * to all available columns, we generate each partition of n+1, but with repeats (A058884). - Jon Perry, Feb 06 2004
Also the number of 1-transitions among all integer partitions of n. A 1-transition is the removal of a digit "1" from a partition containing at least one "1" and subsequent addition of that "1" to another digit in that partition. This other digit may be a "1" also, but all digits of equal amount are considered as undistinquishable (unlabeled). E.g., for n=6 one has the partition [1113] for which the following two 1-transitions are possible: [1113] --> [123] and [1113] --> [114]. The 1-transitions of n form a partial order (poset). For n=6 one has 12 1-transitions: [111111] --> [11112], [11112] --> [1113], [11112] --> [1122], [1113] --> [114], [1113] --> [123], [1122] --> [123], [1122] --> [222], [123] --> [33], [123] --> [24], [114] --> [15], [114] --> [24], [15] --> [6]. - Thomas Wieder, Mar 08 2005
Also number of partitions of 2n+1 where one of the parts is greater than n (also where there are more than n parts) and of 2n+2 where one of the parts is greater than n+1 (or with more than n+1 parts). - Henry Bottomley, Aug 01 2005
Equals left border of triangle A137633 - Gary W. Adamson, Jan 31 2008
Equals row sums of triangle A027293. - Gary W. Adamson, Oct 26 2008
Convolved with A010815 = [1,1,1,...]. n-th partial sum of A000041 convolved with A010815 = the binomial sequence starting (1, n, ...). - Gary W. Adamson, Nov 09 2008
Equals A036469 convolved with A035363. - Gary W. Adamson, Jun 09 2009
a(A004526(n)) = A025065(n). - Reinhard Zumkeller, Jan 23 2010
a(n) = if n <= 1 then A054225(1,n) else A054225(n,1). - Reinhard Zumkeller, Nov 30 2011
Also the total number of 1's among all hook-lengths in all partitions of n. E.g., a(4)=7 because hooks of the partitions of n = 4 comprise the multisets {4,3,2,1}, {4,2,1,1}, {3,2,2,1}, {4,1,2,1}, {4,3,2,1} and their total number of 1's is 7. - T. Amdeberhan, Jun 03 2012
With offset 1, a(n) is also the difference between the sum of largest and the sum of second largest elements in all partitions of n. More generally, the number of occurrences of k in all partitions of n equals the difference between the sum of k-th largest and the sum of (k+1)st largest elements in all partitions of n. And more generally, the sum of the number of occurrences of k, k+1, k+2..k+m in all partitions of n equals the difference between the sum of k-th largest and the sum of (k+m+1)st largest elements in all partitions of n. - Omar E. Pol, Oct 25 2012
a(0) = 1 and 2*a(n-1) >= a(n) for all n > 0. Hence a(n) is a complete sequence. - Frank M Jackson, Apr 08 2013
a(n) is the number of conjugacy classes in the order-preserving, order-decreasing and (order-preserving and order-decreasing) injective transformation semigroups. - Ugbene Ifeanyichukwu, Jun 03 2015
a(n) is also the number of unlabeled subgraphs of the n-cycle C_n. For example, for n = 3, there are 3 unlabeled subgraphs of the triangle C_3 with 0 edges, 2 with 1 edge, 1 with 2 edges, and 1 with 3 edges (C_3 itself), so a(3) = 3 + 2 + 1 + 1 = 7. - John P. McSorley, Nov 21 2016
a(n) is also the number of partitions of 2n with all parts either even or equal to 1. Proof: the number of such partitions of 2n with exactly 2k 1's is p(n-k), for k = 0,..,n. Summing over k gives the formula. - Leonard Chastkofsky, Jul 24 2018
a(n) is the total number of polygamma functions that appear in the expansion of the (n+1)st derivative of x! with respect to x. More specifically, a(n) is the number of times the string "PolyGamma" appears in the expansion of D[x!, {x, n + 1}] in Mathematica. For example, D[x!, {x, 3 + 1}] = Gamma[1 + x] PolyGamma[0, 1 + x]^4 + 6 Gamma[1 + x] PolyGamma[0, 1 + x]^2 PolyGamma[1, 1 + x] + 3 Gamma[1 + x] PolyGamma[1, 1 + x]^2 + 4 Gamma[1 + x] PolyGamma[0, 1 + x] PolyGamma[2, 1 + x] + Gamma[1 + x] PolyGamma[3, 1 + x], and we see that the string "PolyGamma" appears a total of a(3) = 7 times in this expansion. - John M. Campbell, Aug 11 2018
With offset 1, also the number of integer partitions of 2n that do not comprise the multiset of vertex-degrees of any multigraph (i.e., non-multigraphical partitions); see A209816 for multigraphical partitions. - Gus Wiseman, Oct 26 2018
Also a(n) is the number of partitions of 2n+1 with exactly one odd part.
Delete the odd part 2k+1, k=0, ..., n, to get a partition of 2n-2k into even parts. There are as many unrestricted partitions of n-k; now sum those numbers from 0 to n to get a(n). - George Beck, Jul 22 2019
In the Young's lattice, a(n) is the number of branches that connect the (n-1)-th layer to the n-th layer. - Shouvik Datta, Sep 19 2021
a(n) is the number of multiset partitions of the multiset {r^n, s^1}, equivalently, factorization patterns of any number m=p^n*q^1 where p and q are primes. - Joerg Arndt, Jan 01 2024
a(n) is the number of positive integers whose divisors are the parts of the partitions of n + 1. - Omar E. Pol, Nov 07 2024

Examples

			G.f. = 1 + 2*x + 4*x^2 + 7*x^3 + 12*x^4 + 19*x^5 + 30*x^6 + 45*x^7 + 67*x^8 + ...
From _Omar E. Pol_, Oct 25 2012: (Start)
For n = 5 consider the partitions of n+1:
--------------------------------------
.                         Number
Partitions of 6           of 1's
--------------------------------------
6 .......................... 0
3 + 3 ...................... 0
4 + 2 ...................... 0
2 + 2 + 2 .................. 0
5 + 1 ...................... 1
3 + 2 + 1 .................. 1
4 + 1 + 1 .................. 2
2 + 2 + 1 + 1 .............. 2
3 + 1 + 1 + 1 .............. 3
2 + 1 + 1 + 1 + 1 .......... 4
1 + 1 + 1 + 1 + 1 + 1 ...... 6
------------------------------------
35-16 =                     19
.
The difference between the sum of the first column and the sum of the second column of the set of partitions of 6 is 35 - 16 = 19 and equals the number of 1's in all partitions of 6, so the 6th term of this sequence is a(5) = 19.
(End)
From _Gus Wiseman_, Oct 26 2018: (Start)
With offset 1, the a(1) = 1 through a(6) = 19 partitions of 2*n whose greatest part is > n:
  (2)  (4)   (6)    (8)     (A)      (C)
       (31)  (42)   (53)    (64)     (75)
             (51)   (62)    (73)     (84)
             (411)  (71)    (82)     (93)
                    (521)   (91)     (A2)
                    (611)   (622)    (B1)
                    (5111)  (631)    (732)
                            (721)    (741)
                            (811)    (822)
                            (6211)   (831)
                            (7111)   (921)
                            (61111)  (A11)
                                     (7221)
                                     (7311)
                                     (8211)
                                     (9111)
                                     (72111)
                                     (81111)
                                     (711111)
With offset 1, the a(1) = 1 through a(6) = 19 partitions of 2*n whose number of parts is > n:
  (11)  (211)   (2211)    (22211)     (222211)      (2222211)
        (1111)  (3111)    (32111)     (322111)      (3222111)
                (21111)   (41111)     (331111)      (3321111)
                (111111)  (221111)    (421111)      (4221111)
                          (311111)    (511111)      (4311111)
                          (2111111)   (2221111)     (5211111)
                          (11111111)  (3211111)     (6111111)
                                      (4111111)     (22221111)
                                      (22111111)    (32211111)
                                      (31111111)    (33111111)
                                      (211111111)   (42111111)
                                      (1111111111)  (51111111)
                                                    (222111111)
                                                    (321111111)
                                                    (411111111)
                                                    (2211111111)
                                                    (3111111111)
                                                    (21111111111)
                                                    (111111111111)
(End)
From _Joerg Arndt_, Jan 01 2024: (Start)
The a(5) = 19 multiset partitions of the multiset {1^5, 2^1} are:
   1:  {{1, 1, 1, 1, 1, 2}}
   2:  {{1, 1, 1, 1, 1}, {2}}
   3:  {{1, 1, 1, 1, 2}, {1}}
   4:  {{1, 1, 1, 1}, {1, 2}}
   5:  {{1, 1, 1, 1}, {1}, {2}}
   6:  {{1, 1, 1, 2}, {1, 1}}
   7:  {{1, 1, 1, 2}, {1}, {1}}
   8:  {{1, 1, 1}, {1, 1, 2}}
   9:  {{1, 1, 1}, {1, 1}, {2}}
  10:  {{1, 1, 1}, {1, 2}, {1}}
  11:  {{1, 1, 1}, {1}, {1}, {2}}
  12:  {{1, 1, 2}, {1, 1}, {1}}
  13:  {{1, 1, 2}, {1}, {1}, {1}}
  14:  {{1, 1}, {1, 1}, {1, 2}}
  15:  {{1, 1}, {1, 1}, {1}, {2}}
  16:  {{1, 1}, {1, 2}, {1}, {1}}
  17:  {{1, 1}, {1}, {1}, {1}, {2}}
  18:  {{1, 2}, {1}, {1}, {1}, {1}}
  19:  {{1}, {1}, {1}, {1}, {1}, {2}}
(End)
		

References

  • H. Gupta, An asymptotic formula in partitions. J. Indian Math. Soc., (N. S.) 10 (1946), 73-76.
  • H. Gupta et al., Tables of Partitions. Royal Society Mathematical Tables, Vol. 4, Cambridge Univ. Press, 1958, p. 90.
  • R. Honsberger, Mathematical Gems III, M.A.A., 1985, p. 6.
  • D. E. Knuth, The Art of Computer Programming, Vol. 4A, Table A-1, page 778. - N. J. A. Sloane, Dec 30 2018
  • A. M. Odlyzko, Asymptotic Enumeration Methods, p. 19
  • J. Riordan, Combinatorial Identities, Wiley, 1968, p. 199.
  • 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).
  • Stanley, R. P., Exercise 1.26 in Enumerative Combinatorics, Vol. 1. Cambridge, England: Cambridge University Press, p. 59, 1999.

Crossrefs

A diagonal of A066633.
Also second column of A126442. - George Beck, May 07 2011
Row sums of triangle A092905.
Also row sums of triangle A261555. - Omar E. Pol, Sep 14 2016
Also row sums of triangle A278427. - John P. McSorley, Nov 25 2016
Column k=2 of A292508.

Programs

  • GAP
    List([0..45],n->Sum([0..n],k->NrPartitions(k))); # Muniru A Asiru, Jul 25 2018
    
  • Haskell
    a000070 = p a028310_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, Nov 06 2012
    
  • Maple
    with(combinat): a:=n->add(numbpart(j),j=0..n): seq(a(n), n=0..44); # Zerinvary Lajos, Aug 26 2008
  • Mathematica
    CoefficientList[ Series[1/(1 - x)*Product[1/(1 - x^k), {k, 75}], {x, 0, 45}], x] (* Robert G. Wilson v, Jul 13 2004 *)
    Table[ Count[ Flatten@ IntegerPartitions@ n, 1], {n, 45}] (* Robert G. Wilson v, Aug 06 2008 *)
    Join[{1}, Accumulate[PartitionsP[Range[50]]]+1] (* _Harvey P. Dale, Mar 12 2013 *)
    a[ n_] := SeriesCoefficient[ 1 / (1 - x) / QPochhammer[ x], {x, 0, n}]; (* Michael Somos, Nov 09 2013 *)
    Accumulate[PartitionsP[Range[0,49]]] (* George Beck, Oct 23 2014; typo fixed by Virgile Andreani, Jul 10 2016 *)
  • PARI
    {a(n) = if( n<0, 0, polcoeff( 1 / prod(m=1, n, 1 - x^m, 1 + x * O(x^n)) / (1 - x), n))}; /* Michael Somos, Nov 08 2002 */
    
  • PARI
    x='x+O('x^66); Vec(1/((1-x)*eta(x))) /* Joerg Arndt, May 15 2011 */
    
  • PARI
    a(n) = sum(k=0, n, numbpart(k)); \\ Michel Marcus, Sep 16 2016
    
  • Python
    from itertools import accumulate
    def A000070iter(n):
        L = [0]*n; L[0] = 1
        def numpart(n):
            S = 0; J = n-1; k = 2
            while 0 <= J:
                T = L[J]
                S = S+T if (k//2)%2 else S-T
                J -= k  if (k)%2 else k//2
                k += 1
            return S
        for j in range(1, n): L[j] = numpart(j)
        return accumulate(L)
    print(list(A000070iter(100))) # Peter Luschny, Aug 30 2019
    
  • Python
    # Using function A365676Row. Compare also A365675.
    from itertools import accumulate
    def A000070List(size: int) -> list[int]:
        return [sum(accumulate(reversed(A365676Row(n)))) for n in range(size)]
    print(A000070List(45))  # Peter Luschny, Sep 16 2023
  • Sage
    def A000070_list(leng):
        p = [number_of_partitions(n) for n in range(leng)]
        return [add(p[:k+1]) for k in range(leng)]
    A000070_list(45) # Peter Luschny, Sep 15 2014
    

Formula

Euler transform of [ 2, 1, 1, 1, 1, 1, 1, ...].
log(a(n)) ~ -3.3959 + 2.44613*sqrt(n). - Robert G. Wilson v, Jan 11 2002
a(n) = (1/n)*Sum_{k=1..n} (sigma(k)+1)*a(n-k), n > 1, a(0) = 1. - Vladeta Jovovic, Aug 22 2002
G.f.: (1/(1 - x))*Product_{m >= 1} 1/(1 - x^m).
a(n) seems to have the same parity as A027349(n+1). Comment from James Sellers, Mar 08 2006: that is true.
a(n) = A000041(2n+1) - A110618(2n+1) = A000041(2n+2) - A110618(2n+2). - Henry Bottomley, Aug 01 2005
Row sums of triangle A133735. - Gary W. Adamson, Sep 22 2007
a(n) = A092269(n+1) - A195820(n+1). - Omar E. Pol, Oct 20 2011
a(n) = A181187(n+1,1) - A181187(n+1,2). - Omar E. Pol, Oct 25 2012
From Peter Bala, Dec 23 2013: (Start)
Gupta gives the asymptotic result a(n-1) ~ sqrt(6/Pi^2)* sqrt(n)*p(n), where p(n) is the partition function A000041(n).
Let P(2,n) denote the set of partitions of n into parts k >= 2.
a(n-2) = Sum_{parts k in all partitions in P(2,n)} phi(k), where phi(k) is the Euler totient function (see A000010). Using this result and Mertens's theorem on the average order of the phi function, leads to the asymptotic result
a(n-2) ~ (6/Pi^2)*n*(p(n) - p(n-1)) = (6/Pi^2)*A138880(n) as n -> infinity. (End)
a(n) ~ exp(Pi*sqrt(2*n/3)) / (2^(3/2)*Pi*sqrt(n)) * (1 + 11*Pi/(24*sqrt(6*n)) + (73*Pi^2 - 1584)/(6912*n)). - Vaclav Kotesovec, Oct 26 2016
a(n) = A024786(n+2) + A024786(n+1). - Vaclav Kotesovec, Nov 05 2016
G.f.: exp(Sum_{k>=1} (sigma_1(k) + 1)*x^k/k). - Ilya Gutkovskiy, Aug 21 2018
a(n) = A025065(2n). - Gus Wiseman, Oct 26 2018
a(n - 1) = A000041(2n) - A209816(n). - Gus Wiseman, Oct 26 2018

A008483 Number of partitions of n into parts >= 3.

Original entry on oeis.org

1, 0, 0, 1, 1, 1, 2, 2, 3, 4, 5, 6, 9, 10, 13, 17, 21, 25, 33, 39, 49, 60, 73, 88, 110, 130, 158, 191, 230, 273, 331, 391, 468, 556, 660, 779, 927, 1087, 1284, 1510, 1775, 2075, 2438, 2842, 3323, 3872, 4510
Offset: 0

Views

Author

T. Forbes (anthony.d.forbes(AT)googlemail.com)

Keywords

Comments

a(0) = 1 because the empty partition vacuously has each part >= 3. - Jason Kimberley, Jan 11 2011
Number of partitions where the largest part occurs at least three times. - Joerg Arndt, Apr 17 2011
By removing a single part of size 3, an A026796 partition of n becomes an A008483 partition of n - 3.
For n >= 3 the sequence counts the isomorphism classes of authentication codes AC(2,n,n) with perfect secrecy and with largest probability 0.5 that an interceptor could deceive with a substituted message. - E. Keith Lloyd (ekl(AT)soton.ac.uk).
For n >= 1, also the number of regular graphs of degree 2. - Mitch Harris, Jun 22 2005
(1 + 0*x + 0*x^2 + x^3 + x^4 + x^5 + 2*x^6 + ...) = (1 + x + 2*x^2 + 3*x^3 + 5*x^4 + ...) * 1 / (1 + x + 2*x^2 + 2*x^3 + 3*x^4 + 3*x^5 + 4*x^6 + 4*x^7 + ...). - Gary W. Adamson, Jun 30 2009
Because the triangle A051031 is symmetric, a(n) is also the number of (n-3)-regular graphs on n vertices. Since the disconnected (n-3)-regular graph with minimum order is 2K_{n-2}, then for n > 4 there are no disconnected (n-3)-regular graphs on n vertices. Therefore for n > 4, a(n) is also the number of connected (n-3)-regular graphs on n vertices. - Jason Kimberley, Oct 05 2009
Number of partitions of n+2 such that 2*(number of parts) is a part. - Clark Kimberling, Feb 27 2014
For n >= 1, a(n) is the number of (1,1)-separable partitions of n, as defined at A239482. For example, the (1,1)-separable partitions of 11 are [10,1], [7,1,2,1], [6,1,3,1], [5,1,4,1], [4,1,2,1,2,1], [3,1,3,1,2,1], so that a(11) = 6. - Clark Kimberling, Mar 21 2014
From Peter Bala, Dec 01 2024: (Start)
Let P(3, n) denote the set of partitions of n into parts k >= 3. Then A000041(n) = (1/2) * Sum_{parts k in all partitions in P(3, n+3)} phi(k), where phi(k) is the Euler totient function (see A000010). For example, with n = 5, there are 3 partitions of n + 3 = 8 into parts greater then 3, namely, 8, 5 + 3 and 4 + 4, and (1/2)*(phi(8) + phi(5) + phi(3) + 2*phi(4)) = 7 = A000041(5). (End)

Crossrefs

Essentially the same sequence as A026796 and A281356.
From Jason Kimberley, Nov 07 2009 and Jan 05 2011 and Feb 03 2011: (Start)
Not necessarily connected simple regular graphs: A005176 (any degree), A051031 (triangular array), specified degree k: A000012 (k=0), A059841 (k=1), this sequence (k=2), A005638 (k=3), A033301 (k=4), A165626 (k=5), A165627 (k=6), A165628 (k=7).
2-regular simple graphs: A179184 (connected), A165652 (disconnected), this sequence (not necessarily connected).
2-regular not necessarily connected graphs without multiple edges [partitions without 2 as a part]: this sequence (no loops allowed [without 1 as a part]), A027336 (loops allowed [parts may be 1]).
Not necessarily connected 2-regular graphs with girth at least g [partitions into parts >= g]: A026807 (triangle); chosen g: A000041 (g=1 -- multigraphs with loops allowed), A002865 (g=2 -- multigraphs with loops forbidden), this sequence (g=3), A008484 (g=4), A185325 (g=5), A185326 (g=6), A185327 (g=7), A185328 (g=8), A185329 (g=9).
Not necessarily connected 2-regular graphs with girth exactly g [partitions with smallest part g]: A026794 (triangle); chosen g: A002865 (g=2), A026796 (g=3), A026797 (g=4), A026798 (g=5), A026799 (g=6), A026800(g=7), A026801 (g=8), A026802 (g=9), A026803 (g=10), ... (End)
Cf. A008284.

Programs

  • Magma
    p := NumberOfPartitions; A008483 :=  func< n | n eq 0 select 1 else n le 2 select 0 else p(n) - p(n-1) - p(n-2) + p(n-3)>; // Jason Kimberley, Jan 11 2011
    
  • Maple
    series(1/product((1-x^i),i=3..50),x,51);
    ZL := [ B,{B=Set(Set(Z, card>=3))}, unlabeled ]: seq(combstruct[count](ZL, size=n), n=0..46); # Zerinvary Lajos, Mar 13 2007
    with(combstruct):ZL2:=[S,{S=Set(Cycle(Z,card>2))}, unlabeled]:seq(count(ZL2,size=n),n=0..46); # Zerinvary Lajos, Sep 24 2007
    with(combstruct):a:=proc(m) [A,{A=Set(Cycle(Z,card>m))},unlabeled]; end: A008483:=a(2):seq(count(A008483,size=n),n=0..46); # Zerinvary Lajos, Oct 02 2007
  • Mathematica
    f[1, 1] = 1; f[n_, k_] := f[n, k] = If[n < 0, 0, If[k > n, 0, If[k == n, 1, f[n, k + 1] + f[n - k, k]]]]; Table[ f[n, 3], {n, 49}] (* Robert G. Wilson v, Jan 31 2011 *)
    Rest[Table[Count[IntegerPartitions[n], p_ /; MemberQ[p, 2*Length[p]]], {n, 50}]]  (* Clark Kimberling, Feb 27 2014 *)
  • PARI
    a(n) = numbpart(n)-numbpart(n-1)-numbpart(n-2)+numbpart(n-3) \\ Charles R Greathouse IV, Jul 19 2011
    
  • Python
    from sympy import partition
    def A008483(n): return partition(n)-partition(n-1)-partition(n-2)+partition(n-3) # Chai Wah Wu, Jun 10 2025

Formula

a(n) = p(n) - p(n - 1) - p(n - 2) + p(n - 3) where p(n) is the number of unrestricted partitions of n into positive parts (A000041).
G.f.: Product_{m>=3} 1/(1-x^m).
G.f.: (Sum_{n>=0} x^(3*n)) / (Product_{k=1..n} (1 - x^k)). - Joerg Arndt, Apr 17 2011
a(n) = A121081(n+3) - A121659(n+3). - Reinhard Zumkeller, Aug 14 2006
Euler transformation of A179184. a(n) = A179184(n) + A165652(n). - Jason Kimberley, Jan 05 2011
a(n) ~ Pi^2 * exp(Pi*sqrt(2*n/3)) / (12*sqrt(3)*n^2). - Vaclav Kotesovec, Feb 26 2015
G.f.: exp(Sum_{k>=1} x^(3*k)/(k*(1 - x^k))). - Ilya Gutkovskiy, Aug 21 2018
a(n) = Sum_{j=0..floor(n/2)} A008284(n-2*j,j). - Gregory L. Simay, Apr 27 2023
G.f.: 1 + Sum_{n >= 1} x^(n+2)/Product_{k = 0..n-1} (1 - x^(k+3)). - Peter Bala, Dec 01 2024

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).

A257993 Least gap in the partition having Heinz number n; index of the least prime not dividing n.

Original entry on oeis.org

1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 3
Offset: 1

Views

Author

Emeric Deutsch, May 18 2015

Keywords

Comments

The "least gap" of a partition is the least positive integer that is not a part of the partition. For example, the least gap of the partition [7,4,2,2,1] is 3.
We define the Heinz number of a partition p = [p_1, p_2, ..., p_r] as Product(p_j-th prime, j=1...r) (concept used by Alois P. Heinz in A215366 as an "encoding" of a partition). For example, for the partition [1, 1, 2, 4, 10] we get 2*2*3*7*29 = 2436.
In the Maple program the subprogram B yields the partition with Heinz number n.
Sum of least gaps of all partitions of m = A022567(m).
From Antti Karttunen, Aug 22 2016: (Start)
Index of the least prime not dividing n. (After a formula given by Heinz.)
Least k such that A002110(k) does not divide n.
One more than the number of trailing zeros in primorial base representation of n, A049345.
(End)
The least gap is also called the mex (minimal excludant) of the partition. - Gus Wiseman, Apr 20 2021

Examples

			a(18) = 3 because the partition having Heinz number 18 = 2*3*3 is [1,2,2], having least gap equal to 3.
		

References

  • G. E. Andrews and K. Eriksson, Integer Partitions, Cambridge Univ. Press, 2004, Cambridge.
  • Miklós Bóna, A Walk Through Combinatorics, World Scientific Publishing Co., 2002.

Crossrefs

Positions of 1's are A005408.
Positions of 2's are A047235.
The number of gaps is A079067.
The version for crank is A257989.
The triangle counting partitions by this statistic is A264401.
One more than A276084.
The version for greatest difference is A286469 or A286470.
A maximal instead of minimal version is A339662.
Positions of even terms are A342050.
Positions of odd terms are A342051.
A000070 counts partitions with a selected part.
A006128 counts partitions with a selected position.
A056239 adds up prime indices, row sums of A112798.
A073491 lists numbers with gap-free prime indices.
A238709 counts partitions by sum and least difference.
A333214 lists positions of adjacent unequal prime gaps.
A339737 counts partitions by sum and greatest gap.

Programs

  • Maple
    with(numtheory): a := proc (n) local B, q: B := proc (n) local nn, j, m: nn := op(2, ifactors(n)): for j to nops(nn) do m[j] := op(j, nn) end do: [seq(seq(pi(op(1, m[i])), q = 1 .. op(2, m[i])), i = 1 .. nops(nn))] end proc: for q while member(q, B(n)) = true do  end do: q end proc: seq(a(n), n = 1 .. 150);
    # second Maple program:
    a:= n-> `if`(n=1, 1, (s-> min({$1..(max(s)+1)} minus s))(
            {map(x-> numtheory[pi](x[1]), ifactors(n)[2])[]})):
    seq(a(n), n=1..100);  # Alois P. Heinz, May 09 2016
    # faster:
    A257993 := proc(n) local p, c; c := 1; p := 2;
    while n mod p = 0 do p := nextprime(p); c := c + 1 od: c end:
    seq(A257993(n), n=1..100); # Peter Luschny, Jun 04 2017
  • Mathematica
    A053669[n_] := For[p = 2, True, p = NextPrime[p], If[CoprimeQ[p, n], Return[p]]]; a[n_] := PrimePi[A053669[n]]; Array[a, 100] (* Jean-François Alcover, Nov 28 2016 *)
    Table[k = 1; While[! CoprimeQ[Prime@ k, n], k++]; k, {n, 100}] (* Michael De Vlieger, Jun 22 2017 *)
  • PARI
    a(n) = forprime(p=2,, if (n % p, return(primepi(p)))); \\ Michel Marcus, Jun 22 2017
  • Python
    from sympy import nextprime, primepi
    def a053669(n):
        p = 2
        while True:
            if n%p!=0: return p
            else: p=nextprime(p)
    def a(n): return primepi(a053669(n)) # Indranil Ghosh, May 12 2017
    
  • Scheme
    (define (A257993 n) (let loop ((n n) (i 1)) (let* ((p (A000040 i)) (d (modulo n p))) (if (not (zero? d)) i (loop (/ (- n d) p) (+ 1 i))))))
    ;; Antti Karttunen, Aug 22 2016
    

Formula

a(n) = A000720(A053669(n)). - Alois P. Heinz, May 18 2015
From Antti Karttunen, Aug 22-30 2016: (Start)
a(n) = 1 + A276084(n).
a(n) = A055396(A276086(n)).
A276152(n) = A002110(a(n)).
(End)
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 1 + Sum_{k>=1} 1/A002110(k) = 1.705230... (1 + A064648). - Amiram Eldar, Jul 23 2022
a(n) << log n/log log n. - Charles R Greathouse IV, Dec 03 2022

Extensions

A simpler description added to the name by Antti Karttunen, Aug 22 2016

A362614 Irregular triangle read by rows where T(n,k) is the number of integer partitions of n with k modes.

Original entry on oeis.org

1, 0, 1, 0, 2, 0, 2, 1, 0, 4, 1, 0, 5, 2, 0, 7, 3, 1, 0, 11, 3, 1, 0, 16, 4, 2, 0, 21, 6, 3, 0, 29, 8, 4, 1, 0, 43, 7, 5, 1, 0, 54, 13, 8, 2, 0, 78, 12, 8, 3, 0, 102, 17, 11, 5, 0, 131, 26, 12, 6, 1, 0, 175, 29, 17, 9, 1, 0, 233, 33, 18, 11, 2, 0, 295, 47, 25
Offset: 0

Views

Author

Gus Wiseman, May 04 2023

Keywords

Comments

A mode in a multiset is an element that appears at least as many times as each of the others. For example, the modes of {a,a,b,b,b,c,d,d,d} are {b,d}.

Examples

			Triangle begins:
   1
   0   1
   0   2
   0   2   1
   0   4   1
   0   5   2
   0   7   3   1
   0  11   3   1
   0  16   4   2
   0  21   6   3
   0  29   8   4   1
   0  43   7   5   1
   0  54  13   8   2
   0  78  12   8   3
   0 102  17  11   5
   0 131  26  12   6   1
   0 175  29  17   9   1
Row n = 8 counts the following partitions:
  (8)         (53)    (431)
  (44)        (62)    (521)
  (332)       (71)
  (422)       (3311)
  (611)
  (2222)
  (3221)
  (4211)
  (5111)
  (22211)
  (32111)
  (41111)
  (221111)
  (311111)
  (2111111)
  (11111111)
		

Crossrefs

Row sums are A000041.
Row lengths are A002024.
Removing columns 0 and 1 and taking sums gives A362607, ranks A362605.
Column k = 1 is A362608, ranks A356862.
This statistic (mode-count) is ranked by A362611.
For co-modes we have A362615, ranked by A362613.
A008284 counts partitions by length.
A096144 counts partitions by number of minima, A026794 by maxima.
A238342 counts compositions by number of minima, A238341 by maxima.
A275870 counts collapsible partitions.

Programs

  • Mathematica
    msi[ms_]:=Select[Union[ms],Count[ms,#]>=Max@@Length/@Split[ms]&];
    Table[Length[Select[IntegerPartitions[n],Length[msi[#]]==k&]],{n,0,15},{k,0,Floor[(Sqrt[1+8n]-1)/2]}]

Formula

Sum_{k=0..A003056(n)} k * T(n,k) = A372542. - Alois P. Heinz, May 05 2024

A046746 Sum of smallest parts of all partitions of n.

Original entry on oeis.org

0, 1, 3, 5, 9, 12, 20, 25, 38, 49, 69, 87, 123, 152, 205, 260, 341, 425, 555, 687, 882, 1094, 1380, 1702, 2140, 2620, 3254, 3982, 4907, 5967, 7318, 8856, 10787, 13019, 15759, 18943, 22840, 27334, 32794, 39139, 46758, 55595, 66182, 78433, 93021, 109935, 129922
Offset: 0

Views

Author

Keywords

Comments

Also total number of largest parts in all partitions of n. - Vladeta Jovovic, Feb 16 2004
To see this, consider the properties of a partition related through conjugation, such as the total number of parts and the size of the largest parts. The sums over all of the partitions of n of these two properties are equal. The size of the smallest part and the number of largest parts are two such properties (this is immediate when looking at the Ferrers diagram). - Michael Donatz, Apr 17 2011
Starting with offset 1, = the partition triangle A026794 * [1, 2, 3, ...]. - Gary W. Adamson, Feb 13 2008
For n >= 1, a(n) = T(n+1,1) + T(n+2,2) + T(n+3,3)+ ... (sum along a falling diagonal) of the partition triangle A026794. - Bob Selcoe, Jun 22 2013

Examples

			For n = 4 the five partitions of 4 are 4, 2+2, 3+1, 2+1+1, 1+1+1+1, therefore the smallest parts of all partitions of 4 are 4, 2, 1, 1, 1 and the sum is 4+2+1+1+1 = 9, so a(4) = 9. - _Omar E. Pol_, Aug 02 2013
		

Crossrefs

Row sums of A026807.

Programs

  • Maple
    b:= proc(n, i) option remember;
          `if`(n=i, n, 0) +`if`(i<1, 0, b(n, i-1) +`if`(n b(n, n):
    seq(a(n), n=0..100);  # Alois P. Heinz, Mar 28 2012
  • Mathematica
    f[n_] := Plus @@ Min /@ IntegerPartitions@ n; Array[f, 45, 0] (* Robert G. Wilson v, Apr 12 2011 *)
    b[n_, i_] := b[n, i] = If[n==i, n, 0] + If[i<1, 0, b[n, i-1] + If[nJean-François Alcover, Aug 31 2015, after Alois P. Heinz *)
    Join[{0},Table[Total[IntegerPartitions[n][[;;,-1]]],{n,50}]] (* Harvey P. Dale, Aug 24 2025 *)
  • PARI
    N=66; z='z+O('z^N);  gf=sum(k=1,N, k * z^k / prod(j=k,N, 1-z^j ) ); concat([0], Vec(gf)) \\ Joerg Arndt, Apr 17 2011

Formula

G.f.: Sum_{k>=1} k*z^k/Product_{i>=0} (1-z^(k+i)). - Vladeta Jovovic, Jun 22 2003
G.f.: Sum_{k>=1} (-1 + 1/Product_{i>=0} (1-z^(k+i))). - Vladeta Jovovic, Jun 22 2003 [Cannot verify, Joerg Arndt, Apr 17 2011]
G.f.: Sum_{j>=1} (x^j/(1-x^j))/Product_{i=1..j} (1-x^i). - Vladeta Jovovic, Aug 11 2004 [Cannot verify, Joerg Arndt, Apr 17 2011]
G.f.: Sum_{k >= 1} (-1 + z^k/(1-z^k)(1-z^{k+1})(1-z^{k+2})...). - Don Knuth, Aug 08 2002 [Cannot verify, Joerg Arndt, Apr 17 2011]
G.f.: Sum_{n>=1} (x^n/(1-x^n)) / Product_{k=1..n} (1-x^k). - Joerg Arndt, May 26 2012
a(n) = A066186(n) - A066186(n-1) - A182709(n), n >= 1. - Omar E. Pol, Aug 01 2013
a(n) ~ exp(Pi*sqrt(2*n/3)) / (4*sqrt(3)*n) * (1 + (23*Pi/(24*sqrt(6)) - sqrt(3/2)/Pi)/sqrt(n) + (1681*Pi^2/6912 - 23/16)/n). - Vaclav Kotesovec, Jul 06 2019

A362615 Irregular triangle read by rows where T(n,k) is the number of integer partitions of n with k co-modes.

Original entry on oeis.org

1, 0, 1, 0, 2, 0, 2, 1, 0, 4, 1, 0, 5, 2, 0, 7, 3, 1, 0, 10, 4, 1, 0, 13, 7, 2, 0, 16, 11, 3, 0, 23, 14, 4, 1, 0, 30, 19, 6, 1, 0, 35, 29, 11, 2, 0, 50, 34, 14, 3, 0, 61, 46, 23, 5, 0, 73, 69, 27, 6, 1, 0, 95, 81, 44, 10, 1, 0, 123, 105, 53, 14, 2
Offset: 0

Views

Author

Gus Wiseman, May 04 2023

Keywords

Comments

We define a co-mode in a multiset to be an element that appears at most as many times as each of the others. For example, the co-modes of {a,a,b,b,b,c,c} are {a,c}.

Examples

			Triangle begins:
   1
   0   1
   0   2
   0   2   1
   0   4   1
   0   5   2
   0   7   3   1
   0  10   4   1
   0  13   7   2
   0  16  11   3
   0  23  14   4   1
   0  30  19   6   1
   0  35  29  11   2
   0  50  34  14   3
   0  61  46  23   5
   0  73  69  27   6   1
   0  95  81  44  10   1
Row n = 8 counts the following partitions:
  (8)         (53)     (431)
  (44)        (62)     (521)
  (332)       (71)
  (422)       (3221)
  (611)       (3311)
  (2222)      (4211)
  (5111)      (32111)
  (22211)
  (41111)
  (221111)
  (311111)
  (2111111)
  (11111111)
		

Crossrefs

Row sums are A000041.
Row lengths are A002024.
Removing columns 0 and 1 and taking sums gives A362609, ranks A362606.
Column k = 1 is A362610, ranks A359178.
This statistic (co-mode count) is ranked by A362613.
For mode instead of co-mode we have A362614, ranked by A362611.
A008284 counts partitions by length.
A096144 counts partitions by number of minima, A026794 by maxima.
A238342 counts compositions by number of minima, A238341 by maxima.
A275870 counts collapsible partitions.

Programs

  • Mathematica
    comsi[ms_]:=Select[Union[ms],Count[ms,#]<=Min@@Length/@Split[ms]&];
    Table[Length[Select[IntegerPartitions[n],Length[comsi[#]]==k&]],{n,0,15},{k,0,Floor[(Sqrt[1+8n]-1)/2]}]

Formula

Sum_{k=0..A003056(n)} k * T(n,k) = A372632(n). - Alois P. Heinz, May 07 2024

A008484 Number of partitions of n into parts >= 4.

Original entry on oeis.org

1, 0, 0, 0, 1, 1, 1, 1, 2, 2, 3, 3, 5, 5, 7, 8, 11, 12, 16, 18, 24, 27, 34, 39, 50, 57, 70, 81, 100, 115, 140, 161, 195, 225, 269, 311, 371, 427, 505, 583, 688, 791, 928, 1067, 1248, 1434, 1668, 1914, 2223, 2546, 2945, 3370, 3889, 4443, 5113, 5834, 6698
Offset: 0

Views

Author

T. Forbes (anthony.d.forbes(AT)googlemail.com)

Keywords

Comments

a(n) is also the number of not necessarily connected 2-regular graphs on n-vertices with girth at least 4 (all such graphs are simple). The integer i corresponds to the i-cycle; addition of integers corresponds to disconnected union of cycles. - Jason Kimberley, Jan 2011 and Feb 2012
By removing a single part of size 4, an A026797 partition of n becomes an A008484 partition of n - 4. Hence this sequence is essentially the same as A026797. - Jason Kimberley, Feb 2012
Number of partitions of n+3 such that 3*(number of parts) is a part. - Clark Kimberling, Feb 27 2014
Let c(n) be the number of partitions of n such that both (number of parts) and 2*(number of parts) are parts; then c(n) = a(n-6) for n >= 6 and c(n) = 0 for n < 6. - Clark Kimberling, Mar 01 2014
a(n) is also the number of partitions of n for which three times the number of ones is twice the number of parts (conjectured). - George Beck, Aug 19 2017
Proof: Above definition is equivalent to 2 out of 3 parts being equal to 1. Arrange in triples 1, 1, >= 2, etc. Sum of each triple corresponds to sequence definition. - Martin Fuller, Aug 21 2023

Crossrefs

2-regular graphs with girth at least 4: A185114 (connected), A185224 (disconnected), this sequence (not necessarily connected).
Not necessarily connected 2-regular graphs with girth at least g [partitions into parts >= g]: A026807 (triangle); chosen g: A000041 (g=1 -- multigraphs with loops allowed), A002865 (g=2 -- multigraphs with loops forbidden), A008483 (g=3), this sequence (g=4), A185325(g=5), A185326 (g=6), A185327 (g=7), A185328 (g=8), A185329 (g=9).
Not necessarily connected 2-regular graphs with girth exactly g [partitions with smallest part g]: A026794 (triangle); chosen g: A002865 (g=2), A026796 (g=3), A026797 (g=4), A026798 (g=5), A026799 (g=6), A026800(g=7), A026801 (g=8), A026802 (g=9), A026803 (g=10).
Not necessarily connected k-regular simple graphs with girth at least 4: A185314 (any k), A185304 (triangle); specified degree k: this sequence (k=2), A185334 (k=3), A185344 (k=4), A185354 (k=5), A185364 (k=6).

Programs

  • Magma
    a:= func< n | NumberOfPartitions(n)-NumberOfPartitions(n-1)-NumberOfPartitions(n-2)+ NumberOfPartitions(n-4)+NumberOfPartitions(n-5)- NumberOfPartitions(n-6) >; [1,0,0,0,1,1,1] cat [ a(n) : n in [7..60]]; // Vincenzo Librandi, Aug 20 2017
    
  • Magma
    R:=PowerSeriesRing(Integers(), 60); Coefficients(R!( 1/(&*[1-x^(m+4): m in [0..70]]) )); // G. C. Greubel, Nov 03 2019
    
  • Maple
    series(1/product((1-x^i),i=4..65),x,60); # end of program
    ZL := [ B,{B=Set(Set(Z, card>=4))}, unlabeled ]: seq(combstruct[count](ZL, size=n), n=0..60); # Zerinvary Lajos, Mar 13 2007
  • Mathematica
    f[1, 1]=1; f[n_, k_]:= f[n, k] = If[n<0, 0, If[k>n, 0, If[k==n, 1, f[n, k +1] + f[n-k, k]]]]; Table[f[n, 4], {n, 60}] (* end of program *)
    Drop[Table[Count[IntegerPartitions[n], p_ /; MemberQ[p, 3*Length[p]]], {n, 60}],2]  (* Clark Kimberling, Feb 27 2014 *)
    Table[Count[IntegerPartitions[n],
      p_ /; 3 Count[p, 1] == 2 Length[p]], {n, 0, 60}] (* George Beck Aug 19 2017 *)
    CoefficientList[Series[1/QPochhammer[x^4, x], {x,0,60}], x] (* G. C. Greubel, Nov 03 2019 *)
  • PARI
    my(x='x+O('x^60)); Vec(1/prod(m=0,70, 1-x^(m+4))) \\ G. C. Greubel, Nov 03 2019
    
  • Sage
    def A008484_list(prec):
        P. = PowerSeriesRing(ZZ, prec)
        return P( 1/product((1-x^(m+4)) for m in (0..70)) ).list()
    A008484_list(60) # G. C. Greubel, Nov 03 2019

Formula

G.f.: 1 / Product_{m>=4} (1 - x^m).
Euler transformation of A185114. - Jason Kimberley, Jan 30 2011
Given by p(n) - p(n-1) - p(n-2) + p(n-4) + p(n-5) - p(n-6) where p(n) = A000041(n). Generally, 1/Product_{i>=K} (1 - x^i) is given by p({A}), where {A} is defined over the coefficients of Product_{i=1..K-1} (1 - x^i). In this case, K=4, so (1-x)(1-x^2)(1-x^3) = 1 - x - x^2 + x^4 + x^5 - x^6, defining {A} as above. G.f.: 1 + Sum_{i>=1} (x^4i)/Product_{j=1..i}(1 - x^j). - Jon Perry, Jul 04 2004
a(n) ~ exp(Pi*sqrt(2*n/3)) * Pi^3 / (12*sqrt(2)*n^(5/2)). - Vaclav Kotesovec, Jun 02 2018
G.f.: exp(Sum_{k>=1} x^(4*k)/(k*(1 - x^k))). - Ilya Gutkovskiy, Aug 21 2018
G.f.: 1 + Sum_{n >= 1} x^(n+3)/Product_{k = 0..n-1} (1 - x^(k+4)). - Peter Bala, Dec 01 2024

A026796 Number of partitions of n in which the least part is 3.

Original entry on oeis.org

0, 0, 0, 1, 0, 0, 1, 1, 1, 2, 2, 3, 4, 5, 6, 9, 10, 13, 17, 21, 25, 33, 39, 49, 60, 73, 88, 110, 130, 158, 191, 230, 273, 331, 391, 468, 556, 660, 779, 927, 1087, 1284, 1510, 1775, 2075, 2438, 2842, 3323, 3872, 4510, 5237, 6095, 7056, 8182, 9465, 10945, 12625
Offset: 0

Views

Author

Keywords

Comments

Let b(k) be the number of partitions of k for which twice the number of ones is the number of parts, k = 0, 1, 2, ... . Then a(n+4) = b(n), n = 0, 1, 2, ... (conjectured). - George Beck, Aug 19 2017

Crossrefs

Essentially the same sequence as A008483.
Not necessarily connected 2-regular graphs with girth at least g [partitions into parts >= g]: A026807 (triangle); chosen g: A000041 (g=1 -- multigraphs with loops allowed), A002865 (g=2 -- multigraphs with loops forbidden), A008483 (g=3), A008484 (g=4), A185325(g=5), A185326 (g=6), A185327 (g=7), A185328 (g=8), A185329 (g=9).
Not necessarily connected 2-regular graphs with girth exactly g [partitions with smallest part g]: A026794 (triangle); chosen g: A002865 (g=2 -- multigraphs with at least one pair of parallel edges, but loops forbidden), this sequence (g=3), A026797 (g=4), A026798 (g=5), A026799 (g=6), A026800 (g=7), A026801 (g=8), A026802 (g=9), A026803 (g=10).
Not necessarily connected k-regular simple graphs girth exactly 3: A198313 (any k), A185643 (triangle); fixed k: this sequence (k=2), A185133 (k=3), A185143 (k=4), A185153 (k=5), A185163 (k=6).

Programs

  • Magma
    R:=PowerSeriesRing(Integers(), 60); [0,0,0] cat Coefficients(R!( x^3/(&*[1-x^(m+3): m in [0..70]]) )); // G. C. Greubel, Nov 02 2019
    
  • Maple
    seq(coeff(series(x^3/mul(1-x^(m+3), m=0..65), x, n+1), x, n), n = 0 .. 60); # G. C. Greubel, Nov 02 2019
  • Mathematica
    Table[Count[IntegerPartitions[n], p_ /; Min@p==3], {n, 0, 60}] (* George Beck Aug 19 2017 *)
    CoefficientList[Series[x^3/QPochhammer[x^3, x], {x,0,60}], x] (* G. C. Greubel, Nov 02 2019 *)
  • PARI
    a(n) = numbpart(n-3) - numbpart(n-4) - numbpart(n-5) + numbpart(n-6); \\ Michel Marcus, Aug 20 2014
    
  • PARI
    x='x+O('x^66); Vecrev(Pol(x^3*(1-x)*(1-x^2)/eta(x))) \\ Joerg Arndt, Aug 22 2014
    
  • Sage
    def A026796_list(prec):
        P. = PowerSeriesRing(ZZ, prec)
        return P( x^3/product((1-x^(m+3)) for m in (0..65)) ).list()
    A026796_list(60) # G. C. Greubel, Nov 02 2019

Formula

G.f.: x^3 / Product_{m>=3} (1 - x^m).
a(n) = p(n-3) - p(n-4) - p(n-5) + p(n-6), where p(n) = A000041(n). - Bob Selcoe, Aug 07 2014
a(n) ~ exp(Pi*sqrt(2*n/3)) * Pi^2 / (12*sqrt(3)*n^2). - Vaclav Kotesovec, Jun 02 2018
G.f.: Sum_{k>=1} x^(3*k) / Product_{j=1..k-1} (1 - x^j). - Ilya Gutkovskiy, Nov 25 2020

Extensions

More terms from Michel Marcus, Aug 20 2014
a(0) = 0 prepended by Joerg Arndt, Aug 22 2014
Previous Showing 11-20 of 49 results. Next