On Fri, 02 Dec 2011 01:36:45 -0500, quasi wrote:
> On Thu, 1 Dec 2011 19:47:24 -0800, William Elliot wrote:
>>On Thu, 1 Dec 2011, quasi wrote:
>>> Conjecture:
>>>
>>> For all primes q > 2, there exists a prime p < q such that
>>> ord_p(q)=p-1 and ord_q(p)=q-1.
>>
>>What's ord_p mean?
> <
http://en.wikipedia.org/wiki/Multiplicative_order>
>
>>> Question:
>>>
>>> For the above conjecture, which of the following, if any,
>>> characterizes its status?
>>>
>>> (a) It's known to be true.
>>>
>>> (b) It's known to be false.
>>>
>>> (c) It's a known conjecture, but as yet unsettled.
I don't know which of those applies; most likely (c).
Computationally it's easy to verify as true for p, q < 500000,
via pari/gp code between -- lines below . For each q from
2...qmax, the code counts the number of mutual primitive roots
for primes p, q, with p < q and p < pmax. Eg, when pmax=523
and qmax=99991, there are at least q/18 solutions for each q<600,
and from 25 to 55 solutions for each q in {600..100000}; when
pmax=1979 and qmax=499979, there are at least q/21 solutions for
each q<1979, and from 87 to 144 solutions for q in {1979..500000}.
--
/* 4 Dec 2011 jiw: File "mutual-roots", re: mutual primitive
roots for primes p, q, with p < q, ie p primitive mod q and q
primitive mod p.
For production, use bigger numbers, eg:
s=299; pmax=prime(s); qmax=precprime(500000);
Shell code to filter for gnuplot:
gp < mutual-roots | grep '^[0-9]*.[0-9]*$' > t
In gnuplot use commands like:
plot [0:100000][0:100] 't', x/18, 25,55
plot [0:500000][0:300] 't', x/21, 87,144
plot [0:500000][0:300] 't', x/21, 87,144
*/
s=99; pmax=prime(s); qmax=precprime(100000);
b=vectorsmall(pmax); i=t=1;
forprime(p=2, pmax, b[p]=t; t+=p);
r=vectorsmall(t);
forprime(p=2, pmax, c=b[p]; for(i=1,p-1,if(znorder(Mod(i,p))==p-1,r[c+i]=1)));
forprime(q=3, qmax, g=0; qx=min(pmax,q-1); forprime(p=2, qx, g+=r[b[p]+(q%p)]); print(q,"\t",g));
--
jiw