OFFSET
1,2
PROG
(PARI) a(n) = my(p=prod(j=1, n, x+j/6), 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 Poly, rf
from sympy.abc import x
def A076642(n):
y = Poly(rf(6*x+1, n)).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