|
%I #17 Sep 04 2017 23:28:25
%S 0,0,2,3,5,5,10,15,14,10,17,23,28,32,33,39,41,53,58,70,67,64,77,95,95,
%T 85,97,91,100,119,129,158,146,141,148,155,160,176,181,190,197,226,238,
%U 268,273,256,281,323,321,321,308,313,328,376,365,372,359,350,381,430
%N Sum of primes < n that do not divide n.
%C Equals row sums of triangle A143655. - _Gary W. Adamson_, Aug 28 2008
%H Michael De Vlieger, <a href="/A066911/b066911.txt">Table of n, a(n) for n = 1..10000</a>
%F a(prime(n)) = A007504(n-1). - _Michel Marcus_, May 20 2014
%e a(8) = 3 + 5 + 7 = 15 because 3, 5 and 7 are the primes < 8 that do not divide 8.
%t Table[Total@ Complement[Prime@ Range@ PrimePi@ n, FactorInteger[n][[All, 1]]], {n, 60}] (* _Michael De Vlieger_, Sep 04 2017 *)
%o (PARI) a(n) = sum(i=1, n-1, if (isprime(i) && (n%i), i)); \\ _Michel Marcus_, May 20 2014
%Y Cf. A007504, A143655.
%K nonn
%O 1,3
%A _Leroy Quet_, Jan 22 2002
|
