Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
A373984
a(n) = A108951(n) - A373158(n), where A108951 is fully multiplicative and A373158 is fully additive with a(p) = p# for prime p, where x# is the primorial A034386(x).
3
1, 0, 0, 0, 0, 4, 0, 2, 24, 28, 0, 14, 0, 208, 144, 8, 0, 58, 0, 86, 1044, 2308, 0, 36, 840, 30028, 198, 626, 0, 322, 0, 22, 11544, 510508, 6060, 128, 0, 9699688, 150144, 204, 0, 2302, 0, 6926, 1038, 223092868, 0, 82, 43680, 1738, 2552544, 90086, 0, 412, 66960, 1464, 48498444, 6469693228, 0, 680, 0, 200560490128
OFFSET
1,6
PROG
(PARI) A373984(n) = { my(f=factor(n), m=1, s=0); for(i=1, #f~, my(x=prod(i=1, primepi(f[i, 1]), prime(i))); s += f[i, 2]*x; m *= x^f[i, 2]); (m-s); };
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jun 25 2024
STATUS
approved