The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A158846 Primes which are removed with the algorithm of A156284, starting the selection with the interval (2^4, 2^5). 4

%I

%S 19,29,41,47,53,59,61,97,149,167,173,233,239,251,271,283,313,331,349,

%T 373,409,433,439,499,509,521,557,563,593,641,677,743,761,797,827,887,

%U 911,941,953,1013,1019,1021,1039,1051,1129,1171,1237,1279,1291

%N Primes which are removed with the algorithm of A156284, starting the selection with the interval (2^4, 2^5).

%C We iteratively scan integer intervals (2^(m-1)..2^m), first the one with m=5, then m=6, m=7, etc., and start with the set S={3,5,7,11,...} of all odd primes. For each prime p = 2^m-k, 2^(m-1) < p < 2^m, p is removed from S if k is in S. Basically, all the upper primes of primes pairs are removed when the prime pair sums to a power of 2 which are larger than 2^4. The sequence shows all p that are removed from S at any stage m.

%C Powers 2^m, m >= 5, are not expressible as sums of two primes which are not in the sequence.

%p A158846 := proc()

%p local mmax,prrem,m,prm,pi,p,q ;

%p mmax := 12 ; prrem := {} ;

%p for m from 5 to mmax do

%p prm := {} ;

%p for pi from 1 do

%p k := ithprime(pi) ;

%p p := 2^m-k ;

%p if p <= 2^(m-1) then break; end if;

%p if isprime(p) and not k in prrem then prm := prm union {p} ;

%p end if ;

%p end do:

%p prrem := prrem union prm ;

%p end do: print( sort(prrem)) ; return ;

%p end proc:

%p A158846() ; # _R. J. Mathar_, Dec 07 2010

%Y Cf. A156284, A158756, A156759.

%K nonn

%O 1,1

%A _Vladimir Shevelev_, Mar 28 2009

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified April 6 06:21 EDT 2020. Contains 333267 sequences. (Running on oeis4.)