login
Number of partitions of n with zero crank.
29

%I #36 May 10 2023 06:18:16

%S 0,0,1,1,1,1,1,2,2,4,4,7,7,11,12,17,19,27,30,41,48,62,73,95,110,140,

%T 166,206,243,302,354,435,513,622,733,887,1039,1249,1467,1750,2049,

%U 2438,2847,3371,3934,4634,5398,6343,7367,8626,10009,11677,13521,15737,18184

%N Number of partitions of n with zero crank.

%C For a partition p, let l(p) = largest part of p, w(p) = number of 1's in p, m(p) = number of parts of p larger than w(p). The crank of p is given by l(p) if w(p) = 0, otherwise m(p)-w(p).

%H Alois P. Heinz and Vaclav Kotesovec, <a href="/A064410/b064410.txt">Table of n, a(n) for n = 1..10000</a> (terms 1..1000 from Alois P. Heinz)

%H Brian Hopkins and James A. Sellers, <a href="https://arxiv.org/abs/2305.05096">On Blecher and Knopfmacher's Fixed Points for Integer Partitions</a>, arXiv:2305.05096 [math.CO], 2023. Mentions this sequence.

%H Brian Hopkins, James A. Sellers, and Dennis Stanton, <a href="https://arxiv.org/abs/2009.10873">Dyson's Crank and the Mex of Integer Partitions</a>, arXiv:2009.10873 [math.CO], 2020. Mentions this sequence.

%F a(n) = A000041(n)-2*A001522(n). a(n) = A064391(n, 0).

%F a(n) ~ exp(Pi*sqrt(2*n/3)) * Pi / (3 * 2^(9/2) * n^(3/2)). - _Vaclav Kotesovec_, May 06 2018

%F a(n > 1) = A064428(n) - A001522(n), where A001522/A064428 count odd/even-length compositions with alternating parts strictly decreasing. - _Gus Wiseman_, Apr 02 2021

%e a(10)=4 because there are 4 partitions of 10 with zero crank: 1+1+2+3+3, 1+1+4+4, 1+1+3+5 and 1+9.

%e From _Gus Wiseman_, Apr 02 2021: (Start)

%e The a(3) = 1 through a(14) = 11 partitions (A..D = 10..13):

%e 21 31 41 51 61 71 81 91 A1 B1 C1 D1

%e 3311 4311 4411 5411 5511 6511 6611

%e 5311 6311 6411 7411 7511

%e 33211 43211 7311 8311 8411

%e 44211 54211 9311

%e 53211 63211 55211

%e 332211 432211 64211

%e 73211

%e 442211

%e 532211

%e 3322211

%e (End)

%t nmax = 60; Rest[CoefficientList[Series[x - 1 + Sum[(-1)^k*(x^(k*(k + 1)/2) - x^(k*(k - 1)/2)), {k, 1, nmax}] / Product[1 - x^k, {k, 1, nmax}], {x, 0, nmax}], x]] (* _Vaclav Kotesovec_, Sep 26 2016 *)

%t Flatten[{0, Table[PartitionsP[n] - 2*Sum[(-1)^(j+1)*PartitionsP[n - j*((j+1)/2)], {j, 1, Floor[(Sqrt[8*n + 1] - 1)/2]}], {n, 2, 60}]}] (* _Vaclav Kotesovec_, Sep 26 2016 *)

%t ck[y_]:=With[{w=Count[y,1]},If[w==0,Max@@y,Count[y,_?(#>w&)]-w]];

%t Table[Length[Select[IntegerPartitions[n],ck[#]==0&]],{n,0,30}] (* _Gus Wiseman_, Apr 02 2021 *)

%o (Sage)

%o [[p.crank() for p in Partitions(n)].count(0) for n in (1..20)] # _Peter Luschny_, Sep 15 2014

%Y The version for positive crank is A001522.

%Y Central column of A064391.

%Y The version for nonnegative crank is A064428.

%Y The Heinz numbers of these partitions are A342192.

%Y A003242 counts anti-run compositions.

%Y A224958 counts compositions with alternating parts unequal.

%Y A257989 gives the crank of the partition with Heinz number n.

%Y Cf. A000009, A000041, A000726, A027193, A034008, A062968, A114921, A342194, A342343, A342528.

%K nonn

%O 1,8

%A _Vladeta Jovovic_, Sep 29 2001

%E More terms from _Reiner Martin_, Dec 26 2001