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