OFFSET
3,3
COMMENTS
See A332916 for more information.
PROG
(PARI) f(n)=(40+6*n-4*n^2)/2^n-83/2+331/12*n-6*n^2+2/3*n^3;
for(k=3, 70, print1(round(log(denominator(f(k)))/log(2)), ", "))
(Python)
from fractions import Fraction
def f(n): return Fraction(40+6*n-4*n**2, 2**n) - Fraction(83, 2) + Fraction(331*n, 12) - 6*n**2 + Fraction(2*n**3, 3)
def a(n): return (f(n).denominator).bit_length() - 1
print([a(n) for n in range(3, 71)]) # Michael S. Branicky, Aug 31 2021
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Hugo Pfoertner, Mar 05 2020
STATUS
approved