OFFSET
1,1
MATHEMATICA
Select[Range[3600], PrimeQ[#] && OddQ[MultiplicativeOrder[2, (# - 1)/2^IntegerExponent[# - 1, 2]]] &] (* Amiram Eldar, Nov 26 2021 *)
PROG
(PARI) isok(p) = isprime(p) && znorder(Mod(2, (p-1)/2^valuation(p-1, 2)))%2;
CROSSREFS
KEYWORD
nonn
AUTHOR
Arkadiusz Wesolowski, Nov 24 2021
STATUS
approved