src/HOL/MetisExamples/BigO.thy
author wenzelm
Thu, 21 Jun 2007 20:07:26 +0200
changeset 23464 bc2563c37b1a
parent 23449 dd874e6a3282
child 23477 f4b83f03cac9
permissions -rw-r--r--
tuned proofs -- avoid implicit prems;
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
23449
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
     1
(*  Title:      HOL/MetisExamples/BigO.thy
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
     2
    ID:         $Id$
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
     3
    Author:     Lawrence C Paulson, Cambridge University Computer Laboratory
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
     4
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
     5
Testing the metis method
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
     6
*)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
     7
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
     8
header {* Big O notation *}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
     9
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    10
theory BigO
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    11
imports SetsAndFunctions 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    12
begin
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    13
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    14
subsection {* Definitions *}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    15
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    16
constdefs 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    17
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    18
  bigo :: "('a => 'b::ordered_idom) => ('a => 'b) set"    ("(1O'(_'))")
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    19
  "O(f::('a => 'b)) ==   {h. EX c. ALL x. abs (h x) <= c * abs (f x)}"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    20
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    21
ML{*ResAtp.problem_name := "BigO__bigo_pos_const"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    22
lemma bigo_pos_const: "(EX (c::'a::ordered_idom). 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    23
    ALL x. (abs (h x)) <= (c * (abs (f x))))
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    24
      = (EX c. 0 < c & (ALL x. (abs(h x)) <= (c * (abs (f x)))))"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    25
  apply auto
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    26
  apply (case_tac "c = 0", simp)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    27
  apply (rule_tac x = "1" in exI, simp)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    28
  apply (rule_tac x = "abs c" in exI, auto);
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    29
txt{*Version 1: one-shot proof. MUCH SLOWER with types: 24 versus 6.7 seconds*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    30
  apply (metis abs_ge_minus_self abs_ge_zero abs_minus_cancel abs_of_nonneg equation_minus_iff Orderings.xt1(6) abs_le_mult)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    31
  done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    32
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    33
(*** Now various verions with an increasing modulus ***)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    34
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    35
ML{*ResReconstruct.modulus := 1*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    36
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    37
lemma bigo_pos_const: "(EX (c::'a::ordered_idom). 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    38
    ALL x. (abs (h x)) <= (c * (abs (f x))))
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    39
      = (EX c. 0 < c & (ALL x. (abs(h x)) <= (c * (abs (f x)))))"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    40
  apply auto
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    41
  apply (case_tac "c = 0", simp)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    42
  apply (rule_tac x = "1" in exI, simp)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    43
  apply (rule_tac x = "abs c" in exI, auto)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    44
(*hand-modified to give 'a sort ordered_idom and X3 type 'a*)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    45
proof (neg_clausify)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    46
fix c x
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    47
assume 0: "\<And>A. \<bar>h A\<bar> \<le> c * \<bar>f A\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    48
assume 1: "c \<noteq> (0\<Colon>'a::ordered_idom)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    49
assume 2: "\<not> \<bar>h x\<bar> \<le> \<bar>c\<bar> * \<bar>f x\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    50
have 3: "\<And>X1 X3. \<bar>h X3\<bar> < X1 \<or> \<not> c * \<bar>f X3\<bar> < X1"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    51
  by (metis order_le_less_trans 0)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    52
have 4: "\<And>X3. (1\<Colon>'a) * X3 \<le> X3 \<or> \<not> (1\<Colon>'a) \<le> (1\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    53
  by (metis mult_le_cancel_right2 order_refl)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    54
have 5: "\<And>X3. (1\<Colon>'a) * X3 \<le> X3"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    55
  by (metis 4 order_refl)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    56
have 6: "\<And>X3. \<bar>0\<Colon>'a\<bar> = \<bar>X3\<bar> * (0\<Colon>'a) \<or> \<not> (0\<Colon>'a) \<le> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    57
  by (metis abs_mult_pos mult_cancel_right1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    58
have 7: "\<bar>0\<Colon>'a\<bar> = (0\<Colon>'a) \<or> \<not> (0\<Colon>'a) \<le> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    59
  by (metis 6 mult_cancel_right1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    60
have 8: "\<bar>0\<Colon>'a\<bar> = (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    61
  by (metis 7 order_refl)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    62
have 9: "\<not> (0\<Colon>'a) < (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    63
  by (metis abs_not_less_zero 8)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    64
have 10: "\<bar>(1\<Colon>'a) * (0\<Colon>'a)\<bar> = - ((1\<Colon>'a) * (0\<Colon>'a))"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    65
  by (metis abs_of_nonpos 5)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    66
have 11: "(0\<Colon>'a) = - ((1\<Colon>'a) * (0\<Colon>'a))"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    67
  by (metis 10 mult_cancel_right1 8)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    68
have 12: "(0\<Colon>'a) = - (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    69
  by (metis 11 mult_cancel_right1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    70
have 13: "\<And>X3. \<bar>X3\<bar> = X3 \<or> X3 \<le> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    71
  by (metis abs_of_nonneg linorder_linear)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    72
have 14: "c \<le> (0\<Colon>'a) \<or> \<not> \<bar>h x\<bar> \<le> c * \<bar>f x\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    73
  by (metis 2 13)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    74
have 15: "c \<le> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    75
  by (metis 14 0)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    76
have 16: "c = (0\<Colon>'a) \<or> c < (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    77
  by (metis linorder_antisym_conv2 15)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    78
have 17: "\<bar>c\<bar> = - c"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    79
  by (metis abs_of_nonpos 15)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    80
have 18: "c < (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    81
  by (metis 16 1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    82
have 19: "\<not> \<bar>h x\<bar> \<le> - c * \<bar>f x\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    83
  by (metis 2 17)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    84
have 20: "\<And>X3. X3 * (1\<Colon>'a) = X3"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    85
  by (metis mult_cancel_right1 AC_mult.f.commute)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    86
have 21: "\<And>X3. (0\<Colon>'a) \<le> X3 * X3"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    87
  by (metis zero_le_square AC_mult.f.commute)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    88
have 22: "(0\<Colon>'a) \<le> (1\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    89
  by (metis 21 mult_cancel_left1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    90
have 23: "\<And>X3. (0\<Colon>'a) = X3 \<or> (0\<Colon>'a) \<noteq> - X3"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    91
  by (metis neg_equal_iff_equal 12)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    92
have 24: "\<And>X3. (0\<Colon>'a) = - X3 \<or> X3 \<noteq> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    93
  by (metis 23 minus_equation_iff)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    94
have 25: "\<And>X3. \<bar>0\<Colon>'a\<bar> = \<bar>X3\<bar> \<or> X3 \<noteq> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    95
  by (metis abs_minus_cancel 24)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    96
have 26: "\<And>X3. (0\<Colon>'a) = \<bar>X3\<bar> \<or> X3 \<noteq> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    97
  by (metis 25 8)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    98
have 27: "\<And>X1 X3. (0\<Colon>'a) * \<bar>X1\<bar> = \<bar>X3 * X1\<bar> \<or> X3 \<noteq> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
    99
  by (metis abs_mult 26)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   100
have 28: "\<And>X1 X3. (0\<Colon>'a) = \<bar>X3 * X1\<bar> \<or> X3 \<noteq> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   101
  by (metis 27 mult_cancel_left1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   102
have 29: "\<And>X1 X3. (0\<Colon>'a) = X3 * X1 \<or> (0\<Colon>'a) < (0\<Colon>'a) \<or> X3 \<noteq> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   103
  by (metis zero_less_abs_iff 28)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   104
have 30: "\<And>X1 X3. X3 * X1 = (0\<Colon>'a) \<or> X3 \<noteq> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   105
  by (metis 29 9)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   106
have 31: "\<And>X1 X3. (0\<Colon>'a) = X1 * X3 \<or> X3 \<noteq> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   107
  by (metis AC_mult.f.commute 30)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   108
have 32: "\<And>X1 X3. (0\<Colon>'a) = \<bar>X3 * X1\<bar> \<or> \<bar>X1\<bar> \<noteq> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   109
  by (metis abs_mult 31)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   110
have 33: "\<And>X3::'a. \<bar>X3 * X3\<bar> = X3 * X3"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   111
  by (metis abs_mult_self abs_mult AC_mult.f.commute)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   112
have 34: "\<And>X3. (0\<Colon>'a) \<le> \<bar>X3\<bar> \<or> \<not> (0\<Colon>'a) \<le> (1\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   113
  by (metis abs_ge_zero abs_mult_pos 20)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   114
have 35: "\<And>X3. (0\<Colon>'a) \<le> \<bar>X3\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   115
  by (metis 34 22)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   116
have 36: "\<And>X3. X3 * (1\<Colon>'a) = (0\<Colon>'a) \<or> \<bar>X3\<bar> \<noteq> (0\<Colon>'a) \<or> \<not> (0\<Colon>'a) \<le> (1\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   117
  by (metis abs_eq_0 abs_mult_pos 20)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   118
have 37: "\<And>X3. X3 = (0\<Colon>'a) \<or> \<bar>X3\<bar> \<noteq> (0\<Colon>'a) \<or> \<not> (0\<Colon>'a) \<le> (1\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   119
  by (metis 36 20)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   120
have 38: "\<And>X3. X3 = (0\<Colon>'a) \<or> \<bar>X3\<bar> \<noteq> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   121
  by (metis 37 22)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   122
have 39: "\<And>X1 X3. X3 * X1 = (0\<Colon>'a) \<or> \<bar>X1\<bar> \<noteq> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   123
  by (metis 38 32)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   124
have 40: "\<And>X3::'a. \<bar>\<bar>X3\<bar>\<bar> = \<bar>X3\<bar> \<or> \<not> (0\<Colon>'a) \<le> (1\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   125
  by (metis abs_idempotent abs_mult_pos 20)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   126
have 41: "\<And>X3::'a. \<bar>\<bar>X3\<bar>\<bar> = \<bar>X3\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   127
  by (metis 40 22)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   128
have 42: "\<And>X3. \<not> \<bar>X3\<bar> < (0\<Colon>'a) \<or> \<not> (0\<Colon>'a) \<le> (1\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   129
  by (metis abs_not_less_zero abs_mult_pos 20)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   130
have 43: "\<And>X3. \<not> \<bar>X3\<bar> < (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   131
  by (metis 42 22)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   132
have 44: "\<And>X3. X3 * (1\<Colon>'a) = (0\<Colon>'a) \<or> \<not> \<bar>X3\<bar> \<le> (0\<Colon>'a) \<or> \<not> (0\<Colon>'a) \<le> (1\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   133
  by (metis abs_le_zero_iff abs_mult_pos 20)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   134
have 45: "\<And>X3. X3 = (0\<Colon>'a) \<or> \<not> \<bar>X3\<bar> \<le> (0\<Colon>'a) \<or> \<not> (0\<Colon>'a) \<le> (1\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   135
  by (metis 44 20)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   136
have 46: "\<And>X3. X3 = (0\<Colon>'a) \<or> \<not> \<bar>X3\<bar> \<le> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   137
  by (metis 45 22)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   138
have 47: "\<And>X3. X3 * X3 = (0\<Colon>'a) \<or> \<not> X3 * X3 \<le> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   139
  by (metis 46 33)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   140
have 48: "\<And>X3. X3 * X3 = (0\<Colon>'a) \<or> \<not> X3 \<le> (0\<Colon>'a) \<or> \<not> (0\<Colon>'a) \<le> X3"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   141
  by (metis 47 mult_le_0_iff)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   142
have 49: "\<And>X3. \<bar>X3\<bar> = (0\<Colon>'a) \<or> \<not> X3 \<le> (0\<Colon>'a) \<or> \<not> (0\<Colon>'a) \<le> X3"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   143
  by (metis mult_eq_0_iff abs_mult_self 48)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   144
have 50: "\<And>X1 X3.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   145
   (0\<Colon>'a) * \<bar>X1\<bar> = \<bar>\<bar>X3 * X1\<bar>\<bar> \<or>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   146
   \<not> (0\<Colon>'a) \<le> \<bar>X1\<bar> \<or> \<not> \<bar>X3\<bar> \<le> (0\<Colon>'a) \<or> \<not> (0\<Colon>'a) \<le> \<bar>X3\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   147
  by (metis abs_mult_pos abs_mult 49)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   148
have 51: "\<And>X1 X3. X3 * X1 = (0\<Colon>'a) \<or> \<not> X1 \<le> (0\<Colon>'a) \<or> \<not> (0\<Colon>'a) \<le> X1"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   149
  by (metis 39 49)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   150
have 52: "\<And>X1 X3.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   151
   (0\<Colon>'a) = \<bar>\<bar>X3 * X1\<bar>\<bar> \<or>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   152
   \<not> (0\<Colon>'a) \<le> \<bar>X1\<bar> \<or> \<not> \<bar>X3\<bar> \<le> (0\<Colon>'a) \<or> \<not> (0\<Colon>'a) \<le> \<bar>X3\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   153
  by (metis 50 mult_cancel_left1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   154
have 53: "\<And>X1 X3.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   155
   (0\<Colon>'a) = \<bar>X3 * X1\<bar> \<or> \<not> (0\<Colon>'a) \<le> \<bar>X1\<bar> \<or> \<not> \<bar>X3\<bar> \<le> (0\<Colon>'a) \<or> \<not> (0\<Colon>'a) \<le> \<bar>X3\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   156
  by (metis 52 41)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   157
have 54: "\<And>X1 X3. (0\<Colon>'a) = \<bar>X3 * X1\<bar> \<or> \<not> \<bar>X3\<bar> \<le> (0\<Colon>'a) \<or> \<not> (0\<Colon>'a) \<le> \<bar>X3\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   158
  by (metis 53 35)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   159
have 55: "\<And>X1 X3. (0\<Colon>'a) = \<bar>X3 * X1\<bar> \<or> \<not> \<bar>X3\<bar> \<le> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   160
  by (metis 54 35)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   161
have 56: "\<And>X1 X3. \<bar>X1 * X3\<bar> = (0\<Colon>'a) \<or> \<not> \<bar>X3\<bar> \<le> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   162
  by (metis 55 AC_mult.f.commute)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   163
have 57: "\<And>X1 X3. X3 * X1 = (0\<Colon>'a) \<or> \<not> \<bar>X1\<bar> \<le> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   164
  by (metis 38 56)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   165
have 58: "\<And>X3. \<bar>h X3\<bar> \<le> (0\<Colon>'a) \<or> \<not> \<bar>f X3\<bar> \<le> (0\<Colon>'a) \<or> \<not> (0\<Colon>'a) \<le> \<bar>f X3\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   166
  by (metis 0 51)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   167
have 59: "\<And>X3. \<bar>h X3\<bar> \<le> (0\<Colon>'a) \<or> \<not> \<bar>f X3\<bar> \<le> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   168
  by (metis 58 35)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   169
have 60: "\<And>X3. \<bar>h X3\<bar> \<le> (0\<Colon>'a) \<or> (0\<Colon>'a) < \<bar>f X3\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   170
  by (metis 59 linorder_not_le)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   171
have 61: "\<And>X1 X3. X3 * X1 = (0\<Colon>'a) \<or> (0\<Colon>'a) < \<bar>X1\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   172
  by (metis 57 linorder_not_le)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   173
have 62: "(0\<Colon>'a) < \<bar>\<bar>f x\<bar>\<bar> \<or> \<not> \<bar>h x\<bar> \<le> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   174
  by (metis 19 61)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   175
have 63: "(0\<Colon>'a) < \<bar>f x\<bar> \<or> \<not> \<bar>h x\<bar> \<le> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   176
  by (metis 62 41)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   177
have 64: "(0\<Colon>'a) < \<bar>f x\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   178
  by (metis 63 60)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   179
have 65: "\<And>X3. \<bar>h X3\<bar> < (0\<Colon>'a) \<or> \<not> c < (0\<Colon>'a) \<or> \<not> (0\<Colon>'a) < \<bar>f X3\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   180
  by (metis 3 mult_less_0_iff)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   181
have 66: "\<And>X3. \<bar>h X3\<bar> < (0\<Colon>'a) \<or> \<not> (0\<Colon>'a) < \<bar>f X3\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   182
  by (metis 65 18)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   183
have 67: "\<And>X3. \<not> (0\<Colon>'a) < \<bar>f X3\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   184
  by (metis 66 43)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   185
show "False"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   186
  by (metis 67 64)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   187
qed
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   188
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   189
lemma (*bigo_pos_const:*) "(EX (c::'a::ordered_idom). 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   190
    ALL x. (abs (h x)) <= (c * (abs (f x))))
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   191
      = (EX c. 0 < c & (ALL x. (abs(h x)) <= (c * (abs (f x)))))"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   192
  apply auto
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   193
  apply (case_tac "c = 0", simp)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   194
  apply (rule_tac x = "1" in exI, simp)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   195
  apply (rule_tac x = "abs c" in exI, auto);
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   196
ML{*ResReconstruct.modulus:=2*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   197
proof (neg_clausify)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   198
fix c x
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   199
assume 0: "\<And>A. \<bar>h A\<bar> \<le> c * \<bar>f A\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   200
assume 1: "c \<noteq> (0\<Colon>'a::ordered_idom)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   201
assume 2: "\<not> \<bar>h x\<bar> \<le> \<bar>c\<bar> * \<bar>f x\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   202
have 3: "\<And>X3. (1\<Colon>'a) * X3 \<le> X3"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   203
  by (metis mult_le_cancel_right2 order_refl order_refl)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   204
have 4: "\<bar>0\<Colon>'a\<bar> = (0\<Colon>'a) \<or> \<not> (0\<Colon>'a) \<le> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   205
  by (metis abs_mult_pos mult_cancel_right1 mult_cancel_right1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   206
have 5: "\<not> (0\<Colon>'a) < (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   207
  by (metis abs_not_less_zero 4 order_refl)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   208
have 6: "(0\<Colon>'a) = - ((1\<Colon>'a) * (0\<Colon>'a))"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   209
  by (metis abs_of_nonpos 3 mult_cancel_right1 4 order_refl)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   210
have 7: "\<And>X3. \<bar>X3\<bar> = X3 \<or> X3 \<le> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   211
  by (metis abs_of_nonneg linorder_linear)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   212
have 8: "c \<le> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   213
  by (metis 2 7 0)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   214
have 9: "\<bar>c\<bar> = - c"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   215
  by (metis abs_of_nonpos 8)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   216
have 10: "\<not> \<bar>h x\<bar> \<le> - c * \<bar>f x\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   217
  by (metis 2 9)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   218
have 11: "\<And>X3. X3 * (1\<Colon>'a) = X3"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   219
  by (metis mult_cancel_right1 AC_mult.f.commute)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   220
have 12: "(0\<Colon>'a) \<le> (1\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   221
  by (metis zero_le_square AC_mult.f.commute mult_cancel_left1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   222
have 13: "\<And>X3. (0\<Colon>'a) = - X3 \<or> X3 \<noteq> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   223
  by (metis neg_equal_iff_equal 6 mult_cancel_right1 minus_equation_iff)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   224
have 14: "\<And>X3. (0\<Colon>'a) = \<bar>X3\<bar> \<or> X3 \<noteq> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   225
  by (metis abs_minus_cancel 13 4 order_refl)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   226
have 15: "\<And>X1 X3. (0\<Colon>'a) = \<bar>X3 * X1\<bar> \<or> X3 \<noteq> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   227
  by (metis abs_mult 14 mult_cancel_left1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   228
have 16: "\<And>X1 X3. X3 * X1 = (0\<Colon>'a) \<or> X3 \<noteq> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   229
  by (metis zero_less_abs_iff 15 5)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   230
have 17: "\<And>X1 X3. (0\<Colon>'a) = \<bar>X3 * X1\<bar> \<or> \<bar>X1\<bar> \<noteq> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   231
  by (metis abs_mult AC_mult.f.commute 16)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   232
have 18: "\<And>X3. (0\<Colon>'a) \<le> \<bar>X3\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   233
  by (metis abs_ge_zero abs_mult_pos 11 12)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   234
have 19: "\<And>X3. X3 * (1\<Colon>'a) = (0\<Colon>'a) \<or> \<bar>X3\<bar> \<noteq> (0\<Colon>'a) \<or> \<not> (0\<Colon>'a) \<le> (1\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   235
  by (metis abs_eq_0 abs_mult_pos 11)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   236
have 20: "\<And>X3. X3 = (0\<Colon>'a) \<or> \<bar>X3\<bar> \<noteq> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   237
  by (metis 19 11 12)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   238
have 21: "\<And>X3::'a. \<bar>\<bar>X3\<bar>\<bar> = \<bar>X3\<bar> \<or> \<not> (0\<Colon>'a) \<le> (1\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   239
  by (metis abs_idempotent abs_mult_pos 11)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   240
have 22: "\<And>X3. \<not> \<bar>X3\<bar> < (0\<Colon>'a) \<or> \<not> (0\<Colon>'a) \<le> (1\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   241
  by (metis abs_not_less_zero abs_mult_pos 11)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   242
have 23: "\<And>X3. X3 = (0\<Colon>'a) \<or> \<not> \<bar>X3\<bar> \<le> (0\<Colon>'a) \<or> \<not> (0\<Colon>'a) \<le> (1\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   243
  by (metis abs_le_zero_iff abs_mult_pos 11 11)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   244
have 24: "\<And>X3. X3 * X3 = (0\<Colon>'a) \<or> \<not> X3 * X3 \<le> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   245
  by (metis 23 12 abs_mult_self abs_mult AC_mult.f.commute)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   246
have 25: "\<And>X3. \<bar>X3\<bar> = (0\<Colon>'a) \<or> \<not> X3 \<le> (0\<Colon>'a) \<or> \<not> (0\<Colon>'a) \<le> X3"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   247
  by (metis mult_eq_0_iff abs_mult_self 24 mult_le_0_iff)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   248
have 26: "\<And>X1 X3. X3 * X1 = (0\<Colon>'a) \<or> \<not> X1 \<le> (0\<Colon>'a) \<or> \<not> (0\<Colon>'a) \<le> X1"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   249
  by (metis 20 17 25)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   250
have 27: "\<And>X1 X3.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   251
   (0\<Colon>'a) = \<bar>X3 * X1\<bar> \<or> \<not> (0\<Colon>'a) \<le> \<bar>X1\<bar> \<or> \<not> \<bar>X3\<bar> \<le> (0\<Colon>'a) \<or> \<not> (0\<Colon>'a) \<le> \<bar>X3\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   252
  by (metis abs_mult_pos abs_mult 25 mult_cancel_left1 21 12)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   253
have 28: "\<And>X1 X3. (0\<Colon>'a) = \<bar>X3 * X1\<bar> \<or> \<not> \<bar>X3\<bar> \<le> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   254
  by (metis 27 18 18)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   255
have 29: "\<And>X1 X3. X3 * X1 = (0\<Colon>'a) \<or> \<not> \<bar>X1\<bar> \<le> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   256
  by (metis 20 28 AC_mult.f.commute)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   257
have 30: "\<And>X3. \<bar>h X3\<bar> \<le> (0\<Colon>'a) \<or> \<not> \<bar>f X3\<bar> \<le> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   258
  by (metis 0 26 18)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   259
have 31: "\<And>X1 X3. X3 * X1 = (0\<Colon>'a) \<or> (0\<Colon>'a) < \<bar>X1\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   260
  by (metis 29 linorder_not_le)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   261
have 32: "(0\<Colon>'a) < \<bar>f x\<bar> \<or> \<not> \<bar>h x\<bar> \<le> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   262
  by (metis 10 31 21 12)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   263
have 33: "\<And>X3. \<bar>h X3\<bar> < (0\<Colon>'a) \<or> \<not> c < (0\<Colon>'a) \<or> \<not> (0\<Colon>'a) < \<bar>f X3\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   264
  by (metis order_le_less_trans 0 mult_less_0_iff)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   265
have 34: "\<And>X3. \<not> (0\<Colon>'a) < \<bar>f X3\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   266
  by (metis 33 linorder_antisym_conv2 8 1 22 12)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   267
show "False"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   268
  by (metis 34 32 30 linorder_not_le)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   269
qed
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   270
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   271
lemma (*bigo_pos_const:*) "(EX (c::'a::ordered_idom). 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   272
    ALL x. (abs (h x)) <= (c * (abs (f x))))
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   273
      = (EX c. 0 < c & (ALL x. (abs(h x)) <= (c * (abs (f x)))))"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   274
  apply auto
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   275
  apply (case_tac "c = 0", simp)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   276
  apply (rule_tac x = "1" in exI, simp)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   277
  apply (rule_tac x = "abs c" in exI, auto);
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   278
ML{*ResReconstruct.modulus:=3*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   279
proof (neg_clausify)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   280
fix c x
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   281
assume 0: "\<And>A\<Colon>'b\<Colon>type.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   282
   \<bar>(h\<Colon>'b\<Colon>type \<Rightarrow> 'a\<Colon>ordered_idom) A\<bar>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   283
   \<le> (c\<Colon>'a\<Colon>ordered_idom) * \<bar>(f\<Colon>'b\<Colon>type \<Rightarrow> 'a\<Colon>ordered_idom) A\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   284
assume 1: "(c\<Colon>'a\<Colon>ordered_idom) \<noteq> (0\<Colon>'a\<Colon>ordered_idom)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   285
assume 2: "\<not> \<bar>(h\<Colon>'b\<Colon>type \<Rightarrow> 'a\<Colon>ordered_idom) (x\<Colon>'b\<Colon>type)\<bar>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   286
  \<le> \<bar>c\<Colon>'a\<Colon>ordered_idom\<bar> * \<bar>(f\<Colon>'b\<Colon>type \<Rightarrow> 'a\<Colon>ordered_idom) x\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   287
have 3: "\<And>X3\<Colon>'a\<Colon>ordered_idom. (1\<Colon>'a\<Colon>ordered_idom) * X3 \<le> X3"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   288
  by (metis mult_le_cancel_right2 order_refl order_refl)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   289
have 4: "\<bar>0\<Colon>'a\<Colon>ordered_idom\<bar> = (0\<Colon>'a\<Colon>ordered_idom)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   290
  by (metis abs_mult_pos mult_cancel_right1 mult_cancel_right1 order_refl)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   291
have 5: "(0\<Colon>'a\<Colon>ordered_idom) = - ((1\<Colon>'a\<Colon>ordered_idom) * (0\<Colon>'a\<Colon>ordered_idom))"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   292
  by (metis abs_of_nonpos 3 mult_cancel_right1 4)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   293
have 6: "(c\<Colon>'a\<Colon>ordered_idom) \<le> (0\<Colon>'a\<Colon>ordered_idom)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   294
  by (metis 2 abs_of_nonneg linorder_linear 0)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   295
have 7: "(c\<Colon>'a\<Colon>ordered_idom) < (0\<Colon>'a\<Colon>ordered_idom)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   296
  by (metis linorder_antisym_conv2 6 1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   297
have 8: "\<And>X3\<Colon>'a\<Colon>ordered_idom. X3 * (1\<Colon>'a\<Colon>ordered_idom) = X3"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   298
  by (metis mult_cancel_right1 AC_mult.f.commute)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   299
have 9: "\<And>X3\<Colon>'a\<Colon>ordered_idom. (0\<Colon>'a\<Colon>ordered_idom) = X3 \<or> (0\<Colon>'a\<Colon>ordered_idom) \<noteq> - X3"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   300
  by (metis neg_equal_iff_equal 5 mult_cancel_right1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   301
have 10: "\<And>X3\<Colon>'a\<Colon>ordered_idom. (0\<Colon>'a\<Colon>ordered_idom) = \<bar>X3\<bar> \<or> X3 \<noteq> (0\<Colon>'a\<Colon>ordered_idom)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   302
  by (metis abs_minus_cancel 9 minus_equation_iff 4)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   303
have 11: "\<And>(X1\<Colon>'a\<Colon>ordered_idom) X3\<Colon>'a\<Colon>ordered_idom.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   304
   (0\<Colon>'a\<Colon>ordered_idom) * \<bar>X1\<bar> = \<bar>X3 * X1\<bar> \<or> X3 \<noteq> (0\<Colon>'a\<Colon>ordered_idom)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   305
  by (metis abs_mult 10)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   306
have 12: "\<And>(X1\<Colon>'a\<Colon>ordered_idom) X3\<Colon>'a\<Colon>ordered_idom.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   307
   X3 * X1 = (0\<Colon>'a\<Colon>ordered_idom) \<or> X3 \<noteq> (0\<Colon>'a\<Colon>ordered_idom)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   308
  by (metis zero_less_abs_iff 11 mult_cancel_left1 abs_not_less_zero 4)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   309
have 13: "\<And>X3\<Colon>'a\<Colon>ordered_idom. \<bar>X3 * X3\<bar> = X3 * X3"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   310
  by (metis abs_mult_self abs_mult AC_mult.f.commute)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   311
have 14: "\<And>X3\<Colon>'a\<Colon>ordered_idom. (0\<Colon>'a\<Colon>ordered_idom) \<le> \<bar>X3\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   312
  by (metis abs_ge_zero abs_mult_pos 8 zero_le_square AC_mult.f.commute mult_cancel_left1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   313
have 15: "\<And>X3\<Colon>'a\<Colon>ordered_idom.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   314
   X3 = (0\<Colon>'a\<Colon>ordered_idom) \<or>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   315
   \<bar>X3\<bar> \<noteq> (0\<Colon>'a\<Colon>ordered_idom) \<or> \<not> (0\<Colon>'a\<Colon>ordered_idom) \<le> (1\<Colon>'a\<Colon>ordered_idom)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   316
  by (metis abs_eq_0 abs_mult_pos 8 8)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   317
have 16: "\<And>X3\<Colon>'a\<Colon>ordered_idom.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   318
   \<bar>\<bar>X3\<bar>\<bar> = \<bar>X3\<bar> \<or> \<not> (0\<Colon>'a\<Colon>ordered_idom) \<le> (1\<Colon>'a\<Colon>ordered_idom)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   319
  by (metis abs_idempotent abs_mult_pos 8)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   320
have 17: "\<And>X3\<Colon>'a\<Colon>ordered_idom. \<not> \<bar>X3\<bar> < (0\<Colon>'a\<Colon>ordered_idom)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   321
  by (metis abs_not_less_zero abs_mult_pos 8 zero_le_square AC_mult.f.commute mult_cancel_left1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   322
have 18: "\<And>X3\<Colon>'a\<Colon>ordered_idom.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   323
   X3 = (0\<Colon>'a\<Colon>ordered_idom) \<or>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   324
   \<not> \<bar>X3\<bar> \<le> (0\<Colon>'a\<Colon>ordered_idom) \<or>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   325
   \<not> (0\<Colon>'a\<Colon>ordered_idom) \<le> (1\<Colon>'a\<Colon>ordered_idom)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   326
  by (metis abs_le_zero_iff abs_mult_pos 8 8)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   327
have 19: "\<And>X3\<Colon>'a\<Colon>ordered_idom.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   328
   X3 * X3 = (0\<Colon>'a\<Colon>ordered_idom) \<or>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   329
   \<not> X3 \<le> (0\<Colon>'a\<Colon>ordered_idom) \<or> \<not> (0\<Colon>'a\<Colon>ordered_idom) \<le> X3"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   330
  by (metis 18 zero_le_square AC_mult.f.commute mult_cancel_left1 13 mult_le_0_iff)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   331
have 20: "\<And>(X1\<Colon>'a\<Colon>ordered_idom) X3\<Colon>'a\<Colon>ordered_idom.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   332
   X3 * X1 = (0\<Colon>'a\<Colon>ordered_idom) \<or>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   333
   \<not> X1 \<le> (0\<Colon>'a\<Colon>ordered_idom) \<or> \<not> (0\<Colon>'a\<Colon>ordered_idom) \<le> X1"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   334
  by (metis 15 zero_le_square AC_mult.f.commute mult_cancel_left1 abs_mult AC_mult.f.commute 12 mult_eq_0_iff abs_mult_self 19)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   335
have 21: "\<And>(X1\<Colon>'a\<Colon>ordered_idom) X3\<Colon>'a\<Colon>ordered_idom.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   336
   (0\<Colon>'a\<Colon>ordered_idom) = \<bar>X3 * X1\<bar> \<or>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   337
   \<not> \<bar>X3\<bar> \<le> (0\<Colon>'a\<Colon>ordered_idom) \<or> \<not> (0\<Colon>'a\<Colon>ordered_idom) \<le> \<bar>X3\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   338
  by (metis abs_mult_pos abs_mult mult_eq_0_iff abs_mult_self 19 mult_cancel_left1 16 zero_le_square AC_mult.f.commute mult_cancel_left1 14)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   339
have 22: "\<And>(X1\<Colon>'a\<Colon>ordered_idom) X3\<Colon>'a\<Colon>ordered_idom.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   340
   X3 * X1 = (0\<Colon>'a\<Colon>ordered_idom) \<or> \<not> \<bar>X1\<bar> \<le> (0\<Colon>'a\<Colon>ordered_idom)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   341
  by (metis 15 zero_le_square AC_mult.f.commute mult_cancel_left1 21 14 AC_mult.f.commute)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   342
have 23: "\<And>X3\<Colon>'b\<Colon>type.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   343
   \<bar>(h\<Colon>'b\<Colon>type \<Rightarrow> 'a\<Colon>ordered_idom) X3\<bar> \<le> (0\<Colon>'a\<Colon>ordered_idom) \<or>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   344
   (0\<Colon>'a\<Colon>ordered_idom) < \<bar>(f\<Colon>'b\<Colon>type \<Rightarrow> 'a\<Colon>ordered_idom) X3\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   345
  by (metis 0 20 14 linorder_not_le)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   346
have 24: "(0\<Colon>'a\<Colon>ordered_idom) < \<bar>(f\<Colon>'b\<Colon>type \<Rightarrow> 'a\<Colon>ordered_idom) (x\<Colon>'b\<Colon>type)\<bar> \<or>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   347
\<not> \<bar>(h\<Colon>'b\<Colon>type \<Rightarrow> 'a\<Colon>ordered_idom) x\<bar> \<le> (0\<Colon>'a\<Colon>ordered_idom)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   348
  by (metis 2 abs_of_nonpos 6 22 linorder_not_le 16 zero_le_square AC_mult.f.commute mult_cancel_left1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   349
have 25: "\<And>X3\<Colon>'b\<Colon>type.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   350
   \<bar>(h\<Colon>'b\<Colon>type \<Rightarrow> 'a\<Colon>ordered_idom) X3\<bar> < (0\<Colon>'a\<Colon>ordered_idom) \<or>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   351
   \<not> (0\<Colon>'a\<Colon>ordered_idom) < \<bar>(f\<Colon>'b\<Colon>type \<Rightarrow> 'a\<Colon>ordered_idom) X3\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   352
  by (metis order_le_less_trans 0 mult_less_0_iff 7)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   353
show "False"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   354
  by (metis 25 17 24 23)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   355
qed
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   356
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   357
lemma (*bigo_pos_const:*) "(EX (c::'a::ordered_idom). 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   358
    ALL x. (abs (h x)) <= (c * (abs (f x))))
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   359
      = (EX c. 0 < c & (ALL x. (abs(h x)) <= (c * (abs (f x)))))"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   360
  apply auto
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   361
  apply (case_tac "c = 0", simp)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   362
  apply (rule_tac x = "1" in exI, simp)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   363
  apply (rule_tac x = "abs c" in exI, auto);
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   364
ML{*ResReconstruct.modulus:=4*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   365
ML{*ResReconstruct.recon_sorts:=false*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   366
proof (neg_clausify)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   367
fix c x
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   368
assume 0: "\<And>A. \<bar>h A\<bar> \<le> c * \<bar>f A\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   369
assume 1: "c \<noteq> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   370
assume 2: "\<not> \<bar>h x\<bar> \<le> \<bar>c\<bar> * \<bar>f x\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   371
have 3: "\<And>X3. (1\<Colon>'a) * X3 \<le> X3"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   372
  by (metis mult_le_cancel_right2 order_refl order_refl)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   373
have 4: "\<not> (0\<Colon>'a) < (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   374
  by (metis abs_not_less_zero abs_mult_pos mult_cancel_right1 mult_cancel_right1 order_refl)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   375
have 5: "c \<le> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   376
  by (metis 2 abs_of_nonneg linorder_linear 0)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   377
have 6: "\<not> \<bar>h x\<bar> \<le> - c * \<bar>f x\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   378
  by (metis 2 abs_of_nonpos 5)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   379
have 7: "(0\<Colon>'a) \<le> (1\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   380
  by (metis zero_le_square AC_mult.f.commute mult_cancel_left1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   381
have 8: "\<And>X3. (0\<Colon>'a) = \<bar>X3\<bar> \<or> X3 \<noteq> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   382
  by (metis abs_minus_cancel neg_equal_iff_equal abs_of_nonpos 3 mult_cancel_right1 abs_mult_pos mult_cancel_right1 mult_cancel_right1 order_refl mult_cancel_right1 minus_equation_iff abs_mult_pos mult_cancel_right1 mult_cancel_right1 order_refl)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   383
have 9: "\<And>X1 X3. (0\<Colon>'a) = \<bar>X3 * X1\<bar> \<or> X3 \<noteq> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   384
  by (metis abs_mult 8 mult_cancel_left1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   385
have 10: "\<And>X1 X3. (0\<Colon>'a) = \<bar>X3 * X1\<bar> \<or> \<bar>X1\<bar> \<noteq> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   386
  by (metis abs_mult AC_mult.f.commute zero_less_abs_iff 9 4)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   387
have 11: "\<And>X3. (0\<Colon>'a) \<le> \<bar>X3\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   388
  by (metis abs_ge_zero abs_mult_pos mult_cancel_right1 AC_mult.f.commute 7)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   389
have 12: "\<And>X3. X3 = (0\<Colon>'a) \<or> \<bar>X3\<bar> \<noteq> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   390
  by (metis abs_eq_0 abs_mult_pos mult_cancel_right1 AC_mult.f.commute mult_cancel_right1 AC_mult.f.commute 7)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   391
have 13: "\<And>X3. \<not> \<bar>X3\<bar> < (0\<Colon>'a) \<or> \<not> (0\<Colon>'a) \<le> (1\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   392
  by (metis abs_not_less_zero abs_mult_pos mult_cancel_right1 AC_mult.f.commute)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   393
have 14: "\<And>X3. X3 = (0\<Colon>'a) \<or> \<not> \<bar>X3\<bar> \<le> (0\<Colon>'a) \<or> \<not> (0\<Colon>'a) \<le> (1\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   394
  by (metis abs_le_zero_iff abs_mult_pos mult_cancel_right1 AC_mult.f.commute mult_cancel_right1 AC_mult.f.commute)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   395
have 15: "\<And>X3. \<bar>X3\<bar> = (0\<Colon>'a) \<or> \<not> X3 \<le> (0\<Colon>'a) \<or> \<not> (0\<Colon>'a) \<le> X3"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   396
  by (metis mult_eq_0_iff abs_mult_self 14 7 abs_mult_self abs_mult AC_mult.f.commute mult_le_0_iff)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   397
have 16: "\<And>X1 X3.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   398
   (0\<Colon>'a) = \<bar>X3 * X1\<bar> \<or> \<not> (0\<Colon>'a) \<le> \<bar>X1\<bar> \<or> \<not> \<bar>X3\<bar> \<le> (0\<Colon>'a) \<or> \<not> (0\<Colon>'a) \<le> \<bar>X3\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   399
  by (metis abs_mult_pos abs_mult 15 mult_cancel_left1 abs_idempotent abs_mult_pos mult_cancel_right1 AC_mult.f.commute 7)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   400
have 17: "\<And>X1 X3. X3 * X1 = (0\<Colon>'a) \<or> \<not> \<bar>X1\<bar> \<le> (0\<Colon>'a)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   401
  by (metis 12 16 11 11 AC_mult.f.commute)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   402
have 18: "\<And>X1 X3. X3 * X1 = (0\<Colon>'a) \<or> (0\<Colon>'a) < \<bar>X1\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   403
  by (metis 17 linorder_not_le)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   404
have 19: "\<And>X3. \<bar>h X3\<bar> < (0\<Colon>'a) \<or> \<not> c < (0\<Colon>'a) \<or> \<not> (0\<Colon>'a) < \<bar>f X3\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   405
  by (metis order_le_less_trans 0 mult_less_0_iff)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   406
show "False"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   407
  by (metis 19 linorder_antisym_conv2 5 1 13 7 6 18 abs_idempotent abs_mult_pos mult_cancel_right1 AC_mult.f.commute 7 0 12 10 15 11 linorder_not_le)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   408
qed
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   409
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   410
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   411
ML{*ResReconstruct.modulus:=1*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   412
ML{*ResReconstruct.recon_sorts:=true*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   413
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   414
lemma bigo_alt_def: "O(f) = 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   415
    {h. EX c. (0 < c & (ALL x. abs (h x) <= c * abs (f x)))}"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   416
by (auto simp add: bigo_def bigo_pos_const)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   417
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   418
ML{*ResAtp.problem_name := "BigO__bigo_elt_subset"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   419
lemma bigo_elt_subset [intro]: "f : O(g) ==> O(f) <= O(g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   420
  apply (auto simp add: bigo_alt_def)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   421
  apply (rule_tac x = "ca * c" in exI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   422
  apply (rule conjI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   423
  apply (rule mult_pos_pos)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   424
  apply (assumption)+ 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   425
(*sledgehammer*);
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   426
  apply (rule allI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   427
  apply (drule_tac x = "xa" in spec)+
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   428
  apply (subgoal_tac "ca * abs(f xa) <= ca * (c * abs(g xa))");
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   429
  apply (erule order_trans)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   430
  apply (simp add: mult_ac)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   431
  apply (rule mult_left_mono, assumption)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   432
  apply (rule order_less_imp_le, assumption);
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   433
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   434
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   435
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   436
ML{*ResAtp.problem_name := "BigO__bigo_refl"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   437
lemma bigo_refl [intro]: "f : O(f)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   438
  apply(auto simp add: bigo_def)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   439
proof (neg_clausify)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   440
fix x
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   441
assume 0: "\<And>mes_pSG\<Colon>'b\<Colon>ordered_idom.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   442
   \<not> \<bar>(f\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) ((x\<Colon>'b\<Colon>ordered_idom \<Rightarrow> 'a) mes_pSG)\<bar>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   443
     \<le> mes_pSG * \<bar>f (x mes_pSG)\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   444
have 1: "\<And>X3\<Colon>'b. X3 \<le> (1\<Colon>'b) * X3 \<or> \<not> (1\<Colon>'b) \<le> (1\<Colon>'b)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   445
  by (metis Ring_and_Field.mult_le_cancel_right1 order_refl)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   446
have 2: "\<And>X3\<Colon>'b. X3 \<le> (1\<Colon>'b) * X3"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   447
  by (metis 1 order_refl)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   448
show 3: "False"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   449
  by (metis 0 2)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   450
qed
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   451
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   452
ML{*ResAtp.problem_name := "BigO__bigo_zero"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   453
lemma bigo_zero: "0 : O(g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   454
  apply (auto simp add: bigo_def func_zero)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   455
proof (neg_clausify)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   456
fix x
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   457
assume 0: "\<And>mes_mVM\<Colon>'b\<Colon>ordered_idom.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   458
   \<not> (0\<Colon>'b\<Colon>ordered_idom)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   459
     \<le> mes_mVM *
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   460
       \<bar>(g\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   461
         ((x\<Colon>'b\<Colon>ordered_idom \<Rightarrow> 'a) mes_mVM)\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   462
have 1: "(0\<Colon>'b\<Colon>ordered_idom) < (0\<Colon>'b\<Colon>ordered_idom)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   463
  by (metis 0 Ring_and_Field.mult_le_cancel_left1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   464
show 2: "False"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   465
  by (metis Orderings.linorder_class.neq_iff 1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   466
qed
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   467
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   468
lemma bigo_zero2: "O(%x.0) = {%x.0}"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   469
  apply (auto simp add: bigo_def) 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   470
  apply (rule ext)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   471
  apply auto
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   472
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   473
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   474
lemma bigo_plus_self_subset [intro]: 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   475
  "O(f) + O(f) <= O(f)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   476
  apply (auto simp add: bigo_alt_def set_plus)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   477
  apply (rule_tac x = "c + ca" in exI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   478
  apply auto
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   479
  apply (simp add: ring_distrib func_plus)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   480
  apply (blast intro:order_trans abs_triangle_ineq add_mono elim:) 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   481
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   482
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   483
lemma bigo_plus_idemp [simp]: "O(f) + O(f) = O(f)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   484
  apply (rule equalityI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   485
  apply (rule bigo_plus_self_subset)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   486
  apply (rule set_zero_plus2) 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   487
  apply (rule bigo_zero)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   488
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   489
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   490
lemma bigo_plus_subset [intro]: "O(f + g) <= O(f) + O(g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   491
  apply (rule subsetI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   492
  apply (auto simp add: bigo_def bigo_pos_const func_plus set_plus)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   493
  apply (subst bigo_pos_const [symmetric])+
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   494
  apply (rule_tac x = 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   495
    "%n. if abs (g n) <= (abs (f n)) then x n else 0" in exI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   496
  apply (rule conjI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   497
  apply (rule_tac x = "c + c" in exI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   498
  apply (clarsimp)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   499
  apply (auto)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   500
  apply (subgoal_tac "c * abs (f xa + g xa) <= (c + c) * abs (f xa)")
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   501
  apply (erule_tac x = xa in allE)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   502
  apply (erule order_trans)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   503
  apply (simp)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   504
  apply (subgoal_tac "c * abs (f xa + g xa) <= c * (abs (f xa) + abs (g xa))")
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   505
  apply (erule order_trans)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   506
  apply (simp add: ring_distrib)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   507
  apply (rule mult_left_mono)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   508
  apply assumption
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   509
  apply (simp add: order_less_le)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   510
  apply (rule mult_left_mono)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   511
  apply (simp add: abs_triangle_ineq)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   512
  apply (simp add: order_less_le)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   513
  apply (rule mult_nonneg_nonneg)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   514
  apply (rule add_nonneg_nonneg)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   515
  apply auto
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   516
  apply (rule_tac x = "%n. if (abs (f n)) <  abs (g n) then x n else 0" 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   517
     in exI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   518
  apply (rule conjI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   519
  apply (rule_tac x = "c + c" in exI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   520
  apply auto
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   521
  apply (subgoal_tac "c * abs (f xa + g xa) <= (c + c) * abs (g xa)")
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   522
  apply (erule_tac x = xa in allE)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   523
  apply (erule order_trans)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   524
  apply (simp)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   525
  apply (subgoal_tac "c * abs (f xa + g xa) <= c * (abs (f xa) + abs (g xa))")
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   526
  apply (erule order_trans)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   527
  apply (simp add: ring_distrib)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   528
  apply (rule mult_left_mono)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   529
  apply (simp add: order_less_le)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   530
  apply (simp add: order_less_le)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   531
  apply (rule mult_left_mono)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   532
  apply (rule abs_triangle_ineq)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   533
  apply (simp add: order_less_le)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   534
  apply (rule mult_nonneg_nonneg)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   535
  apply (rule add_nonneg_nonneg)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   536
  apply (erule order_less_imp_le)+
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   537
  apply simp
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   538
  apply (rule ext)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   539
  apply (auto simp add: if_splits linorder_not_le)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   540
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   541
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   542
lemma bigo_plus_subset2 [intro]: "A <= O(f) ==> B <= O(f) ==> A + B <= O(f)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   543
  apply (subgoal_tac "A + B <= O(f) + O(f)")
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   544
  apply (erule order_trans)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   545
  apply simp
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   546
  apply (auto del: subsetI simp del: bigo_plus_idemp)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   547
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   548
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   549
ML{*ResAtp.problem_name := "BigO__bigo_plus_eq"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   550
lemma bigo_plus_eq: "ALL x. 0 <= f x ==> ALL x. 0 <= g x ==> 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   551
  O(f + g) = O(f) + O(g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   552
  apply (rule equalityI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   553
  apply (rule bigo_plus_subset)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   554
  apply (simp add: bigo_alt_def set_plus func_plus)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   555
  apply clarify 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   556
(*sledgehammer*); 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   557
  apply (rule_tac x = "max c ca" in exI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   558
  apply (rule conjI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   559
  apply (subgoal_tac "c <= max c ca")
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   560
  apply (erule order_less_le_trans)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   561
  apply assumption
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   562
  apply (rule le_maxI1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   563
  apply clarify
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   564
  apply (drule_tac x = "xa" in spec)+
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   565
  apply (subgoal_tac "0 <= f xa + g xa")
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   566
  apply (simp add: ring_distrib)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   567
  apply (subgoal_tac "abs(a xa + b xa) <= abs(a xa) + abs(b xa)")
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   568
  apply (subgoal_tac "abs(a xa) + abs(b xa) <= 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   569
      max c ca * f xa + max c ca * g xa")
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   570
  apply (blast intro: order_trans)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   571
  defer 1
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   572
  apply (rule abs_triangle_ineq)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   573
  apply (rule add_nonneg_nonneg)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   574
  apply assumption+
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   575
  apply (rule add_mono)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   576
ML{*ResAtp.problem_name := "BigO__bigo_plus_eq_simpler"*} 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   577
(*sledgehammer...fails*); 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   578
  apply (subgoal_tac "c * f xa <= max c ca * f xa")
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   579
  apply (blast intro: order_trans)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   580
  apply (rule mult_right_mono)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   581
  apply (rule le_maxI1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   582
  apply assumption
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   583
  apply (subgoal_tac "ca * g xa <= max c ca * g xa")
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   584
  apply (blast intro: order_trans)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   585
  apply (rule mult_right_mono)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   586
  apply (rule le_maxI2)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   587
  apply assumption
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   588
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   589
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   590
ML{*ResAtp.problem_name := "BigO__bigo_bounded_alt"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   591
lemma bigo_bounded_alt: "ALL x. 0 <= f x ==> ALL x. f x <= c * g x ==> 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   592
    f : O(g)" 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   593
  apply (auto simp add: bigo_def)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   594
(*Version 1: one-shot proof*)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   595
  apply (metis OrderedGroup.abs_ge_self  OrderedGroup.abs_le_D1  OrderedGroup.abs_of_nonneg  Orderings.linorder_class.not_less  order_less_le  Orderings.xt1(12)  Ring_and_Field.abs_mult)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   596
  done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   597
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   598
lemma bigo_bounded_alt: "ALL x. 0 <= f x ==> ALL x. f x <= c * g x ==> 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   599
    f : O(g)" 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   600
  apply (auto simp add: bigo_def)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   601
(*Version 2: single-step proof*)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   602
proof (neg_clausify)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   603
fix x
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   604
assume 0: "\<And>mes_mbt\<Colon>'a.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   605
   (f\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) mes_mbt
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   606
   \<le> (c\<Colon>'b\<Colon>ordered_idom) * (g\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) mes_mbt"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   607
assume 1: "\<And>mes_mbs\<Colon>'b\<Colon>ordered_idom.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   608
   \<not> (f\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) ((x\<Colon>'b\<Colon>ordered_idom \<Rightarrow> 'a) mes_mbs)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   609
     \<le> mes_mbs * \<bar>(g\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) (x mes_mbs)\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   610
have 2: "\<And>X3\<Colon>'a.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   611
   (c\<Colon>'b\<Colon>ordered_idom) * (g\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) X3 =
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   612
   (f\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) X3 \<or>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   613
   \<not> c * g X3 \<le> f X3"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   614
  by (metis Lattices.min_max.less_eq_less_inf.antisym_intro 0)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   615
have 3: "\<And>X3\<Colon>'b\<Colon>ordered_idom.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   616
   \<not> (f\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) ((x\<Colon>'b\<Colon>ordered_idom \<Rightarrow> 'a) \<bar>X3\<bar>)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   617
     \<le> \<bar>X3 * (g\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) (x \<bar>X3\<bar>)\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   618
  by (metis 1 Ring_and_Field.abs_mult)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   619
have 4: "\<And>X3\<Colon>'b\<Colon>ordered_idom. (1\<Colon>'b\<Colon>ordered_idom) * X3 = X3"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   620
  by (metis Ring_and_Field.mult_cancel_left2 Finite_Set.AC_mult.f.commute)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   621
have 5: "\<And>X3\<Colon>'b\<Colon>ordered_idom. X3 * (1\<Colon>'b\<Colon>ordered_idom) = X3"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   622
  by (metis Ring_and_Field.mult_cancel_right2 Finite_Set.AC_mult.f.commute)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   623
have 6: "\<And>X3\<Colon>'b\<Colon>ordered_idom. \<bar>X3\<bar> * \<bar>X3\<bar> = X3 * X3"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   624
  by (metis Ring_and_Field.abs_mult_self Finite_Set.AC_mult.f.commute)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   625
have 7: "\<And>X3\<Colon>'b\<Colon>ordered_idom. (0\<Colon>'b\<Colon>ordered_idom) \<le> X3 * X3"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   626
  by (metis Ring_and_Field.zero_le_square Finite_Set.AC_mult.f.commute)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   627
have 8: "(0\<Colon>'b\<Colon>ordered_idom) \<le> (1\<Colon>'b\<Colon>ordered_idom)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   628
  by (metis 7 Ring_and_Field.mult_cancel_left2)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   629
have 9: "\<And>X3\<Colon>'b\<Colon>ordered_idom. X3 * X3 = \<bar>X3 * X3\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   630
  by (metis Ring_and_Field.abs_mult 6)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   631
have 10: "\<bar>1\<Colon>'b\<Colon>ordered_idom\<bar> = (1\<Colon>'b\<Colon>ordered_idom)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   632
  by (metis 9 4)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   633
have 11: "\<And>X3\<Colon>'b\<Colon>ordered_idom. \<bar>\<bar>X3\<bar>\<bar> = \<bar>X3\<bar> * \<bar>1\<Colon>'b\<Colon>ordered_idom\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   634
  by (metis Ring_and_Field.abs_mult OrderedGroup.abs_idempotent 5)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   635
have 12: "\<And>X3\<Colon>'b\<Colon>ordered_idom. \<bar>\<bar>X3\<bar>\<bar> = \<bar>X3\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   636
  by (metis 11 10 5)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   637
have 13: "\<And>(X1\<Colon>'b\<Colon>ordered_idom) X3\<Colon>'b\<Colon>ordered_idom.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   638
   X3 * (1\<Colon>'b\<Colon>ordered_idom) \<le> X1 \<or>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   639
   \<not> \<bar>X3\<bar> \<le> X1 \<or> \<not> (0\<Colon>'b\<Colon>ordered_idom) \<le> (1\<Colon>'b\<Colon>ordered_idom)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   640
  by (metis OrderedGroup.abs_le_D1 Ring_and_Field.abs_mult_pos 5)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   641
have 14: "\<And>(X1\<Colon>'b\<Colon>ordered_idom) X3\<Colon>'b\<Colon>ordered_idom.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   642
   X3 \<le> X1 \<or> \<not> \<bar>X3\<bar> \<le> X1 \<or> \<not> (0\<Colon>'b\<Colon>ordered_idom) \<le> (1\<Colon>'b\<Colon>ordered_idom)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   643
  by (metis 13 5)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   644
have 15: "\<And>(X1\<Colon>'b\<Colon>ordered_idom) X3\<Colon>'b\<Colon>ordered_idom. X3 \<le> X1 \<or> \<not> \<bar>X3\<bar> \<le> X1"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   645
  by (metis 14 8)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   646
have 16: "\<And>(X1\<Colon>'b\<Colon>ordered_idom) X3\<Colon>'b\<Colon>ordered_idom. X3 \<le> X1 \<or> X1 \<le> \<bar>X3\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   647
  by (metis 15 Orderings.linorder_class.less_eq_less.linear)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   648
have 17: "\<And>X3\<Colon>'b\<Colon>ordered_idom.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   649
   X3 * (g\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) ((x\<Colon>'b\<Colon>ordered_idom \<Rightarrow> 'a) \<bar>X3\<bar>)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   650
   \<le> (f\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) (x \<bar>X3\<bar>)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   651
  by (metis 3 16)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   652
have 18: "(c\<Colon>'b\<Colon>ordered_idom) *
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   653
(g\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) ((x\<Colon>'b\<Colon>ordered_idom \<Rightarrow> 'a) \<bar>c\<bar>) =
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   654
(f\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) (x \<bar>c\<bar>)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   655
  by (metis 2 17)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   656
have 19: "\<And>(X1\<Colon>'b\<Colon>ordered_idom) X3\<Colon>'b\<Colon>ordered_idom. \<bar>X3 * X1\<bar> \<le> \<bar>\<bar>X3\<bar>\<bar> * \<bar>\<bar>X1\<bar>\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   657
  by (metis 15 Ring_and_Field.abs_le_mult Ring_and_Field.abs_mult)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   658
have 20: "\<And>(X1\<Colon>'b\<Colon>ordered_idom) X3\<Colon>'b\<Colon>ordered_idom. \<bar>X3 * X1\<bar> \<le> \<bar>X3\<bar> * \<bar>X1\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   659
  by (metis 19 12 12)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   660
have 21: "\<And>(X1\<Colon>'b\<Colon>ordered_idom) X3\<Colon>'b\<Colon>ordered_idom. X3 * X1 \<le> \<bar>X3\<bar> * \<bar>X1\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   661
  by (metis 15 20)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   662
have 22: "(f\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   663
 ((x\<Colon>'b\<Colon>ordered_idom \<Rightarrow> 'a) \<bar>c\<Colon>'b\<Colon>ordered_idom\<bar>)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   664
\<le> \<bar>c\<bar> * \<bar>(g\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) (x \<bar>c\<bar>)\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   665
  by (metis 21 18)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   666
show 23: "False"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   667
  by (metis 22 1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   668
qed
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   669
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   670
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   671
text{*So here is the easier (and more natural) problem using transitivity*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   672
ML{*ResAtp.problem_name := "BigO__bigo_bounded_alt_trans"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   673
lemma "ALL x. 0 <= f x ==> ALL x. f x <= c * g x ==> f : O(g)" 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   674
  apply (auto simp add: bigo_def)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   675
  (*Version 1: one-shot proof*) 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   676
apply (metis Orderings.leD Orderings.leI abs_ge_self abs_le_D1 abs_mult abs_of_nonneg order_le_less xt1(12));
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   677
  done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   678
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   679
text{*So here is the easier (and more natural) problem using transitivity*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   680
ML{*ResAtp.problem_name := "BigO__bigo_bounded_alt_trans"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   681
lemma "ALL x. 0 <= f x ==> ALL x. f x <= c * g x ==> f : O(g)" 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   682
  apply (auto simp add: bigo_def)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   683
(*Version 2: single-step proof*)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   684
proof (neg_clausify)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   685
fix x
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   686
assume 0: "\<And>mes_mb9\<Colon>'a.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   687
   (f\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) mes_mb9
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   688
   \<le> (c\<Colon>'b\<Colon>ordered_idom) * (g\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) mes_mb9"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   689
assume 1: "\<And>mes_mb8\<Colon>'b\<Colon>ordered_idom.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   690
   \<not> (f\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) ((x\<Colon>'b\<Colon>ordered_idom \<Rightarrow> 'a) mes_mb8)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   691
     \<le> mes_mb8 * \<bar>(g\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) (x mes_mb8)\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   692
have 2: "\<And>X3\<Colon>'a.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   693
   (c\<Colon>'b\<Colon>ordered_idom) * (g\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) X3 =
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   694
   (f\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) X3 \<or>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   695
   \<not> c * g X3 \<le> f X3"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   696
  by (metis Lattices.min_max.less_eq_less_inf.antisym_intro 0)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   697
have 3: "\<And>X3\<Colon>'b\<Colon>ordered_idom.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   698
   \<not> (f\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) ((x\<Colon>'b\<Colon>ordered_idom \<Rightarrow> 'a) \<bar>X3\<bar>)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   699
     \<le> \<bar>X3 * (g\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) (x \<bar>X3\<bar>)\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   700
  by (metis 1 Ring_and_Field.abs_mult)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   701
have 4: "\<And>X3\<Colon>'b\<Colon>ordered_idom. (1\<Colon>'b\<Colon>ordered_idom) * X3 = X3"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   702
  by (metis Ring_and_Field.mult_cancel_left2 Finite_Set.AC_mult.f.commute)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   703
have 5: "\<And>X3\<Colon>'b\<Colon>ordered_idom. X3 * (1\<Colon>'b\<Colon>ordered_idom) = X3"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   704
  by (metis Ring_and_Field.mult_cancel_right2 Finite_Set.AC_mult.f.commute)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   705
have 6: "\<And>X3\<Colon>'b\<Colon>ordered_idom. \<bar>X3\<bar> * \<bar>X3\<bar> = X3 * X3"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   706
  by (metis Ring_and_Field.abs_mult_self Finite_Set.AC_mult.f.commute)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   707
have 7: "\<And>X3\<Colon>'b\<Colon>ordered_idom. (0\<Colon>'b\<Colon>ordered_idom) \<le> X3 * X3"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   708
  by (metis Ring_and_Field.zero_le_square Finite_Set.AC_mult.f.commute)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   709
have 8: "(0\<Colon>'b\<Colon>ordered_idom) \<le> (1\<Colon>'b\<Colon>ordered_idom)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   710
  by (metis 7 Ring_and_Field.mult_cancel_left2)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   711
have 9: "\<And>X3\<Colon>'b\<Colon>ordered_idom. X3 * X3 = \<bar>X3 * X3\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   712
  by (metis Ring_and_Field.abs_mult 6)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   713
have 10: "\<bar>1\<Colon>'b\<Colon>ordered_idom\<bar> = (1\<Colon>'b\<Colon>ordered_idom)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   714
  by (metis 9 4)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   715
have 11: "\<And>X3\<Colon>'b\<Colon>ordered_idom. \<bar>\<bar>X3\<bar>\<bar> = \<bar>X3\<bar> * \<bar>1\<Colon>'b\<Colon>ordered_idom\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   716
  by (metis Ring_and_Field.abs_mult OrderedGroup.abs_idempotent 5)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   717
have 12: "\<And>X3\<Colon>'b\<Colon>ordered_idom. \<bar>\<bar>X3\<bar>\<bar> = \<bar>X3\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   718
  by (metis 11 10 5)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   719
have 13: "\<And>(X1\<Colon>'b\<Colon>ordered_idom) X3\<Colon>'b\<Colon>ordered_idom.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   720
   X3 * (1\<Colon>'b\<Colon>ordered_idom) \<le> X1 \<or>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   721
   \<not> \<bar>X3\<bar> \<le> X1 \<or> \<not> (0\<Colon>'b\<Colon>ordered_idom) \<le> (1\<Colon>'b\<Colon>ordered_idom)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   722
  by (metis OrderedGroup.abs_le_D1 Ring_and_Field.abs_mult_pos 5)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   723
have 14: "\<And>(X1\<Colon>'b\<Colon>ordered_idom) X3\<Colon>'b\<Colon>ordered_idom.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   724
   X3 \<le> X1 \<or> \<not> \<bar>X3\<bar> \<le> X1 \<or> \<not> (0\<Colon>'b\<Colon>ordered_idom) \<le> (1\<Colon>'b\<Colon>ordered_idom)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   725
  by (metis 13 5)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   726
have 15: "\<And>(X1\<Colon>'b\<Colon>ordered_idom) X3\<Colon>'b\<Colon>ordered_idom. X3 \<le> X1 \<or> \<not> \<bar>X3\<bar> \<le> X1"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   727
  by (metis 14 8)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   728
have 16: "\<And>(X1\<Colon>'b\<Colon>ordered_idom) X3\<Colon>'b\<Colon>ordered_idom. X3 \<le> X1 \<or> X1 \<le> \<bar>X3\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   729
  by (metis 15 Orderings.linorder_class.less_eq_less.linear)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   730
have 17: "\<And>X3\<Colon>'b\<Colon>ordered_idom.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   731
   X3 * (g\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) ((x\<Colon>'b\<Colon>ordered_idom \<Rightarrow> 'a) \<bar>X3\<bar>)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   732
   \<le> (f\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) (x \<bar>X3\<bar>)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   733
  by (metis 3 16)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   734
have 18: "(c\<Colon>'b\<Colon>ordered_idom) *
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   735
(g\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) ((x\<Colon>'b\<Colon>ordered_idom \<Rightarrow> 'a) \<bar>c\<bar>) =
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   736
(f\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) (x \<bar>c\<bar>)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   737
  by (metis 2 17)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   738
have 19: "\<And>(X1\<Colon>'b\<Colon>ordered_idom) X3\<Colon>'b\<Colon>ordered_idom. \<bar>X3 * X1\<bar> \<le> \<bar>\<bar>X3\<bar>\<bar> * \<bar>\<bar>X1\<bar>\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   739
  by (metis 15 Ring_and_Field.abs_le_mult Ring_and_Field.abs_mult)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   740
have 20: "\<And>(X1\<Colon>'b\<Colon>ordered_idom) X3\<Colon>'b\<Colon>ordered_idom. \<bar>X3 * X1\<bar> \<le> \<bar>X3\<bar> * \<bar>X1\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   741
  by (metis 19 12 12)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   742
have 21: "\<And>(X1\<Colon>'b\<Colon>ordered_idom) X3\<Colon>'b\<Colon>ordered_idom. X3 * X1 \<le> \<bar>X3\<bar> * \<bar>X1\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   743
  by (metis 15 20)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   744
have 22: "(f\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   745
 ((x\<Colon>'b\<Colon>ordered_idom \<Rightarrow> 'a) \<bar>c\<Colon>'b\<Colon>ordered_idom\<bar>)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   746
\<le> \<bar>c\<bar> * \<bar>(g\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) (x \<bar>c\<bar>)\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   747
  by (metis 21 18)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   748
show 23: "False"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   749
  by (metis 22 1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   750
qed
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   751
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   752
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   753
lemma bigo_bounded: "ALL x. 0 <= f x ==> ALL x. f x <= g x ==> 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   754
    f : O(g)" 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   755
  apply (erule bigo_bounded_alt [of f 1 g])
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   756
  apply simp
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   757
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   758
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   759
ML{*ResAtp.problem_name := "BigO__bigo_bounded2"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   760
lemma bigo_bounded2: "ALL x. lb x <= f x ==> ALL x. f x <= lb x + g x ==>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   761
    f : lb +o O(g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   762
  apply (rule set_minus_imp_plus)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   763
  apply (rule bigo_bounded)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   764
  apply (auto simp add: diff_minus func_minus func_plus)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   765
  prefer 2
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   766
  apply (drule_tac x = x in spec)+ 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   767
  apply arith (*not clear that it's provable otherwise*) 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   768
proof (neg_clausify)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   769
fix x
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   770
assume 0: "\<And>y. lb y \<le> f y"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   771
assume 1: "\<not> (0\<Colon>'b) \<le> f x + - lb x"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   772
have 2: "\<And>X3. (0\<Colon>'b) + X3 = X3"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   773
  by (metis diff_eq_eq right_minus_eq)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   774
have 3: "\<not> (0\<Colon>'b) \<le> f x - lb x"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   775
  by (metis 1 compare_rls(1))
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   776
have 4: "\<not> (0\<Colon>'b) + lb x \<le> f x"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   777
  by (metis 3 le_diff_eq)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   778
show "False"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   779
  by (metis 4 2 0)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   780
qed
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   781
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   782
ML{*ResAtp.problem_name := "BigO__bigo_abs"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   783
lemma bigo_abs: "(%x. abs(f x)) =o O(f)" 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   784
  apply (unfold bigo_def)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   785
  apply auto
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   786
proof (neg_clausify)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   787
fix x
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   788
assume 0: "!!mes_o43::'b::ordered_idom.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   789
   ~ abs ((f::'a::type => 'b::ordered_idom)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   790
           ((x::'b::ordered_idom => 'a::type) mes_o43))
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   791
     <= mes_o43 * abs (f (x mes_o43))"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   792
have 1: "!!X3::'b::ordered_idom.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   793
   X3 <= (1::'b::ordered_idom) * X3 |
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   794
   ~ (1::'b::ordered_idom) <= (1::'b::ordered_idom)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   795
  by (metis mult_le_cancel_right1 order_refl)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   796
have 2: "!!X3::'b::ordered_idom. X3 <= (1::'b::ordered_idom) * X3"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   797
  by (metis 1 order_refl)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   798
show "False"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   799
  by (metis 0 2)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   800
qed
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   801
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   802
ML{*ResAtp.problem_name := "BigO__bigo_abs2"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   803
lemma bigo_abs2: "f =o O(%x. abs(f x))"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   804
  apply (unfold bigo_def)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   805
  apply auto
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   806
proof (neg_clausify)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   807
fix x
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   808
assume 0: "\<And>mes_o4C\<Colon>'b\<Colon>ordered_idom.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   809
   \<not> \<bar>(f\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) ((x\<Colon>'b\<Colon>ordered_idom \<Rightarrow> 'a) mes_o4C)\<bar>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   810
     \<le> mes_o4C * \<bar>f (x mes_o4C)\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   811
have 1: "\<And>X3\<Colon>'b\<Colon>ordered_idom.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   812
   X3 \<le> (1\<Colon>'b\<Colon>ordered_idom) * X3 \<or>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   813
   \<not> (1\<Colon>'b\<Colon>ordered_idom) \<le> (1\<Colon>'b\<Colon>ordered_idom)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   814
  by (metis mult_le_cancel_right1 order_refl)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   815
have 2: "\<And>X3\<Colon>'b\<Colon>ordered_idom. X3 \<le> (1\<Colon>'b\<Colon>ordered_idom) * X3"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   816
  by (metis 1 order_refl)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   817
show "False"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   818
  by (metis 0 2)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   819
qed
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   820
 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   821
lemma bigo_abs3: "O(f) = O(%x. abs(f x))"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   822
  apply (rule equalityI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   823
  apply (rule bigo_elt_subset)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   824
  apply (rule bigo_abs2)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   825
  apply (rule bigo_elt_subset)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   826
  apply (rule bigo_abs)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   827
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   828
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   829
lemma bigo_abs4: "f =o g +o O(h) ==> 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   830
    (%x. abs (f x)) =o (%x. abs (g x)) +o O(h)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   831
  apply (drule set_plus_imp_minus)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   832
  apply (rule set_minus_imp_plus)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   833
  apply (subst func_diff)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   834
proof -
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   835
  assume a: "f - g : O(h)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   836
  have "(%x. abs (f x) - abs (g x)) =o O(%x. abs(abs (f x) - abs (g x)))"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   837
    by (rule bigo_abs2)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   838
  also have "... <= O(%x. abs (f x - g x))"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   839
    apply (rule bigo_elt_subset)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   840
    apply (rule bigo_bounded)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   841
    apply force
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   842
    apply (rule allI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   843
    apply (rule abs_triangle_ineq3)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   844
    done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   845
  also have "... <= O(f - g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   846
    apply (rule bigo_elt_subset)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   847
    apply (subst func_diff)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   848
    apply (rule bigo_abs)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   849
    done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   850
  also have "... <= O(h)"
23464
bc2563c37b1a tuned proofs -- avoid implicit prems;
wenzelm
parents: 23449
diff changeset
   851
    using a by (rule bigo_elt_subset)
23449
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   852
  finally show "(%x. abs (f x) - abs (g x)) : O(h)".
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   853
qed
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   854
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   855
lemma bigo_abs5: "f =o O(g) ==> (%x. abs(f x)) =o O(g)" 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   856
by (unfold bigo_def, auto)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   857
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   858
lemma bigo_elt_subset2 [intro]: "f : g +o O(h) ==> O(f) <= O(g) + O(h)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   859
proof -
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   860
  assume "f : g +o O(h)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   861
  also have "... <= O(g) + O(h)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   862
    by (auto del: subsetI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   863
  also have "... = O(%x. abs(g x)) + O(%x. abs(h x))"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   864
    apply (subst bigo_abs3 [symmetric])+
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   865
    apply (rule refl)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   866
    done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   867
  also have "... = O((%x. abs(g x)) + (%x. abs(h x)))"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   868
    by (rule bigo_plus_eq [symmetric], auto)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   869
  finally have "f : ...".
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   870
  then have "O(f) <= ..."
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   871
    by (elim bigo_elt_subset)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   872
  also have "... = O(%x. abs(g x)) + O(%x. abs(h x))"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   873
    by (rule bigo_plus_eq, auto)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   874
  finally show ?thesis
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   875
    by (simp add: bigo_abs3 [symmetric])
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   876
qed
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   877
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   878
ML{*ResAtp.problem_name := "BigO__bigo_mult"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   879
lemma bigo_mult [intro]: "O(f)*O(g) <= O(f * g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   880
  apply (rule subsetI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   881
  apply (subst bigo_def)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   882
  apply (auto simp del: abs_mult mult_ac
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   883
              simp add: bigo_alt_def set_times func_times)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   884
(*sledgehammer*); 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   885
  apply (rule_tac x = "c * ca" in exI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   886
  apply(rule allI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   887
  apply(erule_tac x = x in allE)+
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   888
  apply(subgoal_tac "c * ca * abs(f x * g x) = 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   889
      (c * abs(f x)) * (ca * abs(g x))")
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   890
ML{*ResAtp.problem_name := "BigO__bigo_mult_simpler"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   891
prefer 2 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   892
apply (metis  Finite_Set.AC_mult.f.assoc  Finite_Set.AC_mult.f.left_commute  OrderedGroup.abs_of_pos  OrderedGroup.mult_left_commute  Ring_and_Field.abs_mult  Ring_and_Field.mult_pos_pos)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   893
  apply(erule ssubst) 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   894
  apply (subst abs_mult)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   895
(*not qute BigO__bigo_mult_simpler_1 (a hard problem!) as abs_mult has
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   896
  just been done*)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   897
proof (neg_clausify)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   898
fix a c b ca x
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   899
assume 0: "(0\<Colon>'b\<Colon>ordered_idom) < (c\<Colon>'b\<Colon>ordered_idom)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   900
assume 1: "\<bar>(a\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) (x\<Colon>'a)\<bar>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   901
\<le> (c\<Colon>'b\<Colon>ordered_idom) * \<bar>(f\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) x\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   902
assume 2: "\<bar>(b\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) (x\<Colon>'a)\<bar>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   903
\<le> (ca\<Colon>'b\<Colon>ordered_idom) * \<bar>(g\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) x\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   904
assume 3: "\<not> \<bar>(a\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) (x\<Colon>'a)\<bar> *
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   905
  \<bar>(b\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) x\<bar>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   906
  \<le> (c\<Colon>'b\<Colon>ordered_idom) * \<bar>(f\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) x\<bar> *
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   907
    ((ca\<Colon>'b\<Colon>ordered_idom) * \<bar>(g\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) x\<bar>)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   908
have 4: "\<bar>c\<Colon>'b\<Colon>ordered_idom\<bar> = c"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   909
  by (metis OrderedGroup.abs_of_pos 0)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   910
have 5: "\<And>X1\<Colon>'b\<Colon>ordered_idom. (c\<Colon>'b\<Colon>ordered_idom) * \<bar>X1\<bar> = \<bar>c * X1\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   911
  by (metis Ring_and_Field.abs_mult 4)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   912
have 6: "(0\<Colon>'b\<Colon>ordered_idom) = (1\<Colon>'b\<Colon>ordered_idom) \<or>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   913
(0\<Colon>'b\<Colon>ordered_idom) < (1\<Colon>'b\<Colon>ordered_idom)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   914
  by (metis OrderedGroup.abs_not_less_zero Ring_and_Field.abs_one Ring_and_Field.linorder_neqE_ordered_idom)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   915
have 7: "(0\<Colon>'b\<Colon>ordered_idom) < (1\<Colon>'b\<Colon>ordered_idom)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   916
  by (metis 6 Ring_and_Field.one_neq_zero)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   917
have 8: "\<bar>1\<Colon>'b\<Colon>ordered_idom\<bar> = (1\<Colon>'b\<Colon>ordered_idom)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   918
  by (metis OrderedGroup.abs_of_pos 7)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   919
have 9: "\<And>X1\<Colon>'b\<Colon>ordered_idom. (0\<Colon>'b\<Colon>ordered_idom) \<le> (c\<Colon>'b\<Colon>ordered_idom) * \<bar>X1\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   920
  by (metis OrderedGroup.abs_ge_zero 5)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   921
have 10: "\<And>X1\<Colon>'b\<Colon>ordered_idom. X1 * (1\<Colon>'b\<Colon>ordered_idom) = X1"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   922
  by (metis Ring_and_Field.mult_cancel_right2 Finite_Set.AC_mult.f.commute)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   923
have 11: "\<And>X1\<Colon>'b\<Colon>ordered_idom. \<bar>\<bar>X1\<bar>\<bar> = \<bar>X1\<bar> * \<bar>1\<Colon>'b\<Colon>ordered_idom\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   924
  by (metis Ring_and_Field.abs_mult OrderedGroup.abs_idempotent 10)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   925
have 12: "\<And>X1\<Colon>'b\<Colon>ordered_idom. \<bar>\<bar>X1\<bar>\<bar> = \<bar>X1\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   926
  by (metis 11 8 10)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   927
have 13: "\<And>X1\<Colon>'b\<Colon>ordered_idom. (0\<Colon>'b\<Colon>ordered_idom) \<le> \<bar>X1\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   928
  by (metis OrderedGroup.abs_ge_zero 12)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   929
have 14: "\<not> (0\<Colon>'b\<Colon>ordered_idom)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   930
  \<le> (c\<Colon>'b\<Colon>ordered_idom) * \<bar>(f\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) (x\<Colon>'a)\<bar> \<or>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   931
\<not> (0\<Colon>'b\<Colon>ordered_idom) \<le> \<bar>(b\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) x\<bar> \<or>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   932
\<not> \<bar>b x\<bar> \<le> (ca\<Colon>'b\<Colon>ordered_idom) * \<bar>(g\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) x\<bar> \<or>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   933
\<not> \<bar>(a\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) x\<bar> \<le> c * \<bar>f x\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   934
  by (metis 3 Ring_and_Field.mult_mono)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   935
have 15: "\<not> (0\<Colon>'b\<Colon>ordered_idom) \<le> \<bar>(b\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) (x\<Colon>'a)\<bar> \<or>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   936
\<not> \<bar>b x\<bar> \<le> (ca\<Colon>'b\<Colon>ordered_idom) * \<bar>(g\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) x\<bar> \<or>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   937
\<not> \<bar>(a\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) x\<bar>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   938
  \<le> (c\<Colon>'b\<Colon>ordered_idom) * \<bar>(f\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) x\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   939
  by (metis 14 9)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   940
have 16: "\<not> \<bar>(b\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) (x\<Colon>'a)\<bar>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   941
  \<le> (ca\<Colon>'b\<Colon>ordered_idom) * \<bar>(g\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) x\<bar> \<or>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   942
\<not> \<bar>(a\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) x\<bar>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   943
  \<le> (c\<Colon>'b\<Colon>ordered_idom) * \<bar>(f\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) x\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   944
  by (metis 15 13)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   945
have 17: "\<not> \<bar>(a\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) (x\<Colon>'a)\<bar>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   946
  \<le> (c\<Colon>'b\<Colon>ordered_idom) * \<bar>(f\<Colon>'a \<Rightarrow> 'b\<Colon>ordered_idom) x\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   947
  by (metis 16 2)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   948
show 18: "False"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   949
  by (metis 17 1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   950
qed
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   951
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   952
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   953
ML{*ResAtp.problem_name := "BigO__bigo_mult2"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   954
lemma bigo_mult2 [intro]: "f *o O(g) <= O(f * g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   955
  apply (auto simp add: bigo_def elt_set_times_def func_times abs_mult)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   956
(*sledgehammer*); 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   957
  apply (rule_tac x = c in exI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   958
  apply clarify
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   959
  apply (drule_tac x = x in spec)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   960
ML{*ResAtp.problem_name := "BigO__bigo_mult2_simpler"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   961
(*sledgehammer*); 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   962
  apply (subgoal_tac "abs(f x) * abs(b x) <= abs(f x) * (c * abs(g x))")
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   963
  apply (simp add: mult_ac)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   964
  apply (rule mult_left_mono, assumption)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   965
  apply (rule abs_ge_zero)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   966
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   967
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   968
ML{*ResAtp.problem_name:="BigO__bigo_mult3"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   969
lemma bigo_mult3: "f : O(h) ==> g : O(j) ==> f * g : O(h * j)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   970
by (metis bigo_mult set_times_intro subset_iff)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   971
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   972
ML{*ResAtp.problem_name:="BigO__bigo_mult4"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   973
lemma bigo_mult4 [intro]:"f : k +o O(h) ==> g * f : (g * k) +o O(g * h)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   974
by (metis bigo_mult2 set_plus_mono_b set_times_intro2 set_times_plus_distrib)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   975
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   976
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   977
lemma bigo_mult5: "ALL x. f x ~= 0 ==>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   978
    O(f * g) <= (f::'a => ('b::ordered_field)) *o O(g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   979
proof -
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   980
  assume "ALL x. f x ~= 0"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   981
  show "O(f * g) <= f *o O(g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   982
  proof
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   983
    fix h
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   984
    assume "h : O(f * g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   985
    then have "(%x. 1 / (f x)) * h : (%x. 1 / f x) *o O(f * g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   986
      by auto
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   987
    also have "... <= O((%x. 1 / f x) * (f * g))"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   988
      by (rule bigo_mult2)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   989
    also have "(%x. 1 / f x) * (f * g) = g"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   990
      apply (simp add: func_times) 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   991
      apply (rule ext)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   992
      apply (simp add: prems nonzero_divide_eq_eq mult_ac)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   993
      done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   994
    finally have "(%x. (1::'b) / f x) * h : O(g)".
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   995
    then have "f * ((%x. (1::'b) / f x) * h) : f *o O(g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   996
      by auto
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   997
    also have "f * ((%x. (1::'b) / f x) * h) = h"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   998
      apply (simp add: func_times) 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
   999
      apply (rule ext)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1000
      apply (simp add: prems nonzero_divide_eq_eq mult_ac)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1001
      done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1002
    finally show "h : f *o O(g)".
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1003
  qed
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1004
qed
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1005
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1006
ML{*ResAtp.problem_name := "BigO__bigo_mult6"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1007
lemma bigo_mult6: "ALL x. f x ~= 0 ==>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1008
    O(f * g) = (f::'a => ('b::ordered_field)) *o O(g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1009
by (metis bigo_mult2 bigo_mult5 order_antisym)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1010
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1011
(*proof requires relaxing relevance: 2007-01-25*)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1012
ML{*ResAtp.problem_name := "BigO__bigo_mult7"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1013
  declare bigo_mult6 [simp]
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1014
lemma bigo_mult7: "ALL x. f x ~= 0 ==>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1015
    O(f * g) <= O(f::'a => ('b::ordered_field)) * O(g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1016
(*sledgehammer*)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1017
  apply (subst bigo_mult6)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1018
  apply assumption
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1019
  apply (rule set_times_mono3) 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1020
  apply (rule bigo_refl)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1021
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1022
  declare bigo_mult6 [simp del]
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1023
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1024
ML{*ResAtp.problem_name := "BigO__bigo_mult8"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1025
  declare bigo_mult7[intro!]
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1026
lemma bigo_mult8: "ALL x. f x ~= 0 ==>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1027
    O(f * g) = O(f::'a => ('b::ordered_field)) * O(g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1028
by (metis bigo_mult bigo_mult7 order_antisym_conv)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1029
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1030
lemma bigo_minus [intro]: "f : O(g) ==> - f : O(g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1031
  by (auto simp add: bigo_def func_minus)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1032
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1033
lemma bigo_minus2: "f : g +o O(h) ==> -f : -g +o O(h)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1034
  apply (rule set_minus_imp_plus)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1035
  apply (drule set_plus_imp_minus)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1036
  apply (drule bigo_minus)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1037
  apply (simp add: diff_minus)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1038
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1039
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1040
lemma bigo_minus3: "O(-f) = O(f)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1041
  by (auto simp add: bigo_def func_minus abs_minus_cancel)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1042
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1043
lemma bigo_plus_absorb_lemma1: "f : O(g) ==> f +o O(g) <= O(g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1044
proof -
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1045
  assume a: "f : O(g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1046
  show "f +o O(g) <= O(g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1047
  proof -
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1048
    have "f : O(f)" by auto
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1049
    then have "f +o O(g) <= O(f) + O(g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1050
      by (auto del: subsetI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1051
    also have "... <= O(g) + O(g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1052
    proof -
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1053
      from a have "O(f) <= O(g)" by (auto del: subsetI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1054
      thus ?thesis by (auto del: subsetI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1055
    qed
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1056
    also have "... <= O(g)" by (simp add: bigo_plus_idemp)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1057
    finally show ?thesis .
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1058
  qed
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1059
qed
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1060
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1061
lemma bigo_plus_absorb_lemma2: "f : O(g) ==> O(g) <= f +o O(g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1062
proof -
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1063
  assume a: "f : O(g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1064
  show "O(g) <= f +o O(g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1065
  proof -
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1066
    from a have "-f : O(g)" by auto
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1067
    then have "-f +o O(g) <= O(g)" by (elim bigo_plus_absorb_lemma1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1068
    then have "f +o (-f +o O(g)) <= f +o O(g)" by auto
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1069
    also have "f +o (-f +o O(g)) = O(g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1070
      by (simp add: set_plus_rearranges)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1071
    finally show ?thesis .
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1072
  qed
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1073
qed
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1074
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1075
ML{*ResAtp.problem_name:="BigO__bigo_plus_absorb"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1076
lemma bigo_plus_absorb [simp]: "f : O(g) ==> f +o O(g) = O(g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1077
by (metis bigo_plus_absorb_lemma1 bigo_plus_absorb_lemma2 order_eq_iff);
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1078
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1079
lemma bigo_plus_absorb2 [intro]: "f : O(g) ==> A <= O(g) ==> f +o A <= O(g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1080
  apply (subgoal_tac "f +o A <= f +o O(g)")
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1081
  apply force+
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1082
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1083
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1084
lemma bigo_add_commute_imp: "f : g +o O(h) ==> g : f +o O(h)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1085
  apply (subst set_minus_plus [symmetric])
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1086
  apply (subgoal_tac "g - f = - (f - g)")
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1087
  apply (erule ssubst)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1088
  apply (rule bigo_minus)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1089
  apply (subst set_minus_plus)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1090
  apply assumption
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1091
  apply  (simp add: diff_minus add_ac)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1092
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1093
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1094
lemma bigo_add_commute: "(f : g +o O(h)) = (g : f +o O(h))"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1095
  apply (rule iffI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1096
  apply (erule bigo_add_commute_imp)+
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1097
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1098
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1099
lemma bigo_const1: "(%x. c) : O(%x. 1)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1100
by (auto simp add: bigo_def mult_ac)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1101
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1102
declare bigo_const1 [skolem]
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1103
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1104
ML{*ResAtp.problem_name:="BigO__bigo_const2"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1105
lemma (*bigo_const2 [intro]:*) "O(%x. c) <= O(%x. 1)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1106
by (metis bigo_const1 bigo_elt_subset);
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1107
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1108
lemma bigo_const2 [intro]: "O(%x. c) <= O(%x. 1)"; 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1109
(*??FAILS because the two occurrences of COMBK have different polymorphic types
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1110
proof (neg_clausify)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1111
assume 0: "\<not> O(COMBK (c\<Colon>'b\<Colon>ordered_idom)) \<subseteq> O(COMBK (1\<Colon>'b\<Colon>ordered_idom))"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1112
have 1: "COMBK (c\<Colon>'b\<Colon>ordered_idom) \<notin> O(COMBK (1\<Colon>'b\<Colon>ordered_idom))"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1113
apply (rule notI) 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1114
apply (rule 0 [THEN notE]) 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1115
apply (rule bigo_elt_subset) 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1116
apply assumption; 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1117
sorry
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1118
  by (metis 0 bigo_elt_subset)  loops??
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1119
show "False"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1120
  by (metis 1 bigo_const1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1121
qed
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1122
*)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1123
  apply (rule bigo_elt_subset)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1124
  apply (rule bigo_const1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1125
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1126
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1127
declare bigo_const2 [skolem]
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1128
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1129
ML{*ResAtp.problem_name := "BigO__bigo_const3"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1130
lemma bigo_const3: "(c::'a::ordered_field) ~= 0 ==> (%x. 1) : O(%x. c)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1131
apply (simp add: bigo_def)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1132
proof (neg_clausify)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1133
assume 0: "(c\<Colon>'a\<Colon>ordered_field) \<noteq> (0\<Colon>'a\<Colon>ordered_field)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1134
assume 1: "\<And>mes_md\<Colon>'a\<Colon>ordered_field. \<not> (1\<Colon>'a\<Colon>ordered_field) \<le> mes_md * \<bar>c\<Colon>'a\<Colon>ordered_field\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1135
have 2: "(0\<Colon>'a\<Colon>ordered_field) = \<bar>c\<Colon>'a\<Colon>ordered_field\<bar> \<or>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1136
\<not> (1\<Colon>'a\<Colon>ordered_field) \<le> (1\<Colon>'a\<Colon>ordered_field)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1137
  by (metis 1 field_inverse)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1138
have 3: "\<bar>c\<Colon>'a\<Colon>ordered_field\<bar> = (0\<Colon>'a\<Colon>ordered_field)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1139
  by (metis 2 order_refl)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1140
have 4: "(0\<Colon>'a\<Colon>ordered_field) = (c\<Colon>'a\<Colon>ordered_field)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1141
  by (metis OrderedGroup.abs_eq_0 3)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1142
show 5: "False"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1143
  by (metis 4 0)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1144
qed
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1145
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1146
lemma bigo_const4: "(c::'a::ordered_field) ~= 0 ==> O(%x. 1) <= O(%x. c)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1147
by (rule bigo_elt_subset, rule bigo_const3, assumption)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1148
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1149
lemma bigo_const [simp]: "(c::'a::ordered_field) ~= 0 ==> 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1150
    O(%x. c) = O(%x. 1)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1151
by (rule equalityI, rule bigo_const2, rule bigo_const4, assumption)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1152
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1153
ML{*ResAtp.problem_name := "BigO__bigo_const_mult1"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1154
lemma bigo_const_mult1: "(%x. c * f x) : O(f)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1155
  apply (simp add: bigo_def abs_mult) 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1156
proof (neg_clausify)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1157
fix x
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1158
assume 0: "\<And>mes_vAL\<Colon>'b.
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1159
   \<not> \<bar>c\<Colon>'b\<bar> *
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1160
     \<bar>(f\<Colon>'a \<Rightarrow> 'b) ((x\<Colon>'b \<Rightarrow> 'a) mes_vAL)\<bar>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1161
     \<le> mes_vAL * \<bar>f (x mes_vAL)\<bar>"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1162
have 1: "\<And>Y\<Colon>'b. Y \<le> Y"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1163
  by (metis order_refl)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1164
show 2: "False"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1165
  by (metis 0 1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1166
qed
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1167
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1168
lemma bigo_const_mult2: "O(%x. c * f x) <= O(f)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1169
by (rule bigo_elt_subset, rule bigo_const_mult1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1170
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1171
ML{*ResAtp.problem_name := "BigO__bigo_const_mult3"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1172
lemma bigo_const_mult3: "(c::'a::ordered_field) ~= 0 ==> f : O(%x. c * f x)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1173
  apply (simp add: bigo_def)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1174
(*sledgehammer*); 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1175
  apply (rule_tac x = "abs(inverse c)" in exI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1176
  apply (simp only: abs_mult [symmetric] mult_assoc [symmetric])
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1177
apply (subst left_inverse) 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1178
apply (auto ); 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1179
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1180
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1181
lemma bigo_const_mult4: "(c::'a::ordered_field) ~= 0 ==> 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1182
    O(f) <= O(%x. c * f x)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1183
by (rule bigo_elt_subset, rule bigo_const_mult3, assumption)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1184
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1185
lemma bigo_const_mult [simp]: "(c::'a::ordered_field) ~= 0 ==> 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1186
    O(%x. c * f x) = O(f)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1187
by (rule equalityI, rule bigo_const_mult2, erule bigo_const_mult4)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1188
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1189
ML{*ResAtp.problem_name := "BigO__bigo_const_mult5"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1190
lemma bigo_const_mult5 [simp]: "(c::'a::ordered_field) ~= 0 ==> 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1191
    (%x. c) *o O(f) = O(f)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1192
  apply (auto del: subsetI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1193
  apply (rule order_trans)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1194
  apply (rule bigo_mult2)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1195
  apply (simp add: func_times)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1196
  apply (auto intro!: subsetI simp add: bigo_def elt_set_times_def func_times)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1197
  apply (rule_tac x = "%y. inverse c * x y" in exI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1198
apply (rename_tac g d) 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1199
apply safe;
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1200
apply (rule_tac [2] ext) 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1201
(*sledgehammer*); 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1202
  apply (simp_all del: mult_assoc add: mult_assoc [symmetric] abs_mult)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1203
  apply (rule_tac x = "abs (inverse c) * d" in exI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1204
  apply (rule allI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1205
  apply (subst mult_assoc)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1206
  apply (rule mult_left_mono)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1207
  apply (erule spec)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1208
apply (simp add: ); 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1209
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1210
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1211
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1212
ML{*ResAtp.problem_name := "BigO__bigo_const_mult6"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1213
lemma bigo_const_mult6 [intro]: "(%x. c) *o O(f) <= O(f)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1214
  apply (auto intro!: subsetI
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1215
    simp add: bigo_def elt_set_times_def func_times
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1216
    simp del: abs_mult mult_ac)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1217
(*sledgehammer*); 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1218
  apply (rule_tac x = "ca * (abs c)" in exI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1219
  apply (rule allI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1220
  apply (subgoal_tac "ca * abs(c) * abs(f x) = abs(c) * (ca * abs(f x))")
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1221
  apply (erule ssubst)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1222
  apply (subst abs_mult)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1223
  apply (rule mult_left_mono)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1224
  apply (erule spec)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1225
  apply simp
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1226
  apply(simp add: mult_ac)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1227
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1228
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1229
lemma bigo_const_mult7 [intro]: "f =o O(g) ==> (%x. c * f x) =o O(g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1230
proof -
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1231
  assume "f =o O(g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1232
  then have "(%x. c) * f =o (%x. c) *o O(g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1233
    by auto
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1234
  also have "(%x. c) * f = (%x. c * f x)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1235
    by (simp add: func_times)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1236
  also have "(%x. c) *o O(g) <= O(g)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1237
    by (auto del: subsetI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1238
  finally show ?thesis .
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1239
qed
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1240
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1241
lemma bigo_compose1: "f =o O(g) ==> (%x. f(k x)) =o O(%x. g(k x))"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1242
by (unfold bigo_def, auto)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1243
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1244
lemma bigo_compose2: "f =o g +o O(h) ==> (%x. f(k x)) =o (%x. g(k x)) +o 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1245
    O(%x. h(k x))"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1246
  apply (simp only: set_minus_plus [symmetric] diff_minus func_minus
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1247
      func_plus)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1248
  apply (erule bigo_compose1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1249
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1250
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1251
subsection {* Setsum *}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1252
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1253
lemma bigo_setsum_main: "ALL x. ALL y : A x. 0 <= h x y ==> 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1254
    EX c. ALL x. ALL y : A x. abs(f x y) <= c * (h x y) ==>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1255
      (%x. SUM y : A x. f x y) =o O(%x. SUM y : A x. h x y)"  
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1256
  apply (auto simp add: bigo_def)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1257
  apply (rule_tac x = "abs c" in exI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1258
  apply (subst abs_of_nonneg) back back
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1259
  apply (rule setsum_nonneg)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1260
  apply force
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1261
  apply (subst setsum_right_distrib)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1262
  apply (rule allI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1263
  apply (rule order_trans)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1264
  apply (rule setsum_abs)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1265
  apply (rule setsum_mono)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1266
apply (blast intro: order_trans mult_right_mono abs_ge_self) 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1267
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1268
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1269
ML{*ResAtp.problem_name := "BigO__bigo_setsum1"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1270
lemma bigo_setsum1: "ALL x y. 0 <= h x y ==> 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1271
    EX c. ALL x y. abs(f x y) <= c * (h x y) ==>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1272
      (%x. SUM y : A x. f x y) =o O(%x. SUM y : A x. h x y)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1273
  apply (rule bigo_setsum_main)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1274
(*sledgehammer*); 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1275
  apply force
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1276
  apply clarsimp
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1277
  apply (rule_tac x = c in exI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1278
  apply force
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1279
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1280
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1281
lemma bigo_setsum2: "ALL y. 0 <= h y ==> 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1282
    EX c. ALL y. abs(f y) <= c * (h y) ==>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1283
      (%x. SUM y : A x. f y) =o O(%x. SUM y : A x. h y)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1284
by (rule bigo_setsum1, auto)  
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1285
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1286
ML{*ResAtp.problem_name := "BigO__bigo_setsum3"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1287
lemma bigo_setsum3: "f =o O(h) ==>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1288
    (%x. SUM y : A x. (l x y) * f(k x y)) =o
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1289
      O(%x. SUM y : A x. abs(l x y * h(k x y)))"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1290
  apply (rule bigo_setsum1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1291
  apply (rule allI)+
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1292
  apply (rule abs_ge_zero)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1293
  apply (unfold bigo_def)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1294
  apply (auto simp add: abs_mult);
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1295
(*sledgehammer*); 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1296
  apply (rule_tac x = c in exI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1297
  apply (rule allI)+
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1298
  apply (subst mult_left_commute)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1299
  apply (rule mult_left_mono)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1300
  apply (erule spec)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1301
  apply (rule abs_ge_zero)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1302
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1303
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1304
lemma bigo_setsum4: "f =o g +o O(h) ==>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1305
    (%x. SUM y : A x. l x y * f(k x y)) =o
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1306
      (%x. SUM y : A x. l x y * g(k x y)) +o
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1307
        O(%x. SUM y : A x. abs(l x y * h(k x y)))"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1308
  apply (rule set_minus_imp_plus)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1309
  apply (subst func_diff)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1310
  apply (subst setsum_subtractf [symmetric])
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1311
  apply (subst right_diff_distrib [symmetric])
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1312
  apply (rule bigo_setsum3)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1313
  apply (subst func_diff [symmetric])
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1314
  apply (erule set_plus_imp_minus)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1315
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1316
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1317
ML{*ResAtp.problem_name := "BigO__bigo_setsum5"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1318
lemma bigo_setsum5: "f =o O(h) ==> ALL x y. 0 <= l x y ==> 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1319
    ALL x. 0 <= h x ==>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1320
      (%x. SUM y : A x. (l x y) * f(k x y)) =o
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1321
        O(%x. SUM y : A x. (l x y) * h(k x y))" 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1322
  apply (subgoal_tac "(%x. SUM y : A x. (l x y) * h(k x y)) = 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1323
      (%x. SUM y : A x. abs((l x y) * h(k x y)))")
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1324
  apply (erule ssubst)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1325
  apply (erule bigo_setsum3)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1326
  apply (rule ext)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1327
  apply (rule setsum_cong2)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1328
  apply (thin_tac "f \<in> O(h)") 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1329
(*sledgehammer*); 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1330
  apply (subst abs_of_nonneg)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1331
  apply (rule mult_nonneg_nonneg)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1332
  apply auto
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1333
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1334
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1335
lemma bigo_setsum6: "f =o g +o O(h) ==> ALL x y. 0 <= l x y ==>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1336
    ALL x. 0 <= h x ==>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1337
      (%x. SUM y : A x. (l x y) * f(k x y)) =o
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1338
        (%x. SUM y : A x. (l x y) * g(k x y)) +o
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1339
          O(%x. SUM y : A x. (l x y) * h(k x y))" 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1340
  apply (rule set_minus_imp_plus)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1341
  apply (subst func_diff)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1342
  apply (subst setsum_subtractf [symmetric])
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1343
  apply (subst right_diff_distrib [symmetric])
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1344
  apply (rule bigo_setsum5)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1345
  apply (subst func_diff [symmetric])
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1346
  apply (drule set_plus_imp_minus)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1347
  apply auto
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1348
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1349
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1350
subsection {* Misc useful stuff *}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1351
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1352
lemma bigo_useful_intro: "A <= O(f) ==> B <= O(f) ==>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1353
  A + B <= O(f)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1354
  apply (subst bigo_plus_idemp [symmetric])
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1355
  apply (rule set_plus_mono2)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1356
  apply assumption+
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1357
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1358
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1359
lemma bigo_useful_add: "f =o O(h) ==> g =o O(h) ==> f + g =o O(h)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1360
  apply (subst bigo_plus_idemp [symmetric])
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1361
  apply (rule set_plus_intro)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1362
  apply assumption+
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1363
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1364
  
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1365
lemma bigo_useful_const_mult: "(c::'a::ordered_field) ~= 0 ==> 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1366
    (%x. c) * f =o O(h) ==> f =o O(h)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1367
  apply (rule subsetD)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1368
  apply (subgoal_tac "(%x. 1 / c) *o O(h) <= O(h)")
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1369
  apply assumption
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1370
  apply (rule bigo_const_mult6)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1371
  apply (subgoal_tac "f = (%x. 1 / c) * ((%x. c) * f)")
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1372
  apply (erule ssubst)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1373
  apply (erule set_times_intro2)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1374
  apply (simp add: func_times) 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1375
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1376
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1377
ML{*ResAtp.problem_name := "BigO__bigo_fix"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1378
lemma bigo_fix: "(%x. f ((x::nat) + 1)) =o O(%x. h(x + 1)) ==> f 0 = 0 ==>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1379
    f =o O(h)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1380
  apply (simp add: bigo_alt_def)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1381
(*sledgehammer*); 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1382
  apply clarify
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1383
  apply (rule_tac x = c in exI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1384
  apply safe
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1385
  apply (case_tac "x = 0")
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1386
prefer 2
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1387
  apply (subgoal_tac "x = Suc (x - 1)")
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1388
  apply (erule ssubst) back
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1389
  apply (erule spec)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1390
  apply (rule Suc_pred') 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1391
  apply simp
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1392
apply (metis OrderedGroup.abs_ge_zero  OrderedGroup.abs_zero  order_less_le  Ring_and_Field.split_mult_pos_le) 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1393
  done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1394
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1395
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1396
lemma bigo_fix2: 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1397
    "(%x. f ((x::nat) + 1)) =o (%x. g(x + 1)) +o O(%x. h(x + 1)) ==> 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1398
       f 0 = g 0 ==> f =o g +o O(h)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1399
  apply (rule set_minus_imp_plus)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1400
  apply (rule bigo_fix)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1401
  apply (subst func_diff)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1402
  apply (subst func_diff [symmetric])
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1403
  apply (rule set_plus_imp_minus)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1404
  apply simp
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1405
  apply (simp add: func_diff)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1406
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1407
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1408
subsection {* Less than or equal to *}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1409
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1410
constdefs 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1411
  lesso :: "('a => 'b::ordered_idom) => ('a => 'b) => ('a => 'b)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1412
      (infixl "<o" 70)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1413
  "f <o g == (%x. max (f x - g x) 0)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1414
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1415
lemma bigo_lesseq1: "f =o O(h) ==> ALL x. abs (g x) <= abs (f x) ==>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1416
    g =o O(h)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1417
  apply (unfold bigo_def)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1418
  apply clarsimp
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1419
apply (blast intro: order_trans) 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1420
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1421
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1422
lemma bigo_lesseq2: "f =o O(h) ==> ALL x. abs (g x) <= f x ==>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1423
      g =o O(h)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1424
  apply (erule bigo_lesseq1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1425
apply (blast intro: abs_ge_self order_trans) 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1426
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1427
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1428
lemma bigo_lesseq3: "f =o O(h) ==> ALL x. 0 <= g x ==> ALL x. g x <= f x ==>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1429
      g =o O(h)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1430
  apply (erule bigo_lesseq2)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1431
  apply (rule allI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1432
  apply (subst abs_of_nonneg)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1433
  apply (erule spec)+
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1434
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1435
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1436
lemma bigo_lesseq4: "f =o O(h) ==>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1437
    ALL x. 0 <= g x ==> ALL x. g x <= abs (f x) ==>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1438
      g =o O(h)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1439
  apply (erule bigo_lesseq1)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1440
  apply (rule allI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1441
  apply (subst abs_of_nonneg)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1442
  apply (erule spec)+
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1443
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1444
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1445
ML{*ResAtp.problem_name:="BigO__bigo_lesso1"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1446
lemma bigo_lesso1: "ALL x. f x <= g x ==> f <o g =o O(h)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1447
  apply (unfold lesso_def)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1448
  apply (subgoal_tac "(%x. max (f x - g x) 0) = 0")
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1449
(*
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1450
?? abstractions don't work: abstraction function gets the wrong type?
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1451
proof (neg_clausify)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1452
assume 0: "llabs_subgoal_1 f g = 0"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1453
assume 1: "llabs_subgoal_1 f g \<notin> O(h)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1454
show "False"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1455
  by (metis 1 0 bigo_zero)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1456
*)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1457
  apply (erule ssubst)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1458
  apply (rule bigo_zero)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1459
  apply (unfold func_zero)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1460
  apply (rule ext)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1461
  apply (simp split: split_max)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1462
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1463
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1464
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1465
ML{*ResAtp.problem_name := "BigO__bigo_lesso2"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1466
lemma bigo_lesso2: "f =o g +o O(h) ==>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1467
    ALL x. 0 <= k x ==> ALL x. k x <= f x ==>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1468
      k <o g =o O(h)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1469
  apply (unfold lesso_def)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1470
  apply (rule bigo_lesseq4)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1471
  apply (erule set_plus_imp_minus)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1472
  apply (rule allI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1473
  apply (rule le_maxI2)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1474
  apply (rule allI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1475
  apply (subst func_diff)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1476
apply (erule thin_rl)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1477
(*sledgehammer*);  
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1478
  apply (case_tac "0 <= k x - g x")
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1479
  apply (simp del: compare_rls diff_minus);
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1480
  apply (subst abs_of_nonneg)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1481
  apply (drule_tac x = x in spec) back
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1482
ML{*ResAtp.problem_name := "BigO__bigo_lesso2_simpler"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1483
(*sledgehammer*);  
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1484
  apply (simp add: compare_rls del: diff_minus)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1485
  apply (subst diff_minus)+
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1486
  apply (rule add_right_mono)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1487
  apply (erule spec)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1488
  apply (rule order_trans) 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1489
  prefer 2
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1490
  apply (rule abs_ge_zero)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1491
(*
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1492
  apply (simp only: compare_rls min_max.below_sup.above_sup_conv 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1493
             linorder_not_le order_less_imp_le)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1494
*)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1495
  apply (simp add: compare_rls del: diff_minus)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1496
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1497
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1498
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1499
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1500
ML{*ResAtp.problem_name := "BigO__bigo_lesso3"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1501
lemma bigo_lesso3: "f =o g +o O(h) ==>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1502
    ALL x. 0 <= k x ==> ALL x. g x <= k x ==>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1503
      f <o k =o O(h)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1504
  apply (unfold lesso_def)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1505
  apply (rule bigo_lesseq4)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1506
  apply (erule set_plus_imp_minus)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1507
  apply (rule allI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1508
  apply (rule le_maxI2)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1509
  apply (rule allI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1510
  apply (subst func_diff)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1511
apply (erule thin_rl) 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1512
(*sledgehammer*); 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1513
  apply (case_tac "0 <= f x - k x")
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1514
  apply (simp del: compare_rls diff_minus);
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1515
  apply (subst abs_of_nonneg)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1516
  apply (drule_tac x = x in spec) back
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1517
ML{*ResAtp.problem_name := "BigO__bigo_lesso3_simpler"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1518
(*sledgehammer*);  
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1519
  apply (simp del: diff_minus)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1520
  apply (subst diff_minus)+
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1521
  apply (rule add_left_mono)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1522
  apply (rule le_imp_neg_le)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1523
  apply (erule spec)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1524
  apply (rule order_trans) 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1525
  prefer 2
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1526
  apply (rule abs_ge_zero)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1527
  apply (simp del: diff_minus)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1528
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1529
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1530
lemma bigo_lesso4: "f <o g =o O(k::'a=>'b::ordered_field) ==>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1531
    g =o h +o O(k) ==> f <o h =o O(k)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1532
  apply (unfold lesso_def)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1533
  apply (drule set_plus_imp_minus)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1534
  apply (drule bigo_abs5) back
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1535
  apply (simp add: func_diff)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1536
  apply (drule bigo_useful_add)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1537
  apply assumption
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1538
  apply (erule bigo_lesseq2) back
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1539
  apply (rule allI)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1540
  apply (auto simp add: func_plus func_diff compare_rls 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1541
    split: split_max abs_split)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1542
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1543
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1544
ML{*ResAtp.problem_name := "BigO__bigo_lesso5"*}
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1545
lemma bigo_lesso5: "f <o g =o O(h) ==>
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1546
    EX C. ALL x. f x <= g x + C * abs(h x)"
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1547
  apply (simp only: lesso_def bigo_alt_def)
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1548
  apply clarsimp
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1549
(*sledgehammer*);  
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1550
apply (auto simp add: compare_rls add_ac) 
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1551
done
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1552
dd874e6a3282 integration of Metis prover
paulson
parents:
diff changeset
  1553
end