OFFSET
1,1
COMMENTS
For n>1, a(n) = n iff mu(n) = 0; 0 if n is prime; and (n+1) iff mu(n) = 1.
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..20000
EXAMPLE
a(4) = 4 = 4 + A140664(4) = 4 + 0.
MATHEMATICA
Table[Exp[MangoldtLambda[n]]*MoebiusMu[n] + n, {n, 100}] (* G. C. Greubel, Feb 13 2019 *)
PROG
(PARI)
(Sage)
def A140668(n): return simplify(exp(add(moebius(d)*log(n/d) for d in divisors(n))))*moebius(n) + n
[A140668(n) for n in (1..100)] # G. C. Greubel, Feb 13 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Gary W. Adamson and Mats Granvik, May 20 2008
EXTENSIONS
More terms from Antti Karttunen, Feb 13 2019
STATUS
approved