login
A359431
a(n) = A325973(n) - A326043(n).
2
0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 4, 0, 0, 2, 0, 0, 0, 0, 5, 0, 0, 0, 3, 0, 0, 0, 6, 0, 0, 0, 0, 0, 0, 0, 8, 0, 0, 0, 0, 0, 6, 0, 8, 0, 0, 0, 0, 0, 0, 0, 12, 0, 0, 0, 0, 0, 0, 0, 16, 0, 0, 0, 0, 0, 0, 0, 12, 12, 0, 0, 0, 0, 0, 0, 12, 0, 0, 0, 0, 0, 0, 0, 20, 0, 0, 0, 10, 0, 0, 0, 14, 0
OFFSET
1,16
LINKS
PROG
(PARI)
A325973(n) = (1/2)*sumdiv(n, d, d*(issquarefree(d) + (1==gcd(d, n/d))));
A326043(n) = if(1==n, n, my(f = factor(n)); prod(i=1, #f~, floor((1/f[i, 2]) * ((f[i, 2]-1) + (((f[i, 1]^(1+f[i, 2])) - 1)/(f[i, 1]-1))))));
A359431(n) = (A325973(n)-A326043(n));
CROSSREFS
Cf. A048107 (positions of 0's), A325973, A326043, A359471.
Cf. also comments in A325981.
Sequence in context: A083929 A182660 A286100 * A338210 A122698 A002483
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jan 04 2023
STATUS
approved