|
%I #57 Jun 08 2022 14:31:14
%S 1,2,3,11,17,23,27,29,35,37,41,47,51,53,57,59,65,67,71,77,79,83,87,89,
%T 93,95,97,101,107,113,117,119,121,123,125,127,131,135,137,143,145,147,
%U 149,155,157,161,163,167,171,173,177,179,185,187,189,191,197,203,205,207,209
%N Numbers that are not the sum of 2 primes.
%C Suggested by the Goldbach conjecture that every even number larger than 2 is the sum of 2 primes.
%C Since (if we believe the Goldbach conjecture) all the entries > 2 in this sequence are odd, they are equal to 2 + an odd composite number (or 1).
%C Otherwise said, the sequence consists of 2 and odd numbers k such that k-2 is not prime. In particular there is no element from A006512, greater of a twin prime pair. - _M. F. Hasler_, Sep 18 2012
%C Values of k such that A061358(k) = 0. - _Emeric Deutsch_, Apr 03 2006
%C Values of k such that A073610(k) = 0. - _Graeme McRae_, Jul 18 2006
%D G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, 5th ed., Oxford Univ. Press, 1979, Section 2.8 (for Goldbach conjecture).
%H Reinhard Zumkeller, <a href="/A014092/b014092.txt">Table of n, a(n) for n = 1..10000</a>
%H <a href="/index/Go#Goldbach">Index entries for sequences related to Goldbach conjecture</a>
%F Odd composite numbers + 2 (essentially A014076(n) + 2 ).
%F Equals {2} union A005408 \ A052147, i.e., essentially the complement of A052147 (or rather A048974) within the odd numbers A005408. - _M. F. Hasler_, Sep 18 2012
%p g:=sum(sum(x^(ithprime(i)+ithprime(j)),i=1..j),j=1..50): gser:=series(g,x=0,230): a:=proc(n) if coeff(gser,x^n)=0 then n else fi end: seq(a(n),n=1..225); # _Emeric Deutsch_, Apr 03 2006
%t s1falsifiziertQ[s_]:= Module[{ip=IntegerPartitions[s, {2}], widerlegt=False},Do[If[PrimeQ[ip[[i,1]] ] ~And~ PrimeQ[ip[[i,2]] ], widerlegt = True; Break[]],{i,1,Length[ip]}];widerlegt]; Select[Range[250],s1falsifiziertQ[ # ]==False&] (* _Michael Taktikos_, Dec 30 2007 *)
%t Join[{1,2},Select[Range[3,300,2],!PrimeQ[#-2]&]] (* _Zak Seidov_, Nov 27 2010 *)
%t Select[Range[250],Count[IntegerPartitions[#,{2}],_?(AllTrue[#,PrimeQ]&)]==0&] (* _Harvey P. Dale_, Jun 08 2022 *)
%o (PARI) isA014092(n)=local(p,i) ; i=1 ; p=prime(i); while(p<n, if( isprime(n-p), return(0)); i++; p=prime(i)); 1
%o n=1; for(a=1,200, if(isA014092(a), print(n," ",a); n++)) \\ _R. J. Mathar_, Aug 20 2006
%o (Haskell)
%o a014092 n = a014092_list !! (n-1)
%o a014092_list = filter (\x ->
%o all ((== 0) . a010051) $ map (x -) $ takeWhile (< x) a000040_list) [1..]
%o -- _Reinhard Zumkeller_, Sep 28 2011
%o (Python)
%o from sympy import prime, isprime
%o def ok(n):
%o i=1
%o x=prime(i)
%o while x<n:
%o if isprime(n - x): return False
%o i+=1
%o x=prime(i)
%o return True
%o print([n for n in range(1, 301) if ok(n)]) # _Indranil Ghosh_, Apr 29 2017
%Y Cf. A002372, A002373, A002374, A048974, A061358.
%Y Cf. A010051, A000040, A051035 (composites).
%Y Equivalent sequence for prime powers: A071331.
%Y Numbers that can be expressed as the sum of two primes in k ways for k=0..10: this sequence (k=0), A067187 (k=1), A067188 (k=2), A067189 (k=3), A067190 (k=4), A067191 (k=5), A066722 (k=6), A352229 (k=7), A352230 (k=8), A352231 (k=9), A352233 (k=10).
%K nonn,nice,easy
%O 1,2
%A _N. J. A. Sloane_
|
