OFFSET
1,3
PROG
(PARI) a(n) = my(p=prod(j=1, n, x+j/2), m=vecmax(Vec(p))); for (i=0, poldegree(p), if (polcoef(p, i)==m, return(i))); \\ Michel Marcus, Feb 19 2021
(Python)
from sympy import prod, Poly
from sympy.abc import x
def A076634(n):
y = Poly(prod(2*x+i for i in range(1, n+1))).all_coeffs()[::-1]
return y.index(max(y)) # Chai Wah Wu, Mar 07 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Benoit Cloitre, Nov 10 2002
STATUS
approved