login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A301776
Prime numbers p with the property that all even numbers n (2 < n <= 2p) are the sum of two primes <= p.
3
2, 3, 5, 7, 13, 19, 109
OFFSET
1,1
COMMENTS
Conjecture: this sequence is finite - it has 7 terms only. Conjecture verified up to first 10^5 primes.
This sequence is related to the Goldbach Strong Conjecture.
LINKS
Marcin Barylski, C++ program
Marcin Barylski, Sum building from first primes vs. theoretical maximum - the first 200 rounds of the algorithm. If red and green lines ever touch each other, we have a new term.
EXAMPLE
a(1)=2 because all even numbers 2 < n <= 2*2 (there is just one such number: 4) can be expressed as a sum of 2 only: 4=2+2.
a(2)=3 because 4=2+2, 6=3+3.
a(3)=5 because 4=2+2, 6=3+3, 8=5+3, 10=5+5.
a(4)=7 because 4=2+2, 6=3+3, 8=5+3, 10=5+5, 12=5+7, 14=7+7.
a(5)=13 (and is not 11) because 20 cannot be expressed as a sum of two primes from a set {2,3,5,7,11} but all even numbers 2 < n <= 26 can be expressed as a sum of two primes from a set {2,3,5,7,11,13}.
MATHEMATICA
Select[Prime@ Range[500], Function[p, SameQ[Select[Union@ Map[Total, Tuples[Prime@ Range@ PrimePi@ p, 2]], And[EvenQ@ #, # > p] &], Range[p + 1 + Boole@ EvenQ@ p, 2 p, 2]]]] (* Michael De Vlieger, Apr 10 2018 *)
PROG
(C++) See Barylski link.
(PARI) isok(p) = {vp = primes(primepi(p)); slist = List(); for (i=1, #vp, for (j=1, i, if (!((vp[i]+vp[j]) % 2), listput(slist, vp[i]+vp[j])); ); ); #Set(slist) == (p-1); }
lista(nn) = forprime(p=2, nn, if (isok(p), print1(p, ", "))); \\ Michel Marcus, Apr 09 2018
CROSSREFS
Cf. A002372 (number of ordered Goldbach partitions).
Sequence in context: A341650 A341640 A104189 * A178570 A294443 A293160
KEYWORD
nonn,more
AUTHOR
Marcin Barylski, Mar 26 2018
STATUS
approved