OFFSET
1,1
COMMENTS
Among the initial 15275 terms, in 110 cases A276085(k) is a multiple of 5^5, and for no cases a multiple of 7^7.
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..15275
EXAMPLE
A276085(55) = 216 = 2^3 * 3^3, which although it has a divisor of the form p^p, with p an odd prime, it is also a multiple of 4, and therefore 55 is NOT included in this sequence.
A276085(174) = 223092873 = 3^3 * 3 * 1063 * 2591, which has a divisor of the form p^p, with p an odd prime, thus 174 is included in this sequence.
A276085(4823) = 614889782588493750 = 2 * 3 * 5^5 * 13 * 2522624749081, thus 4823 is included.
A276085(1104299) = 11231250 = 2 * 3 * 5^5 * 599, thus 1104299 is included.
PROG
(PARI)
A083345(n) = { my(f=factor(n)); numerator(vecsum(vector(#f~, i, f[i, 2]/f[i, 1]))); };
A276085(n) = { my(f = factor(n), pr=1, i=1, s=0); for(k=1, #f~, while(i <= primepi(f[k, 1])-1, pr *= prime(i); i++); s += f[k, 2]*pr); (s); };
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Nov 11 2024
STATUS
approved