OFFSET
0,7
LINKS
Rémy Sigrist, Table of n, a(n) for n = 0..10000
FORMULA
a(n) = Sum_{k = 1..n} A360710(k).
PROG
(PARI) { my (s=0, f); for (n=1, #a=vector(78), print1 (s", "); f=factor(n); s+=a[n]=if (#f~==1, if (s, -sign(s), a[n-1]), prod(k=1, #f~, a[f[k, 1]^f[k, 2]]))) }
CROSSREFS
KEYWORD
sign,look
AUTHOR
Rémy Sigrist, Feb 17 2023
STATUS
approved