login
a(n) = sum of nonprimes <= n.
6

%I #24 Jan 08 2024 02:24:30

%S 0,1,1,1,5,5,11,11,19,28,38,38,50,50,64,79,95,95,113,113,133,154,176,

%T 176,200,225,251,278,306,306,336,336,368,401,435,470,506,506,544,583,

%U 623,623,665,665,709,754,800,800,848,897,947,998,1050,1050,1104,1159,1215

%N a(n) = sum of nonprimes <= n.

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

%F a(n) = A000217(n) - A034387(n) = A101256(n) + 1.

%t Accumulate[Table[If[PrimeQ[n],0,n],{n,0,60}]] (* _Harvey P. Dale_, Oct 02 2020 *)

%o (PARI) my(s=0); for(k=0,100,if(!isprime(k),s+=k);print1(s", ")); \\ _Cino Hilliard_, Feb 04 2006

%o (Haskell)

%o a101203 n = a101203_list !! (n-1)

%o a101203_list = scanl (+) 0 $ zipWith (*) [1..] $ map (1 -) a010051_list

%o -- _Reinhard Zumkeller_, Oct 10 2013

%Y Cf. A062298, A018252.

%Y Partial sums of A191558.

%Y Cf. A000217, A034387, A101256.

%K nonn

%O 0,5

%A _Reinhard Zumkeller_, Jan 23 2005