next up previous
Next: Modular Computation Up: Interpolation and Rational Reconstruction Previous: The Chinese Remaindering Algorithm

Rational Function Reconstruction

Let $ \bf k$ mathend000# be a field, let m $ \in$ $ \bf k$[x] mathend000# be a polynomial of degree n > 0 mathend000#. Given a polynomial f $ \in$ $ \bf k$[x] mathend000# of degree less than n mathend000# and an integer d $ \in$ {1,..., n} mathend000#, we want to find a rational function r/t $ \in$ $ \bf k$(x) mathend000# with r, t $ \in$ $ \bf k$[x] mathend000# satisfying

gcd(t, m) = 1  and  rt-1 $\displaystyle \equiv$ f modm,  degr < d,  degt $\displaystyle \leq$ n - d. (86)

Let us denote this problem by RFR(m, n, f, d ) mathend000#. A solution to this problem is also a solution to the following weaker problem: compute a rational function r/t $ \in$ $ \bf k$(x) mathend000# with r, t $ \in$ $ \bf k$[x] mathend000# satisfying

r $\displaystyle \equiv$ tf modm,  degr < d,  degt $\displaystyle \leq$ n - d. (87)

We denote this second problem by WRFR(m, n, f, d ) mathend000#.

We recall that a rational function r/t $ \in$ $ \bf k$(x) mathend000# is said to be in canonical form if t mathend000# is monic and if gcd(r, t) = 1 mathend000#. Clearly, every rational function has a unique canonical form.

The following Theorem 5:

As a result, we obtain an algorithm for solving both problems WRFR(m, n, f, d ) mathend000# and RFR(m, n, f, d ) mathend000#.

Theorem 5   Let rj, sj, tj $ \in$ $ \bf k$[x] mathend000# be the j mathend000#-th row of the Extended Euclidean Algorithm applied to (a, b) = (m, f ) mathend000# where j mathend000# is minimal such that degrj < d mathend000#.
(i) mathend000#
The couple (r, t) = (rj, tj) mathend000# is a solution to problem WRFR(m, n, f, d ) mathend000#.
(ii) mathend000#
If gcd(rj, tj) = 1 mathend000#, then the couple (r, t) = (rj, tj) mathend000# is a solution to problem RFR(m, n, f, d ) mathend000#.
(iii) mathend000#
If (r, t) mathend000# is a solution to problem RFR(m, n, f, d ) mathend000# and if r/t mathend000# is a canonical form, then we have (r, t) = (wj-1rj, wj-1tj) mathend000# where wj = lc(tj) mathend000#.
(iv) mathend000#
Problem RFR(m, n, f, d ) mathend000# admits a solution if and only if gcd(rj, tj) = 1 mathend000#.

Proof. First, we observe that rj, sj, tj mathend000# are uniquely defined. Indeed, we have 1 $ \leq$ d $ \leq$ n mathend000# and,

- $\displaystyle \infty$ = degrk+1 < degrk < ... < degri+1 < degri < ... < dega = n. (88)

Hence, there exists a unique j $ \in$ {1,..., k + 1} mathend000# such that we have

degrj < d $\displaystyle \leq$ degrj-1. (89)

From Point (iv) mathend000# of Proposition 1 we have

rj = sjm + tjf (90)

leading to rj $ \equiv$ tjf modm mathend000#. From Proposition 2, we deduce

degtj = degm - degrj-1 $\displaystyle \leq$ n - d. (91)

Since degrj < d mathend000# holds, we have proved that (rj, tj) mathend000# is a solution to problem WRFR(m, n, f, d ) mathend000#. This shows (i) mathend000#. Now, we prove (ii) mathend000#. We assume that gcd(rj, tj) = 1 mathend000# holds. From Point (vii) mathend000# of Proposition 1 we have

gcd(m, tj) = gcd(rj, tj) = 1. (92)

This implies that tj mathend000# is invertible modulo m mathend000#, leading to rt-1 $ \equiv$ f modm mathend000#. Hence, under the assumption gcd(rj, tj) = 1 mathend000#, the couple (rj, tj) mathend000# is also a solution to RFR(m, n, f, d ) mathend000#.

Next, we prove (iii) mathend000#. So we consider (r, t) mathend000# a solution to RFR(m, n, f, d ) mathend000# such that the fraction r/t mathend000# is in canonical form. There exists s $ \in$ $ \bf k$[x] mathend000# such that we have r = sm + tf mathend000#. Since degr + degt < degm mathend000# holds, we can apply Proposition 3. So, let $ \ell$ $ \in$ {1,..., k + 1} mathend000# be such that

degr$\scriptstyle \ell$ $\displaystyle \leq$ degr < degr$\scriptstyle \ell$-1. (93)

Then, there exists $ \alpha$ $ \in$ $ \bf k$[x] mathend000# such that r = $ \alpha$r$\scriptstyle \ell$ mathend000# and t = $ \alpha$t$\scriptstyle \ell$ mathend000# both hold.

If degrj $ \leq$ degr mathend000# then we have $ \ell$ = j mathend000#. If degr < degrj mathend000# then we have $ \ell$ > j mathend000#.

We prove that $ \ell$ = j mathend000# holds. If t = tj mathend000# then r - rj $ \equiv$ 0 modm mathend000#, and thus r = rj mathend000#, since deg(r - rj) < n mathend000# holds, which implies $ \ell$ = j mathend000#. If t $ \neq$ tj mathend000#, we can apply Proposition 3 again. Then, there exists $ \beta$ $ \in$ $ \bf k$[x] mathend000# such that r - rj = $ \beta$r$\scriptstyle \ell$ mathend000# and t - tj = $ \beta$t$\scriptstyle \ell$ mathend000# both hold. It follows that we have rj = ($ \alpha$ - $ \beta$)r$\scriptstyle \ell$ mathend000# and tj = ($ \alpha$ - $ \beta$)t$\scriptstyle \ell$ mathend000#. Hence, we deduce

msj = rj - tjf = ($\displaystyle \alpha$ - $\displaystyle \beta$)ms$\scriptstyle \ell$ (94)

leading to sj = ($ \alpha$ - $ \beta$)s$\scriptstyle \ell$ mathend000#. Since gcd(sj, tj) = 1 mathend000# holds we must have $ \alpha$ - $ \beta$ $ \in$ $ \bf k$ mathend000# and thus rj = r$\scriptstyle \ell$ mathend000#. Therefore, we have proved that $ \ell$ = j mathend000# holds. This implies

r = $\displaystyle \alpha$rj  and  t = $\displaystyle \alpha$tj (95)

Since r/t mathend000# is in canonical form, we deduce that gcd(tj, rj) = 1 mathend000# and $ \alpha$ = (lc(tj))-1 mathend000# hold. This proves (iii) mathend000#. Finally, (iv) mathend000# follows from (ii) mathend000# and (iii) mathend000#. $ \qedsymbol$


next up previous
Next: Modular Computation Up: Interpolation and Rational Reconstruction Previous: The Chinese Remaindering Algorithm
Marc Moreno Maza
2007-01-10