19.07.2014 Views

Contents - Student subdomain for University of Bath

Contents - Student subdomain for University of Bath

Contents - Student subdomain for University of Bath

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

126 CHAPTER 4. MODULAR METHODS<br />

This notation is analogous to the previous section, where A p was the remainder<br />

on dividing A by p.<br />

Observation 8 Clearly R[x, y] ≡ R[x][y] also, and the definition <strong>of</strong> g.c.d. is<br />

independent <strong>of</strong> this choice. Algorithmically, though, it seems as if we must make<br />

such a choice. Some systems may already have imposed a default choice, but if<br />

we have a free hand it is usual to choose as the main variable (x in Notation<br />

19) the one which minimises min(deg(A), deg(B)).<br />

4.2.2 The evaluation–interpolation relationship<br />

In this sub-section, we answer a question analogous to that in section 4.1.2:<br />

what do we do if the g.c.d. <strong>of</strong> the evaluations is not the image under evaluation<br />

<strong>of</strong> the g.c.d. calculated be<strong>for</strong>e evaluation?<br />

Lemma 7 If y−v does not divide the leading coefficient <strong>of</strong> gcd(A, B), the degree<br />

<strong>of</strong> gcd(A y−v , B y−v ) is greater than or equal to that <strong>of</strong> gcd(A, B).<br />

Pro<strong>of</strong>. Since gcd(A, B) divides A, then (gcd(A, B)) y−v divides A y−v . Similarly,<br />

it divides B y−v , and there<strong>for</strong>e it divides gcd(A y−v , B y−v ). This implies that the<br />

degree <strong>of</strong> gcd(A y−v , B y−v ) is greater than or equal to that <strong>of</strong> gcd(A, B) y−v .<br />

But the degree <strong>of</strong> gcd(A, B) y−v is equal to that <strong>of</strong> gcd(A, B), <strong>for</strong> the leading<br />

coefficient <strong>of</strong> gcd(A, B) does not cancel when it is evaluated at v.<br />

This lemma is not very easy to use on its own, <strong>for</strong> it supposes that we know<br />

the g.c.d. (or at least its leading coefficient) be<strong>for</strong>e we are able to check whether<br />

the modular reduction has the same degree. But this leading coefficient has<br />

to divide the two leading coefficients <strong>of</strong> A and B, and this gives a <strong>for</strong>mulation<br />

which is easier to use.<br />

Corollary 11 If y − v does not divide the leading coefficients <strong>of</strong> A and <strong>of</strong> B<br />

(it may divide one, but not both), then the degree <strong>of</strong> gcd(A y−v , B y−v ) is greater<br />

than or equal to that <strong>of</strong> gcd(A, B).<br />

As the g.c.d. is the only polynomial (to within a multiple from R[y]) <strong>of</strong> its degree<br />

(in x) which divides A and B, we can test the correctness <strong>of</strong> our calculations <strong>of</strong><br />

the g.c.d.: if the result has the degree <strong>of</strong> gcd(A y−v , B y−v ) (where v satisfies the<br />

hypothesis <strong>of</strong> this corollary) and if it divides A and B, then it is the g.c.d. (to<br />

within a multiple from R[y]).<br />

As in section 4.1.2, it is quite possible that we could find a gcd(A y−v , B y−v )<br />

<strong>of</strong> too high a degree: consider A = x − 1, B = x − y and the evaluation y ↦→ 1.<br />

The following lemma shows that this possibility can only arise <strong>for</strong> a finite number<br />

<strong>of</strong> v.<br />

Lemma 8 Let C = gcd(A, B). If v satisfies the condition <strong>of</strong> the corollary above,<br />

and if y − v does not divide Res x (A/C, B/C), then gcd(A y−v , B y−v ) = C y−v .<br />

Pro<strong>of</strong>. A/C and B/C are relatively prime, <strong>for</strong> otherwise C would not be the<br />

g.c.d. <strong>of</strong> A and B. By the corollary, C y−v does not vanish. There<strong>for</strong>e<br />

gcd(A y−v , B y−v ) = C y−v gcd(A y−v /C y−v , B y−v /C y−v ).

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!