login
A388978
a(n) = A326128(n) * A389077(n).
4
0, 0, 0, 9, 0, 0, 0, 60, 40, 0, 0, 252, 0, 0, 0, 585, 0, 736, 0, 510, 0, 0, 0, 1620, 168, 0, 624, 966, 0, 0, 0, 3540, 0, 0, 0, 6195, 0, 0, 0, 3300, 0, 0, 0, 1716, 2640, 0, 0, 14220, 528, 3744, 0, 2496, 0, 11232, 0, 6300, 0, 0, 0, 12330, 0, 0, 4816, 22869, 0, 0, 0, 3570, 0, 0, 0, 39340, 0, 0, 6624, 4674, 0, 0, 0
OFFSET
1,4
FORMULA
a(n) = (n-A007913(n)) * (A000203(A003961(n))-A007913(A003961(n))-A003961(n)).
a(n) = A341528(n) + A191002(A007913(n)) - A191002(n) - n*A007913(A003961(n)) - A003973(n)*A007913(n) + A007913(n)*A003961(n).
a(n) <= A388979(n).
PROG
(PARI)
A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
A388978(n) = { my(u=A003961(n)); ((n-core(n)) * (sigma(u)-core(u)-u)); };
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Sep 23 2025
STATUS
approved