OFFSET
1,1
LINKS
Wikipedia, Primitive root modulo n.
EXAMPLE
For m = 41 the roots of x^2 - x - 1 (mod 41) are 7 and 35. 7 and 35 are both primitive roots modulo 41.
MATHEMATICA
test[p_]:=Module[{inv2, sqr}, If[JacobiSymbol[5, p]==1, inv2=ModularInverse[2, p]; sqr=PowerMod[5, 1/2, p]; {MultiplicativeOrder[Mod[inv2*(sqr-1), p], p], MultiplicativeOrder[Mod[inv2*(-sqr-1), p], p]} == {p-1, p-1}, False]]; Cases[Prime[Range[4, 5000]], _?(test[#] &)] (* Shenghui Yang, Jun 01 2025 *)
PROG
(PARI) { forprime(p=2, 3600, s=polrootsmod(x^2 - x - 1, p);
if( #s==2 && p-1==znorder(Mod(s[1], p)) && p-1==znorder(Mod(s[2], p)),
print1(p, ", "); ); ); } \\ Joerg Arndt, May 31 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Jay Anderson, May 31 2025
STATUS
approved
