login
Sum of first n composite numbers.
23

%I #49 Feb 17 2018 18:32:27

%S 0,4,10,18,27,37,49,63,78,94,112,132,153,175,199,224,250,277,305,335,

%T 367,400,434,469,505,543,582,622,664,708,753,799,847,896,946,997,1049,

%U 1103,1158,1214,1271,1329,1389,1451,1514,1578,1643,1709,1777,1846,1916,1988

%N Sum of first n composite numbers.

%C a(A196415(n)) = A036691(A196415(n)) / A141092(n). - _Reinhard Zumkeller_, Oct 03 2011

%H Reinhard Zumkeller, <a href="/A053767/b053767.txt">Table of n, a(n) for n = 0..10000</a>

%H C. Rivera, <a href="http://www.primepuzzles.net/puzzles/puzz_089.htm">See also</a>

%F a(n) = A000217(A002808(n)) - A034387(A002808(n)) - 1 . - _Robert Israel_, Jan 09 2015

%F a(n) = A051349(n+1) - 1. - _Michel Marcus_, Feb 16 2018

%p A053767 := proc(n)

%p add(A002808(i),i=1..n) ;

%p end proc: # _R. J. Mathar_, Oct 03 2011

%p ListTools[PartialSums](remove(isprime,[$2..1000])); # _Robert Israel_, Jan 09 2015

%t lst={};s=0;Do[If[ !PrimeQ[n], s=s+n;AppendTo[lst, s]], {n, 2, 10^2}];lst (* _Vladimir Joseph Stephan Orlovsky_, Aug 14 2008 *)

%t Accumulate[Complement[Range[2,100],Prime[Range[PrimePi[100]]]]] (* _Harvey P. Dale_, Dec 28 2010 *)

%t Accumulate[Select[Range[2, 100], ! PrimeQ[#] &]]

%o (Haskell)

%o a053767 n = a053767_list !! (n-1)

%o a053767_list = scanl1 (+) a002808_list -- _Reinhard Zumkeller_, Oct 03 2011

%o (PARI) lista(nn) = {my(s=0); forcomposite(n=0, nn, print1(s, ", "); s += n;);} \\ _Michel Marcus_, Jan 09 2015

%Y Cf. A002808, A051349.

%Y First differences of A023539.

%K nonn

%O 0,2

%A _G. L. Honaker, Jr._, Mar 29 2000

%E a(0)=0 prepended by _Max Alekseyev_, Feb 10 2018