OFFSET
1,1
FORMULA
EXAMPLE
a(12) = 4 = 1+3 because (1 -1/2)(1 -1/3) = 1/3.
MAPLE
a:=proc(n) local b, ct, f; with(numtheory): b:=convert(factorset(n), list): ct:=nops(b): f:=simplify(product(1-1/b[j], j=1..ct)):numer(f)+denom(f) end: seq(a(n), n=1..100); # Emeric Deutsch, May 28 2005
CROSSREFS
KEYWORD
nonn
AUTHOR
Leroy Quet, May 27 2005
EXTENSIONS
More terms from Emeric Deutsch, May 28 2005
STATUS
approved