OFFSET
1,1
COMMENTS
FORMULA
a(n) ~ (3/2)*n*log((3/2)*n).
PROG
(Magma) [p: p in [3..419 by 2] | IsPrime(p) and (p-1)/Modorder(2, p) le 2];
(PARI) isok(p) = p%2 && isprime(p) && (p-1)/znorder(Mod(2, p))<=2;
(Python)
from itertools import islice
from sympy import nextprime, n_order
def A363286_gen(startvalue=3): # generator of terms >= startvalue
p = max(startvalue, 3)-1
while (p:=nextprime(p)):
if n_order(2, p)<<1 >= p-1:
yield p
CROSSREFS
KEYWORD
nonn
AUTHOR
Arkadiusz Wesolowski, May 25 2023
STATUS
approved