cp's OEIS Frontend

A070187 Primes p such that x^11 = 2 has a solution mod p, but x^(11^2) = 2 has no solution mod p.

%I A070187 #8 Sep 21 2012 02:44:51
%S A070187 12101,13553,30493,32429,44771,66067,103577,128987,180533,182711,
%T A070187 187793,201829,242243,257489,264749,299113,314359,330331,337349,
%U A070187 341947,356467,371471,431729,442619,475289,484243,505781,513767,540871,558053,564103,573299,581527,582011,586367,593869,596047,630169
%N A070187 Primes p such that x^11 = 2 has a solution mod p, but x^(11^2) = 2 has no solution mod p.
%o A070187 (PARI) forprime(p=2,550000,x=0; while(x<p&&x^11%p!=2%p,x++); if(x<p,y=0; while(y<p&&y^(11^2)%p!=2%p,y++); if(y==p,print1(p,","))))
%o A070187 (PARI)
%o A070187 N=10^6;  default(primelimit,N);
%o A070187 ok(p, r, k1, k2)={
%o A070187     if (  Mod(r,p)^((p-1)/gcd(k1,p-1))!=1, return(0) );
%o A070187     if (  Mod(r,p)^((p-1)/gcd(k2,p-1))==1, return(0) );
%o A070187     return(1);
%o A070187 }
%o A070187 forprime(p=2,N, if (ok(p,2,11,11^2),print1(p,", ")));
%o A070187 /* _Joerg Arndt_, Sep 21 2012 */
%Y A070187 Cf. A049543, A059667, A070179 - A070186, A070188.
%K A070187 nonn
%O A070187 1,1
%A A070187 _Klaus Brockhaus_, Apr 29 2002