OFFSET
1,36
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..65537
MATHEMATICA
a[n_]:=a[n]=If[n==1, 1, -Sum[(-1)^PrimeOmega[d]*a[n/d], {d, Select[Rest[Divisors[n]], GCD@@FactorInteger[#][[All, 2]]==1&]}]];
Table[Sum[a[d]*MoebiusMu[n/d], {d, Divisors[n]}], {n, 100}]
PROG
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, May 19 2018
EXTENSIONS
More terms from Antti Karttunen, Jul 29 2018
STATUS
approved