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”).
%I #43 Oct 13 2019 11:21:47
%S 3,5,7,13,23,29,31,37,47,61,101,127,229,269,431,503,509,997,1021,1319,
%T 2039,3853,4093,7949,8111,8191,14197,16141,16381,32687,45853,65293,
%U 130343,130829,131063,131071,347141,502829,524261,524287,1028893,1046389,1048549
%N Primes of the form 2^j - 3^k, for j >= 0, k >= 0.
%C The numbers in A007643 are not in this sequence.
%C For n > 1, a(n) is of the form 8k - 1 or 8k - 3.
%C In this sequence, only 3 and 7 make both j and k even numbers.
%C Generally, the way to prove that a number is not in this sequence is to successively take residues modulo 3, 8, 5, and 16 on both sides of the equation 2^j - 3^k = x.
%H H. Gauchman and I. Rosenholtz (Proposers), R. Martin (Solver), <a href="https://www.jstor.org/stable/2691457">Difference of prime powers, Problem 1404</a>, Math. Mag., 65 (No. 4, 1992), 265; <a href="https://www.jstor.org/stable/2690747">Solution</a>, Math. Mag., 66 (No. 4, 1993), 269.
%F Intersection of A000040 and A192110.
%e 7 = 2^3 - 3^0, so 7 is a term.
%o (PARI) forprime(p=1,1000,k=0;x=2;y=1;while(k<p+1,while(x<y+p,x=2*x);if(x-y==p,print1(p,", ");k=p);k++;y=3*y))
%Y Cf. A000040, A007643, A192110.
%Y Cf. A004051 (primes of the form 2^a + 3^b).
%Y Cf. A063005.
%K nonn
%O 1,1
%A _Jinyuan Wang_, Nov 16 2018
%E More terms from _Alois P. Heinz_, Nov 16 2018