OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
filter:= proc(p) local k;
if not isprime(p) then return false fi;
k:= (p-1)/2^padic:-ordp(p-1, 2);
numtheory:-order(2, k)::odd
end proc:
select(filter, [2, seq(i, i=3..10000, 2)]); # Robert Israel, Feb 02 2025
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