login
Search: a089353 -id:a089353
     Sort: relevance | references | number | modified | created      Format: long | short | data
The positive integers. Also called the natural numbers, the whole numbers or the counting numbers, but these terms are ambiguous.
(Formerly M0472 N0173)
+10
2160
1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77
OFFSET
1,2
COMMENTS
For some authors, the terms "natural numbers" and "counting numbers" include 0, i.e., refer to the nonnegative integers A001477; the term "whole numbers" frequently also designates the whole set of (signed) integers A001057.
a(n) is smallest positive integer which is consistent with sequence being monotonically increasing and satisfying a(a(n)) = n (cf. A007378).
Inverse Euler transform of A000219.
The rectangular array having A000027 as antidiagonals is the dispersion of the complement of the triangular numbers, A000217 (which triangularly form column 1 of this array). The array is also the transpose of A038722. - Clark Kimberling, Apr 05 2003
For nonzero x, define f(n) = floor(nx) - floor(n/x). Then f=A000027 if and only if x=tau or x=-tau. - Clark Kimberling, Jan 09 2005
Numbers of form (2^i)*k for odd k (i.e., n = A006519(n)*A000265(n)); thus n corresponds uniquely to an ordered pair (i,k) where i=A007814, k=A000265 (with A007814(2n)=A001511(n), A007814(2n+1)=0). - Lekraj Beedassy, Apr 22 2006
If the offset were changed to 0, we would have the following pattern: a(n)=binomial(n,0) + binomial(n,1) for the present sequence (number of regions in 1-space defined by n points), A000124 (number of regions in 2-space defined by n straight lines), A000125 (number of regions in 3-space defined by n planes), A000127 (number of regions in 4-space defined by n hyperplanes), A006261, A008859, A008860, A008861, A008862 and A008863, where the last six sequences are interpreted analogously and in each "... by n ..." clause an offset of 0 has been assumed, resulting in a(0)=1 for all of them, which corresponds to the case of not cutting with a hyperplane at all and therefore having one region. - Peter C. Heinig (algorithms(AT)gmx.de), Oct 19 2006
Define a number of points on a straight line to be in general arrangement when no two points coincide. Then these are the numbers of regions defined by n points in general arrangement on a straight line, when an offset of 0 is assumed. For instance, a(0)=1, since using no point at all leaves one region. The sequence satisfies the recursion a(n) = a(n-1) + 1. This has the following geometrical interpretation: Suppose there are already n-1 points in general arrangement, thus defining the maximal number of regions on a straight line obtainable by n-1 points, and now one more point is added in general arrangement. Then it will coincide with no other point and act as a dividing wall thereby creating one new region in addition to the a(n-1)=(n-1)+1=n regions already there, hence a(n)=a(n-1)+1. Cf. the comments on A000124 for an analogous interpretation. - Peter C. Heinig (algorithms(AT)gmx.de), Oct 19 2006
The sequence a(n)=n (for n=1,2,3) and a(n)=n+1 (for n=4,5,...) gives to the rank (minimal cardinality of a generating set) for the semigroup I_n\S_n, where I_n and S_n denote the symmetric inverse semigroup and symmetric group on [n]. - James East, May 03 2007
The sequence a(n)=n (for n=1,2), a(n)=n+1 (for n=3) and a(n)=n+2 (for n=4,5,...) gives the rank (minimal cardinality of a generating set) for the semigroup PT_n\T_n, where PT_n and T_n denote the partial transformation semigroup and transformation semigroup on [n]. - James East, May 03 2007
"God made the integers; all else is the work of man." This famous quotation is a translation of "Die ganzen Zahlen hat der liebe Gott gemacht, alles andere ist Menschenwerk," spoken by Leopold Kronecker in a lecture at the Berliner Naturforscher-Versammlung in 1886. Possibly the first publication of the statement is in Heinrich Weber's "Leopold Kronecker," Jahresberichte D.M.V. 2 (1893) 5-31. - Clark Kimberling, Jul 07 2007
Binomial transform of A019590, inverse binomial transform of A001792. - Philippe Deléham, Oct 24 2007
Writing A000027 as N, perhaps the simplest one-to-one correspondence between N X N and N is this: f(m,n) = ((m+n)^2 - m - 3n + 2)/2. Its inverse is given by I(k)=(g,h), where g = k - J(J-1)/2, h = J + 1 - g, J = floor((1 + sqrt(8k - 7))/2). Thus I(1)=(1,1), I(2)=(1,2), I(3)=(2,1) and so on; the mapping I fills the first-quadrant lattice by successive antidiagonals. - Clark Kimberling, Sep 11 2008
a(n) is also the mean of the first n odd integers. - Ian Kent, Dec 23 2008
Equals INVERTi transform of A001906, the even-indexed Fibonacci numbers starting (1, 3, 8, 21, 55, ...). - Gary W. Adamson, Jun 05 2009
These are also the 2-rough numbers: positive integers that have no prime factors less than 2. - Michael B. Porter, Oct 08 2009
Totally multiplicative sequence with a(p) = p for prime p. Totally multiplicative sequence with a(p) = a(p-1) + 1 for prime p. - Jaroslav Krizek, Oct 18 2009
Triangle T(k,j) of natural numbers, read by rows, with T(k,j) = binomial(k,2) + j = (k^2-k)/2 + j where 1 <= j <= k. In other words, a(n) = n = binomial(k,2) + j where k is the largest integer such that binomial(k,2) < n and j = n - binomial(k,2). For example, T(4,1)=7, T(4,2)=8, T(4,3)=9, and T(4,4)=10. Note that T(n,n)=A000217(n), the n-th triangular number. - Dennis P. Walsh, Nov 19 2009
Hofstadter-Conway-like sequence (see A004001): a(n) = a(a(n-1)) + a(n-a(n-1)) with a(1) = 1, a(2) = 2. - Jaroslav Krizek, Dec 11 2009
a(n) is also the dimension of the irreducible representations of the Lie algebra sl(2). - Leonid Bedratyuk, Jan 04 2010
Floyd's triangle read by rows. - Paul Muljadi, Jan 25 2010
Number of numbers between k and 2k where k is an integer. - Giovanni Teofilatto, Mar 26 2010
Generated from a(2n) = r*a(n), a(2n+1) = a(n) + a(n+1), r = 2; in an infinite set, row 2 of the array shown in A178568. - Gary W. Adamson, May 29 2010
1/n = continued fraction [n]. Let barover[n] = [n,n,n,...] = 1/k. Then k - 1/k = n. Example: [2,2,2,...] = (sqrt(2) - 1) = 1/k, with k = (sqrt(2) + 1). Then 2 = k - 1/k. - Gary W. Adamson, Jul 15 2010
Number of n-digit numbers the binary expansion of which contains one run of 1's. - Vladimir Shevelev, Jul 30 2010
From Clark Kimberling, Jan 29 2011: (Start)
Let T denote the "natural number array A000027":
1 2 4 7 ...
3 5 8 12 ...
6 9 13 18 ...
10 14 19 25 ...
T(n,k) = n+(n+k-2)*(n+k-1)/2. See A185787 for a list of sequences based on T, such as rows, columns, diagonals, and sub-arrays. (End)
The Stern polynomial B(n,x) evaluated at x=2. See A125184. - T. D. Noe, Feb 28 2011
The denominator in the Maclaurin series of log(2), which is 1 - 1/2 + 1/3 - 1/4 + .... - Mohammad K. Azarian, Oct 13 2011
As a function of Bernoulli numbers B_n (cf. A027641: (1, -1/2, 1/6, 0, -1/30, 0, 1/42, ...)): let V = a variant of B_n changing the (-1/2) to (1/2). Then triangle A074909 (the beheaded Pascal's triangle) * [1, 1/2, 1/6, 0, -1/30, ...] = the vector [1, 2, 3, 4, 5, ...]. - Gary W. Adamson, Mar 05 2012
Number of partitions of 2n+1 into exactly two parts. - Wesley Ivan Hurt, Jul 15 2013
Integers n dividing u(n) = 2u(n-1) - u(n-2); u(0)=0, u(1)=1 (Lucas sequence A001477). - Thomas M. Bridge, Nov 03 2013
For this sequence, the generalized continued fraction a(1)+a(1)/(a(2)+a(2)/(a(3)+a(3)/(a(4)+...))), evaluates to 1/(e-2) = A194807. - Stanislav Sykora, Jan 20 2014
Engel expansion of e-1 (A091131 = 1.71828...). - Jaroslav Krizek, Jan 23 2014
a(n) is the number of permutations of length n simultaneously avoiding 213, 231 and 321 in the classical sense which are breadth-first search reading words of increasing unary-binary trees. For more details, see the entry for permutations avoiding 231 at A245898. - Manda Riehl, Aug 05 2014
a(n) is also the number of permutations simultaneously avoiding 213, 231 and 321 in the classical sense which can be realized as labels on an increasing strict binary tree with 2n-1 nodes. See A245904 for more information on increasing strict binary trees. - Manda Riehl, Aug 07 2014
a(n) = least k such that 2*Pi - Sum_{h=1..k} 1/(h^2 - h + 3/16) < 1/n. - Clark Kimberling, Sep 28 2014
a(n) = least k such that Pi^2/6 - Sum_{h=1..k} 1/h^2 < 1/n. - Clark Kimberling, Oct 02 2014
Determinants of the spiral knots S(2,k,(1)). a(k) = det(S(2,k,(1))). These knots are also the torus knots T(2,k). - Ryan Stees, Dec 15 2014
As a function, the restriction of the identity map on the nonnegative integers {0,1,2,3...}, A001477, to the positive integers {1,2,3,...}. - M. F. Hasler, Jan 18 2015
See also A131685(k) = smallest positive number m such that c(i) = m (i^1 + 1) (i^2 + 2) ... (i^k+ k) / k! takes integral values for all i>=0: For k=1, A131685(k)=1, which implies that this is a well defined integer sequence. - Alexander R. Povolotsky, Apr 24 2015
a(n) is the number of compositions of n+2 into n parts avoiding the part 2. - Milan Janjic, Jan 07 2016
Does not satisfy Benford's law [Berger-Hill, 2017] - N. J. A. Sloane, Feb 07 2017
Parametrization for the finite multisubsets of the positive integers, where, for p_j the j-th prime, n = Product_{j} p_j^(e_j) corresponds to the multiset containing e_j copies of j ('Heinz encoding' -- see A056239, A003963, A289506, A289507, A289508, A289509). - Christopher J. Smyth, Jul 31 2017
The arithmetic function v_1(n,1) as defined in A289197. - Robert Price, Aug 22 2017
For n >= 3, a(n)=n is the least area that can be obtained for an irregular octagon drawn in a square of n units side, whose sides are parallel to the axes, with 4 vertices that coincide with the 4 vertices of the square, and the 4 remaining vertices having integer coordinates. See Affaire de Logique link. - Michel Marcus, Apr 28 2018
a(n+1) is the order of rowmotion on a poset defined by a disjoint union of chains of length n. - Nick Mayers, Jun 08 2018
Number of 1's in n-th generation of 1-D Cellular Automata using Rules 50, 58, 114, 122, 178, 186, 206, 220, 238, 242, 250 or 252 in the Wolfram numbering scheme, started with a single 1. - Frank Hollstein, Mar 25 2019
(1, 2, 3, 4, 5, ...) is the fourth INVERT transform of (1, -2, 3, -4, 5, ...). - Gary W. Adamson, Jul 15 2019
REFERENCES
T. M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1976, page 1.
T. M. Apostol, Modular Functions and Dirichlet Series in Number Theory, Springer-Verlag, 1990, page 25.
John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See p. 22.
W. Fulton and J. Harris, Representation theory: a first course, (1991), page 149. [From Leonid Bedratyuk, Jan 04 2010]
I. S. Gradstein and I. M. Ryshik, Tables of series, products, and integrals, Volume 1, Verlag Harri Deutsch, 1981.
R. E. Schwartz, You Can Count on Monsters: The First 100 numbers and Their Characters, A. K. Peters and MAA, 2010.
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).
LINKS
N. J. A. Sloane, Table of n, a(n) for n = 1..500000 [a large file]
Archimedes Laboratory, What's special about this number?
Affaire de Logique, Pick et Pick et Colegram (in French), No. 1051, 18-04-2018.
Paul Barry, A Catalan Transform and Related Transformations on Integer Sequences, Journal of Integer Sequences, Vol. 8 (2005), Article 05.4.5.
James Barton, The Numbers
A. Berger and T. P. Hill, What is Benford's Law?, Notices, Amer. Math. Soc., 64:2 (2017), 132-134.
A. Breiland, L. Oesper, and L. Taalman, p-Coloring classes of torus knots, Online Missouri J. Math. Sci., 21 (2009), 120-126.
N. Brothers, S. Evans, L. Taalman, L. Van Wyk, D. Witczak, and C. Yarnall, Spiral knots, Missouri J. of Math. Sci., 22 (2010).
C. K. Caldwell, Prime Curios
Case and Abiessu, interesting number
M. DeLong, M. Russell, and J. Schrock, Colorability and determinants of T(m,n,r,s) twisted torus knots for n equiv. +/-1(mod m), Involve, Vol. 8 (2015), No. 3, 361-384.
P. Flajolet and R. Sedgewick, Analytic Combinatorics, 2009; see page 371.
Robert R. Forslund, A logical alternative to the existing positional number system, Southwest Journal of Pure and Applied Mathematics, Vol. 1 1995 pp. 27-29.
Kival Ngaokrajang, Illustration about relation to many other sequences, when the sequence is considered as a triangular table read by its antidiagonals. Additional illustrations when the sequence is considered as a centered triangular table read by rows.
Leonardo of Pisa [Leonardo Pisano], Illustration of initial terms, from Liber Abaci [The Book of Calculation], 1202 (photo by David Singmaster).
Feihu Liu, Guoce Xin, and Chen Zhang, Ehrhart Polynomials of Order Polytopes: Interpreting Combinatorial Sequences on the OEIS, arXiv:2412.18744 [math.CO], 2024. See pp. 15, 24.
G. Pfeiffer, Counting Transitive Relations, Journal of Integer Sequences, Vol. 7 (2004), Article 04.3.2.
J. Striker, Dynamical Algebraic Combinatorics: Promotion, Rowmotion, and Resonance, Notices of the AMS, June/July 2017, pp. 543-549.
G. Villemin's Almanac of Numbers, NOMBRES en BREF (in French)
FORMULA
a(2k+1) = A005408(k), k >= 0, a(2k) = A005843(k), k >= 1.
Multiplicative with a(p^e) = p^e. - David W. Wilson, Aug 01 2001
Another g.f.: Sum_{n>0} phi(n)*x^n/(1-x^n) (Apostol).
When seen as an array: T(k, n) = n+1 + (k+n)*(k+n+1)/2. Main diagonal is 2n*(n+1)+1 (A001844), antidiagonal sums are n*(n^2+1)/2 (A006003). - Ralf Stephan, Oct 17 2004
Dirichlet generating function: zeta(s-1). - Franklin T. Adams-Watters, Sep 11 2005
G.f.: x/(1-x)^2. E.g.f.: x*exp(x). a(n)=n. a(-n)=-a(n).
Series reversion of g.f. A(x) is x*C(-x)^2 where C(x) is the g.f. of A000108. - Michael Somos, Sep 04 2006
G.f. A(x) satisfies 0 = f(A(x), A(x^2)) where f(u, v) = u^2 - v - 4*u*v. - Michael Somos, Oct 03 2006
Convolution of A000012 (the all-ones sequence) with itself. - Tanya Khovanova, Jun 22 2007
a(n) = 2*a(n-1)-a(n-2); a(1)=1, a(2)=2. a(n) = 1+a(n-1). - Philippe Deléham, Nov 03 2008
a(n) = A000720(A000040(n)). - Juri-Stepan Gerasimov, Nov 29 2009
a(n+1) = Sum_{k=0..n} A101950(n,k). - Philippe Deléham, Feb 10 2012
a(n) = Sum_{d | n} phi(d) = Sum_{d | n} A000010(d). - Jaroslav Krizek, Apr 20 2012
G.f.: x * Product_{j>=0} (1+x^(2^j))^2 = x * (1+2*x+x^2) * (1+2*x^2+x^4) * (1+2*x^4+x^8) * ... = x + 2x^2 + 3x^3 + ... . - Gary W. Adamson, Jun 26 2012
a(n) = det(binomial(i+1,j), 1 <= i,j <= n). - Mircea Merca, Apr 06 2013
E.g.f.: x*E(0), where E(k) = 1 + 1/(x - x^3/(x^2 + (k+1)/E(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Aug 03 2013
From Wolfdieter Lang, Oct 09 2013: (Start)
a(n) = Product_{k=1..n-1} 2*sin(Pi*k/n), n > 1.
a(n) = Product_{k=1..n-1} (2*sin(Pi*k/(2*n)))^2, n > 1.
These identities are used in the calculation of products of ratios of lengths of certain lines in a regular n-gon. For the first identity see the Gradstein-Ryshik reference, p. 62, 1.392 1., bringing the first factor there to the left hand side and taking the limit x -> 0 (L'Hôpital). The second line follows from the first one. Thanks to Seppo Mustonen who led me to consider n-gon lengths products. (End)
a(n) = Sum_{j=0..k} (-1)^(j-1)*j*binomial(n,j)*binomial(n-1+k-j,k-j), k>=0. - Mircea Merca, Jan 25 2014
a(n) = A052410(n)^A052409(n). - Reinhard Zumkeller, Apr 06 2014
a(n) = Sum_{k=1..n^2+2*n} 1/(sqrt(k)+sqrt(k+1)). - Pierre CAMI, Apr 25 2014
a(n) = floor(1/sin(1/n)) = floor(cot(1/(n+1))) = ceiling(cot(1/n)). - Clark Kimberling, Oct 08 2014
a(n) = floor(1/(log(n+1)-log(n))). - Thomas Ordowski, Oct 10 2014
a(k) = det(S(2,k,1)). - Ryan Stees, Dec 15 2014
a(n) = 1/(1/(n+1) + 1/(n+1)^2 + 1/(n+1)^3 + ...). - Pierre CAMI, Jan 22 2015
a(n) = Sum_{m=0..n-1} Stirling1(n-1,m)*Bell(m+1), for n >= 1. This corresponds to Bell(m+1) = Sum_{k=0..m} Stirling2(m, k)*(k+1), for m >= 0, from the fact that Stirling2*Stirling1 = identity matrix. See A048993, A048994 and A000110. - Wolfdieter Lang, Feb 03 2015
a(n) = Sum_{k=1..2n-1}(-1)^(k+1)*k*(2n-k). In addition, surprisingly, a(n) = Sum_{k=1..2n-1}(-1)^(k+1)*k^2*(2n-k)^2. - Charlie Marion, Jan 05 2016
G.f.: x/(1-x)^2 = (x * r(x) *r(x^3) * r(x^9) * r(x^27) * ...), where r(x) = (1 + x + x^2)^2 = (1 + 2x + 3x^2 + 2x^3 + x^4). - Gary W. Adamson, Jan 11 2017
a(n) = floor(1/(Pi/2-arctan(n))). - Clark Kimberling, Mar 11 2020
a(n) = Sum_{d|n} mu(n/d)*sigma(d). - Ridouane Oudra, Oct 03 2020
a(n) = Sum_{k=1..n} phi(gcd(n,k))/phi(n/gcd(n,k)). - Richard L. Ollerton, May 09 2021
a(n) = S(n-1, 2), with the Chebyshev S-polynomials A049310. - Wolfdieter Lang, Mar 09 2023
From Peter Bala, Nov 02 2024: (Start)
For positive integer m, a(n) = (1/m)* Sum_{k = 1..2*m*n-1} (-1)^(k+1) * k * (2*m*n - k) = (1/m) * Sum_{k = 1..2*m*n-1} (-1)^(k+1) * k^2 * (2*m*n - k)^2 (the case m = 1 is given above).
a(n) = Sum_{k = 0..3*n} (-1)^(n+k+1) * k * binomial(3*n+k, 2*k). (End)
MAPLE
A000027 := n->n; seq(A000027(n), n=1..100);
MATHEMATICA
Range@ 77 (* Robert G. Wilson v, Mar 31 2015 *)
PROG
(Magma) [ n : n in [1..100]];
(PARI) {a(n) = n};
(R) 1:100
(Shell) seq 1 100
(Haskell)
a000027 = id
a000027_list = [1..] -- Reinhard Zumkeller, May 07 2012
(Maxima) makelist(n, n, 1, 30); /* Martin Ettl, Nov 07 2012 */
(Python)
def A000027(n): return n # Chai Wah Wu, May 09 2022
(Julia) print([n for n in 1:280]) # Paul Muljadi, Apr 09 2024
(Perl) print join(", ", 1..280) # Paul Muljadi, May 29 2024
CROSSREFS
A001477 = nonnegative numbers.
Partial sums of A000012.
Cf. A026081 = integers in reverse alphabetical order in U.S. English, A107322 = English name for number and its reverse have the same number of letters, A119796 = zero through ten in alphabetical order of English reverse spelling, A005589, etc. Cf. A185787 (includes a list of sequences based on the natural number array A000027).
Cf. Boustrophedon transforms: A000737, A231179;
Cf. A038722 (mirrored when seen as triangle), A056011 (boustrophedon).
Cf. A048993, A048994, A000110 (see the Feb 03 2015 formula).
KEYWORD
core,nonn,easy,mult,tabl,changed
EXTENSIONS
Links edited by Daniel Forgues, Oct 07 2009.
STATUS
approved
Number of plane partitions (or planar partitions) of n.
(Formerly M2566 N1016)
+10
279
1, 1, 3, 6, 13, 24, 48, 86, 160, 282, 500, 859, 1479, 2485, 4167, 6879, 11297, 18334, 29601, 47330, 75278, 118794, 186475, 290783, 451194, 696033, 1068745, 1632658, 2483234, 3759612, 5668963, 8512309, 12733429, 18974973, 28175955, 41691046, 61484961, 90379784, 132441995, 193487501, 281846923
OFFSET
0,3
COMMENTS
Two-dimensional partitions of n in which no row or column is longer than the one before it (compare A001970). E.g., a(4) = 13:
4.31.3.22.2.211.21..2.1111.111.11.11.1 but not 2
.....1....2.....1...1......1...11.1..1........ 11
....................1.............1..1
.....................................1
In the above, one also must require that rows & columns are nondecreasing, e.g., [1,1; 2] is also forbidden (which implies that row and column lengths are nondecreasing, if empty cells are identified with cells filled with 0's). - M. F. Hasler, Sep 22 2018
Can also be regarded as number of "safe pilings" of cubes in the corner of a room: the height should not increase away from the corner. - Wouter Meeussen
Also number of partitions of n objects of 2 colors, each part containing at least one black object; see example. - Christian G. Bower, Jan 08 2004
Number of partitions of n into 1 type of part 1, 2 types of part 2, ..., k types of part k. E.g., n=3 gives 111, 12, 12', 3, 3', 3''. - Jon Perry, May 27 2004
The bijection between the partitions in the two preceding comments goes by identifying a part with k black objects with a part of type k. - David Scambler and Joerg Arndt, May 01 2013
Can also be regarded as the number of Jordan canonical forms for an n X n matrix. (I.e., a 5 X 5 matrix has 24 distinct Jordan canonical forms, dependent on the algebraic and geometric multiplicity of each eigenvalue.) - Aaron Gable (agable(AT)hmc.edu), May 26 2009
(1/n) * convolution product of n terms * A001157 (sum of squares of divisors of n): (1, 5, 10, 21, 26, 50, 50, 85, ...) = a(n). As shown by [Bressoud, p. 12]: 1/6 * [1*24 + 5*13 + 10*6 + 21*3 + 26*1 + 50*1] = 288/6 = 48. - Gary W. Adamson, Jun 13 2009
Convolved with the aerated version (1, 0, 1, 0, 3, 0, 6, 0, 13, ...) = A026007: (1, 1, 2, 5, 8, 16, 28, 49, 83, ...). - Gary W. Adamson, Jun 13 2009
Starting with offset 1 = row sums of triangle A162453. - Gary W. Adamson, Jul 03 2009
Unfortunately, Wright's formula is also incomplete in the paper by G. Almkvist: "Asymptotic formulas and generalized Dedekind sums", p. 344, (the denominator should have sqrt(3*Pi) not sqrt(Pi)). This error was already corrected in the paper by Steven Finch: "Integer Partitions". - Vaclav Kotesovec, Aug 17 2015
Also the number of non-isomorphic weight-n chains of multisets whose dual is also a chain of multisets. The dual of a multiset partition has, for each vertex, one block consisting of the indices (or positions) of the blocks containing that vertex, counted with multiplicity. The weight of a multiset partition is the sum of sizes of its parts. - Gus Wiseman, Sep 25 2018
REFERENCES
G. Almkvist, The differences of the number of plane partitions, Manuscript, circa 1991.
G. E. Andrews, The Theory of Partitions, Addison-Wesley, 1976, p. 241.
D. M. Bressoud, Proofs and Confirmations, Camb. Univ. Press, 1999; pp(n) on p. 10.
Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 575.
L. Carlitz, Generating functions and partition problems, pp. 144-169 of A. L. Whiteman, ed., Theory of Numbers, Proc. Sympos. Pure Math., 8 (1965). Amer. Math. Soc., see p. 145, eq. (1.6).
I. P. Goulden and D. M. Jackson, Combinatorial Enumeration, Wiley, N.Y., 1983, (5.4.5).
P. A. MacMahon, Memoir on the theory of partitions of numbers - Part VI, Phil. Trans. Royal Soc., 211 (1912), 345-373.
P. A. MacMahon, Combinatory Analysis. Cambridge Univ. Press, London and New York, Vol. 1, 1915 and Vol. 2, 1916; see vol. 2, p 332.
P. A. MacMahon, The connexion between the sum of the squares of the divisors and the number of partitions of a given number, Messenger Math., 54 (1924), 113-116. Collected Papers, MIT Press, 1978, Vol. I, pp. 1364-1367. See Table II. - N. J. A. Sloane, May 21 2014
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).
LINKS
Suresh Govindarajan, Table of n, a(n) for n = 0..6500 (first 401 terms from T. D. Noe)
G. Almkvist, Asymptotic formulas and generalized Dedekind sums, Exper. Math., 7 (No. 4, 1998), pp. 343-359.
G. E. Andrews and P. Paule, MacMahon's partition analysis XII: Plane Partitions, J. Lond. Math. Soc., 76 (2007), 647-666.
A. O. L. Atkin, P. Bratley, I. G. McDonald and J. K. S. McKay, Some computations for m-dimensional partitions, Proc. Camb. Phil. Soc., 63 (1967), 1097-1100.
A. O. L. Atkin, P. Bratley, I. G. McDonald and J. K. S. McKay, Some computations for m-dimensional partitions, Proc. Camb. Phil. Soc., 63 (1967), 1097-1100. [Annotated scanned copy]
Michael Beeler, R. William Gosper and Richard C. Schroeppel, HAKMEM, ITEM 18, Memo AIM-239, Artificial Intelligence Laboratory, Massachusetts Institute of Technology, Cambridge, Mass., 1972.
Edward A. Bender, Asymptotic methods in enumeration, SIAM Review 16 (1974), no. 4, p. 509.
E. A. Bender and D. E. Knuth, Enumeration of Plane Partitions, J. Combin. Theory A. 13, 40-54, 1972.
S. Benvenuti, B. Feng, A. Hanany and Y. H. He, Counting BPS operators in gauge theories: Quivers, syzygies and plethystics, arXiv:hep-th/0608050, p. 41-42.
D. M. Bressoud and J. Propp, How the alternating sign matrix conjecture was solved, Notices Amer. Math. Soc., 46 (No. 6, 1999), 637-646.
Shouvik Datta, M. R. Gaberdiel, W. Li, and C. Peng, Twisted sectors from plane partitions, arXiv preprint arXiv:1606.07070 [hep-th], 2016.
Wenjie Fang, Hsien-Kuei Hwang, and Mihyun Kang, Phase transitions from exp(n^(1/2)) to exp(n^(2/3)) in the asymptotics of banded plane partitions, arXiv:2004.08901 [math.CO], 2020.
Steven Finch, Integer Partitions, September 22, 2004. [Cached copy, with permission of the author]
P. Flajolet and R. Sedgewick, Analytic Combinatorics, 2009; see page 580.
Bernhard Heim, Markus Neuhauser and Robert Tröger, Inequalities for Plane Partitions, arXiv:2109.15145 [math.CO], 2021.
Vaclav Kotesovec, A method of finding the asymptotics of q-series based on the convolution of generating functions, arXiv:1509.08708 [math.CO], 2015-2016, p. 18.
D. E. Knuth, A Note on Solid Partitions, Math. Comp. 24, 955-961, 1970.
Oleg Lazarev, Matt Mizuhara and Ben Reid, Some Results in Partitions, Plane Partitions, and Multipartitions, 13 August 2010.
P. A. MacMahon, Combinatory analysis.
J. Mangual, McMahon's Formula via Free Fermions, arXiv preprint arXiv:1210.7109 [math.CO], 2012. - From N. J. A. Sloane, Jan 01 2013
Ville Mustonen and R. Rajesh, Numerical Estimation of the Asymptotic Behaviour of Solid Partitions ..., arXiv:cond-mat/0303607 [cond-mat.stat-mech], 2003.
L. Mutafchiev and E. Kamenov, On The Asymptotic Formula for the Number of Plane Partitions..., arXiv:math/0601253 [math.CO], 2006; C. R. Acad. Bulgare Sci. 59(2006), No. 4, 361-366.
Ken Ono, Sudhir Pujahari and Larry Rolen, Turán inequalities for the plane partition function, arXiv:2201.01352 [math.NT], 2022.
I. Pak, Partition bijections, a survey, Ramanujan J. 12 (2006) 5-75.
A. Rovenchak, Enumeration of plane partitions with a restricted number of parts, arXiv preprint arXiv:1401.4367 [math-ph], 2014.
Raphael Schumacher, The self-counting identity, Fib. Quart., 55 (No. 2 2017), 157-167.
N. J. A. Sloane, Transforms
J. Stienstra, Mahler measure, Eisenstein series and dimers, arXiv:math/0502197 [math.NT], 2005.
Balázs Szendrői, Non-commutative Donaldson-Thomas invariants and the conifold, Geometry & Topology 12.2 (2008): 1171-1202.
Eric Weisstein's World of Mathematics, Plane Partition
E. M. Wright, Rotatable partitions, J. London Math. Soc., 43 (1968), 501-505.
FORMULA
G.f.: Product_{k >= 1} 1/(1 - x^k)^k. - MacMahon, 1912.
Euler transform of sequence [1, 2, 3, ...].
a(n) ~ (c_2 / n^(25/36)) * exp( c_1 * n^(2/3) ), where c_1 = A249387 = 2.00945... and c_2 = A249386 = 0.23151... - Wright, 1931. Corrected Jun 01 2010 by Rod Canfield - see Mutafchiev and Kamenov. The exact value of c_2 is e^(2c)*2^(-11/36)*zeta(3)^(7/36)*(3*Pi)^(-1/2), where c = Integral_{y=0..inf} (y*log(y)/(e^(2*Pi*y)-1))dy = (1/2)*zeta'(-1).
The exact value of c_1 is 3*2^(-2/3)*Zeta(3)^(1/3) = 2.0094456608770137530649... - Vaclav Kotesovec, Sep 14 2014
a(n) = (1/n) * Sum_{k=1..n} a(n-k)*sigma_2(k), n > 0, a(0)=1, where sigma_2(n) = A001157(n) = sum of squares of divisors of n. - Vladeta Jovovic, Jan 20 2002
G.f.: exp(Sum_{n>0} sigma_2(n)*x^n/n). a(n) = Sum_{pi} Product_{i=1..n} binomial(k(i)+i-1, k(i)) where pi runs through all nonnegative solutions of k(1)+2*k(2)+..+n*k(n)=n. - Vladeta Jovovic, Jan 10 2003
From Vaclav Kotesovec, Nov 07 2016: (Start)
More precise asymptotics: a(n) ~ Zeta(3)^(7/36) * exp(3 * Zeta(3)^(1/3) * (n/2)^(2/3) + 1/12) / (A * sqrt(3*Pi) * 2^(11/36) * n^(25/36))
* (1 + c1/n^(2/3) + c2/n^(4/3) + c3/n^2), where
c1 = -0.23994424421250649114273759... = -277/(864*(2*Zeta(3))^(1/3)) - Zeta(3)^(2/3)/(1440*2^(1/3))
c2 = -0.02576771365117401620018082... = 353*Zeta(3)^(1/3)/(248832*2^(2/3)) - 17*Zeta(3)^(4/3)/(3225600*2^(2/3)) - 71575/(1492992*(2*Zeta(3))^(2/3))
c3 = -0.00533195302658826100834286... = -629557/859963392 - 42944125/(7739670528*Zeta(3)) + 14977*Zeta(3)/1114767360 - 22567*Zeta(3)^2/250822656000
and A = A074962 is the Glaisher-Kinkelin constant.
(End)
EXAMPLE
A planar partition of 13:
4 3 1 1
2 1
1
a(5) = (1/5!)*(sigma_2(1)^5+10*sigma_2(2)*sigma_2(1)^3+20*sigma_2(3)*sigma_2(1)^2+ 15*sigma_2(1)*sigma_2(2)^2+30*sigma_2(4)*sigma_2(1)+20*sigma_2(2)*sigma_2(3)+24*sigma_2(5)) = 24. - Vladeta Jovovic, Jan 10 2003
From David Scambler and Joerg Arndt, May 01 2013: (Start)
There are a(4) = 13 partitions of 4 objects of 2 colors ('b' and 'w'), each part containing at least one black object:
1 black part:
[ bwww ]
2 black parts:
[ bbww ]
[ bww, b ]
[ bw, bw ]
3 black parts:
[ bbbw ]
[ bbw, b ]
[ bb, bw ]
(but not: [bw, bb ] )
[ bw, b, b ]
4 black parts:
[ bbbb ]
[ bbb, b ]
[ bb, bb ]
[ bb, b, b ]
[ b, b, b, b ]
(End)
From Geoffrey Critzer, Nov 29 2014: (Start)
The corresponding partitions of the integer 4 are:
4'''
4''
3'' + 1
2' + 2'
4'
3' + 1
2 + 2'
2' + 1 + 1
4
3 + 1
2 + 2
2 + 1 + 1
1 + 1 + 1 + 1.
(End)
From Gus Wiseman, Sep 25 2018: (Start)
Non-isomorphic representatives of the a(4) = 13 chains of multisets whose dual is also a chain of multisets:
{{1,1,1,1}}
{{1,1,2,2}}
{{1,2,2,2}}
{{1,2,3,3}}
{{1,2,3,4}}
{{1},{1,1,1}}
{{2},{1,2,2}}
{{3},{1,2,3}}
{{1,1},{1,1}}
{{1,2},{1,2}}
{{1},{1},{1,1}}
{{2},{2},{1,2}}
{{1},{1},{1},{1}}
(End)
G.f. = 1 + x + 3*x^2 + 6*x^3 + 13*x^4 + 24*x^5 + 48*x^6 + 86*x^7 + 160*x^8 + ...
MAPLE
series(mul((1-x^k)^(-k), k=1..64), x, 63);
# second Maple program:
a:= proc(n) option remember; `if`(n=0, 1, add(
a(n-j)*numtheory[sigma][2](j), j=1..n)/n)
end:
seq(a(n), n=0..50); # Alois P. Heinz, Aug 17 2015
MATHEMATICA
CoefficientList[Series[Product[(1 - x^k)^-k, {k, 64}], {x, 0, 64}], x]
Zeta[3]^(7/36)/2^(11/36)/Sqrt[3 Pi]/Glaisher E^(3 Zeta[3]^(1/3) (n/2)^(2/3) + 1/12)/n^(25/36) (* asymptotic formula after Wright; Vaclav Kotesovec, Jun 23 2014 *)
a[0] = 1; a[n_] := a[n] = Sum[a[n - j] DivisorSigma[2, j], {j, n}]/n; Table[a[n], {n, 0, 50}] (* Jean-François Alcover, Sep 21 2015, after Alois P. Heinz *)
CoefficientList[Series[Exp[Sum[DivisorSigma[2, n] x^n/n, {n, 50}]], {x, 0, 50}], x] (* Eric W. Weisstein, Feb 01 2018 *)
PROG
(PARI) {a(n) = if( n<0, 0, polcoeff( exp( sum( k=1, n, x^k / (1 - x^k)^2 / k, x * O(x^n))), n))}; /* Michael Somos, Jan 29 2005 */
(PARI) {a(n) = if( n<0, 0, polcoeff( prod( k=1, n, (1 - x^k + x * O(x^n))^-k), n))}; /* Michael Somos, Jan 29 2005 */
(PARI) my(N=66, x='x+O('x^N)); Vec( prod(n=1, N, (1-x^n)^-n) ) \\ Joerg Arndt, Mar 25 2014
(PARI) A000219(n)=#PlanePartitions(n) \\ See A091298 for PlanePartitions(). For illustrative use: much slower than the above. - M. F. Hasler, Sep 24 2018
(Python)
from sympy import cacheit
from sympy.ntheory import divisor_sigma
@cacheit
def A000219(n):
if n <= 1:
return 1
return sum(A000219(n - k) * divisor_sigma(k, 2) for k in range(1, n + 1)) // n
print([A000219(n) for n in range(20)])
# R. J. Mathar, Oct 18 2009
(Julia)
using Nemo, Memoize
@memoize function a(n)
if n == 0 return 1 end
s = sum(a(n - j) * divisor_sigma(j, 2) for j in 1:n)
return div(s, n)
end
[a(n) for n in 0:20] # Peter Luschny, May 03 2020
(SageMath) # uses[EulerTransform from A166861]
b = EulerTransform(lambda n: n)
print([b(n) for n in range(37)]) # Peter Luschny, Nov 11 2020
CROSSREFS
Differences: A191659, A191660, A191661.
Row sums of A089353 and A091438 and A091298.
Column k=1 of A144048. - Alois P. Heinz, Nov 02 2012
Sequences "number of r-line partitions": A000041 (r=1), A000990 (r=2), A000991 (r=3), A002799 (r=4), A001452 (r=5), A225196 (r=6), A225197 (r=7), A225198 (r=8), A225199 (r=9).
KEYWORD
nonn,nice,easy,core
EXTENSIONS
Corrected by N. J. A. Sloane, Jul 29 2006
Minor edits by Vaclav Kotesovec, Oct 27 2014
STATUS
approved
Expansion of 1 / Product_{k>=1} (1-x^k)^(k+1).
(Formerly M1601)
+10
37
1, 2, 6, 14, 33, 70, 149, 298, 591, 1132, 2139, 3948, 7199, 12894, 22836, 39894, 68982, 117948, 199852, 335426, 558429, 922112, 1511610, 2460208, 3977963, 6390942, 10206862, 16207444, 25596941, 40214896, 62868772, 97814358
OFFSET
0,2
COMMENTS
Also, a(n) = number of partitions of the integer n where there are k+1 different kinds of part k for k = 1, 2, 3, ....
Also, a(n) = number of partitions of n objects of 2 colors. These are set partitions, the n objects are not labeled but colored, using two colors. For each subset of size k there are k+1 different possibilities, i=0..k white and k-i black objects.
Also, a(n) = number of simple unlabeled graphs with n nodes of 2 colors whose components are complete graphs. - Geoffrey Critzer, Sep 27 2012
REFERENCES
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see Exercise 7.99, p. 484 and pp. 548-549.
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..10000 (first 1001 terms from T. D. Noe)
P. J. Cameron, Some sequences of integers, Discrete Math., 75 (1989), 89-102; also in "Graph Theory and Combinatorics 1988", ed. B. Bollobas, Annals of Discrete Math., 43 (1989), 89-102.
Carlos A. A. Florentino, Plethystic exponential calculus and permutation polynomials, arXiv:2105.13049 [math.CO], 2021. Mentions this sequence.
P. A. MacMahon, Memoir on symmetric functions of the roots of systems of equations, Phil. Trans. Royal Soc. London, 181 (1890), 481-536; Coll. Papers II, 32-87.
N. J. A. Sloane, Transforms
R. P. Stanley, Theory and Applications of Plane Partitions: Part 2, Studies in Appl. Math., 1 (1971), 259-279.
R. P. Stanley, The conjugate trace and trace of a plane partition, J. Combin. Theory, vol. A14 53-65 1973, esp. p. 64.
FORMULA
EULER transform of b(n) = n+1.
a(n) ~ Zeta(3)^(13/36) * exp(1/12 - Pi^4/(432*Zeta(3)) + Pi^2 * n^(1/3) / (3*2^(4/3)*Zeta(3)^(1/3)) + 3*Zeta(3)^(1/3) * n^(2/3) / 2^(2/3)) / (A * 2^(23/36) * 3^(1/2) * Pi * n^(31/36)), where A = A074962 = 1.2824271291... is the Glaisher-Kinkelin constant and Zeta(3) = A002117 = 1.202056903... . - Vaclav Kotesovec, Mar 07 2015
a(n) = A089353(n+m, m), n >= 1, for each m >= n. a(0) =1. See the Stanley reference, Exercise 7.99. - Wolfdieter Lang, Mar 09 2015
G.f.: exp(Sum_{k>=1} (sigma_1(k) + sigma_2(k))*x^k/k). - Ilya Gutkovskiy, Aug 11 2018
EXAMPLE
We represent each summand, k, in a partition of n as k identical objects. Then we color each object. We have no regard for the order of the colored objects.
a(3) = 14 because we have: www; wwb; wbb; bbb; ww + w; ww + b; wb + w; wb + b; bb + w; bb + b; w + w + w; w + w + b; w + b + b; b + b + b, where the 2 colors are black b and white w. - Geoffrey Critzer, Sep 27 2012
a(3) = 14 because we have: 3; 3'; 3''; 3'''; 2 + 1; 2 + 1'; 2' + 1; 2' + 1'; 2'' + 1; 2'' + 1'; 1 + 1 + 1; 1 + 1 + 1'; 1 + 1' + 1'; 1' + 1' + 1', where a part k of different sorts is given as k, k', k'', etc. - Joerg Arndt, Mar 09 2015
From Alois P. Heinz, Mar 09 2015: (Start)
The a(4) = 33 = 5 + 9 + 6 + 8 + 5 partitions of 4 objects of 2 colors are:
5 partitions for the integer partition of 4 = 1 + 1 + 1 + 1:
01: {{b}, {b}, {b}, {b}}
02: {{b}, {b}, {b}, {w}}
03: {{b}, {b}, {w}, {w}}
04: {{b}, {w}, {w}, {w}}
05: {{w}, {w}, {w}, {w}}
9 partitions for the integer partition of 4 = 1 + 1 + 2:
06: {{b}, {b}, {b,b}}
07: {{b}, {w}, {b,b}}
08: {{w}, {w}, {b,b}}
09: {{b}, {b}, {w,b}}
10: {{b}, {w}, {w,b}}
11: {{w}, {w}, {w,b}}
12: {{b}, {b}, {w,w}}
13: {{b}, {w}, {w,w}}
14: {{w}, {w}, {w,w}}
6 partitions for the integer partition of 4 = 2 + 2:
15: {{b,b}, {b,b}}
16: {{b,b}, {w,b}}
17: {{b,b}, {w,w}}
18: {{w,b}, {w,b}}
19: {{w,b}, {w,w}}
20: {{w,w}, {w,w}}
8 partitions for the integer partition of 4 = 1 + 3:
21: {{b}, {b,b,b}}
22: {{w}, {b,b,b}}
23: {{b}, {w,b,b}}
24: {{w}, {w,b,b}}
25: {{b}, {w,w,b}}
26: {{w}, {w,w,b}}
27: {{b}, {w,w,w}}
28: {{w}, {w,w,w}}
5 partitions for the integer partition of 4 = 4:
29: {{b,b,b,b}}
30: {{w,b,b,b}}
31: {{w,w,b,b}}
32: {{w,w,w,b}}
33: {{w,w,w,w}}
Some see number partitions, others see set partitions, ...
(End)
It is obvious from the example of Alois P. Heinz that a(n) enumerates multi-set partitions of a multi-set of n elements of two kinds. In the case that there is only one kind, this reduces to the usual case of numerical partitions. In the case that all the n elements are distinct, then this reduces to the case of set partitions. - Michael Somos, Mar 09 2015
There are a(3) = 14 plane partitions of 6 with trace 3; of 7 with trace 4; of 8 with trace 5; etc. See a formula above with the Stanley Exercise 7.99. - Wolfdieter Lang, Mar 09 2015
From Daniel Forgues, Mar 09 2015: (Start)
The a(3) = 14 = 4 + 6 + 4 partitions of 3 objects of 2 colors are:
4 partitions for the integer partition of 3 = 1 + 1 + 1:
01: {{b}, {b}, {b}}
02: {{b}, {b}, {w}}
03: {{b}, {w}, {w}}
04: {{w}, {w}, {w}}
6 partitions for the integer partition of 3 = 1 + 2:
05: {{b}, {b,b}}
06: {{w}, {b,b}}
07: {{b}, {w,b}}
08: {{w}, {w,b}}
09: {{b}, {w,w}}
10: {{w}, {w,w}}
4 partitions for the integer partition of 3 = 3:
11: {{b,b,b}}
12: {{w,b,b}}
13: {{w,w,b}}
14: {{w,w,w}}
The a(2) = 6 = 3 + 3 partitions of 2 objects of 2 colors are:
3 partitions for the integer partition of 2 = 1 + 1:
01: {{b}, {b}}
02: {{b}, {w}}
03: {{w}, {w}}
3 partitions for the integer partition of 2 = 2:
04: {{b,b}}
05: {{w,b}}
06: {{w,w}}
The a(1) = 2 partitions of 1 object of 2 colors are:
2 partitions for the integer partition of 1 = 1:
01: {{b}}
02: {{w}}
a(0) = 1: the empty partition, since empty sum is 0.
Triangle(sort of, since n_th row has p(n) = A000041 terms):
1: 2
2: 3, 3
3: 4, 6, 4
4: 5, 9, 6, 8, 5
5: 6, ?, ?, ?, ?, ?, 6
6: 7, ?, ?, ?, ?, ?, ?, ?, ?, ?, 7
Can we find a recurrence relation? (End)
MAPLE
mul( (1-x^i)^(-i-1), i=1..80); series(%, x, 80); seriestolist(%);
# second Maple program:
with(numtheory): etr:= proc(p) local b; b:=proc(n) option remember; local d, j; if n=0 then 1 else add(add(d*p(d), d=divisors(j)) *b(n-j), j=1..n)/n fi end end: a:=etr(n-> n+1): seq(a(n), n=0..40); # Alois P. Heinz, Sep 08 2008
MATHEMATICA
max = 31; f[x_] = Product[ 1/(1-x^k)^(k+1), {k, 1, max}]; CoefficientList[ Series[ f[x], {x, 0, max}], x] (* Jean-François Alcover, Nov 08 2011, after g.f. *)
etr[p_] := Module[{b}, b[n_] := b[n] = Module[{d, j}, If[n==0, 1, Sum[ Sum[ d*p[d], {d, Divisors[j]}]*b[n-j], {j, 1, n}]/n]]; b]; a = etr[#+1&]; Table[a[n], {n, 0, 40}] (* Jean-François Alcover, Nov 23 2015, after Alois P. Heinz *)
PROG
(PARI) a(n)=polcoeff(prod(i=1, n, (1-x^i+x*O(x^n))^-(i+1)), n)
CROSSREFS
Row sums of A054225.
Column k=2 of A075196.
KEYWORD
nonn,easy,nice
EXTENSIONS
Edited by Christian G. Bower, Sep 07 2002
New name from Joerg Arndt, Mar 09 2015
Restored 1995 name. - N. J. A. Sloane, Mar 09 2015
STATUS
approved
Triangle read by rows: T(n,k) is the number of multisets of k odd numbers with a cap of the total sum set to n.
+10
3
1, 1, 1, 2, 1, 1, 2, 3, 1, 1, 3, 4, 3, 1, 1, 3, 8, 5, 3, 1, 1, 4, 10, 10, 5, 3, 1, 1, 4, 16, 15, 11, 5, 3, 1, 1, 5, 20, 27, 17, 11, 5, 3, 1, 1, 5, 29, 38, 32, 18, 11, 5, 3, 1, 1, 6, 35, 60, 49, 34, 18, 11, 5, 3, 1, 1, 6, 47, 84, 83, 54, 35, 18, 11, 5, 3, 1, 1, 7, 56, 122, 123
OFFSET
1,4
COMMENTS
By considering the partitions of n into k parts we set a cap on the odd numbers of each part and count the multisets (ordered k-tuples) of odd numbers where each number is not larger than the cap of its part.
Multiset transformation of A110654 or A065033.
FORMULA
T(n,1) = A110654(n).
T(n,k) = Sum_{c_i*N_i=n,i=1..k} binomial(T(N_i,1)+c_i-1,c_i) for 1 < k <= n.
G.f.: Product_{j>=1} (1-y*x^j)^(-ceiling(j/2)). - Alois P. Heinz, Apr 13 2017
EXAMPLE
T(6,2) = 3+2+3 = 8 counts {1,1} {1,3}, and {3,3} from taking two odd numbers <= 3; it counts {1,1} and {1,3} from taking an odd number <= 2 and an odd number <= 4; and it counts {1,1}, {1,3} and {1,5} from taking an odd number <= 1 and an odd number <= 5.
T(6,3) = 1+2+2 = 5 counts {1,1,1} from taking three odd numbers <= 2; it counts {1,1,1} and {1,1,3} from taking an odd number <= 1 and an odd number <= 2 and an odd number <= 3; and it counts {1,1,1} and {1,1,3} from taking two odd numbers <= 1 and an odd number <= 4.
1
1 1
2 1 1
2 3 1 1
3 4 3 1 1
3 8 5 3 1 1
4 10 10 5 3 1 1
4 16 15 11 5 3 1 1
5 20 27 17 11 5 3 1 1
5 29 38 32 18 11 5 3 1 1
6 35 60 49 34 18 11 5 3 1 1
6 47 84 83 54 35 18 11 5 3 1 1
7 56 122 123 94 56 35 18 11 5 3 1 1
7 72 164 192 146 99 57 35 18 11 5 3 1 1
MAPLE
b:= proc(n, i, p) option remember; `if`(p>n, 0, `if`(n=0, 1,
`if`(min(i, p)<1, 0, add(b(n-i*j, i-1, p-j)*
binomial(ceil(i/2)+j-1, j), j=0..min(n/i, p)))))
end:
T:= (n, k)-> b(n$2, k):
seq(seq(T(n, k), k=1..n), n=1..16); # Alois P. Heinz, Apr 13 2017
MATHEMATICA
b[n_, i_, p_] := b[n, i, p] = If[p > n, 0, If[n == 0, 1, If[Min[i, p] < 1, 0, Sum[b[n - i*j, i - 1, p - j]*Binomial[Ceiling[i/2] + j - 1, j], {j, 0, Min[n/i, p]}]]]];
T[n_, k_] := b[n, n, k];
Table[T[n, k], {n, 1, 16}, {k, 1, n}] // Flatten (* Jean-François Alcover, May 19 2018, after Alois P. Heinz *)
CROSSREFS
Cf. A110654 (column 1), A003293 (row sums?), A089353 (equivalent Multiset transformation of A000027), A005232 (2nd column?), A097513 (3rd column?).
T(2n,n) gives A269628.
KEYWORD
nonn,tabl
AUTHOR
R. J. Mathar, Jul 27 2016
STATUS
approved
Sum of the traces of all plane partitions of n.
+10
3
1, 4, 10, 26, 56, 126, 252, 512, 980, 1866, 3427, 6258, 11121, 19618, 33975, 58328, 98732, 165804, 275246, 453544, 740338, 1200088, 1929897, 3083898, 4893775, 7720826, 12106814, 18883104, 29291740, 45215386, 69451631, 106197524, 161656759, 245050410, 369935066
OFFSET
1,2
COMMENTS
Convolution of A000203 and A000219. - Vaclav Kotesovec, Sep 25 2016
Convolution of A340793 and A091360. - Omar E. Pol, Feb 16 2021
REFERENCES
G. E. Andrews, The Theory of Partitions, Addison-Wesley, 1976, pp. 179-201.
LINKS
FORMULA
G.f.: g(x) = Sum_{j>=1} (j*x^j/(1-x^j))/Product_{k>=1} (1-x^k)^k.
a(n) = Sum(k*A089353(n,k), k>=1).
EXAMPLE
a(3) = 10 because the 6 (=A000219(3)) planar partitions of 3 are [3], [2,1], [2;1], [1,1,1], [1;1;1], [1,1;1] (; indicates a new row); the sum of their traces is 3+2+2+1+1+1 = 10.
MAPLE
g:= (sum(j*x^j/(1-x^j), j = 1..100))/(product((1-x^k)^k, k = 1..100)): gser := series(g, x = 0, 40): seq(coeff(gser, x, m), m = 1 .. 35);
# second Maple program:
b:= proc(n, i) option remember; `if`(n=0, [1, 0], `if`(i<1, 0, add((p
->p+[0, j*p[1]])(b(n-i*j, i-1))*binomial(i+j-1, j), j=0..n/i)))
end:
a:= n-> b(n$2)[2]:
seq(a(n), n=1..50); # Alois P. Heinz, Sep 24 2018
MATHEMATICA
nmax = 50; Rest[CoefficientList[Series[Sum[j*x^j/(1-x^j), {j, 1, nmax}]*Product[1/(1-x^k)^k, {k, 1, nmax}], {x, 0, nmax}], x]] (* Vaclav Kotesovec, Sep 25 2016 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Emeric Deutsch, Sep 24 2016
STATUS
approved
Number of planar partitions of n with trace 4.
+10
1
1, 2, 6, 14, 33, 64, 127, 228, 404, 672, 1100, 1724, 2661, 3974, 5849, 8402, 11911, 16556, 22751, 30772, 41198, 54436, 71283, 92316, 118609, 150950, 190753, 239090, 297783, 368236, 452782, 553240, 672532, 812980, 978211, 1171144, 1396235
OFFSET
4,2
COMMENTS
Also number of partitions of n objects of 2 colors into 4 parts, each part containing at least one black object.
REFERENCES
G. E. Andrews, The Theory of Partitions, Addison-Wesley, 1976 (Ch. XI, exercise 5 and Ch. XII, exercise 5).
FORMULA
G.f.: q^4*(q^12+q^10+2*q^9+4*q^8+2*q^7+4*q^6+2*q^5+4*q^4+2*q^3+q^2+1) / ((-1+q^4)^2*(-1+q^3)^2*(-1+q^2)^2*(-1+q)^2).
CROSSREFS
Column 4 of A089353. Cf. A000219, A005380, A005993 (trace 2), A050531 (trace 3).
KEYWORD
easy,nonn
AUTHOR
EXTENSIONS
Edited and extended by Christian G. Bower, Jan 08 2004
STATUS
approved

Search completed in 0.019 seconds