login
A333178
a(n) = Sum_{d|n, gcd(d, n/d) = 1} prime(d).
1
2, 5, 7, 9, 13, 23, 19, 21, 25, 45, 33, 51, 43, 65, 65, 55, 61, 89, 69, 91, 97, 115, 85, 115, 99, 147, 105, 133, 111, 223, 129, 133, 175, 203, 179, 183, 159, 235, 215, 205, 181, 337, 193, 233, 233, 287, 213, 283, 229, 331, 299, 289, 243, 359, 301, 301, 343, 385, 279, 461
OFFSET
1,1
FORMULA
Sum_{d|n, gcd(d, n/d) = 1} (-1)^omega(n/d) * a(d) = prime(n).
MATHEMATICA
a[n_] := Sum[If[GCD[n/d, d] == 1, Prime[d], 0], {d, Divisors[n]}]; Table[a[n], {n, 1, 60}]
PROG
(PARI) a(n) = sumdiv(n, d, if (gcd(d, n/d) ==1, prime(d))); \\ Michel Marcus, Mar 10 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Mar 10 2020
STATUS
approved