OFFSET
1,2
COMMENTS
Dirichlet inverse of sequence b(n) = 1+A095112(n).
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..16384
FORMULA
a(1) = 1, and for n > 1, a(n) = -Sum_{d|n, d<n} (1+A095112(n/d)) * a(d).
MATHEMATICA
A095112[n_] := n/Flatten[#[[1]]^Range[#[[2]]]& /@ FactorInteger[n]] // Total;
a[n_] := a[n] = If[n == 1, 1, -Sum[(1 + A095112[n/d]) a[d], {d, Most@ Divisors[n]}]];
Table[a[n], {n, 1, 100}] (* Jean-François Alcover, Nov 26 2023 *)
PROG
CROSSREFS
KEYWORD
sign
AUTHOR
Antti Karttunen, Nov 22 2023
STATUS
approved