A078319 - OEIS

%I #13 Sep 26 2015 17:25:26

%S 2,5,7,3,15,37,7,17,9,101,63,73,24,197,115,14,36,109,183,70,32,102,60,

%T 34,12,677,43,134,423,70,52,18,116,102,615,38,144,22,763,401,31,358,

%U 44,39,15,102,37,17,45,170,1303,55,288,18,108,162,20,678,1743,1801,1863

%N Sum of distinct prime factors of n*rad(n)+1, where rad=A007947 (squarefree kernel).

%C a(n)=A008472(A078310(n)).

%H Reinhard Zumkeller, <a href="/A078319/b078319.txt">Table of n, a(n) for n = 1..10000</a>

%o (PARI) vecprod(v)=prod(i=1,#v,v[i])

%o rad(n)=vecprod(factor(n)[,1])

%o a(n)=vecsum(factor(n*rad(n)+1)[,1]) \\ _Charles R Greathouse IV_, Jul 09 2013

%o (Haskell)

%o a078319 = a008472 . a078310  -- _Reinhard Zumkeller_, Jul 23 2013

%Y Cf. A078320, A078318, A078313.

%K nonn

%O 1,1

%A _Reinhard Zumkeller_, Nov 23 2002