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”).

A238735
Number of prime pairs {2^n + (2k + 1), (2k + 1)*2^n + 1}, k < n.
1
1, 2, 1, 2, 0, 3, 2, 2, 0, 1, 0, 2, 0, 0, 0, 3, 0, 1, 0, 3, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
OFFSET
1,2
COMMENTS
If k = 0, then the two numbers in the "prime pair" are actually the same number, 2^n + 1 (which is either 2 or a Fermat prime; see A019434, A092506).
LINKS
EXAMPLE
a(1) = 1 because 2^1+(2*0+1)=3 and (2*0+1)*2^1+1=3 is prime pair for k=0,
a(2) = 2 because 2^2+(2*0+1)=5 and (2*0+1)*2^2+1=5 is prime pair for k=0, 2^2+(2*1+1)=7 and (2*1+1)*2^2+1=13 is prime pair for k=1,
a(3) = 1 because 2^3+(2*2+1)=13 and (2*2+1)*2^3+1=41 is prime pair for k=2.
MATHEMATICA
a[n_] := Length@Select[Range[0, n-1], PrimeQ[2^n + (2*# + 1)] && PrimeQ[(2*# + 1)*2^n + 1] &]; Array[a, 100] (* Giovanni Resta, Mar 04 2014 *)
PROG
(PARI) a(n)=sum(k=0, n-1, isprime(2^n+2*k+1)&&isprime((2*k+1)<<n+1)) \\ Charles R Greathouse IV, Mar 06 2014
CROSSREFS
Sequence in context: A199920 A177995 A332104 * A356006 A258120 A147786
KEYWORD
nonn
AUTHOR
EXTENSIONS
a(47)-a(87) from Giovanni Resta, Mar 04 2014
STATUS
approved