login
A339566
Primes p such that A007088(p) == 1 (mod p).
1
5, 137, 3967, 25087, 242899421
OFFSET
1,1
EXAMPLE
a(3) = 3967 is in the sequence because 3967 = 111101111111_2 and 111101111111 == 1 (mod 3967).
MAPLE
p:= 1: R:= NULL:
while p < 3*10^8 do
p:= nextprime(p);
if convert(p, binary) mod p = 1 then R:= R, p fi
od:
R;
PROG
(Python)
from sympy import nextprime
A339566_list, p = [], 2
while p < 10**10:
if int(bin(p)[2:]) % p == 1:
A339566_list.append(p)
p = nextprime(p) # Chai Wah Wu, Dec 14 2020
CROSSREFS
Primes in A339567.
Sequence in context: A265875 A248523 A159546 * A012215 A012142 A012207
KEYWORD
nonn,more
AUTHOR
Robert Israel, Dec 09 2020
STATUS
approved