7998
|
1 |
(*
|
|
2 |
Long division of polynomials
|
|
3 |
$Id$
|
|
4 |
Author: Clemens Ballarin, started 23 June 1999
|
|
5 |
*)
|
|
6 |
|
13735
|
7 |
(* legacy bindings and theorems *)
|
|
8 |
|
|
9 |
val deg_aboveI = thm "deg_aboveI";
|
|
10 |
val smult_l_minus = thm "smult_l_minus";
|
|
11 |
val deg_monom_ring = thm "deg_monom_ring";
|
|
12 |
val deg_smult_ring = thm "deg_smult_ring";
|
|
13 |
val smult_l_distr = thm "smult_l_distr";
|
|
14 |
val smult_r_distr = thm "smult_r_distr";
|
|
15 |
val smult_r_minus = thm "smult_r_minus";
|
|
16 |
val smult_assoc2 = thm "smult_assoc2";
|
|
17 |
val smult_assoc1 = thm "smult_assoc1";
|
|
18 |
val monom_mult_smult = thm "monom_mult_smult";
|
|
19 |
val field_ax = thm "field_ax";
|
|
20 |
val lcoeff_nonzero = thm "lcoeff_nonzero";
|
|
21 |
|
|
22 |
Goal
|
|
23 |
"!! f::(nat=>'a::ring). \
|
|
24 |
\ (ALL i. i < m --> f i = 0) --> \
|
|
25 |
\ setsum (%i. f (i+m)) {..d} = setsum f {..m+d}";
|
|
26 |
by (induct_tac "d" 1);
|
|
27 |
(* Base case *)
|
|
28 |
by (induct_tac "m" 1);
|
|
29 |
by (Simp_tac 1);
|
|
30 |
by (Force_tac 1);
|
|
31 |
(* Induction step *)
|
|
32 |
by (asm_simp_tac (simpset() addsimps add_ac) 1);
|
|
33 |
val SUM_shrink_below_lemma = result();
|
|
34 |
|
|
35 |
Goal
|
|
36 |
"!! f::(nat=>'a::ring). \
|
|
37 |
\ [| m <= n; !!i. i < m ==> f i = 0; P (setsum (%i. f (i+m)) {..n-m}) |] \
|
|
38 |
\ ==> P (setsum f {..n})";
|
|
39 |
by (asm_full_simp_tac
|
|
40 |
(simpset() addsimps [SUM_shrink_below_lemma, add_diff_inverse, leD]) 1);
|
|
41 |
qed "SUM_extend_below";
|
|
42 |
|
|
43 |
Goal
|
|
44 |
"!! p::'a::ring up. \
|
|
45 |
\ [| deg p <= n; P (setsum (%i. monom (coeff p i) i) {..n}) |] \
|
|
46 |
\ ==> P p";
|
|
47 |
by (asm_full_simp_tac (simpset() addsimps [thm "up_repr_le"]) 1);
|
|
48 |
qed "up_repr2D";
|
|
49 |
|
|
50 |
(* Start of LongDiv *)
|
|
51 |
|
7998
|
52 |
Goal
|
|
53 |
"!!p::('a::ring up). \
|
|
54 |
\ [| deg p <= deg r; deg q <= deg r; \
|
13735
|
55 |
\ coeff p (deg r) = - (coeff q (deg r)); deg r ~= 0 |] ==> \
|
7998
|
56 |
\ deg (p + q) < deg r";
|
|
57 |
by (res_inst_tac [("j", "deg r - 1")] le_less_trans 1);
|
10962
|
58 |
by (arith_tac 2);
|
7998
|
59 |
by (rtac deg_aboveI 1);
|
|
60 |
by (strip_tac 1);
|
|
61 |
by (case_tac "deg r = m" 1);
|
|
62 |
by (Clarify_tac 1);
|
|
63 |
by (Asm_full_simp_tac 1);
|
|
64 |
(* case "deg q ~= m" *)
|
10959
|
65 |
by (subgoal_tac "deg p < m & deg q < m" 1);
|
7998
|
66 |
by (asm_simp_tac (simpset() addsimps [deg_aboveD]) 1);
|
10959
|
67 |
by (arith_tac 1);
|
7998
|
68 |
qed "deg_lcoeff_cancel";
|
|
69 |
|
|
70 |
Goal
|
|
71 |
"!!p::('a::ring up). \
|
|
72 |
\ [| deg p <= deg r; deg q <= deg r; \
|
13735
|
73 |
\ p ~= -q; coeff p (deg r) = - (coeff q (deg r)) |] ==> \
|
7998
|
74 |
\ deg (p + q) < deg r";
|
|
75 |
by (rtac deg_lcoeff_cancel 1);
|
|
76 |
by (REPEAT (atac 1));
|
|
77 |
by (rtac classical 1);
|
|
78 |
by (Clarify_tac 1);
|
|
79 |
by (etac notE 1);
|
|
80 |
by (res_inst_tac [("p", "p")] up_repr2D 1 THEN atac 1);
|
|
81 |
by (res_inst_tac [("p", "q")] up_repr2D 1 THEN atac 1);
|
|
82 |
by (rotate_tac ~1 1);
|
|
83 |
by (asm_full_simp_tac (simpset() addsimps [smult_l_minus]) 1);
|
|
84 |
qed "deg_lcoeff_cancel2";
|
|
85 |
|
|
86 |
Goal
|
11093
|
87 |
"!!g::('a::ring up). g ~= 0 ==> \
|
7998
|
88 |
\ Ex (% (q, r, k). \
|
|
89 |
\ (lcoeff g)^k *s f = q * g + r & (eucl_size r < eucl_size g))";
|
|
90 |
by (res_inst_tac [("P", "%f. Ex (% (q, r, k). \
|
|
91 |
\ (lcoeff g)^k *s f = q * g + r & (eucl_size r < eucl_size g))")]
|
|
92 |
wf_induct 1);
|
|
93 |
(* TO DO: replace by measure_induct *)
|
|
94 |
by (res_inst_tac [("f", "eucl_size")] wf_measure 1);
|
|
95 |
by (case_tac "eucl_size x < eucl_size g" 1);
|
11093
|
96 |
by (res_inst_tac [("x", "(0, x, 0)")] exI 1);
|
7998
|
97 |
by (Asm_simp_tac 1);
|
|
98 |
(* case "eucl_size x >= eucl_size g" *)
|
13735
|
99 |
by (dres_inst_tac [("x", "lcoeff g *s x - (monom (lcoeff x) (deg x - deg g)) * g")] spec 1);
|
7998
|
100 |
by (etac impE 1);
|
|
101 |
by (full_simp_tac (simpset() addsimps [inv_image_def, measure_def, lcoeff_def]) 1);
|
11093
|
102 |
by (case_tac "x = 0" 1);
|
7998
|
103 |
by (rotate_tac ~1 1);
|
|
104 |
by (asm_full_simp_tac (simpset() addsimps [eucl_size_def]) 1);
|
11093
|
105 |
(* case "x ~= 0 *)
|
7998
|
106 |
by (rotate_tac ~1 1);
|
|
107 |
by (asm_full_simp_tac (simpset() addsimps [eucl_size_def]) 1);
|
13783
|
108 |
(* by (Simp_tac 1); *)
|
7998
|
109 |
by (rtac impI 1);
|
|
110 |
by (rtac deg_lcoeff_cancel2 1);
|
|
111 |
(* replace by linear arithmetic??? *)
|
13783
|
112 |
by (rtac le_trans 2);
|
|
113 |
by (rtac deg_smult_ring 2);
|
|
114 |
by (Simp_tac 2);
|
7998
|
115 |
by (Simp_tac 1);
|
|
116 |
by (rtac le_trans 1);
|
|
117 |
by (rtac deg_mult_ring 1);
|
|
118 |
by (rtac le_trans 1);
|
13735
|
119 |
(**)
|
|
120 |
by (rtac add_le_mono 1); by (rtac le_refl 1);
|
|
121 |
(* term order forces to use this instead of add_le_mono1 *)
|
|
122 |
by (rtac deg_monom_ring 1);
|
|
123 |
by (Asm_simp_tac 1);
|
|
124 |
(**)
|
|
125 |
(*
|
7998
|
126 |
by (rtac add_le_mono1 1);
|
|
127 |
by (rtac deg_smult_ring 1);
|
11093
|
128 |
(* by (asm_simp_tac (simpset() addsimps [leI]) 1); *)
|
|
129 |
by (Asm_simp_tac 1);
|
|
130 |
by (cut_inst_tac [("m", "deg x - deg g"), ("'a", "'a")] deg_monom_ring 1);
|
|
131 |
by (arith_tac 1);
|
13735
|
132 |
*)
|
8707
|
133 |
by (Force_tac 1);
|
7998
|
134 |
by (Simp_tac 1);
|
13735
|
135 |
(**)
|
|
136 |
(* This change is probably caused by application of commutativity *)
|
|
137 |
by (res_inst_tac [("m", "deg g"), ("n", "deg x")] SUM_extend 1);
|
|
138 |
by (Simp_tac 1);
|
|
139 |
by (Asm_simp_tac 1);
|
|
140 |
by (arith_tac 1);
|
|
141 |
by (res_inst_tac [("m", "deg g"), ("n", "deg g")] SUM_extend_below 1);
|
|
142 |
by (rtac le_refl 1);
|
|
143 |
by (Asm_simp_tac 1);
|
|
144 |
by (arith_tac 1);
|
|
145 |
by (Simp_tac 1);
|
|
146 |
(**)
|
|
147 |
(*
|
7998
|
148 |
by (res_inst_tac [("m", "deg x - deg g"), ("n", "deg x")] SUM_extend 1);
|
|
149 |
by (Simp_tac 1);
|
|
150 |
by (asm_simp_tac (simpset() addsimps [less_not_refl2 RS not_sym]) 1);
|
|
151 |
by (res_inst_tac [("m", "deg x - deg g"), ("n", "deg x - deg g")]
|
|
152 |
SUM_extend_below 1);
|
|
153 |
by (rtac le_refl 1);
|
|
154 |
by (asm_simp_tac (simpset() addsimps [less_not_refl2]) 1);
|
|
155 |
by (asm_simp_tac (simpset() addsimps [diff_diff_right, leI, m_comm]) 1);
|
13735
|
156 |
*)
|
7998
|
157 |
(* end of subproof deg f1 < deg f *)
|
|
158 |
by (etac exE 1);
|
13735
|
159 |
by (res_inst_tac [("x", "((%(q,r,k). (monom (lcoeff g ^ k * lcoeff x) (deg x - deg g) + q)) xa, (%(q,r,k). r) xa, (%(q,r,k). Suc k) xa)")] exI 1);
|
7998
|
160 |
by (Clarify_tac 1);
|
|
161 |
by (dtac sym 1);
|
13735
|
162 |
by (simp_tac (simpset() addsimps [l_distr, a_assoc]
|
|
163 |
delsimprocs [ring_simproc]) 1);
|
|
164 |
by (asm_simp_tac (simpset() delsimprocs [ring_simproc]) 1);
|
|
165 |
by (simp_tac (simpset() addsimps [minus_def, smult_r_distr, smult_r_minus,
|
|
166 |
monom_mult_smult, smult_assoc1, smult_assoc2]
|
|
167 |
delsimprocs [ring_simproc]) 1);
|
|
168 |
by (Simp_tac 1);
|
7998
|
169 |
qed "long_div_eucl_size";
|
|
170 |
|
|
171 |
Goal
|
11093
|
172 |
"!!g::('a::ring up). g ~= 0 ==> \
|
7998
|
173 |
\ Ex (% (q, r, k). \
|
11093
|
174 |
\ (lcoeff g)^k *s f = q * g + r & (r = 0 | deg r < deg g))";
|
7998
|
175 |
by (forw_inst_tac [("f", "f")]
|
13735
|
176 |
(simplify (simpset() addsimps [eucl_size_def]
|
|
177 |
delsimprocs [ring_simproc]) long_div_eucl_size) 1);
|
|
178 |
by (auto_tac (claset(), simpset() delsimprocs [ring_simproc]));
|
11093
|
179 |
by (case_tac "aa = 0" 1);
|
7998
|
180 |
by (Blast_tac 1);
|
11093
|
181 |
(* case "aa ~= 0 *)
|
7998
|
182 |
by (rotate_tac ~1 1);
|
|
183 |
by Auto_tac;
|
|
184 |
qed "long_div_ring";
|
|
185 |
|
13735
|
186 |
(* Next one fails *)
|
7998
|
187 |
Goal
|
13735
|
188 |
"!!g::('a::ring up). [| g ~= 0; (lcoeff g) dvd 1 |] ==> \
|
11093
|
189 |
\ Ex (% (q, r). f = q * g + r & (r = 0 | deg r < deg g))";
|
7998
|
190 |
by (forw_inst_tac [("f", "f")] long_div_ring 1);
|
|
191 |
by (etac exE 1);
|
10448
|
192 |
by (res_inst_tac [("x", "((%(q,r,k). (inverse(lcoeff g ^k) *s q)) x, \
|
|
193 |
\ (%(q,r,k). inverse(lcoeff g ^k) *s r) x)")] exI 1);
|
7998
|
194 |
by (Clarify_tac 1);
|
13735
|
195 |
(* by (Simp_tac 1); *)
|
7998
|
196 |
by (rtac conjI 1);
|
|
197 |
by (dtac sym 1);
|
13735
|
198 |
by (asm_simp_tac (simpset() addsimps [smult_r_distr RS sym, smult_assoc2]
|
|
199 |
delsimprocs [ring_simproc]) 1);
|
7998
|
200 |
by (asm_simp_tac (simpset() addsimps [l_inverse_ring, unit_power, smult_assoc1 RS sym]) 1);
|
|
201 |
(* degree property *)
|
|
202 |
by (etac disjE 1);
|
|
203 |
by (Asm_simp_tac 1);
|
|
204 |
by (rtac disjI2 1);
|
|
205 |
by (rtac le_less_trans 1);
|
|
206 |
by (rtac deg_smult_ring 1);
|
8006
|
207 |
by (Asm_simp_tac 1);
|
7998
|
208 |
qed "long_div_unit";
|
|
209 |
|
|
210 |
Goal
|
11093
|
211 |
"!!g::('a::field up). g ~= 0 ==> \
|
|
212 |
\ Ex (% (q, r). f = q * g + r & (r = 0 | deg r < deg g))";
|
7998
|
213 |
by (rtac long_div_unit 1);
|
|
214 |
by (assume_tac 1);
|
13735
|
215 |
by (asm_simp_tac (simpset() addsimps [lcoeff_def, lcoeff_nonzero, field_ax]) 1);
|
7998
|
216 |
qed "long_div_theorem";
|
|
217 |
|
13735
|
218 |
Goal "- (0::'a::ring) = 0";
|
|
219 |
by (Simp_tac 1);
|
|
220 |
val uminus_zero = result();
|
|
221 |
|
|
222 |
Goal "!!a::'a::ring. a - b = 0 ==> a = b";
|
13783
|
223 |
by (res_inst_tac [("s", "a - (a - b)")] trans 1);
|
|
224 |
by (asm_simp_tac (simpset() delsimprocs [ring_simproc]) 1);
|
|
225 |
by (Simp_tac 1);
|
|
226 |
by (Simp_tac 1);
|
13735
|
227 |
val diff_zero_imp_eq = result();
|
|
228 |
|
|
229 |
Goal "!!a::'a::ring. a = b ==> a + (-b) = 0";
|
|
230 |
by (Asm_simp_tac 1);
|
|
231 |
val eq_imp_diff_zero = result();
|
|
232 |
|
7998
|
233 |
Goal
|
11093
|
234 |
"!!g::('a::field up). [| g ~= 0; \
|
|
235 |
\ f = q1 * g + r1; (r1 = 0 | deg r1 < deg g); \
|
|
236 |
\ f = q2 * g + r2; (r2 = 0 | deg r2 < deg g) |] ==> q1 = q2";
|
13735
|
237 |
by (subgoal_tac "(q1 - q2) * g = r2 - r1" 1); (* 1 *)
|
7998
|
238 |
by (thin_tac "f = ?x" 1);
|
|
239 |
by (thin_tac "f = ?x" 1);
|
|
240 |
by (rtac diff_zero_imp_eq 1);
|
|
241 |
by (rtac classical 1);
|
|
242 |
by (etac disjE 1);
|
11093
|
243 |
(* r1 = 0 *)
|
7998
|
244 |
by (etac disjE 1);
|
11093
|
245 |
(* r2 = 0 *)
|
13735
|
246 |
by (asm_full_simp_tac (simpset()
|
|
247 |
addsimps [thm "integral_iff", minus_def, l_zero, uminus_zero]
|
|
248 |
delsimprocs [ring_simproc]) 1);
|
11093
|
249 |
(* r2 ~= 0 *)
|
13735
|
250 |
by (dres_inst_tac [("f", "deg"), ("y", "r2 - r1")] arg_cong 1);
|
|
251 |
by (asm_full_simp_tac (simpset() addsimps [minus_def, l_zero, uminus_zero]
|
|
252 |
delsimprocs [ring_simproc]) 1);
|
11093
|
253 |
(* r1 ~=0 *)
|
7998
|
254 |
by (etac disjE 1);
|
11093
|
255 |
(* r2 = 0 *)
|
13735
|
256 |
by (dres_inst_tac [("f", "deg"), ("y", "r2 - r1")] arg_cong 1);
|
|
257 |
by (asm_full_simp_tac (simpset() addsimps [minus_def, l_zero, uminus_zero]
|
|
258 |
delsimprocs [ring_simproc]) 1);
|
11093
|
259 |
(* r2 ~= 0 *)
|
13735
|
260 |
by (dres_inst_tac [("f", "deg"), ("y", "r2 - r1")] arg_cong 1);
|
|
261 |
by (asm_full_simp_tac (simpset() addsimps [minus_def]
|
|
262 |
delsimprocs [ring_simproc]) 1);
|
7998
|
263 |
by (dtac (order_eq_refl RS add_leD2) 1);
|
|
264 |
by (dtac leD 1);
|
|
265 |
by (etac notE 1 THEN rtac (deg_add RS le_less_trans) 1);
|
|
266 |
by (Asm_simp_tac 1);
|
|
267 |
(* proof of 1 *)
|
|
268 |
by (rtac diff_zero_imp_eq 1);
|
13601
|
269 |
by (hyp_subst_tac 1);
|
7998
|
270 |
by (dres_inst_tac [("a", "?x+?y")] eq_imp_diff_zero 1);
|
13735
|
271 |
by (Asm_full_simp_tac 1);
|
7998
|
272 |
qed "long_div_quo_unique";
|
|
273 |
|
|
274 |
Goal
|
11093
|
275 |
"!!g::('a::field up). [| g ~= 0; \
|
|
276 |
\ f = q1 * g + r1; (r1 = 0 | deg r1 < deg g); \
|
|
277 |
\ f = q2 * g + r2; (r2 = 0 | deg r2 < deg g) |] ==> r1 = r2";
|
7998
|
278 |
by (subgoal_tac "q1 = q2" 1);
|
|
279 |
by (Clarify_tac 1);
|
13735
|
280 |
by (res_inst_tac [("a", "q2 * g + r1 - q2 * g"), ("b", "q2 * g + r2 - q2 * g")]
|
|
281 |
box_equals 1);
|
|
282 |
by (Asm_full_simp_tac 1);
|
|
283 |
by (Simp_tac 1);
|
|
284 |
by (Simp_tac 1);
|
7998
|
285 |
by (rtac long_div_quo_unique 1);
|
|
286 |
by (REPEAT (atac 1));
|
|
287 |
qed "long_div_rem_unique";
|