OFFSET
1,2
FORMULA
a(n) = 3*(n-1) + 1 - A073089(n).
PROG
(Python)
n, f, i, p, q, base = 1, 1, 0, 0, 1, 2
while i < 10000000:
i, p, q = i+1, p*base, q*base
if i == f:
p, n = p+1, n+1
f = f*n
n, a, j = 0, 0, 0
while p%q > 0:
a, f, p, q = a+1, p//q, q, p%q
if f == 1:
n = n+1
print(n, a-1)
CROSSREFS
KEYWORD
nonn
AUTHOR
A.H.M. Smeets, Jul 30 2018
STATUS
approved