A240302 - OEIS

%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