OFFSET
1,2
EXAMPLE
MATHEMATICA
b[1] = 0;
b[n_] := b[n] = Product[If[n > d^2, 1, 1 - b[d]], {d, Select[Range[n-1], Mod[n, #] == 0&]}];
a[n_] := Select[Divisors[n], b[#] == 1&] // Total;
Array[a, 100] (* Jean-François Alcover, Sep 18 2020 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Sep 23 2004
STATUS
approved