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

A006128 Total number of parts in all partitions of n. Also, sum of largest parts of all partitions of n.

Original entry on oeis.org

0, 1, 3, 6, 12, 20, 35, 54, 86, 128, 192, 275, 399, 556, 780, 1068, 1463, 1965, 2644, 3498, 4630, 6052, 7899, 10206, 13174, 16851, 21522, 27294, 34545, 43453, 54563, 68135, 84927, 105366, 130462, 160876, 198014, 242812, 297201, 362587, 441546, 536104, 649791, 785437, 947812, 1140945, 1371173, 1644136, 1968379, 2351597, 2805218, 3339869, 3970648, 4712040, 5584141, 6606438, 7805507, 9207637
Offset: 0

Views

Author

Keywords

Comments

a(n) = degree of Kac determinant at level n as polynomial in the conformal weight (called h). (Cf. C. Itzykson and J.-M. Drouffe, Statistical Field Theory, Vol. 2, p. 533, eq.(98); reference p. 643, Cambridge University Press, (1989).) - Wolfdieter Lang
Also the number of one-element transitions from the integer partitions of n to the partitions of n-1 for labeled parts with the assumption that from any part z > 1 one can take an element of amount 1 in one way only. That means z is composed of z unlabeled parts of amount 1, i.e. z = 1 + 1 + ... + 1. E.g., for n=3 to n=2 we have a(3) = 6 and [111] --> [11], [111] --> [11], [111] --> [11], [12] --> [11], [12] --> [2], [3] --> [2]. For the case of z composed by labeled elements, z = 1_1 + 1_2 + ... + 1_z, see A066186. - Thomas Wieder, May 20 2004
Number of times a derivative of any order (not 0 of course) appears when expanding the n-th derivative of 1/f(x). For instance (1/f(x))'' = (2 f'(x)^2-f(x) f''(x)) / f(x)^3 which makes a(2) = 3 (by counting k times the k-th power of a derivative). - Thomas Baruchel, Nov 07 2005
Starting with offset 1, = the partition triangle A008284 * [1, 2, 3, ...]. - Gary W. Adamson, Feb 13 2008
Starting with offset 1 equals A000041: (1, 1, 2, 3, 5, 7, 11, ...) convolved with A000005: (1, 2, 2, 3, 2, 4, ...). - Gary W. Adamson, Jun 16 2009
Apart from initial 0 row sums of triangle A066633, also the Möbius transform is A085410. - Gary W. Adamson, Mar 21 2011
More generally, the total number of parts >= k in all partitions of n equals the sum of k-th largest parts of all partitions of n. In this case k = 1. Apart from initial 0 the first column of A181187. - Omar E. Pol, Feb 14 2012
Row sums of triangle A221530. - Omar E. Pol, Jan 21 2013
From Omar E. Pol, Feb 04 2021: (Start)
a(n) is also the total number of divisors of all positive integers in a sequence with n blocks where the m-th block consists of A000041(n-m) copies of m, with 1 <= m <= n. The mentioned divisors are also all parts of all partitions of n.
Apart from initial zero this is also as follows:
Convolution of A000005 and A000041.
Convolution of A006218 and A002865.
Convolution of A341062 and A000070.
Row sums of triangles A221531, A245095, A339258, A340525, A340529. (End)
Number of ways to choose a part index of an integer partition of n, i.e., partitions of n with a selected position. Selecting a part value instead of index gives A000070. - Gus Wiseman, Apr 19 2021

Examples

			For n = 4 the partitions of 4 are [4], [2, 2], [3, 1], [2, 1, 1], [1, 1, 1, 1]. The total number of parts is 12. On the other hand, the sum of the largest parts of all partitions is 4 + 2 + 3 + 2 + 1 = 12, equaling the total number of parts, so a(4) = 12. - _Omar E. Pol_, Oct 12 2018
		

References

  • S. M. Luthra, On the average number of summands in partitions of n, Proc. Nat. Inst. Sci. India Part. A, 23 (1957), p. 483-498.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Main diagonal of A210485.
Column k=1 of A256193.
The version for normal multisets is A001787.
The unordered version is A001792.
The strict case is A015723.
The version for factorizations is A066637.
A000041 counts partitions.
A000070 counts partitions with a selected part.
A336875 counts compositions with a selected part.
A339564 counts factorizations with a selected factor.

Programs

  • GAP
    List([0..60],n->Length(Flat(Partitions(n)))); # Muniru A Asiru, Oct 12 2018
  • Haskell
    a006128 = length . concat . ps 1 where
       ps _ 0 = [[]]
       ps i j = [t:ts | t <- [i..j], ts <- ps t (j - t)]
    -- Reinhard Zumkeller, Jul 13 2013
    
  • Maple
    g:= add(n*x^n*mul(1/(1-x^k), k=1..n), n=1..61):
    a:= n-> coeff(series(g,x,62),x,n):
    seq(a(n), n=0..61);
    # second Maple program:
    a:= n-> add(combinat[numbpart](n-j)*numtheory[tau](j), j=1..n):
    seq(a(n), n=0..61);  # Alois P. Heinz, Aug 23 2019
  • Mathematica
    a[n_] := Sum[DivisorSigma[0, m] PartitionsP[n - m], {m, 1, n}]; Table[ a[n], {n, 0, 41}]
    CoefficientList[ Series[ Sum[n*x^n*Product[1/(1 - x^k), {k, n}], {n, 100}], {x, 0, 100}], x]
    a[n_] := Plus @@ Max /@ IntegerPartitions@ n; Array[a, 45] (* Robert G. Wilson v, Apr 12 2011 *)
    Join[{0}, ((Log[1 - x] + QPolyGamma[1, x])/(Log[x] QPochhammer[x]) + O[x]^60)[[3]]] (* Vladimir Reshetnikov, Nov 17 2016 *)
    Length /@ Table[IntegerPartitions[n] // Flatten, {n, 50}] (* Shouvik Datta, Sep 12 2021 *)
  • PARI
    f(n)= {local(v,i,k,s,t);v=vector(n,k,0);v[n]=2;t=0;while(v[1]1,i--;s+=i*(v[i]=(n-s)\i));t+=sum(k=1,n,v[k]));t } /* Thomas Baruchel, Nov 07 2005 */
    
  • PARI
    a(n) = sum(m=1, n, numdiv(m)*numbpart(n-m)) \\ Michel Marcus, Jul 13 2013
    
  • Python
    from sympy import divisor_count, npartitions
    def a(n): return sum([divisor_count(m)*npartitions(n - m) for m in range(1, n + 1)]) # Indranil Ghosh, Apr 25 2017
    

Formula

G.f.: Sum_{n>=1} n*x^n / Product_{k=1..n} (1-x^k).
G.f.: Sum_{k>=1} x^k/(1-x^k) / Product_{m>=1} (1-x^m).
a(n) = Sum_{k=1..n} k*A008284(n, k).
a(n) = Sum_{m=1..n} of the number of divisors of m * number of partitions of n-m.
Note that the formula for the above comment is a(n) = Sum_{m=1..n} d(m)*p(n-m) = Sum_{m=1..n} A000005(m)*A000041(n-m), if n >= 1. - Omar E. Pol, Jan 21 2013
Erdős and Lehner show that if u(n) denotes the average largest part in a partition of n, then u(n) ~ constant*sqrt(n)*log n.
a(n) = A066897(n) + A066898(n), n>0. - Reinhard Zumkeller, Mar 09 2012
a(n) = A066186(n) - A196087(n), n >= 1. - Omar E. Pol, Apr 22 2012
a(n) = A194452(n) + A024786(n+1). - Omar E. Pol, May 19 2012
a(n) = A000203(n) + A220477(n). - Omar E. Pol, Jan 17 2013
a(n) = Sum_{m=1..p(n)} A194446(m) = Sum_{m=1..p(n)} A141285(m), where p(n) = A000041(n), n >= 1. - Omar E. Pol, May 12 2013
a(n) = A198381(n) + A026905(n), n >= 1. - Omar E. Pol, Aug 10 2013
a(n) = O(sqrt(n)*log(n)*p(n)), where p(n) is the partition function A000041(n). - Peter Bala, Dec 23 2013
a(n) = Sum_{m=1..n} A006218(m)*A002865(n-m), n >= 1. - Omar E. Pol, Jul 14 2014
From Vaclav Kotesovec, Jun 23 2015: (Start)
Asymptotics (Luthra, 1957): a(n) = p(n) * (C*N^(1/2) + C^2/2) * (log(C*N^(1/2)) + gamma) + (1+C^2)/4 + O(N^(-1/2)*log(N)), where N = n - 1/24, C = sqrt(6)/Pi, gamma is the Euler-Mascheroni constant A001620 and p(n) is the partition function A000041(n).
The formula a(n) = p(n) * (sqrt(3*n/(2*Pi)) * (log(n) + 2*gamma - log(Pi/6)) + O(log(n)^3)) in the abstract of the article by Kessler and Livingston (cited also in the book by Sandor, p. 495) is incorrect!
Right is: a(n) = p(n) * (sqrt(3*n/2)/Pi * (log(n) + 2*gamma - log(Pi^2/6)) + O(log(n)^3))
or a(n) ~ exp(Pi*sqrt(2*n/3)) * (log(6*n/Pi^2) + 2*gamma) / (4*Pi*sqrt(2*n)).
(End)
a(n) = Sum_{m=1..n} A341062(m)*A000070(n-m), n >= 1. - Omar E. Pol, Feb 05 2021 2014

A206437 Triangle read by rows: T(j,k) is the k-th part of the j-th region of the set of partitions of n, if 1 <= j <= A000041(n).

Original entry on oeis.org

1, 2, 1, 3, 1, 1, 2, 4, 2, 1, 1, 1, 3, 5, 2, 1, 1, 1, 1, 1, 2, 4, 2, 3, 6, 3, 2, 2, 1, 1, 1, 1, 1, 1, 1, 3, 5, 2, 4, 7, 3, 2, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 4, 2, 3, 6, 3, 2, 2, 5, 4, 8, 4, 3, 2, 2, 2, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1
Offset: 1

Views

Author

Omar E. Pol, Feb 14 2012

Keywords

Comments

Here the j-th "region" of the set of partitions of n (or more simply the j-th "region" of n) is defined to be the first h elements of the sequence formed by the smallest parts in nonincreasing order of the partitions of the largest part of the j-th partition of n, with the list of partitions in colexicographic order, where h = j - i, and i is the index of the previous partition of n whose largest part is greater than the largest part of the j-th partition of n, or i = 0 if such previous largest part does not exist. The largest part of the j-th region of n is A141285(j) and the number of parts is h = A194446(j).
Some properties of the regions of n:
- The number of regions of n equals the number of partitions of n (see A000041).
- The set of regions of n contain the sets of regions of all positive integers previous to n.
- The first j regions of n are also first j regions of all integers greater than n.
- The sums of all largest parts of all regions of n equals the total number of parts of all regions of n. See A006128(n).
- If T(j,1) is a record in the sequence then the leading diagonals of triangle formed by the first j rows give the partitions of n (see example).
- The rank of a region is the largest part minus the number of parts (see A194447).
- The sum of all ranks of the regions of n is equal to zero.
How to make a diagram of the regions and partitions of n: in the first quadrant of the square grid we draw a horizontal line {[0, 0],[n, 0]} of length n. Then we draw a vertical line {[0, 0],[0, p(n)]} of length p(n) where p(n) is the number of partitions of n. Then, for j = 1..p(n), we draw a horizontal line {[0, j],[g, j]} where g = A141285(j) is the largest part of the j-th partition of n, with the list of partitions in colexicographic order. Then, for n = 1 .. p(n), we draw a vertical line from the point [g,j] down to intercept the next segment in a lower row. So we have a number of closed regions. Then we divide each region of n in horizontal rectangles with shorter sides = 1. We can see that in the original rectangle of area n*p(n) each row contains a set of rectangles whose areas are equal to the parts of one of the partitions of n. Then each region of n is labeled according to the position of its largest part on axis "y". Note that each region of n is similar to a mirror version of the Young diagram of one of the partitions of s, where s is the sum of all parts of the region. See the illustrations of the seven regions of 5 in the Links section.
Note that if row j of triangle contains parts of size 1 then the parts of row j are the smallest parts of all partitions of T(j,1), (see A046746), and also T(j,1) is a record in the sequence and also j is the number of partitions of T(j,1), (see A000041). Otherwise, if row j does not contain parts of size 1 then the parts of row j are the emergent parts of the next record in the sequence (see A183152). Row j is also the partition of A186412(j).
Also triangle read by rows in which row r lists the parts of the last section of the set of partitions of r, ordered by regions, such that the previous parts to the part of size r are the emergent parts of the partitions of r (see A138152) and the rest are the smallest parts of the partitions of r (see example). - Omar E. Pol, Apr 28 2012

Examples

			-------------------------------------------
  Region j   Triangle of parts
-------------------------------------------
  1          1;
  2          2,1;
  3          3,1,1;
  4          2;
  5          4,2,1,1,1;
  6          3;
  7          5,2,1,1,1,1,1;
  8          2;
  9          4,2;
  10         3;
  11         6,3,2,2,1,1,1,1,1,1,1;
  12         3;
  13         5,2;
  14         4;
  15         7,3,2,2,1,1,1,1,1,1,1,1,1,1,1;
.
The rotated triangle shows each row as a partition:
                             7
                           4   3
                         5       2
                       3   2       2
                     6               1
                   3   3               1
                 4       2               1
               2   2       2               1
             5               1               1
           3   2               1               1
         4       1               1               1
       2   2       1               1               1
     3       1       1               1               1
   2   1       1       1               1               1
 1   1   1       1       1               1               1
.
Alternative interpretation of this sequence:
Triangle read by rows in which row r lists the parts of the last section of the set of partitions of r ordered by regions (see comments):
   [1];
   [2,1];
   [3,1,1];
   [2],[4,2,1,1,1];
   [3],[5,2,1,1,1,1,1];
   [2],[4,2],[3],[6,3,2,2,1,1,1,1,1,1,1];
   [3],[5,2],[4],[7,3,2,2,1,1,1,1,1,1,1,1,1,1,1];
		

Crossrefs

Positive integers in A193870. Column 1 is A141285. Row j has length A194446(j). Row sums give A186412. Records are A000027.

Programs

  • Mathematica
    lex[n_]:=DeleteCases[Sort@PadRight[Reverse /@ IntegerPartitions@n], x_ /; x==0, 2];
    reg = {}; l = {};
    For[j = 1, j <= 22, j++,
      mx = Max@lex[j][[j]]; AppendTo[l, mx];
      For[i = j, i > 0, i--, If[l[[i]] > mx, Break[]]];
      AppendTo[reg, Take[Reverse[First /@ lex[mx]], j - i]];
      ];
    Flatten@reg  (* Robert Price, Apr 21 2020, revised Jul 24 2020 *)

Extensions

Further edited by Omar E. Pol, Mar 31 2012, Jan 27 2013
Minor edits by Omar E. Pol, Apr 23 2020
Comments corrected (following a suggestion from Peter Munn) by Omar E. Pol, Jul 20 2025

A026794 Triangular array T read by rows: T(n,k) = number of partitions of n in which least part is k, 1<=k<=n.

Original entry on oeis.org

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

Views

Author

Keywords

Comments

At least one part is k and each part is at least k.
From Emeric Deutsch, Feb 19 2006: (Start)
Also number of partitions of n in which the largest part occurs exactly k times. Example: T(6,2)=2 because we have [3,3] and [2,2,1,1].
G.f. of column k is x^k/prod(j>=k, 1-x^j ) (k>=1).
Row sums yield the partition numbers (A000041).
T(n,1) = A000041(n-1) (the partition numbers).
T(n,2) = A002865(n-2) (n>=2).
T(n,3)=A026796(n). T(n,4) = A026797(n). T(n,5) = A026798(n). T(n,6) = A026799(n). T(n,7) = A026800(n). T(n,8) = A026801(n). T(n,9) = A026802(n). T(n,10) = A026803(n).
Sum(k*T(n,k),k=1..n) = A046746(n). (End)
Triangle inverse = A161363. - Gary W. Adamson, Jun 07 2009
T(n,g) is also the number of not necessarily connected 2-regular graphs with girth exactly g: the part i corresponds to the i-cycle; addition of integers corresponds to disconnected union of cycles. - Jason Kimberley, Feb 05 2012
From Bob Selcoe, Jul 24 2014 (Start):
Below is a process to generate equations for column k.
Let P be the partition numbers A000041(n-j) and let f(k) denote equations which generate column k.
To find f(k), start with f(1) = P(n-j), j=1. Thus T(n,1) = f(1) = P(n-1). This is the equation for column 1.
To find f(k) k>1, first sum the terms of f(k-1) replacing the value j with j+1, and then subtract the terms of f(k-1) replacing the value j with j+k. So to find f(2) (i.e., the equation for column 2, where k=2), start with f(1) = P(n-1); first replace j with j+1 (yielding P(n-2)), and then replace j with j+2 (yielding P(n-3)). Subtracting the second term from the first, we get: f(2) = P(n-2) - P(n-3).
To find f(3), start with f(2), replace j with j+1 (yielding P(n-3) - P(n-4)) and then replace j with j+3 (yielding P(n-5) - P(n-6)). Subtracting the second group of terms from the first, we get: f(3) = P(n-3) - P(n-4) - P(n-5) + P(n-6). This is the equation for column 3; also the equation for T(n,3) = A026796(n). So for example, T(13,3) = 5 because P(13-3) - P(13-4) - P(13-5) + P(13-6) = 42 - 30 - 22 + 15 = 5.
Continue as above to find f(k) k={4..inf.}. This will generate equations for T(n,4) = A026797(n), T(n,5) = A026798(n), T(n,6) = A026799(n), ad inf.
(End)

Examples

			T(12,3) = 4 because we have [9,3], [6,3,3], [5,4,3] and [3,3,3,3]. - Edited by _Bob Selcoe_, Sep 03 2016
Triangle starts:
    1;
    1,  1;
    2,  0, 1;
    3,  1, 0, 1;
    5,  1, 0, 0, 1;
    7,  2, 1, 0, 0, 1;
   11,  2, 1, 0, 0, 0, 1;
   15,  4, 1, 1, 0, 0, 0, 1;
   22,  4, 2, 1, 0, 0, 0, 0, 1;
   30,  7, 2, 1, 1, 0, 0, 0, 0, 1;
   42,  8, 3, 1, 1, 0, 0, 0, 0, 0, 1;
   56, 12, 4, 2, 1, 1, 0, 0, 0, 0, 0, 1;
   77, 14, 5, 2, 1, 1, 0, 0, 0, 0, 0, 0, 1;
  101, 21, 6, 3, 1, 1, 1, 0, 0, 0, 0, 0, 0, 1;
  135, 24, 9, 3, 2, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1;
  ...
		

Crossrefs

Row sums give A000041.
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). For g >= 3, girth at least g implies no loops or parallel edges. - Jason Kimberley, Feb 05 2012
Not necessarily connected 2-regular graphs with girth exactly g [partitions with smallest part g]: this sequence (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). - Jason Kimberley, Feb 05 2012

