STATUS
proposed
approved
proposed
approved
editing
proposed
if isprime((k:=(i<<1)+1):
return k-(m<<1)+1 # Chai Wah Wu, Jan 27 2023
proposed
editing
editing
proposed
a(n) == 3 (mod 4). - Chai Wah Wu, Jan 27 2023
proposed
editing
editing
proposed
return k-(m<<1)+1 # Chai Wah Wu, Jan 27 2023
(Python)
from sympy import isprime, nextprime
def A360081(n):
m = 1<<3*(1<<n)-1
i = m
while i:=nextprime(i):
if isprime((k:=i<<1)+1):
return k-(m<<1)+1 # Chai Wah Wu, Jan 27 2023
a(n) == 3 mod 4. - Chai Wah Wu, Jan 27 2023
approved
editing
proposed
approved