A328638 a(n) is the sum of m such that 2 <= m <= n and gpf(m) = gpf(n), where gpf(i) is the greatest prime factor of i (A006530), with a(1) = 1.
1, 2, 3, 6, 5, 9, 7, 14, 18, 15, 11, 30, 13, 21, 30, 30, 17, 48, 19, 50, 42, 33, 23, 72, 75, 39, 99, 70, 29, 105, 31, 62, 66, 51, 105, 135, 37, 57, 78, 145, 41, 147, 43, 110, 190, 69, 47, 183, 196, 240, 102, 130, 53, 237, 165, 252, 114, 87, 59, 300
Offset: 1
Examples
5, 10, 15 and 20 have same gpf as 20, so a(20) = 5 + 10 + 15 + 20 = 50.
Programs
-
PARI
gpf(n)={my(f=factor(n)[,1]); f[#f]} a(n)={if(n<=1, n==1, my(t=gpf(n)); sum(i=2, n, if(gpf(i)==t, i, 0)))} \\ Andrew Howroyd, Oct 28 2019
Comments