login
Primes whose smallest positive primitive root (A001918) is not prime.
4

%I #26 May 18 2013 16:49:26

%S 2,41,109,151,229,251,271,313,337,367,409,439,733,761,971,991,1021,

%T 1031,1069,1289,1297,1303,1429,1471,1489,1759,1783,1789,1811,1871,

%U 1873,1879,2137,2411,2441,2551,2749,2791,2971,3001,3061,3079,3109,3221,3229

%N Primes whose smallest positive primitive root (A001918) is not prime.

%C Subsequence of A222717 = primes whose smallest positive quadratic nonresidue is not a primitive root. (Proof. If p is not in A222717, then the smallest positive quadratic nonresidue of p is a primitive root g. Since the smallest positive quadratic nonresidue is always a prime, g is prime. But since all primitive roots are quadratic nonresidues, g is the smallest positive primitive root of p. Hence p is not in A047936.) - _Jonathan Sondow_, Mar 13 2013.

%H Charles R Greathouse IV, <a href="/A047936/b047936.txt">Table of n, a(n) for n = 1..10000</a>

%H <a href="/index/Pri#primes_root">Index entries for primes by primitive root</a>

%t lst={}; Do[p=Prime[n]; pr=PrimitiveRoot[p]; If[pr>1&&!PrimeQ[pr], AppendTo[lst, p]], {n, 7!}]; lst (* _Vladimir Joseph Stephan Orlovsky_, Oct 24 2009 *)

%t Select[Prime[Range[500]],!PrimeQ[PrimitiveRoot[#]]&] (* _Harvey P. Dale_, Oct 24 2011 *)

%o (PARI) select(p->!isprime(lift(znprimroot(p))),primes(999)) \\ reverse order of arguments if using an old version of GP

%o \\ _Charles R Greathouse_ IV, Oct 24 2011

%Y Cf. A222717, A223036.

%K nonn,easy

%O 1,1

%A _Felice Russo_

%E More terms from _James A. Sellers_, Dec 22 1999