|
|
A370606
|
|
Primes p such that valuation(p-1,2) is a record.
|
|
2
|
|
|
2, 3, 5, 17, 97, 193, 257, 7681, 12289, 40961, 65537, 786433, 5767169, 7340033, 23068673, 104857601, 167772161, 469762049, 2013265921, 3221225473, 75161927681, 77309411329, 206158430209, 2061584302081, 2748779069441, 6597069766657, 39582418599937, 79164837199873
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
LINKS
|
|
|
PROG
|
(PARI) r=-1; forprime(p=2, 10^12, v=valuation(p-1, 2); if(v>r, print1(p, ", "); r=v))
(Python)
from itertools import count, islice
from sympy import isprime
def A370606_gen(): # generator of terms
a = 1
while True:
for q in count(a, a):
if isprime(q+1):
yield q+1
a = (q&-q)<<1
break
|
|
CROSSREFS
|
Cf. A084924 (valuation(p+1,2) is a record).
Cf. A370607 (corresponding 2-valuations).
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|