A359431 - OEIS

A359431

a(n) = A325973(n) - A326043(n).

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

(list; graph; refs; listen; history; text; internal format)

OFFSET

1,16

LINKS

Antti Karttunen, Table of n, a(n) for n = 1..65537

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.