NTIC A Better Way to Primitive Roots

Section 10.2 A Better Way to Primitive Roots

Subsection 10.2.1 A useful lemma

In order to find primitive roots, we might want a better approach than simply trying every single power of

a

for every

a

until we find one. Let’s walk through an example to motivate a new approach, using a small modulus.

🔗

Example 10.2.1. A motivating example.

🔗

Let’s take a number

n

such that

ϕ (n)

has some, but not too many, factors – say,

n = 11,

ϕ (11) = 10 .

Okay, we know that every element

a \in U_{11}

will have

a^{10} \equiv 1 (mod 11),

🔗

but which elements don’t reach the unit before the tenth power?

🔗

We know by Theorem 8.3.12 that the order of an element has to divide

ϕ (11) = 10,

so we could try

a^{2}

and

a^{5};

no other

a^{k}

could yield 1. In fact, if those aren’t

\equiv 1,

there aren’t any other possible orders out there, so that

a

would work as a primitive root.

Let’s try this with $a = 2 .$

$2^{2} \equiv 4 ≢ 1 (mod 11) and 2^{5} = 32 \equiv - 1 ≢ 1 (mod 11),$

so $2$ must be a primitive root.
What about with $a = 3 ?$

$3^{2} = 9 ≢ 1 (mod 11),$

so that seems promising, but

$3^{5} = 9 \cdot 9 \cdot 3 \equiv (- 2)^{2} \cdot 3 \equiv 12 \equiv 1 (mod 11)$

so $3$ cannot be a primitive root modulo eleven (and in fact has order five).

🔗

The moral is that we didn’t have to check all ten possible powers of

a = 2

a = 3

to decide whether

a

was a primitive root modulo eleven. If you aren’t confident of this idea, try using this strategy to determine which of

a = 4, 5, 6

is a primitive root (exactly one of them is).

🔗

Now we formalize and rephrase our strategy slightly more efficiently.

🔗

Sage note 10.2.2. How Sage does primitive roots.

🔗

As far as I understand, the following strategy is how even Sage tests for finding primitive roots, at least for basic cases. You can check for yourself by looking at the code from the component program, PARI

pari.math.u-bordeaux.fr/cgi-bin/gitweb.cgi?p=pari.git;a=blob;f=src/basemath/arith1.c

; look for is_gener_expo and is_gener_Fp.

🔗

Lemma 10.2.3. Testing for Primitive Roots.

🔗

An element

a \in U_{n}

is a primitive root if and only if

a^{ϕ (n) / q} ≢ 1 in U_{n} for each prime q ∣ ϕ (n) .

Proof.

a

is in fact a primitive root, then

ϕ (n)

is the smallest number

k

such that

a^{k} \equiv 1,

so certainly for numbers smaller than

ϕ (n),

ϕ (n) / q,

those powers shouldn’t be

\equiv 1 .

On the other hand, if

a

isn’t a primitive root, then its order

k

must be a proper divisor of

ϕ (n) .

Now look at the prime divisors

q

ϕ (n) / k .

For such a divisor,

q ∣ ϕ (n) / k so q k ℓ = ϕ (n) for some ℓ \in Z .

That means

ϕ (n) / q = k ℓ

and so the power

ϕ (n) / q

in the statement is actually a multiple of the order

k .

Since

a^{k} \equiv 1,

then certainly

a^{k ℓ} = a^{ϕ (n) / q} \equiv 1 (mod n)

as well, which completes the proof.

🔗

This proof is a little terse, so let’s unpack this test. Essentially, we change two things from the initial idea of trying all divisors of

ϕ (n) :

Instead of trying powers which are divisors of $ϕ (n),$ we try powers which are $ϕ (n)$ divided by divisors. So $2^{5}$ becomes $2^{10 / 2}$ and $3^{2}$ becomes $3^{10 / 5} .$ That seems like it’s not doing anything other than rewriting, but at least it organizes things differently.
Then, instead of having to try all $ϕ (n) / d,$ we use a trick to just need prime divisors $d .$ (See the proof.)

🔗

Doing some examples slowly will help it make sense. Once you have done so, try the interact.


    
        
xxxxxxxxxx
 
1
@interact
2
def _(n=(19,[2..100]),a=3):
3
    phi=euler_phi(n)
4
    pds=prime_divisors(phi)
5
    if gcd(a,n)!=1:
6
        pretty_print(html("Make sure $a$ and $n$ are relatively prime!"))
7
    else:
8
        a = mod(a,n)
9
        pretty_print(html("Is $%s$ a primitive root of $%s$?"%(a,n)))
10
        pretty_print(html(r"The prime divisors of $\phi(%s)=%s$ are $%s$"%(n, euler_phi(n), ','.join([str(pd) for pd in pds]))))
11
        pretty_print(html("The powers are "+' and '.join([r'$%s^{%s/%s}\equiv %s$'%(a,phi,pd,a^(phi/pd)) for pd in pds])))
12
        pretty_print(html("And the order of a=$%s$ is <tt>a.multiplicative_order()</tt>=$%s$"%( a , a.multiplicative_order())))

    
    
    
    
        
            
                Language:
                
            
        
    
    




    
    
        
        Messages

🔗

Subsection 10.2.2 Using the test lemma

🔗

If you tried various

n

and various attempts at primitive roots

a

above, you will see that Lemma 10.2.3 really works. Make sure you are trying

a

that are actually coprime to

n,

though! As it turns out, there aren’t very many test powers to try, since

ϕ (n)

in general doesn’t have a lot of prime divisors, even if

n

is a fairly large prime.

🔗

Why not try it by hand for

n = 17 ?

There is only one prime divisor of

ϕ (17),

which makes things easier. Fill in Table 10.2.4, where PR means primitive root.

🔗

Table 10.2.4. Check which numbers are primitive roots

$a$	1	2	3	4	5	6	7	8	9	10	11	12	13	14	15	16
PR?	No															No

🔗

The lemma also makes easy some statements that would otherwise be quite hard. For instance, you should (Exercise 10.6.2) see how to use the test lemma to prove that if

a

is a primitive root of

n,

then so is

a^{- 1}

(modulo

n

🔗

Here’s something harder, to show the power of this approach.

🔗

Proposition 10.2.5.

🔗

a

is a primitive root of

n,

then so is

n - a

4 ∣ ϕ (n) .

Proof.

Let’s think in terms of powers. If

a^{ϕ (n) / q} ≢ 1,

then

(n - a)^{ϕ (n) / q} \equiv (- a)^{ϕ (n) / q} \equiv (- 1)^{ϕ (n) / q} a^{ϕ (n) / q} .

So, as long as

ϕ (n) / q

is even for all prime divisors of

ϕ (n),

the two powers (the one of

a

and the one of

n - a

) come to the same thing.

Since

ϕ (n)

is already assumed to be even, the only possible odd

ϕ (n) / q

comes from

q = 2,

but

ϕ (n)

is assumed to be divisible by four, so

ϕ (n) / q

will be even.

🔗

Example 10.2.6.

🔗

If you did the table at the beginning of this subsection properly, you will note that

3

and

14

are a pair of primitive roots of seventeen. There should be three other such pairs.

🔗

On the other hand, from the proof we can see that if

ϕ (n)

is even, but not divisible by four, then we expect that if

a

is a PR then

n - a

will not be. For example, since two is a primitive root of eleven in Example 10.2.1, we expect that nine is not; try computing this yourself.

Prev Top Next