|
%I #21 May 25 2022 08:54:40
%S 0,2,4,5,7,7,8,8,13,11,10,10,11,12,25,17,31,12,16,14,12,18,43,17,13,
%T 14,22,55,13,15,15,20,23,28,41,15,20,14,85,47,91,15,26,15,18,19,17,42,
%U 115,26,35,26,16,17,16,133,71,16,30,18,52,77,25,38,22,83
%N Sum of prime factors of prime(n)-1 (counted with multiplicity).
%C From an idea of _Michael B. Porter_.
%H Paolo P. Lava, <a href="/A210936/b210936.txt">Table of n, a(n) for n = 1..10000</a>
%F a(n) = A001414(A006093(n)). - _Michel Marcus_, Oct 05 2013
%e prime(10) = 29, and 29-1 = 28 = 2*2*7, so a(10) = 2+2+7 = 11.
%p with(numtheory);
%p P:=proc(i)
%p local a,k,n;
%p for n from 1 to i do
%p a:=ifactors(ithprime(n)-1)[2]; print(add(a[k][1]*a[k][2],k=1..nops(a)));
%p od; end:
%p # alternative
%p A210936 := proc(n)
%p local p,pplus,f ;
%p p := ithprime(n) ;
%p pplus := ifactors(p-1)[2] ;
%p add(op(1,f)*op(2,f),f=pplus) ;
%p end proc:
%p seq(A210936(n),n=1..300) ; # _R. J. Mathar_, May 25 2022
%Y Cf. A023508, A023514, A210934
%K nonn
%O 1,2
%A _Paolo P. Lava_, Mar 30 2012
|
