OFFSET
1,1
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..8192
PROG
(Python)
from sympy.ntheory.factor_ import digits
from sympy import factorint, nextprime
from operator import mul
from functools import reduce
def a006047(n):
d=digits(n, 3)
return n + 1 if n<3 else reduce(mul, [1 + d[i] for i in range(1, len(d))])
def a048673(n):
f = factorint(n)
return 1 if n==1 else (1 + reduce(mul, [nextprime(i)**f[i] for i in f]))//2
def a(n): return a006047(a048673(n))
print([a(n) for n in range(1, 101)]) # Indranil Ghosh, Jun 12 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, May 31 2017
STATUS
approved