|
%I #12 Jun 05 2021 08:43:07
%S 0,0,0,0,1,2,3,7,10,16,23,35,47,70,93,126,169,228,294,391,501,648,827,
%T 1057,1329,1683,2105,2631,3266,4056,4992,6156,7538,9221,11234,13664,
%U 16549,20033,24152,29077,34904,41844,50012,59710,71100,84541,100318,118869
%N Number of partitions of n such that (maximal multiplicity of parts) > (multiplicity of the maximal part).
%H Alois P. Heinz, <a href="/A240302/b240302.txt">Table of n, a(n) for n = 0..1000</a>
%F a(n) + A171979(n) = A000041(n) for n >= 1.
%e a(7) counts these 7 partitions: 511, 4111, 322, 3211, 31111, 22111, 211111.
%p b:= proc(n, i, k) option remember; `if`(n=0, `if`(k=0, 1, 0),
%p `if`(i<1, 0, b(n, i-1, k) +add(b(n-i*j, i-1, `if`(k=-1, j,
%p `if`(k=0, 0, `if`(j>k, 0, k)))), j=1..n/i)))
%p end:
%p a:= n-> b(n$2, -1):
%p seq(a(n), n=0..70); # _Alois P. Heinz_, Apr 12 2014
%t z = 60; f[n_] := f[n] = IntegerPartitions[n]; m[p_] := Max[Map[Length, Split[p]]] (* maximal multiplicity *)
%t Table[Count[f[n], p_ /; m[p] == Count[p, Max[p]]], {n, 0, z}] (* A171979 *)
%t Table[Count[f[n], p_ /; m[p] > Count[p, Max[p]]], {n, 0, z}] (* A240302 *)
%t (* Second program: *)
%t b[n_, i_, k_] := b[n, i, k] = If[n == 0, If[k == 0, 1, 0],
%t If[i < 1, 0, b[n, i - 1, k] + Sum[b[n - i*j, i - 1, If[k == -1, j,
%t If[k == 0, 0, If[j > k, 0, k]]]], {j, 1, n/i}]]];
%t a[n_] := b[n, n, -1];
%t a /@ Range[0, 70] (* _Jean-François Alcover_, Jun 05 2021, after _Alois P. Heinz_ *)
%Y Cf. A240221, A000041.
%K nonn,easy
%O 0,6
%A _Clark Kimberling_, Apr 04 2014
|
