|
%I #16 Nov 21 2023 05:20:47
%S 3659401,8453041,10319761,14155681,16391761,18094561,19616689,
%T 20456329,21677041,22628929,27275161,32051881,34228489,37728601,
%U 38884561,39191881,40101071,40167241,42163969,47931601,48461449,49460161,50389441,54932329,56219281,57590569
%N Primes with 71 as smallest positive primitive root.
%H Robert Price, <a href="/A114686/b114686.txt">Table of n, a(n) for n = 1..1366</a>
%t t={}; Do[ If[ PrimitiveRoot[ Prime@n] == 71, AppendTo[t, n]; Print@ Prime@n], {n, 3280000}]; Prime@t
%o (PARI) is(n)=if(n<72,return(0));for(k=2,70,if(znorder(Mod(k,n))==n-1,return(0)));znorder(Mod(71,n))==n-1&&isprime(n) \\ _Charles R Greathouse IV_, Jul 19 2011
%o (PARI) is(n)=isprime(n)&&lift(znprimroot(n))==71 \\ relies on implementation details, may not always work
%o \\ _Charles R Greathouse IV_, Jul 19 2011
%Y Cf. A001122-A001126, A061323-A061335, A061730-A061741, A114657-A114686.
%K nonn
%O 1,1
%A _Robert G. Wilson v_, Dec 21 2005
%E a(23) and beyond from _Robert Price_, Nov 20 2023
|
