OFFSET
0,4
COMMENTS
The number of submultisets of a partition is the product of its multiplicities, each plus one.
The Heinz numbers of these partitions are given by A325795.
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..500
FORMULA
EXAMPLE
The a(1) = 1 through a(8) = 10 partitions:
(1) (11) (21) (211) (221) (321) (421) (3221)
(111) (1111) (311) (2211) (2221) (3311)
(2111) (3111) (3211) (4211)
(11111) (21111) (4111) (22211)
(111111) (22111) (32111)
(31111) (41111)
(211111) (221111)
(1111111) (311111)
(2111111)
(11111111)
MAPLE
b:= proc(n, i, p) option remember; `if`(n=0 or i=1,
`if`(n=p-1, 1, 0), add(`if`(irem(p, j+1, 'r')=0,
(w-> b(w, min(w, i-1), r))(n-i*j), 0), j=0..n/i))
end:
a:= n-> combinat[numbpart](n)-add(b(n$2, k), k=0..n):
seq(a(n), n=0..55); # Alois P. Heinz, Aug 17 2019
MATHEMATICA
Table[Length[Select[IntegerPartitions[n], Times@@(1+Length/@Split[#])>n&]], {n, 0, 30}]
(* Second program: *)
b[n_, i_, p_] := b[n, i, p] = If[n == 0 || i == 1,
If[n == p - 1, 1, 0], Sum[If[Mod[p, j + 1] == 0,
Function[w, b[w, Min[w, i-1], p/(j+1)]][n-i*j], 0], {j, 0, n/i}]];
a[n_] := PartitionsP[n] - Sum[b[n, n, k], {k, 0, n}];
a /@ Range[0, 55] (* Jean-François Alcover, May 13 2021, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, May 25 2019
STATUS
approved