OFFSET
1,1
COMMENTS
See A371639 for the connection with Voronoi's congruence.
FORMULA
a(n) = valuation(denominator(Voronoi(3, n))) where Voronoi(c, n) = ((c^n - 1) * bernoulli(n)) / (n * c^(n - 1)).
MATHEMATICA
Array[2 # + IntegerExponent[#, 3] &, 64] (* Michael De Vlieger, Mar 31 2024 *)
PROG
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Peter Luschny, Mar 30 2024
STATUS
approved