%I #16 Dec 09 2023 09:00:16
%S 0,0,2,4,6,8,12,14,17,20,23,26,30,33,38,41,44,47,52,55,60,63,68,71,76,
%T 78,84,89,93,97,103,106,111,115,121,124,128,131,138,142,146,151,159,
%U 162,168,171,176,181,187,190,196,201,206,210,218,221,227,232,238
%N Number of multisets of size n with elements from [n] whose element sum is larger than the product of all elements.
%H Alois P. Heinz, <a href="/A360971/b360971.txt">Table of n, a(n) for n = 0..400</a>
%e a(2) = 2: [1,1], [1,2].
%e a(3) = 4: [1,1,1], [1,1,2], [1,2,2], [1,1,3].
%e a(4) = 6: [1,1,1,1], [1,1,1,2], [1,1,2,2], [1,1,1,3], [1,1,2,3], [1,1,1,4].
%e a(8) = 17: [1,1,1,1,1,1,1,1], [1,1,1,1,1,1,1,2], [1,1,1,1,1,1,2,2], [1,1,1,1,1,2,2,2], [1,1,1,1,1,1,1,3], [1,1,1,1,1,1,2,3], [1,1,1,1,1,1,3,3], [1,1,1,1,1,1,1,4], [1,1,1,1,1,1,2,4], [1,1,1,1,1,1,3,4], [1,1,1,1,1,1,1,5], [1,1,1,1,1,1,2,5], [1,1,1,1,1,1,1,6], [1,1,1,1,1,1,2,6], [1,1,1,1,1,1,1,7], [1,1,1,1,1,1,2,7], [1,1,1,1,1,1,1,8].
%p b:= proc(n, i, s, p) `if`(s+n*i<=p, 0, `if`(n=0 or i=1, 1, g(n, i, s, p))) end:
%p g:= proc(n, i, s, p) option remember; add(b(n-1, j, s+j, p*j), j=1..i) end:
%p a:= n-> b(n$2, 0, 1):
%p seq(a(n), n=0..60);
%t b[n_, i_, s_, p_] := If[s + n*i <= p, 0, If[n == 0 || i == 1, 1, g[n, i, s, p]]];
%t g[n_, i_, s_, p_] := g[n, i, s, p] = Sum[b[n-1, j, s+j, p*j], {j, 1, i}];
%t a[n_] := b[n, n, 0, 1];
%t Table[a[n], {n, 0, 60}] (* _Jean-François Alcover_, Dec 09 2023, after _Alois P. Heinz_ *)
%Y Cf. A103586, A360538, A360563.
%K nonn
%O 0,3
%A _Alois P. Heinz_, Feb 27 2023