Programs

  • Maple
    g:=sum(t^i*x^i/product(1-x^j,j=i..30),i=1..30): gser:=simplify(series(g,x=0,19)): for n from 1 to 15 do P[n]:=coeff(gser,x^n) od: for n from 1 to 15 do seq(coeff(P[n],t^j),j=1..n) od; # Emeric Deutsch, Feb 19 2006
    nmax:=13; for n from 1 to nmax do T(n,n):=1 od: for n from 1 to nmax do for k from floor(n/2)+1 to n-1 do T(n,k):=0 od: od: for n from 2 to nmax do for k from 1 to floor(n/2) do T(n,k):=sum(T(n-k,i),i=k..n-k) od: od: seq(seq(T(n,k),k=1..n), n=1..nmax); # Johannes W. Meijer, Jun 21 2010
    nmax:=13; with(combinat): for n from 1 to nmax do for k from n+1 to nmax do T(n,k):=0 od: od: for n from 1 to nmax do T(n,1):=numbpart(n-1) od: for n from 1 to nmax do T(n,n):=1 od: for n from 2 to nmax do for k from 2 to n-1 do T(n,k) := T(n-1,k-1) - T(n-k,k-1) od: od: seq(seq(T(n,k),k=1..n), n=1..nmax); # Johannes W. Meijer, Jun 21 2010
    #
    p:= (f, g)-> zip((x,y)-> x+y, f, g, 0):
    b:= proc(n, i) option remember; local h;
          h:= `if`(n=i and i>0, [0$(i-1), 1], []);
          `if`(i<1, h, p(p(h, b(n, i-1)), `if`(n b(n, n)[]:
    seq(T(n), n=1..14); # Alois P. Heinz, Mar 28 2012
  • Mathematica
    t[n_, k_] /; k<1 || k>n = 0; t[n_, n_] = 1; t[n_, k_] := t[n, k] = Sum[t[n-k, i], {i, k, n-k}]; Flatten[ Table[t[n, k], {n, 1, 14}, {k, 1, n}]] (* Jean-François Alcover, May 11 2012, after PARI *)
  • PARI
    {T(n, k) = if( k<1 || k>n, 0, if( n==k, 1, sum(i=k, n-k, T(n-k, i))))} \\ Michael Somos, Feb 06 2003
    
  • PARI
    A026794(n,k)=#select(p->p[1]==k,partitions(n,[k,n])) \\ For illustration: Creates the list of all partitions of n with smallest part equal to k. - M. F. Hasler, Jun 14 2018

Formula

T(n, k) = sum{T(n-k, i), k<=i<=n-k} for k=1, 2, ..., m, T(n, k)=0 for k=m+1, ..., n-1, where m=floor(n/2); T(n, n)=1 for n >= 1.
G.f.: G(t,x)=sum(t^i*x^i/product(1-x^j, j=i..infinity), i=1..infinity). - Emeric Deutsch, Feb 19 2006
G.f.: Sum_{k>=1} tx^k/(1-tx^k)/product(1-x^j,j=1..k-1). - Emeric Deutsch, Mar 13 2006
T(n,k) = T(n-1,k-1) - T(n-k,k-1) for n>=2 and 2<=k<=(n-1) with T(n,1) = A000041(n-1), T(n,n) = 1 for n>=1 and T(n,k) = 0 for k>n. - Johannes W. Meijer, Jun 21 2010
T(k,k) = 1 and T(n,1) = row sum (n-1); thus Meijer's 2010 formula generates the triangle without a priori reference to A000041 (the partition sequence). - Bob Selcoe, Sep 03 2016

Extensions

More terms from Emeric Deutsch, Feb 19 2006

A193870 Triangle of regions and partitions of integers (see Comments lines for definition).

Original entry on oeis.org

1, 2, 1, 3, 1, 1, 2, 0, 0, 0, 4, 2, 1, 1, 1, 3, 0, 0, 0, 0, 0, 5, 2, 1, 1, 1, 1, 1, 2, 0, 0, 0, 0, 0, 0, 0, 4, 2, 0, 0, 0, 0, 0, 0, 0, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 6, 3, 2, 2, 1, 1, 1, 1, 1, 1, 1, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 5, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 7, 3, 2, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1
Offset: 1

Views

Author

Omar E. Pol, Aug 07 2011

Keywords

Comments

Triangle T(n,k) read by rows in which, from rows 1..n, if r = T(n,k) is a record in the sequence then the set of positive integers in every row (from 1 to n) is called a “region” of r. Note that n, the number of regions of r is also the number of partitions of r. The consecutive records "r" are the natural numbers A000027. The triangle has the property that, for rows n..1, the diagonals (without the zeros) are also the partitions of r, in juxtaposed reverse-lexicographical order. Note that a record "r" is the initial term of a row if such row contains 1’s. If T(n,k) is a record in the sequence then A000041(T(n,k)) = n. Note that if T(n,k) < 2 is not the last term of the row n then T(n,k+1) = T(n,k). The union of the rows that contain 1's gives A182715.

Examples

			Triangle begins:
1,
2, 1,
3, 1, 1,
2, 0, 0, 0,
4, 2, 1, 1, 1,
3, 0, 0, 0, 0, 0,
5, 2, 1, 1, 1, 1, 1,
2, 0, 0, 0, 0, 0, 0, 0,
4, 2, 0, 0, 0, 0, 0, 0, 0,
3, 0, 0, 0, 0, 0, 0, 0, 0, 0,
6, 3, 2, 2, 1, 1, 1, 1, 1, 1, 1,
…
For n = 11 note that the row n contains the 6th record in the sequence: T(11,1) = a(56) = 6, then consider the first 11 rows of triangle. Note that the diagonals d, from d = n..1, without the zeros, are also the partitions of 6 in juxtaposed reverse-lexicographical order: [6], [3, 3], [4, 2], [2, 2, 2], [5, 1], [3, 2, 1], [4, 1, 1], [2, 2, 1, 1], [3, 1, 1, 1], [2, 1, 1, 1, 1], [1, 1, 1, 1, 1, 1]. See A026792.
		

Crossrefs

Mirror of triangle A186114. Column 1 gives A141285. Right diagonal gives A167392.

Programs

Formula

T(n,1) = A141285(n).
T(n,k) = A167392(n), if k = n.

A026804 Number of partitions of n in which the least part is odd.

Original entry on oeis.org

1, 1, 3, 3, 6, 8, 13, 16, 25, 33, 47, 61, 84, 109, 148, 189, 249, 319, 413, 522, 670, 842, 1066, 1330, 1668, 2068, 2574, 3171, 3915, 4800, 5888, 7175, 8753, 10617, 12879, 15552, 18772, 22570, 27125, 32480, 38867, 46372, 55275, 65707, 78047, 92470, 109456
Offset: 1

Views

Author

Keywords

Comments

Also number of partitions of n in which the largest part occurs an odd number of times. Example: a(5)=6 because we have [5],[4,1],[3,2],[3,1,1],[2,1,1,1] and [1,1,1,1,1] ([2,2,1] does not qualify). - Emeric Deutsch, Apr 04 2006

Examples

			a(5)=6 because we have [5],[4,1],[3,1,1],[2,2,1],[2,1,1,1] and [1,1,1,1,1] ([3,2] does not qualify).
		

Crossrefs

Cf. A046746.

Programs

  • Maple
    g:=sum(x^(2*k-1)/product(1-x^j,j=2*k-1..50),k=1..50): gser:=series(g,x=0,45): seq(coeff(gser,x,n),n=1..43); # Emeric Deutsch, Apr 04 2006
    # second Maple program:
    b:= proc(n, i) option remember; `if`(n<1 or i<1, 0, b(n, i-1)+
          `if`(n=i, irem(n, 2), 0)+`if`(i>n, 0, b(n-i, i)))
        end:
    a:= n-> b(n$2):
    seq(a(n), n=1..60);  # Alois P. Heinz, Jul 26 2015
  • Mathematica
    b[n_, i_] := b[n, i] = If[n < 1 || i < 1, 0, b[n, i - 1] + If[n == i, Mod[n, 2], 0] + If[i > n, 0, b[n - i, i]]]; a[n_] :=  b[n, n]; Table[a[n], {n, 1, 60}] (* Jean-François Alcover, Oct 09 2015, after Alois P. Heinz *)
  • PARI
    b(n, i) = if(n<1 || i<1, 0, b(n, i - 1) + if(n==i, n%2 , 0) + if(i>n, 0, b(n - i, i)));
    a(n) = b(n, n); \\ Indranil Ghosh, Jun 22 2017, after Maple code by Alois P. Heinz

Formula

G.f.: Sum_{k>=1}((-1)^(k+1)*(-1+1/Product_{i>=k} (1-x^i))). a(n) = Sum_{k=1..n}(-1)^(k+1)*A026807(n, k). - Vladeta Jovovic, Aug 26 2003
G.f.: Sum_{j>=1}(x^j/(1+x^j)/Product_{i=1..j}(1-x^i)). - Vladeta Jovovic, Aug 11 2004
G.f.: Sum_{k>=1}(x^(2k-1)/Product_{j>=2k-1}(1-x^j)). - Emeric Deutsch, Apr 04 2006
a(n) ~ exp(Pi*sqrt(2*n/3)) / (4*sqrt(3)*n) * (1 - (sqrt(3/2)/Pi + 25*Pi/(24*sqrt(6))) / sqrt(n) + (25/16 + 2929*Pi^2/6912)/n). - Vaclav Kotesovec, Jul 06 2019, extended Nov 02 2019

A186412 Sum of all parts in the n-th region of the set of partitions of j, if 1<=n<=A000041(j).

Original entry on oeis.org

1, 3, 5, 2, 9, 3, 12, 2, 6, 3, 20, 3, 7, 4, 25, 2, 6, 3, 13, 5, 4, 38, 3, 7, 4, 14, 3, 9, 5, 49, 2, 6, 3, 13, 5, 4, 23, 4, 10, 6, 5, 69, 3, 7, 4, 14, 3, 9, 5, 27, 5, 4, 15, 7, 6, 87, 2, 6, 3, 13, 5, 4, 23, 4, 10, 6, 5, 39, 3, 9, 5, 19, 4, 12, 7, 6, 123
Offset: 1

Views

Author

Omar E. Pol, Aug 12 2011

Keywords

Comments

Also triangle read by rows: T(j,k) = sum of all parts in the k-th region of the last section of the set of partitions of j. See Example section. For more information see A135010. - Omar E. Pol, Nov 26 2011
For the definition of "region" see A206437. - Omar E. Pol, Aug 19 2013

Examples

			Contribution from Omar E. Pol, Nov 26 2011 (Start):
Written as a triangle:
1;
3;
5;
2,9;
3,12;
2,6,3,20;
3,7,4,25;
2,6,3,13,5,4,38;
3,7,4,14,3,9,5,49;
2,6,3,13,5,4,23,4,10,6,5,69;
3,7,4,14,3,9,5,27,5,4,15,7,6,87;
2,6,3,13,5,4,23,4,10,6,5,39,3,9,5,19,4,12,7,6,123;
(End)
From _Omar E. Pol_, Aug 18 2013: (Start)
Illustration of initial terms (first seven regions):
.                                             _ _ _ _ _
.                                     _ _ _  |_ _ _ _ _|
.                           _ _ _ _  |_ _ _|       |_ _|
.                     _ _  |_ _ _ _|                 |_|
.             _ _ _  |_ _|     |_ _|                 |_|
.       _ _  |_ _ _|             |_|                 |_|
.   _  |_ _|     |_|             |_|                 |_|
.  |_|   |_|     |_|             |_|                 |_|
.
.   1     3       5     2         9       3          12
.
(End)
		

Crossrefs

Row sums of triangle A186114 and of triangle A193870.
Row j has length A187219(j).
Row sums give A138879.
Right border gives A046746, j >= 1.
Records give A046746, j >= 1.
Partial sums give A182244.

Programs

  • Mathematica
    lex[n_]:=DeleteCases[Sort@PadRight[Reverse /@ IntegerPartitions@n], x_ /; x==0,2];
    A186412 = {}; l = {};
    For[j = 1, j <= 50, j++,
      mx = Max@lex[j][[j]]; AppendTo[l, mx];
      For[i = j, i > 0, i--, If[l[[i]] > mx, Break[]]];
      AppendTo[A186412, Total@Take[Reverse[First /@ lex[mx]], j - i]];
      ];
    A186412  (* Robert Price, Jul 25 2020 *)

Formula

a(A000041(n)) = A046746(n).

A026807 Triangular array T read by rows: T(n,k) = number of partitions of n in which every part is >=k, for k=1,2,...,n.

Original entry on oeis.org

1, 2, 1, 3, 1, 1, 5, 2, 1, 1, 7, 2, 1, 1, 1, 11, 4, 2, 1, 1, 1, 15, 4, 2, 1, 1, 1, 1, 22, 7, 3, 2, 1, 1, 1, 1, 30, 8, 4, 2, 1, 1, 1, 1, 1, 42, 12, 5, 3, 2, 1, 1, 1, 1, 1, 56, 14, 6, 3, 2, 1, 1, 1, 1, 1, 1, 77, 21, 9, 5, 3, 2, 1, 1, 1, 1, 1, 1, 101, 24, 10, 5, 3, 2, 1, 1, 1, 1, 1, 1, 1, 135, 34, 13
Offset: 1

Views

Author

Keywords

Comments

T(n,g) is also the number of not necessarily connected 2-regular graphs with girth at least g: the part i corresponds to the i-cycle; addition of integers corresponds to disconnected union of cycles. - Jason Kimberley, Feb 05 2012

Examples

			Sum_{k>=1} y^k*(-1+1/Product_{i>=0} (1-x^(k+i))) = y*x+(2*y+y^2)*x^2+(3*y+y^2+y^3)*x^3+(5*y+2*y^2+y^3+y^4)*x^4+(7*y+2*y^2+y^3+y^4+y^5)*x^5+...
Triangle starts:  - _Jason Kimberley_, Feb 05 2012
1;
2, 1;
3, 1, 1;
5, 2, 1, 1;
7, 2, 1, 1, 1;
11, 4, 2, 1, 1, 1;
15, 4, 2, 1, 1, 1, 1;
22, 7, 3, 2, 1, 1, 1, 1;
30, 8, 4, 2, 1, 1, 1, 1, 1;
42, 12, 5, 3, 2, 1, 1, 1, 1, 1;
56, 14, 6, 3, 2, 1, 1, 1, 1, 1, 1;
77, 21, 9, 5, 3, 2, 1, 1, 1, 1, 1, 1;
101, 24, 10, 5, 3, 2, 1, 1, 1, 1, 1, 1, 1;
From _Tilman Piesk_, Feb 20 2016: (Start)
n = 12, k = 4, t = A000217(k-1) = 6
vp = A000041(n..n-t) = A000041(12..6) = (77, 56, 42, 30, 22, 15, 11)
vc = A231599(k-1, 0..t) = A231599(3, 0..6) = (1,-1,-1, 0, 1, 1,-1)
T(12, 4) = vp * transpose(vc) = 77-56-42+22+15-11 = 5
(End)
		

Crossrefs

Row sums give A046746.
Cf. A026835.
Cf. A026794.
Cf. A231599.
Not necessarily connected 2-regular graphs with girth at least g [partitions into parts >= g]: this sequence (triangle); columns of this sequence: 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). For g >= 3, girth at least g implies no loops or parallel edges. - Jason Kimberley, Feb 05 2012
Not necessarily connected 2-regular simple 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). - Jason Kimberley, Feb 05 2012
Cf. A002260.

Programs

  • Haskell
    import Data.List (tails)
    a026807 n k = a026807_tabl !! (n-1) !! (k-1)
    a026807_row n = a026807_tabl !! (n-1)
    a026807_tabl = map
       (\row -> map (p $ last row) $ init $ tails row) a002260_tabl
       where p 0  _ = 1
             p _ [] = 0
             p m ks'@(k:ks) = if m < k then 0 else p (m - k) ks' + p m ks
    -- Reinhard Zumkeller, Dec 01 2012
    
  • Maple
    T:= proc(n, k) option remember;
          `if`(k<1 or k>n, 0, `if`(n=k, 1, T(n, k+1) +T(n-k, k)))
        end:
    seq(seq(T(n, k), k=1..n), n=1..14); # Alois P. Heinz, Mar 28 2012
  • Mathematica
    T[n_, k_] := T[n, k] = If[ k<1 || k>n, 0, If[n == k, 1, T[n, k+1] + T[n-k, k]]]; Table [Table[ T[n, k], {k, 1, n}], {n, 1, 14}] // Flatten (* Jean-François Alcover, Jan 28 2015, after Alois P. Heinz *)
  • Python
    from see_there import a231599_row  # A231599
    from sympy.ntheory import npartitions  # A000041
    def a026807(n, k):
        if k > n:
            return 0
        elif k > n/2:
            return 1
        else:
            vc = a231599_row(k-1)
            t = len(vc)
            vp_range = range(n-t, n+1)
            vp_range = vp_range[::-1]  # reverse
            r = 0
            for i in range(0, t):
                r += vc[i] * npartitions(vp_range[i])
            return r
    # Tilman Piesk, Feb 21 2016

Formula

T(n,1)=A000041(n), T(n,2)=A002865(n) for n>1, T(n,3)=A008483(n) for n>2, T(n,4)=A008484(n) for n>3.
G.f.: Sum_{k>=1} y^k*(-1+1/Product_{i>=0} (1-x^(k+i))). - Vladeta Jovovic, Jun 22 2003
T(n, k) = T(n, k+1) + T(n-k, k), T(n, k) = 1 if n/2 < k <= n. - Franklin T. Adams-Watters, Jan 24 2005; Tilman Piesk, Feb 20 2016
T(n, k) = A000041(n..n-t) * transpose(A231599(k-1, 0..t)) with t = A000217(k-1). - Tilman Piesk, Feb 20 2016
Equals A026794 * A000012 as infinite lower triangular matrices. - Gary W. Adamson, Jan 31 2008

A182709 Sum of the emergent parts of the partitions of n.

Original entry on oeis.org

0, 0, 0, 2, 3, 11, 14, 33, 45, 81, 109, 185, 237, 372, 490, 715, 928, 1326, 1693, 2348, 2998, 4032, 5119, 6795, 8530, 11132, 13952, 17927, 22314, 28417, 35126, 44279, 54532, 68062, 83422, 103427, 126063, 155207, 188506, 230547, 278788, 339223, 408482
Offset: 1

Views

Author

Omar E. Pol, Nov 28 2010, Nov 29 2010

Keywords

Comments

Here the "emergent parts" of the partitions of n are defined to be the parts (with multiplicity) of all the partitions that do not contain "1" as a part, removed by one copy of the smallest part of every partition. Note that these parts are located in the head of the last section of the set of partitions of n. For more information see A182699.
Also total sum of parts of the regions that do not contain 1 as a part in the last section of the set of partitions of n (Cf. A083751, A187219). - Omar E. Pol, Mar 04 2012

Examples

			For n=7 the partitions of 7 that do not contain "1" as a part are
7
4 + 3
5 + 2
3 + 2 + 2
Then remove one copy of the smallest part of every partition. The rest are the emergent parts:
.,
4, .
5, .
3, 2, .
The sum of these parts is 4 + 5 + 3 + 2 = 14, so a(7)=14.
For n=10 the illustration in the link shows the location of the emergent parts (colored yellow and green) and the location of the filler parts (colored blue) in the last section of the set of partitions of 10.
		

Crossrefs

Programs

  • Maple
    b:= proc(n, i) option remember;
          if n<0 then 0
        elif n=0 then 1
        elif i<2 then 0
        else b(n, i-1) +b(n-i, i)
          fi
        end:
    c:= proc(n, i, k) option remember;
          if n<0 then 0
        elif n=0 then k
        elif i<2 then 0
        else c(n, i-1, k) +c(n-i, i, i)
          fi
        end:
    a:= n-> n*b(n, n) - c(n, n, 0):
    seq(a(n), n=1..40);  #  Alois P. Heinz, Dec 01 2010
  • Mathematica
    f[n_]:=Total[Flatten[Most/@Select[IntegerPartitions[n],!MemberQ[#,1]&]]]; Table[f[i],{i,50}] (* Harvey P. Dale, Dec 28 2010 *)
    b[n_, i_] := b[n, i] = Which[n<0, 0, n==0, 1, i<2, 0, True, b[n, i-1] + b[n - i, i]]; c[n_, i_, k_] := c[n, i, k] = Which[n<0, 0, n==0, k, i<2, 0, True, c[n, i-1, k] + c[n-i, i, i]]; a[n_] := n*b[n, n] - c[n, n, 0]; Table[a[n], {n, 1, 40}] (* Jean-François Alcover, Oct 08 2015, after Alois P. Heinz *)

Formula

a(n) = A138880(n) - A182708(n).
a(n) = A066186(n) - A066186(n-1) - A046746(n) = A138879(n) - A046746(n). - Omar E. Pol, Aug 01 2013
a(n) ~ Pi * exp(Pi*sqrt(2*n/3)) / (12*sqrt(2*n)) * (1 - (3*sqrt(3/2)/Pi + 13*Pi/(24*sqrt(6)))/sqrt(n)). - Vaclav Kotesovec, Jan 03 2019, extended Jul 05 2019

Extensions

More terms from Alois P. Heinz, Dec 01 2010

A092309 Sum of smallest parts (counted with multiplicity) of all partitions of n.

Original entry on oeis.org

1, 4, 7, 15, 19, 39, 46, 80, 106, 160, 201, 318, 390, 554, 729, 998, 1262, 1727, 2168, 2894, 3670, 4749, 5963, 7737, 9635, 12232, 15257, 19206, 23727, 29723, 36509, 45296, 55512, 68292, 83298, 102079, 123805, 150697, 182254, 220790, 265766
Offset: 1

Views

Author

Vladeta Jovovic, Feb 16 2004

Keywords

Examples

			Partitions of 4 are: [1,1,1,1], [1,1,2], [2,2], [1,3], [4]; thus a(4)=4*1+2*1+2*2+1*1+1*4=15.
		

Crossrefs

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(irem(n, i)=0, n, 0)
           +`if`(i>1, add(b(n-i*j, i-1), j=0..(n-1)/i), 0)
        end:
    a:= n-> b(n$2):
    seq(a(n), n=1..50);  # Alois P. Heinz, Feb 04 2016
  • Mathematica
    ss[n_]:=Module[{m=Min[n]},Select[n,#==m&]]; Table[Total[Flatten[ss/@ IntegerPartitions[n]]],{n,50}] (* Harvey P. Dale, Dec 16 2013 *)
    b[n_, i_] := b[n, i] = If[Mod[n, i] == 0, n, 0] + If[i > 1, Sum[b[n - i*j, i - 1], {j, 0, (n - 1)/i}], 0]; a[n_] := b[n, n]; Table[a[n], {n, 1, 50}] (* Jean-François Alcover, Aug 29 2016, after Alois P. Heinz *)

Formula

G.f.: Sum(n*x^n/(1-x^n)*Product(1/(1-x^k), k = n .. infinity), n = 1 .. infinity).
a(n) ~ sqrt(2) * exp(Pi*sqrt(2*n/3)) / (4*Pi*sqrt(n)). - Vaclav Kotesovec, Jul 06 2019

Extensions

More terms from Pab Ter (pabrlos(AT)yahoo.com), May 25 2004

A097979 Total number of largest parts in all compositions of n.

Original entry on oeis.org

1, 3, 6, 12, 23, 46, 91, 183, 367, 737, 1478, 2962, 5928, 11858, 23707, 47384, 94698, 189260, 378277, 756160, 1511730, 3022672, 6044472, 12088395, 24177600, 48359695, 96732370, 193495606, 387057584, 774248858, 1548754115, 3097980230, 6196797193, 12395022288
Offset: 1

Views

Author

Vladeta Jovovic, Sep 07 2004

Keywords

Comments

Also number of compositions of n+1 with unique largest part. - Vladeta Jovovic, Apr 03 2005

Crossrefs

Column k=1 of A238341.

Programs

  • Mathematica
    nn=32; Drop[CoefficientList[Series[Sum[x^j/(1 - (x - x^(j + 1))/(1 - x))^2, {j, 1, nn}], {x, 0, nn}], x], 1] (* Geoffrey Critzer, Mar 31 2014 *)
    b[n_, p_, i_] := b[n, p, i] = If[n == 0, p!, If[i<1, 0, Sum[b[n-i*j, p+j, i-1]/j!, {j, 0, n/i}]]]; a[n_, k_] := Sum[b[n-i*k, k, i-1]/k!, {i, 1, n/k}]; a[0, 0] = 1; a[, 0] = 0; a[n] := a[n+1, 1]; Table[a[n], {n, 1, 32}] (* Jean-François Alcover, Feb 10 2015, after A238341 *)
  • PARI
    { b(t)=local(r);sum(k=1,t, forstep(s=t%k,t-k,k,u=(t-k-s)\k;r+=binomial(-2,s)*(-2)^(s-u)*binomial(s,u)));r }
    { a(n)=b(n)-2*b(n-1)+b(n-2) } \\ Max Alekseyev, Apr 16 2005

Formula

G.f.: (1-x)^2 * Sum_{k>=1} x^k/(1-2*x+x^(k+1))^2.
a(n) ~ 2^(n-1)/log(2). - Vaclav Kotesovec, Apr 30 2014

Extensions

More terms from Max Alekseyev, Apr 16 2005
Showing 1-10 of 47 results. Next