OFFSET
0,3
MAPLE
b:= proc(n) b(n):= `if`(n=0, 1, expand(b(n-1)*(1+n*x)^n)) end:
a:= n-> max(coeffs(b(n))):
seq(a(n), n=0..11); # Alois P. Heinz, Jan 31 2024
PROG
(PARI) a(n)=vecmax(Vec(prod(k=1, n, (1+k*x)^k)));
vector(20, n, a(n-1))
(Python)
from collections import Counter
from math import comb
def A369770(n):
c = {0:1}
for k in range(1, n+1):
d = Counter(c)
for j in c:
a = c[j]
for i in range(1, k+1):
d[j+i] += comb(k, i)*k**i*a
c = d
return max(c.values()) # Chai Wah Wu, Jan 31 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Joerg Arndt, Jan 31 2024
STATUS
approved