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

A138578
Numbers k such that 2^(2*k-7)-7 is prime.
0
23, 361, 995, 1163, 1253, 1891, 6413, 31795, 60781, 60923
OFFSET
1,1
EXAMPLE
2^(2*23-7) - 7 = 549755813881 is prime.
PROG
(PARI) is(n)=isprime(2^(2*n-7)-7) \\ Charles R Greathouse IV, Jun 13 2017
(Python)
from sympy import isprime
def afind(limit):
k = 5; pow2 = 2**(2*k-7)
while k < limit:
while not isprime(pow2 - 7) and k <= limit: k += 1; pow2 *= 4
if k > limit: return
print(k, end=", "); k += 1; pow2 *= 4
afind(2000) # Michael S. Branicky, Apr 23 2021
CROSSREFS
Sequence in context: A025964 A180361 A021684 * A077515 A025942 A020000
KEYWORD
nonn,more
AUTHOR
EXTENSIONS
a(6)-a(7) from Michael S. Branicky, Apr 23 2021
a(8)-a(10) from Michael S. Branicky, Sep 12 2024
STATUS
approved