cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A213650 Numbers k such that the sum of the first k primes is semiprime.

This page as a plain text file.
%I A213650 #22 Mar 30 2021 12:43:36
%S A213650 3,7,8,10,16,18,22,28,32,34,36,38,44,46,48,54,55,58,59,65,66,72,75,82,
%T A213650 92,93,94,104,106,110,118,120,133,136,137,138,140,141,142,144,148,150,
%U A213650 154,156,164,168,170,174,190,194,202,210,212,218,224,226,232,234
%N A213650 Numbers k such that the sum of the first k primes is semiprime.
%C A213650 Numbers k such that A007504(k) is included in A001358.
%H A213650 Amiram Eldar, <a href="/A213650/b213650.txt">Table of n, a(n) for n = 1..10000</a>
%e A213650 8 is in the sequence because the sum of the first 8 primes is  2 + 3 + 5 + 7 + 11 + 13 + 17 + 19 = 77 = 7*11, which is semiprime.
%p A213650 with(numtheory): for n from 1 to 500 do:s:=sum(‘ithprime(k)’, ’k’=1..n):if bigomega(s)=2 then printf(`%d, `, n):else fi:od:
%t A213650 Flatten[Position[Accumulate[Prime[Range[300]]],_?(PrimeOmega[#]==2&)]]
%o A213650 (PARI) isok(n) = bigomega(vecsum(primes(n))) == 2; \\ _Michel Marcus_, Sep 18 2017
%Y A213650 Cf. A001358, A007504, A013916, A092189 (numbers n such that sum of first n semiprimes is a semiprime), A092190 (semiprimes that are the sum of first n semiprimes for some n), A180152 (numbers n such that the sum of the first n semiprimes is a prime).
%K A213650 nonn
%O A213650 1,1
%A A213650 _Michel Lagneau_, Jun 17 2012