3 Author: Jeremy Dawson, NICTA
6 header {* Useful Numerical Lemmas *}
12 lemma contentsI: "y = {x} ==> contents y = x"
13 unfolding contents_def by auto
15 lemma prod_case_split: "prod_case = split"
18 lemmas split_split = prod.split [unfolded prod_case_split]
19 lemmas split_split_asm = prod.split_asm [unfolded prod_case_split]
20 lemmas "split.splits" = split_split split_split_asm
22 lemmas funpow_0 = funpow.simps(1)
23 lemmas funpow_Suc = funpow.simps(2)
25 lemma nonemptyE: "S ~= {} ==> (!!x. x : S ==> R) ==> R"
26 apply (erule contrapos_np)
31 lemma gt_or_eq_0: "0 < y \<or> 0 = (y::nat)" by auto
34 mod_alt :: "'a => 'a => 'a :: Divides.div"
35 "mod_alt n m == n mod m"
37 -- "alternative way of defining @{text bin_last}, @{text bin_rest}"
38 bin_rl :: "int => int * bit"
39 "bin_rl w == SOME (r, l). w = r BIT l"
41 declare iszero_0 [iff]
43 lemmas xtr1 = xtrans(1)
44 lemmas xtr2 = xtrans(2)
45 lemmas xtr3 = xtrans(3)
46 lemmas xtr4 = xtrans(4)
47 lemmas xtr5 = xtrans(5)
48 lemmas xtr6 = xtrans(6)
49 lemmas xtr7 = xtrans(7)
50 lemmas xtr8 = xtrans(8)
52 lemma Min_ne_Pls [iff]:
54 unfolding Min_def Pls_def by auto
56 lemmas Pls_ne_Min [iff] = Min_ne_Pls [symmetric]
58 lemmas PlsMin_defs [intro!] =
59 Pls_def Min_def Pls_def [symmetric] Min_def [symmetric]
61 lemmas PlsMin_simps [simp] = PlsMin_defs [THEN Eq_TrueI]
63 lemma number_of_False_cong:
64 "False ==> number_of x = number_of y"
67 lemmas nat_simps = diff_add_inverse2 diff_add_inverse
68 lemmas nat_iffs = le_add1 le_add2
70 lemma sum_imp_diff: "j = k + i ==> j - i = (k :: nat)"
71 by (clarsimp simp add: nat_simps)
74 "0 < (number_of w :: nat) ==>
75 number_of w - (1 :: nat) = number_of (Int.pred w)"
76 apply (unfold nat_number_of_def One_nat_def nat_1 [symmetric] pred_def)
77 apply (simp add: number_of_eq nat_diff_distrib [symmetric])
81 "of_int (a ^ n) = (of_int a ^ n :: 'a :: {recpower, comm_ring_1})"
82 by (induct n) (auto simp add: power_Suc)
84 lemma zless2: "0 < (2 :: int)"
87 lemmas zless2p [simp] = zless2 [THEN zero_less_power]
88 lemmas zle2p [simp] = zless2p [THEN order_less_imp_le]
90 lemmas pos_mod_sign2 = zless2 [THEN pos_mod_sign [where b = "2::int"]]
91 lemmas pos_mod_bound2 = zless2 [THEN pos_mod_bound [where b = "2::int"]]
93 -- "the inverse(s) of @{text number_of}"
94 lemma nmod2: "n mod (2::int) = 0 | n mod 2 = 1"
95 using pos_mod_sign2 [of n] pos_mod_bound2 [of n]
96 unfolding mod_alt_def [symmetric] by auto
99 "even n ==> even d ==> 0 <= d ==> (n + 1) mod (d :: int) = (n mod d) + 1"
100 apply (simp add: add_commute)
101 apply (safe dest!: even_equiv_def [THEN iffD1])
102 apply (subst pos_zmod_mult_2)
104 apply (simp add: zmod_zmult_zmult1)
107 lemmas eme1p = emep1 [simplified add_commute]
109 lemma le_diff_eq': "(a \<le> c - b) = (b + a \<le> (c::int))"
110 by (simp add: le_diff_eq add_commute)
112 lemma less_diff_eq': "(a < c - b) = (b + a < (c::int))"
113 by (simp add: less_diff_eq add_commute)
115 lemma diff_le_eq': "(a - b \<le> c) = (a \<le> b + (c::int))"
116 by (simp add: diff_le_eq add_commute)
118 lemma diff_less_eq': "(a - b < c) = (a < b + (c::int))"
119 by (simp add: diff_less_eq add_commute)
122 lemmas m1mod2k = zless2p [THEN zmod_minus1]
123 lemmas m1mod22k = mult_pos_pos [OF zless2 zless2p, THEN zmod_minus1]
124 lemmas p1mod22k' = zless2p [THEN order_less_imp_le, THEN pos_zmod_mult_2]
125 lemmas z1pmod2' = zero_le_one [THEN pos_zmod_mult_2, simplified]
126 lemmas z1pdiv2' = zero_le_one [THEN pos_zdiv_mult_2, simplified]
129 "(2 * b + 1) mod (2 * 2 ^ n) = 2 * (b mod 2 ^ n) + (1::int)"
130 by (simp add: p1mod22k' add_commute)
133 "(2 * b + 1) mod 2 = (1::int)"
134 by (simp add: z1pmod2' add_commute)
137 "(2 * b + 1) div 2 = (b::int)"
138 by (simp add: z1pdiv2' add_commute)
140 lemmas zdiv_le_dividend = xtr3 [OF zdiv_1 [symmetric] zdiv_mono2,
141 simplified int_one_le_iff_zero_less, simplified, standard]
143 (** ways in which type Bin resembles a datatype **)
145 lemma BIT_eq: "u BIT b = v BIT c ==> u = v & b = c"
146 apply (unfold Bit_def)
147 apply (simp (no_asm_use) split: bit.split_asm)
149 apply (drule_tac f=even in arg_cong, clarsimp)+
152 lemmas BIT_eqE [elim!] = BIT_eq [THEN conjE, standard]
154 lemma BIT_eq_iff [simp]:
155 "(u BIT b = v BIT c) = (u = v \<and> b = c)"
158 lemmas BIT_eqI [intro!] = conjI [THEN BIT_eq_iff [THEN iffD2]]
161 "(v BIT b < w BIT c) = (v < w | v <= w & b = bit.B0 & c = bit.B1)"
162 unfolding Bit_def by (auto split: bit.split)
165 "(v BIT b <= w BIT c) = (v < w | v <= w & (b ~= bit.B1 | c ~= bit.B0))"
166 unfolding Bit_def by (auto split: bit.split)
169 assumes ne: "y \<noteq> bit.B1"
170 assumes y: "y = bit.B0 \<Longrightarrow> P"
173 apply (cases y rule: bit.exhaust, simp)
177 lemma bin_ex_rl: "EX w b. w BIT b = bin"
178 apply (unfold Bit_def)
179 apply (cases "even bin")
180 apply (clarsimp simp: even_equiv_def)
181 apply (auto simp: odd_equiv_def split: bit.split)
185 assumes Q: "\<And>x b. bin = x BIT b \<Longrightarrow> Q"
187 apply (insert bin_ex_rl [of bin])
193 lemma bin_rl_char: "(bin_rl w = (r, l)) = (r BIT l = w)"
194 apply (unfold bin_rl_def)
196 apply (cases w rule: bin_exhaust)
200 lemmas bin_rl_simps [THEN bin_rl_char [THEN iffD2], standard, simp] =
201 Pls_0_eq Min_1_eq refl
204 "bin = (w BIT b) ==> ~ bin = Int.Min --> ~ bin = Int.Pls -->
205 nat (abs w) < nat (abs bin)"
206 apply (clarsimp simp add: bin_rl_char)
207 apply (unfold Pls_def Min_def Bit_def)
209 apply (clarsimp, arith)
210 apply (clarsimp, arith)
214 assumes PPls: "P Int.Pls"
215 and PMin: "P Int.Min"
216 and PBit: "!!bin bit. P bin ==> P (bin BIT bit)"
218 apply (rule_tac P=P and a=bin and f1="nat o abs"
219 in wf_measure [THEN wf_induct])
220 apply (simp add: measure_def inv_image_def)
221 apply (case_tac x rule: bin_exhaust)
222 apply (frule bin_abs_lem)
223 apply (auto simp add : PPls PMin PBit)
226 lemma no_no [simp]: "number_of (number_of i) = i"
227 unfolding number_of_eq by simp
230 "k BIT bit.B0 = k + k"
231 by (unfold Bit_def) simp
234 "k BIT bit.B1 = k + k + 1"
235 by (unfold Bit_def) simp
237 lemma Bit_B0_2t: "k BIT bit.B0 = 2 * k"
238 by (rule trans, rule Bit_B0) simp
240 lemma Bit_B1_2t: "k BIT bit.B1 = 2 * k + 1"
241 by (rule trans, rule Bit_B1) simp
244 "X = 2 ==> (w BIT bit.B1) mod X = 1 & (w BIT bit.B0) mod X = 0"
245 apply (simp (no_asm) only: Bit_B0 Bit_B1)
246 apply (simp add: z1pmod2)
249 lemmas B1_mod_2 [simp] = B_mod_2' [OF refl, THEN conjunct1, standard]
250 lemmas B0_mod_2 [simp] = B_mod_2' [OF refl, THEN conjunct2, standard]
253 "a + m + m = b + n + n ==> (a = 0 | a = 1) ==> (b = 0 | b = 1) ==>
254 a = b & m = (n :: int)"
256 apply (drule_tac f="%n. n mod 2" in arg_cong)
257 apply (clarsimp simp: z1pmod2)
258 apply (drule_tac f="%n. n mod 2" in arg_cong)
259 apply (clarsimp simp: z1pmod2)
263 "(1 + x + x) mod 2 = (1 :: int) & (0 + x + x) mod 2 = (0 :: int)"
264 by simp (rule z1pmod2)
267 "(1 + x + x) div 2 = (x :: int) & (0 + x + x) div 2 = (x :: int)"
268 by simp (rule z1pdiv2)
270 lemmas iszero_minus = trans [THEN trans,
271 OF iszero_def neg_equal_0_iff_equal iszero_def [symmetric], standard]
273 lemmas zadd_diff_inverse = trans [OF diff_add_cancel [symmetric] add_commute,
276 lemmas add_diff_cancel2 = add_commute [THEN diff_eq_eq [THEN iffD2], standard]
278 lemma zmod_uminus: "- ((a :: int) mod b) mod b = -a mod b"
279 by (simp add : zmod_zminus1_eq_if)
281 lemma zmod_zsub_distrib: "((a::int) - b) mod c = (a mod c - b mod c) mod c"
282 apply (unfold diff_int_def)
283 apply (rule trans [OF _ zmod_zadd1_eq [symmetric]])
284 apply (simp add: zmod_uminus zmod_zadd1_eq [symmetric])
287 lemma zmod_zsub_right_eq: "((a::int) - b) mod c = (a - b mod c) mod c"
288 apply (unfold diff_int_def)
289 apply (rule trans [OF _ zmod_zadd_right_eq [symmetric]])
290 apply (simp add : zmod_uminus zmod_zadd_right_eq [symmetric])
293 lemma zmod_zsub_left_eq: "((a::int) - b) mod c = (a mod c - b) mod c"
294 by (rule zmod_zadd_left_eq [where b = "- b", simplified diff_int_def [symmetric]])
296 lemma zmod_zsub_self [simp]:
297 "((b :: int) - a) mod a = b mod a"
298 by (simp add: zmod_zsub_right_eq)
300 lemma zmod_zmult1_eq_rev:
301 "b * a mod c = b mod c * a mod (c::int)"
302 apply (simp add: mult_commute)
303 apply (subst zmod_zmult1_eq)
307 lemmas rdmods [symmetric] = zmod_uminus [symmetric]
308 zmod_zsub_left_eq zmod_zsub_right_eq zmod_zadd_left_eq
309 zmod_zadd_right_eq zmod_zmult1_eq zmod_zmult1_eq_rev
311 lemma mod_plus_right:
312 "((a + x) mod m = (b + x) mod m) = (a mod m = b mod (m :: nat))"
314 apply (simp_all add: mod_Suc)
318 lemma nat_minus_mod: "(n - n mod m) mod m = (0 :: nat)"
319 by (induct n) (simp_all add : mod_Suc)
321 lemmas nat_minus_mod_plus_right = trans [OF nat_minus_mod mod_0 [symmetric],
322 THEN mod_plus_right [THEN iffD2], standard, simplified]
324 lemmas push_mods' = zmod_zadd1_eq [standard]
325 zmod_zmult_distrib [standard] zmod_zsub_distrib [standard]
326 zmod_uminus [symmetric, standard]
328 lemmas push_mods = push_mods' [THEN eq_reflection, standard]
329 lemmas pull_mods = push_mods [symmetric] rdmods [THEN eq_reflection, standard]
331 zmod_zmult_self1 [THEN eq_reflection] zmod_zmult_self2 [THEN eq_reflection]
332 mod_mod_trivial [THEN eq_reflection]
335 "!!b. b < n ==> a mod n = b mod n ==> a mod n = (b :: nat)"
338 lemmas nat_mod_eq' = refl [THEN [2] nat_mod_eq]
341 "(0 :: nat) < n ==> b < n = (b mod n = b)"
343 apply (erule nat_mod_eq')
345 apply (erule mod_less_divisor)
349 "(x :: nat) < z ==> y < z ==>
350 (x + y) mod z = (if x + y < z then x + y else x + y - z)"
351 apply (rule nat_mod_eq)
354 apply (rule le_mod_geq)
356 apply (rule nat_mod_eq')
361 "(x :: nat) < z ==> (x - y) mod z = x - y"
362 by (rule nat_mod_eq') arith
365 "(0 :: int) < n ==> (0 <= b & b < n) = (b mod n = b)"
367 apply (erule (1) mod_pos_pos_trivial)
368 apply (erule_tac [!] subst)
373 "(0 :: int) <= b ==> b < n ==> a mod n = b mod n ==> a mod n = b"
374 by clarsimp (rule mod_pos_pos_trivial)
376 lemmas int_mod_eq' = refl [THEN [3] int_mod_eq]
378 lemma int_mod_le: "0 <= a ==> 0 < (n :: int) ==> a mod n <= a"
379 apply (cases "a < n")
380 apply (auto dest: mod_pos_pos_trivial pos_mod_bound [where a=a])
383 lemma int_mod_le': "0 <= b - n ==> 0 < (n :: int) ==> b mod n <= b - n"
384 by (rule int_mod_le [where a = "b - n" and n = n, simplified])
386 lemma int_mod_ge: "a < n ==> 0 < (n :: int) ==> a <= a mod n"
387 apply (cases "0 <= a")
388 apply (drule (1) mod_pos_pos_trivial)
390 apply (rule order_trans [OF _ pos_mod_sign])
395 lemma int_mod_ge': "b < 0 ==> 0 < (n :: int) ==> b + n <= b mod n"
396 by (rule int_mod_ge [where a = "b + n" and n = n, simplified])
399 "(x :: int) < z ==> y < z ==> 0 <= y ==> 0 <= x ==> 0 <= z ==>
400 (x + y) mod z = (if x + y < z then x + y else x + y - z)"
401 by (auto intro: int_mod_eq)
404 "(x :: int) < z ==> y < z ==> 0 <= y ==> 0 <= x ==> 0 <= z ==>
405 (x - y) mod z = (if y <= x then x - y else x - y + z)"
406 by (auto intro: int_mod_eq)
408 lemmas zmde = zmod_zdiv_equality [THEN diff_eq_eq [THEN iffD2], symmetric]
409 lemmas mcl = mult_cancel_left [THEN iffD1, THEN make_pos_rule]
411 (* already have this for naturals, div_mult_self1/2, but not for ints *)
412 lemma zdiv_mult_self: "m ~= (0 :: int) ==> (a + m * n) div m = a div m + n"
416 apply (simp add: zmde ring_distribs)
417 apply (simp add: push_mods)
421 lemma eqne: "equiv A r ==> X : A // r ==> X ~= {}"
422 unfolding equiv_def refl_def quotient_def Image_def by auto
424 lemmas Rep_Integ_ne = Integ.Rep_Integ
425 [THEN equiv_intrel [THEN eqne, simplified Integ_def [symmetric]], standard]
427 lemmas riq = Integ.Rep_Integ [simplified Integ_def]
428 lemmas intrel_refl = refl [THEN equiv_intrel_iff [THEN iffD1], standard]
429 lemmas Rep_Integ_equiv = quotient_eq_iff
430 [OF equiv_intrel riq riq, simplified Integ.Rep_Integ_inject, standard]
431 lemmas Rep_Integ_same =
432 Rep_Integ_equiv [THEN intrel_refl [THEN rev_iffD2], standard]
434 lemma RI_int: "(a, 0) : Rep_Integ (int a)"
435 unfolding int_def by auto
437 lemmas RI_intrel [simp] = UNIV_I [THEN quotientI,
438 THEN Integ.Abs_Integ_inverse [simplified Integ_def], standard]
440 lemma RI_minus: "(a, b) : Rep_Integ x ==> (b, a) : Rep_Integ (- x)"
441 apply (rule_tac z=x in eq_Abs_Integ)
442 apply (clarsimp simp: minus)
446 "(a, b) : Rep_Integ x ==> (c, d) : Rep_Integ y ==>
447 (a + c, b + d) : Rep_Integ (x + y)"
448 apply (rule_tac z=x in eq_Abs_Integ)
449 apply (rule_tac z=y in eq_Abs_Integ)
450 apply (clarsimp simp: add)
453 lemma mem_same: "a : S ==> a = b ==> b : S"
456 (* two alternative proofs of this *)
457 lemma RI_eq_diff': "(a, b) : Rep_Integ (int a - int b)"
458 apply (unfold diff_def)
459 apply (rule mem_same)
460 apply (rule RI_minus RI_add RI_int)+
464 lemma RI_eq_diff: "((a, b) : Rep_Integ x) = (int a - int b = x)"
466 apply (rule Rep_Integ_same)
469 apply (rule RI_eq_diff')+
473 "a > 1 ==> a ^ n mod a ^ m = (if m <= n then 0 else (a :: int) ^ n)"
476 apply (simp add: zdvd_iff_zmod_eq_0 [symmetric])
477 apply (drule le_iff_add [THEN iffD1])
478 apply (force simp: zpower_zadd_distrib)
479 apply (rule mod_pos_pos_trivial)
481 apply (rule power_strict_increasing)
485 lemma min_pm [simp]: "min a b + (a - b) = (a :: nat)"
488 lemmas min_pm1 [simp] = trans [OF add_commute min_pm]
490 lemma rev_min_pm [simp]: "min b a + (a - b) = (a::nat)"
493 lemmas rev_min_pm1 [simp] = trans [OF add_commute rev_min_pm]
497 a >= c & b <= d | a <= c & b >= (d :: nat)"
498 apply (cut_tac n=a and m=c in nat_le_linear)
499 apply (safe dest!: le_iff_add [THEN iffD1])
503 lemmas pl_pl_rels' = add_commute [THEN [2] trans, THEN pl_pl_rels]
505 lemma minus_eq: "(m - k = m) = (k = 0 | m = (0 :: nat))"
508 lemma pl_pl_mm: "(a :: nat) + b = c + d ==> a - c = d - b"
511 lemmas pl_pl_mm' = add_commute [THEN [2] trans, THEN pl_pl_mm]
513 lemma min_minus [simp] : "min m (m - k) = (m - k :: nat)"
516 lemmas min_minus' [simp] = trans [OF min_max.inf_commute min_minus]
519 "(number_of b :: int) >= 0 ==> (number_of c :: int) >= 0 ==>
520 (number_of b = (number_of c :: nat)) = (b = c)"
521 apply (unfold nat_number_of_def)
523 apply (drule (2) eq_nat_nat_iff [THEN iffD1])
524 apply (simp add: number_of_eq)
527 lemmas dme = box_equals [OF div_mod_equality add_0_right add_0_right]
528 lemmas dtle = xtr3 [OF dme [symmetric] le_add1]
529 lemmas th2 = order_trans [OF order_refl [THEN [2] mult_le_mono] dtle]
532 "0 < c ==> (a >= b * c) = (a div c >= (b :: nat))"
534 apply (erule (1) xtr4 [OF div_le_mono div_mult_self_is_m])
538 lemmas td_gal_lt = td_gal [simplified le_def, simplified]
540 lemma div_mult_le: "(a :: nat) div b * b <= a"
543 apply (rule order_refl [THEN th2])
547 lemmas sdl = split_div_lemma [THEN iffD1, symmetric]
549 lemma given_quot: "f > (0 :: nat) ==> (f * l + (f - 1)) div f = l"
550 by (rule sdl, assumption) (simp (no_asm))
552 lemma given_quot_alt: "f > (0 :: nat) ==> (l * f + f - Suc 0) div f = l"
553 apply (frule given_quot)
557 apply (rule_tac f="%n. n div f" in arg_cong)
558 apply (simp add : mult_ac)
561 lemma diff_mod_le: "(a::nat) < d ==> b dvd d ==> a - a mod b <= d - b"
562 apply (unfold dvd_def)
567 apply (cases "b > 0")
568 apply (drule mult_commute [THEN xtr1])
569 apply (frule (1) td_gal_lt [THEN iffD1])
570 apply (clarsimp simp: le_simps)
571 apply (rule mult_div_cancel [THEN [2] xtr4])
572 apply (rule mult_mono)
577 "w * c < b * c ==> 0 \<le> c ==> (w + 1) * c \<le> b * (c::int)"
578 apply (rule mult_right_mono)
579 apply (rule zless_imp_add1_zle)
580 apply (erule (1) mult_right_less_imp_less)
584 lemmas less_le_mult = less_le_mult' [simplified left_distrib, simplified]
586 lemmas less_le_mult_minus = iffD2 [OF le_diff_eq less_le_mult,
587 simplified left_diff_distrib, standard]
590 assumes d: "(i::nat) \<le> j \<or> m < j'"
591 assumes R1: "i * k \<le> j * k \<Longrightarrow> R"
592 assumes R2: "Suc m * k' \<le> j' * k' \<Longrightarrow> R"
595 apply (rule R1, erule mult_le_mono1)
596 apply (rule R2, erule Suc_le_eq [THEN iffD2 [THEN mult_le_mono1]])
599 lemma lrlem: "(0::nat) < sc ==>
600 (sc - n + (n + lb * n) <= m * n) = (sc + lb * n <= m * n)"
603 apply (case_tac "sc >= n")
605 apply (insert linorder_le_less_linear [of m lb])
606 apply (erule_tac k=n and k'=n in lrlem')
611 lemma gen_minus: "0 < n ==> f n = f (Suc (n - 1))"
614 lemma mpl_lem: "j <= (i :: nat) ==> k < j ==> i - j + k < i"
615 apply (induct i, clarsimp)
616 apply (cases j, clarsimp)
620 lemma nonneg_mod_div:
621 "0 <= a ==> 0 <= b ==> 0 <= (a mod b :: int) & 0 <= a div b"
622 apply (cases "b = 0", clarsimp)
623 apply (auto intro: pos_imp_zdiv_nonneg_iff [THEN iffD2])