|
| |
|
|
A064174
|
|
Number of partitions of n with nonnegative rank.
|
|
9
|
|
|
|
1, 1, 2, 3, 4, 6, 9, 12, 17, 23, 31, 42, 56, 73, 96, 125, 161, 207, 265, 336, 426, 536, 672, 840, 1046, 1296, 1603, 1975, 2425, 2970, 3628, 4417, 5367, 6503, 7861, 9482, 11412, 13702, 16423, 19642, 23447, 27938, 33231, 39453, 46767, 55342, 65386, 77135
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
|
OFFSET
|
1,3
|
|
|
COMMENTS
|
The rank of a partition is the largest summand minus the number of summands.
This sequence (up to proof) equals "partitions of 2n with even number of parts, ending in 1, with max descent of 1, where the number of odd parts in odd places equals the number of odd parts in even places. (See link and 2nd Mathematica line.) - Wouter Meeussen, Mar 29 2013
Number of partitions p of n such that max(max(p), number of parts of p) is a part of p. - Clark Kimberling, Feb 28 2014
|
|
|
LINKS
|
Alois P. Heinz, Table of n, a(n) for n = 1..1000
Rekha Biswal bijection between number of partitions of 2n satisfying certain conditions with number of partitions of n (Mathoverflow)
|
|
|
FORMULA
|
a(n) = (A000041(n) + A047993(n))/2.
a(n) = p(n-1) - p(n-5) + p(n-12) - ... -(-1)^k*p(n-(3*k^2-k)/2) + ..., where p() is A000041(). - Vladeta Jovovic, Aug 04 2004
G.f.: Sum_{n>=1} x^n * Product_{k=1..n} (1 - x^(n+k-1))/(1 - x^k). - Paul D. Hanna, Aug 03 2015
|
|
|
EXAMPLE
|
a(20) = p(19) - p(15) + p(8) = 490 - 176 + 22 = 336.
|
|
|
MAPLE
|
f:= n -> add((-1)^(k+1)*combinat:-numbpart(n-(3*k^2-k)/2), k=1..floor((1+sqrt(24*n+1))/6)):
map(f, [$1..100]); # Robert Israel, Aug 03 2015
|
|
|
MATHEMATICA
|
Table[Count[IntegerPartitions[n], q_ /; First[q] >= Length[q]], {n, 16}]
(* also *)
Table[Count[IntegerPartitions[2n], q_/; Last[q]===1 && Max[q-PadRight[Rest[q], Length[q]]]<=1 && Count[First/@Partition[q, 2], _?OddQ]==Count[Last/@Partition[q, 2], _?OddQ]], {n, 16}]
(* also *)
Table[Count[IntegerPartitions[n], p_ /; MemberQ[p, Max[Max[p], Length[p]]]], {n, 50}] (* Clark Kimberling, Feb 28 2014 *)
|
|
|
PROG
|
(PARI) {a(n) = my(A=1); A = sum(m=0, n, x^m*prod(k=1, m, (1-x^(m+k-1))/(1-x^k +x*O(x^n)))); polcoeff(A, n)}
for(n=1, 60, print1(a(n), ", ")) \\ Paul D. Hanna, Aug 03 2015
|
|
|
CROSSREFS
|
Cf. A063995, A064173.
Row sums of triangle A105806.
Sequence in context: A035992 A036003 A027338 * A062121 A094995 A018591
Adjacent sequences: A064171 A064172 A064173 * A064175 A064176 A064177
|
|
|
KEYWORD
|
nonn
|
|
|
AUTHOR
|
Vladeta Jovovic, Sep 20 2001
|
|
|
EXTENSIONS
|
Mathematica programs modified by Clark Kimberling, Feb 12 2014
|
|
|
STATUS
|
approved
|
| |
|
|
