src/HOL/Wellfounded.thy
author krauss
Mon Oct 26 23:26:57 2009 +0100 (2009-10-26)
changeset 33216 7c61bc5d7310
parent 33215 6fd85372981e
child 33217 ab979f6e99f4
permissions -rw-r--r--
point-free characterization of well-foundedness
wenzelm@32960
     1
(*  Title:      HOL/Wellfounded.thy
wenzelm@32960
     2
    Author:     Tobias Nipkow
wenzelm@32960
     3
    Author:     Lawrence C Paulson
wenzelm@32960
     4
    Author:     Konrad Slind
wenzelm@32960
     5
    Author:     Alexander Krauss
krauss@26748
     6
*)
krauss@26748
     7
krauss@26748
     8
header {*Well-founded Recursion*}
krauss@26748
     9
krauss@26748
    10
theory Wellfounded
haftmann@30989
    11
imports Finite_Set Transitive_Closure
haftmann@31775
    12
uses ("Tools/Function/size.ML")
krauss@26748
    13
begin
krauss@26748
    14
krauss@26976
    15
subsection {* Basic Definitions *}
krauss@26976
    16
krauss@26748
    17
constdefs
krauss@26748
    18
  wf         :: "('a * 'a)set => bool"
krauss@26748
    19
  "wf(r) == (!P. (!x. (!y. (y,x):r --> P(y)) --> P(x)) --> (!x. P(x)))"
krauss@26748
    20
krauss@26748
    21
  wfP :: "('a => 'a => bool) => bool"
krauss@26748
    22
  "wfP r == wf {(x, y). r x y}"
krauss@26748
    23
krauss@26748
    24
  acyclic :: "('a*'a)set => bool"
krauss@26748
    25
  "acyclic r == !x. (x,x) ~: r^+"
krauss@26748
    26
krauss@26748
    27
abbreviation acyclicP :: "('a => 'a => bool) => bool" where
krauss@26748
    28
  "acyclicP r == acyclic {(x, y). r x y}"
krauss@26748
    29
krauss@26748
    30
lemma wfP_wf_eq [pred_set_conv]: "wfP (\<lambda>x y. (x, y) \<in> r) = wf r"
krauss@26748
    31
  by (simp add: wfP_def)
krauss@26748
    32
krauss@26748
    33
lemma wfUNIVI: 
krauss@26748
    34
   "(!!P x. (ALL x. (ALL y. (y,x) : r --> P(y)) --> P(x)) ==> P(x)) ==> wf(r)"
krauss@26748
    35
  unfolding wf_def by blast
krauss@26748
    36
krauss@26748
    37
lemmas wfPUNIVI = wfUNIVI [to_pred]
krauss@26748
    38
krauss@26748
    39
text{*Restriction to domain @{term A} and range @{term B}.  If @{term r} is
krauss@26748
    40
    well-founded over their intersection, then @{term "wf r"}*}
krauss@26748
    41
lemma wfI: 
krauss@26748
    42
 "[| r \<subseteq> A <*> B; 
krauss@26748
    43
     !!x P. [|\<forall>x. (\<forall>y. (y,x) : r --> P y) --> P x;  x : A; x : B |] ==> P x |]
krauss@26748
    44
  ==>  wf r"
krauss@26748
    45
  unfolding wf_def by blast
krauss@26748
    46
krauss@26748
    47
lemma wf_induct: 
krauss@26748
    48
    "[| wf(r);           
krauss@26748
    49
        !!x.[| ALL y. (y,x): r --> P(y) |] ==> P(x)  
krauss@26748
    50
     |]  ==>  P(a)"
krauss@26748
    51
  unfolding wf_def by blast
krauss@26748
    52
krauss@26748
    53
lemmas wfP_induct = wf_induct [to_pred]
krauss@26748
    54
krauss@26748
    55
lemmas wf_induct_rule = wf_induct [rule_format, consumes 1, case_names less, induct set: wf]
krauss@26748
    56
krauss@26748
    57
lemmas wfP_induct_rule = wf_induct_rule [to_pred, induct set: wfP]
krauss@26748
    58
krauss@26748
    59
lemma wf_not_sym: "wf r ==> (a, x) : r ==> (x, a) ~: r"
krauss@26748
    60
  by (induct a arbitrary: x set: wf) blast
krauss@26748
    61
krauss@33215
    62
lemma wf_asym:
krauss@33215
    63
  assumes "wf r" "(a, x) \<in> r"
krauss@33215
    64
  obtains "(x, a) \<notin> r"
krauss@33215
    65
  by (drule wf_not_sym[OF assms])
krauss@26748
    66
krauss@26748
    67
lemma wf_not_refl [simp]: "wf r ==> (a, a) ~: r"
krauss@26748
    68
  by (blast elim: wf_asym)
krauss@26748
    69
krauss@33215
    70
lemma wf_irrefl: assumes "wf r" obtains "(a, a) \<notin> r"
krauss@33215
    71
by (drule wf_not_refl[OF assms])
krauss@26748
    72
haftmann@27823
    73
lemma wf_wellorderI:
haftmann@27823
    74
  assumes wf: "wf {(x::'a::ord, y). x < y}"
haftmann@27823
    75
  assumes lin: "OFCLASS('a::ord, linorder_class)"
haftmann@27823
    76
  shows "OFCLASS('a::ord, wellorder_class)"
haftmann@27823
    77
using lin by (rule wellorder_class.intro)
haftmann@27823
    78
  (blast intro: wellorder_axioms.intro wf_induct_rule [OF wf])
haftmann@27823
    79
haftmann@27823
    80
lemma (in wellorder) wf:
haftmann@27823
    81
  "wf {(x, y). x < y}"
haftmann@27823
    82
unfolding wf_def by (blast intro: less_induct)
haftmann@27823
    83
haftmann@27823
    84
krauss@26976
    85
subsection {* Basic Results *}
krauss@26976
    86
krauss@33216
    87
text {* Point-free characterization of well-foundedness *}
krauss@33216
    88
krauss@33216
    89
lemma wfE_pf:
krauss@33216
    90
  assumes wf: "wf R"
krauss@33216
    91
  assumes a: "A \<subseteq> R `` A"
krauss@33216
    92
  shows "A = {}"
krauss@33216
    93
proof -
krauss@33216
    94
  { fix x
krauss@33216
    95
    from wf have "x \<notin> A"
krauss@33216
    96
    proof induct
krauss@33216
    97
      fix x assume "\<And>y. (y, x) \<in> R \<Longrightarrow> y \<notin> A"
krauss@33216
    98
      then have "x \<notin> R `` A" by blast
krauss@33216
    99
      with a show "x \<notin> A" by blast
krauss@33216
   100
    qed
krauss@33216
   101
  } thus ?thesis by auto
krauss@33216
   102
qed
krauss@33216
   103
krauss@33216
   104
lemma wfI_pf:
krauss@33216
   105
  assumes a: "\<And>A. A \<subseteq> R `` A \<Longrightarrow> A = {}"
krauss@33216
   106
  shows "wf R"
krauss@33216
   107
proof (rule wfUNIVI)
krauss@33216
   108
  fix P :: "'a \<Rightarrow> bool" and x
krauss@33216
   109
  let ?A = "{x. \<not> P x}"
krauss@33216
   110
  assume "\<forall>x. (\<forall>y. (y, x) \<in> R \<longrightarrow> P y) \<longrightarrow> P x"
krauss@33216
   111
  then have "?A \<subseteq> R `` ?A" by blast
krauss@33216
   112
  with a show "P x" by blast
krauss@33216
   113
qed
krauss@33216
   114
krauss@33216
   115
text{*Minimal-element characterization of well-foundedness*}
krauss@33216
   116
krauss@33216
   117
lemma wfE_min:
krauss@33216
   118
  assumes wf: "wf R" and Q: "x \<in> Q"
krauss@33216
   119
  obtains z where "z \<in> Q" "\<And>y. (y, z) \<in> R \<Longrightarrow> y \<notin> Q"
krauss@33216
   120
  using Q wfE_pf[OF wf, of Q] by blast
krauss@33216
   121
krauss@33216
   122
lemma wfI_min:
krauss@33216
   123
  assumes a: "\<And>x Q. x \<in> Q \<Longrightarrow> \<exists>z\<in>Q. \<forall>y. (y, z) \<in> R \<longrightarrow> y \<notin> Q"
krauss@33216
   124
  shows "wf R"
krauss@33216
   125
proof (rule wfI_pf)
krauss@33216
   126
  fix A assume b: "A \<subseteq> R `` A"
krauss@33216
   127
  { fix x assume "x \<in> A"
krauss@33216
   128
    from a[OF this] b have "False" by blast
krauss@33216
   129
  }
krauss@33216
   130
  thus "A = {}" by blast
krauss@33216
   131
qed
krauss@33216
   132
krauss@33216
   133
lemma wf_eq_minimal: "wf r = (\<forall>Q x. x\<in>Q --> (\<exists>z\<in>Q. \<forall>y. (y,z)\<in>r --> y\<notin>Q))"
krauss@33216
   134
apply auto
krauss@33216
   135
apply (erule wfE_min, assumption, blast)
krauss@33216
   136
apply (rule wfI_min, auto)
krauss@33216
   137
done
krauss@33216
   138
krauss@33216
   139
lemmas wfP_eq_minimal = wf_eq_minimal [to_pred]
krauss@33216
   140
krauss@33216
   141
text{* Well-foundedness of transitive closure *}
krauss@33216
   142
krauss@26748
   143
lemma wf_trancl:
krauss@26748
   144
  assumes "wf r"
krauss@26748
   145
  shows "wf (r^+)"
krauss@26748
   146
proof -
krauss@26748
   147
  {
krauss@26748
   148
    fix P and x
krauss@26748
   149
    assume induct_step: "!!x. (!!y. (y, x) : r^+ ==> P y) ==> P x"
krauss@26748
   150
    have "P x"
krauss@26748
   151
    proof (rule induct_step)
krauss@26748
   152
      fix y assume "(y, x) : r^+"
krauss@26748
   153
      with `wf r` show "P y"
krauss@26748
   154
      proof (induct x arbitrary: y)
wenzelm@32960
   155
        case (less x)
wenzelm@32960
   156
        note hyp = `\<And>x' y'. (x', x) : r ==> (y', x') : r^+ ==> P y'`
wenzelm@32960
   157
        from `(y, x) : r^+` show "P y"
wenzelm@32960
   158
        proof cases
wenzelm@32960
   159
          case base
wenzelm@32960
   160
          show "P y"
wenzelm@32960
   161
          proof (rule induct_step)
wenzelm@32960
   162
            fix y' assume "(y', y) : r^+"
wenzelm@32960
   163
            with `(y, x) : r` show "P y'" by (rule hyp [of y y'])
wenzelm@32960
   164
          qed
wenzelm@32960
   165
        next
wenzelm@32960
   166
          case step
wenzelm@32960
   167
          then obtain x' where "(x', x) : r" and "(y, x') : r^+" by simp
wenzelm@32960
   168
          then show "P y" by (rule hyp [of x' y])
wenzelm@32960
   169
        qed
krauss@26748
   170
      qed
krauss@26748
   171
    qed
krauss@26748
   172
  } then show ?thesis unfolding wf_def by blast
krauss@26748
   173
qed
krauss@26748
   174
krauss@26748
   175
lemmas wfP_trancl = wf_trancl [to_pred]
krauss@26748
   176
krauss@26748
   177
lemma wf_converse_trancl: "wf (r^-1) ==> wf ((r^+)^-1)"
krauss@26748
   178
  apply (subst trancl_converse [symmetric])
krauss@26748
   179
  apply (erule wf_trancl)
krauss@26748
   180
  done
krauss@26748
   181
krauss@33216
   182
text {* Well-foundedness of subsets *}
krauss@26748
   183
krauss@26748
   184
lemma wf_subset: "[| wf(r);  p<=r |] ==> wf(p)"
krauss@26748
   185
  apply (simp (no_asm_use) add: wf_eq_minimal)
krauss@26748
   186
  apply fast
krauss@26748
   187
  done
krauss@26748
   188
krauss@26748
   189
lemmas wfP_subset = wf_subset [to_pred]
krauss@26748
   190
krauss@26748
   191
text {* Well-foundedness of the empty relation *}
krauss@33216
   192
krauss@33216
   193
lemma wf_empty [iff]: "wf {}"
krauss@26748
   194
  by (simp add: wf_def)
krauss@26748
   195
haftmann@32205
   196
lemma wfP_empty [iff]:
haftmann@32205
   197
  "wfP (\<lambda>x y. False)"
haftmann@32205
   198
proof -
haftmann@32205
   199
  have "wfP bot" by (fact wf_empty [to_pred bot_empty_eq2])
haftmann@32205
   200
  then show ?thesis by (simp add: bot_fun_eq bot_bool_eq)
haftmann@32205
   201
qed
krauss@26748
   202
krauss@26748
   203
lemma wf_Int1: "wf r ==> wf (r Int r')"
krauss@26748
   204
  apply (erule wf_subset)
krauss@26748
   205
  apply (rule Int_lower1)
krauss@26748
   206
  done
krauss@26748
   207
krauss@26748
   208
lemma wf_Int2: "wf r ==> wf (r' Int r)"
krauss@26748
   209
  apply (erule wf_subset)
krauss@26748
   210
  apply (rule Int_lower2)
krauss@26748
   211
  done  
krauss@26748
   212
krauss@33216
   213
text {* Exponentiation *}
krauss@33216
   214
krauss@33216
   215
lemma wf_exp:
krauss@33216
   216
  assumes "wf (R ^^ n)"
krauss@33216
   217
  shows "wf R"
krauss@33216
   218
proof (rule wfI_pf)
krauss@33216
   219
  fix A assume "A \<subseteq> R `` A"
krauss@33216
   220
  then have "A \<subseteq> (R ^^ n) `` A" by (induct n) force+
krauss@33216
   221
  with `wf (R ^^ n)`
krauss@33216
   222
  show "A = {}" by (rule wfE_pf)
krauss@33216
   223
qed
krauss@33216
   224
krauss@33216
   225
text {* Well-foundedness of insert *}
krauss@33216
   226
krauss@26748
   227
lemma wf_insert [iff]: "wf(insert (y,x) r) = (wf(r) & (x,y) ~: r^*)"
krauss@26748
   228
apply (rule iffI)
krauss@26748
   229
 apply (blast elim: wf_trancl [THEN wf_irrefl]
krauss@26748
   230
              intro: rtrancl_into_trancl1 wf_subset 
krauss@26748
   231
                     rtrancl_mono [THEN [2] rev_subsetD])
krauss@26748
   232
apply (simp add: wf_eq_minimal, safe)
krauss@26748
   233
apply (rule allE, assumption, erule impE, blast) 
krauss@26748
   234
apply (erule bexE)
krauss@26748
   235
apply (rename_tac "a", case_tac "a = x")
krauss@26748
   236
 prefer 2
krauss@26748
   237
apply blast 
krauss@26748
   238
apply (case_tac "y:Q")
krauss@26748
   239
 prefer 2 apply blast
krauss@26748
   240
apply (rule_tac x = "{z. z:Q & (z,y) : r^*}" in allE)
krauss@26748
   241
 apply assumption
krauss@26748
   242
apply (erule_tac V = "ALL Q. (EX x. x : Q) --> ?P Q" in thin_rl) 
krauss@26748
   243
  --{*essential for speed*}
krauss@26748
   244
txt{*Blast with new substOccur fails*}
krauss@26748
   245
apply (fast intro: converse_rtrancl_into_rtrancl)
krauss@26748
   246
done
krauss@26748
   247
krauss@26748
   248
text{*Well-foundedness of image*}
krauss@33216
   249
krauss@26748
   250
lemma wf_prod_fun_image: "[| wf r; inj f |] ==> wf(prod_fun f f ` r)"
krauss@26748
   251
apply (simp only: wf_eq_minimal, clarify)
krauss@26748
   252
apply (case_tac "EX p. f p : Q")
krauss@26748
   253
apply (erule_tac x = "{p. f p : Q}" in allE)
krauss@26748
   254
apply (fast dest: inj_onD, blast)
krauss@26748
   255
done
krauss@26748
   256
krauss@26748
   257
krauss@26976
   258
subsection {* Well-Foundedness Results for Unions *}
krauss@26748
   259
krauss@26748
   260
lemma wf_union_compatible:
krauss@26748
   261
  assumes "wf R" "wf S"
krauss@32235
   262
  assumes "R O S \<subseteq> R"
krauss@26748
   263
  shows "wf (R \<union> S)"
krauss@26748
   264
proof (rule wfI_min)
krauss@26748
   265
  fix x :: 'a and Q 
krauss@26748
   266
  let ?Q' = "{x \<in> Q. \<forall>y. (y, x) \<in> R \<longrightarrow> y \<notin> Q}"
krauss@26748
   267
  assume "x \<in> Q"
krauss@26748
   268
  obtain a where "a \<in> ?Q'"
krauss@26748
   269
    by (rule wfE_min [OF `wf R` `x \<in> Q`]) blast
krauss@26748
   270
  with `wf S`
krauss@26748
   271
  obtain z where "z \<in> ?Q'" and zmin: "\<And>y. (y, z) \<in> S \<Longrightarrow> y \<notin> ?Q'" by (erule wfE_min)
krauss@26748
   272
  { 
krauss@26748
   273
    fix y assume "(y, z) \<in> S"
krauss@26748
   274
    then have "y \<notin> ?Q'" by (rule zmin)
krauss@26748
   275
krauss@26748
   276
    have "y \<notin> Q"
krauss@26748
   277
    proof 
krauss@26748
   278
      assume "y \<in> Q"
krauss@26748
   279
      with `y \<notin> ?Q'` 
krauss@26748
   280
      obtain w where "(w, y) \<in> R" and "w \<in> Q" by auto
krauss@32235
   281
      from `(w, y) \<in> R` `(y, z) \<in> S` have "(w, z) \<in> R O S" by (rule rel_compI)
krauss@32235
   282
      with `R O S \<subseteq> R` have "(w, z) \<in> R" ..
krauss@26748
   283
      with `z \<in> ?Q'` have "w \<notin> Q" by blast 
krauss@26748
   284
      with `w \<in> Q` show False by contradiction
krauss@26748
   285
    qed
krauss@26748
   286
  }
krauss@26748
   287
  with `z \<in> ?Q'` show "\<exists>z\<in>Q. \<forall>y. (y, z) \<in> R \<union> S \<longrightarrow> y \<notin> Q" by blast
krauss@26748
   288
qed
krauss@26748
   289
krauss@26748
   290
krauss@26748
   291
text {* Well-foundedness of indexed union with disjoint domains and ranges *}
krauss@26748
   292
krauss@26748
   293
lemma wf_UN: "[| ALL i:I. wf(r i);  
krauss@26748
   294
         ALL i:I. ALL j:I. r i ~= r j --> Domain(r i) Int Range(r j) = {}  
krauss@26748
   295
      |] ==> wf(UN i:I. r i)"
krauss@26748
   296
apply (simp only: wf_eq_minimal, clarify)
krauss@26748
   297
apply (rename_tac A a, case_tac "EX i:I. EX a:A. EX b:A. (b,a) : r i")
krauss@26748
   298
 prefer 2
krauss@26748
   299
 apply force 
krauss@26748
   300
apply clarify
krauss@26748
   301
apply (drule bspec, assumption)  
krauss@26748
   302
apply (erule_tac x="{a. a:A & (EX b:A. (b,a) : r i) }" in allE)
krauss@26748
   303
apply (blast elim!: allE)  
krauss@26748
   304
done
krauss@26748
   305
haftmann@32263
   306
lemma wfP_SUP:
haftmann@32263
   307
  "\<forall>i. wfP (r i) \<Longrightarrow> \<forall>i j. r i \<noteq> r j \<longrightarrow> inf (DomainP (r i)) (RangeP (r j)) = bot \<Longrightarrow> wfP (SUPR UNIV r)"
haftmann@32704
   308
  by (rule wf_UN [where I=UNIV and r="\<lambda>i. {(x, y). r i x y}", to_pred SUP_UN_eq2])
haftmann@32704
   309
    (simp_all add: Collect_def)
krauss@26748
   310
krauss@26748
   311
lemma wf_Union: 
krauss@26748
   312
 "[| ALL r:R. wf r;  
krauss@26748
   313
     ALL r:R. ALL s:R. r ~= s --> Domain r Int Range s = {}  
krauss@26748
   314
  |] ==> wf(Union R)"
krauss@26748
   315
apply (simp add: Union_def)
krauss@26748
   316
apply (blast intro: wf_UN)
krauss@26748
   317
done
krauss@26748
   318
krauss@26748
   319
(*Intuition: we find an (R u S)-min element of a nonempty subset A
krauss@26748
   320
             by case distinction.
krauss@26748
   321
  1. There is a step a -R-> b with a,b : A.
krauss@26748
   322
     Pick an R-min element z of the (nonempty) set {a:A | EX b:A. a -R-> b}.
krauss@26748
   323
     By definition, there is z':A s.t. z -R-> z'. Because z is R-min in the
krauss@26748
   324
     subset, z' must be R-min in A. Because z' has an R-predecessor, it cannot
krauss@26748
   325
     have an S-successor and is thus S-min in A as well.
krauss@26748
   326
  2. There is no such step.
krauss@26748
   327
     Pick an S-min element of A. In this case it must be an R-min
krauss@26748
   328
     element of A as well.
krauss@26748
   329
krauss@26748
   330
*)
krauss@26748
   331
lemma wf_Un:
krauss@26748
   332
     "[| wf r; wf s; Domain r Int Range s = {} |] ==> wf(r Un s)"
krauss@26748
   333
  using wf_union_compatible[of s r] 
krauss@26748
   334
  by (auto simp: Un_ac)
krauss@26748
   335
krauss@26748
   336
lemma wf_union_merge: 
krauss@32235
   337
  "wf (R \<union> S) = wf (R O R \<union> S O R \<union> S)" (is "wf ?A = wf ?B")
krauss@26748
   338
proof
krauss@26748
   339
  assume "wf ?A"
krauss@26748
   340
  with wf_trancl have wfT: "wf (?A^+)" .
krauss@26748
   341
  moreover have "?B \<subseteq> ?A^+"
krauss@26748
   342
    by (subst trancl_unfold, subst trancl_unfold) blast
krauss@26748
   343
  ultimately show "wf ?B" by (rule wf_subset)
krauss@26748
   344
next
krauss@26748
   345
  assume "wf ?B"
krauss@26748
   346
krauss@26748
   347
  show "wf ?A"
krauss@26748
   348
  proof (rule wfI_min)
krauss@26748
   349
    fix Q :: "'a set" and x 
krauss@26748
   350
    assume "x \<in> Q"
krauss@26748
   351
krauss@26748
   352
    with `wf ?B`
krauss@26748
   353
    obtain z where "z \<in> Q" and "\<And>y. (y, z) \<in> ?B \<Longrightarrow> y \<notin> Q" 
krauss@26748
   354
      by (erule wfE_min)
krauss@26748
   355
    then have A1: "\<And>y. (y, z) \<in> R O R \<Longrightarrow> y \<notin> Q"
krauss@32235
   356
      and A2: "\<And>y. (y, z) \<in> S O R \<Longrightarrow> y \<notin> Q"
krauss@26748
   357
      and A3: "\<And>y. (y, z) \<in> S \<Longrightarrow> y \<notin> Q"
krauss@26748
   358
      by auto
krauss@26748
   359
    
krauss@26748
   360
    show "\<exists>z\<in>Q. \<forall>y. (y, z) \<in> ?A \<longrightarrow> y \<notin> Q"
krauss@26748
   361
    proof (cases "\<forall>y. (y, z) \<in> R \<longrightarrow> y \<notin> Q")
krauss@26748
   362
      case True
krauss@26748
   363
      with `z \<in> Q` A3 show ?thesis by blast
krauss@26748
   364
    next
krauss@26748
   365
      case False 
krauss@26748
   366
      then obtain z' where "z'\<in>Q" "(z', z) \<in> R" by blast
krauss@26748
   367
krauss@26748
   368
      have "\<forall>y. (y, z') \<in> ?A \<longrightarrow> y \<notin> Q"
krauss@26748
   369
      proof (intro allI impI)
krauss@26748
   370
        fix y assume "(y, z') \<in> ?A"
krauss@26748
   371
        then show "y \<notin> Q"
krauss@26748
   372
        proof
krauss@26748
   373
          assume "(y, z') \<in> R" 
krauss@26748
   374
          then have "(y, z) \<in> R O R" using `(z', z) \<in> R` ..
krauss@26748
   375
          with A1 show "y \<notin> Q" .
krauss@26748
   376
        next
krauss@26748
   377
          assume "(y, z') \<in> S" 
krauss@32235
   378
          then have "(y, z) \<in> S O R" using  `(z', z) \<in> R` ..
krauss@26748
   379
          with A2 show "y \<notin> Q" .
krauss@26748
   380
        qed
krauss@26748
   381
      qed
krauss@26748
   382
      with `z' \<in> Q` show ?thesis ..
krauss@26748
   383
    qed
krauss@26748
   384
  qed
krauss@26748
   385
qed
krauss@26748
   386
krauss@26748
   387
lemma wf_comp_self: "wf R = wf (R O R)"  -- {* special case *}
krauss@26748
   388
  by (rule wf_union_merge [where S = "{}", simplified])
krauss@26748
   389
krauss@26748
   390
krauss@26748
   391
subsubsection {* acyclic *}
krauss@26748
   392
krauss@26748
   393
lemma acyclicI: "ALL x. (x, x) ~: r^+ ==> acyclic r"
krauss@26748
   394
  by (simp add: acyclic_def)
krauss@26748
   395
krauss@26748
   396
lemma wf_acyclic: "wf r ==> acyclic r"
krauss@26748
   397
apply (simp add: acyclic_def)
krauss@26748
   398
apply (blast elim: wf_trancl [THEN wf_irrefl])
krauss@26748
   399
done
krauss@26748
   400
krauss@26748
   401
lemmas wfP_acyclicP = wf_acyclic [to_pred]
krauss@26748
   402
krauss@26748
   403
lemma acyclic_insert [iff]:
krauss@26748
   404
     "acyclic(insert (y,x) r) = (acyclic r & (x,y) ~: r^*)"
krauss@26748
   405
apply (simp add: acyclic_def trancl_insert)
krauss@26748
   406
apply (blast intro: rtrancl_trans)
krauss@26748
   407
done
krauss@26748
   408
krauss@26748
   409
lemma acyclic_converse [iff]: "acyclic(r^-1) = acyclic r"
krauss@26748
   410
by (simp add: acyclic_def trancl_converse)
krauss@26748
   411
krauss@26748
   412
lemmas acyclicP_converse [iff] = acyclic_converse [to_pred]
krauss@26748
   413
krauss@26748
   414
lemma acyclic_impl_antisym_rtrancl: "acyclic r ==> antisym(r^*)"
krauss@26748
   415
apply (simp add: acyclic_def antisym_def)
krauss@26748
   416
apply (blast elim: rtranclE intro: rtrancl_into_trancl1 rtrancl_trancl_trancl)
krauss@26748
   417
done
krauss@26748
   418
krauss@26748
   419
(* Other direction:
krauss@26748
   420
acyclic = no loops
krauss@26748
   421
antisym = only self loops
krauss@26748
   422
Goalw [acyclic_def,antisym_def] "antisym( r^* ) ==> acyclic(r - Id)
krauss@26748
   423
==> antisym( r^* ) = acyclic(r - Id)";
krauss@26748
   424
*)
krauss@26748
   425
krauss@26748
   426
lemma acyclic_subset: "[| acyclic s; r <= s |] ==> acyclic r"
krauss@26748
   427
apply (simp add: acyclic_def)
krauss@26748
   428
apply (blast intro: trancl_mono)
krauss@26748
   429
done
krauss@26748
   430
krauss@26748
   431
text{* Wellfoundedness of finite acyclic relations*}
krauss@26748
   432
krauss@26748
   433
lemma finite_acyclic_wf [rule_format]: "finite r ==> acyclic r --> wf r"
krauss@26748
   434
apply (erule finite_induct, blast)
krauss@26748
   435
apply (simp (no_asm_simp) only: split_tupled_all)
krauss@26748
   436
apply simp
krauss@26748
   437
done
krauss@26748
   438
krauss@26748
   439
lemma finite_acyclic_wf_converse: "[|finite r; acyclic r|] ==> wf (r^-1)"
krauss@26748
   440
apply (erule finite_converse [THEN iffD2, THEN finite_acyclic_wf])
krauss@26748
   441
apply (erule acyclic_converse [THEN iffD2])
krauss@26748
   442
done
krauss@26748
   443
krauss@26748
   444
lemma wf_iff_acyclic_if_finite: "finite r ==> wf r = acyclic r"
krauss@26748
   445
by (blast intro: finite_acyclic_wf wf_acyclic)
krauss@26748
   446
krauss@26748
   447
krauss@26748
   448
subsection {* @{typ nat} is well-founded *}
krauss@26748
   449
krauss@26748
   450
lemma less_nat_rel: "op < = (\<lambda>m n. n = Suc m)^++"
krauss@26748
   451
proof (rule ext, rule ext, rule iffI)
krauss@26748
   452
  fix n m :: nat
krauss@26748
   453
  assume "m < n"
krauss@26748
   454
  then show "(\<lambda>m n. n = Suc m)^++ m n"
krauss@26748
   455
  proof (induct n)
krauss@26748
   456
    case 0 then show ?case by auto
krauss@26748
   457
  next
krauss@26748
   458
    case (Suc n) then show ?case
krauss@26748
   459
      by (auto simp add: less_Suc_eq_le le_less intro: tranclp.trancl_into_trancl)
krauss@26748
   460
  qed
krauss@26748
   461
next
krauss@26748
   462
  fix n m :: nat
krauss@26748
   463
  assume "(\<lambda>m n. n = Suc m)^++ m n"
krauss@26748
   464
  then show "m < n"
krauss@26748
   465
    by (induct n)
krauss@26748
   466
      (simp_all add: less_Suc_eq_le reflexive le_less)
krauss@26748
   467
qed
krauss@26748
   468
krauss@26748
   469
definition
krauss@26748
   470
  pred_nat :: "(nat * nat) set" where
krauss@26748
   471
  "pred_nat = {(m, n). n = Suc m}"
krauss@26748
   472
krauss@26748
   473
definition
krauss@26748
   474
  less_than :: "(nat * nat) set" where
krauss@26748
   475
  "less_than = pred_nat^+"
krauss@26748
   476
krauss@26748
   477
lemma less_eq: "(m, n) \<in> pred_nat^+ \<longleftrightarrow> m < n"
krauss@26748
   478
  unfolding less_nat_rel pred_nat_def trancl_def by simp
krauss@26748
   479
krauss@26748
   480
lemma pred_nat_trancl_eq_le:
krauss@26748
   481
  "(m, n) \<in> pred_nat^* \<longleftrightarrow> m \<le> n"
krauss@26748
   482
  unfolding less_eq rtrancl_eq_or_trancl by auto
krauss@26748
   483
krauss@26748
   484
lemma wf_pred_nat: "wf pred_nat"
krauss@26748
   485
  apply (unfold wf_def pred_nat_def, clarify)
krauss@26748
   486
  apply (induct_tac x, blast+)
krauss@26748
   487
  done
krauss@26748
   488
krauss@26748
   489
lemma wf_less_than [iff]: "wf less_than"
krauss@26748
   490
  by (simp add: less_than_def wf_pred_nat [THEN wf_trancl])
krauss@26748
   491
krauss@26748
   492
lemma trans_less_than [iff]: "trans less_than"
krauss@26748
   493
  by (simp add: less_than_def trans_trancl)
krauss@26748
   494
krauss@26748
   495
lemma less_than_iff [iff]: "((x,y): less_than) = (x<y)"
krauss@26748
   496
  by (simp add: less_than_def less_eq)
krauss@26748
   497
krauss@26748
   498
lemma wf_less: "wf {(x, y::nat). x < y}"
krauss@26748
   499
  using wf_less_than by (simp add: less_than_def less_eq [symmetric])
krauss@26748
   500
krauss@26748
   501
krauss@26748
   502
subsection {* Accessible Part *}
krauss@26748
   503
krauss@26748
   504
text {*
krauss@26748
   505
 Inductive definition of the accessible part @{term "acc r"} of a
krauss@26748
   506
 relation; see also \cite{paulin-tlca}.
krauss@26748
   507
*}
krauss@26748
   508
krauss@26748
   509
inductive_set
krauss@26748
   510
  acc :: "('a * 'a) set => 'a set"
krauss@26748
   511
  for r :: "('a * 'a) set"
krauss@26748
   512
  where
krauss@26748
   513
    accI: "(!!y. (y, x) : r ==> y : acc r) ==> x : acc r"
krauss@26748
   514
krauss@26748
   515
abbreviation
krauss@26748
   516
  termip :: "('a => 'a => bool) => 'a => bool" where
krauss@26748
   517
  "termip r == accp (r\<inverse>\<inverse>)"
krauss@26748
   518
krauss@26748
   519
abbreviation
krauss@26748
   520
  termi :: "('a * 'a) set => 'a set" where
krauss@26748
   521
  "termi r == acc (r\<inverse>)"
krauss@26748
   522
krauss@26748
   523
lemmas accpI = accp.accI
krauss@26748
   524
krauss@26748
   525
text {* Induction rules *}
krauss@26748
   526
krauss@26748
   527
theorem accp_induct:
krauss@26748
   528
  assumes major: "accp r a"
krauss@26748
   529
  assumes hyp: "!!x. accp r x ==> \<forall>y. r y x --> P y ==> P x"
krauss@26748
   530
  shows "P a"
krauss@26748
   531
  apply (rule major [THEN accp.induct])
krauss@26748
   532
  apply (rule hyp)
krauss@26748
   533
   apply (rule accp.accI)
krauss@26748
   534
   apply fast
krauss@26748
   535
  apply fast
krauss@26748
   536
  done
krauss@26748
   537
krauss@26748
   538
theorems accp_induct_rule = accp_induct [rule_format, induct set: accp]
krauss@26748
   539
krauss@26748
   540
theorem accp_downward: "accp r b ==> r a b ==> accp r a"
krauss@26748
   541
  apply (erule accp.cases)
krauss@26748
   542
  apply fast
krauss@26748
   543
  done
krauss@26748
   544
krauss@26748
   545
lemma not_accp_down:
krauss@26748
   546
  assumes na: "\<not> accp R x"
krauss@26748
   547
  obtains z where "R z x" and "\<not> accp R z"
krauss@26748
   548
proof -
krauss@26748
   549
  assume a: "\<And>z. \<lbrakk>R z x; \<not> accp R z\<rbrakk> \<Longrightarrow> thesis"
krauss@26748
   550
krauss@26748
   551
  show thesis
krauss@26748
   552
  proof (cases "\<forall>z. R z x \<longrightarrow> accp R z")
krauss@26748
   553
    case True
krauss@26748
   554
    hence "\<And>z. R z x \<Longrightarrow> accp R z" by auto
krauss@26748
   555
    hence "accp R x"
krauss@26748
   556
      by (rule accp.accI)
krauss@26748
   557
    with na show thesis ..
krauss@26748
   558
  next
krauss@26748
   559
    case False then obtain z where "R z x" and "\<not> accp R z"
krauss@26748
   560
      by auto
krauss@26748
   561
    with a show thesis .
krauss@26748
   562
  qed
krauss@26748
   563
qed
krauss@26748
   564
krauss@26748
   565
lemma accp_downwards_aux: "r\<^sup>*\<^sup>* b a ==> accp r a --> accp r b"
krauss@26748
   566
  apply (erule rtranclp_induct)
krauss@26748
   567
   apply blast
krauss@26748
   568
  apply (blast dest: accp_downward)
krauss@26748
   569
  done
krauss@26748
   570
krauss@26748
   571
theorem accp_downwards: "accp r a ==> r\<^sup>*\<^sup>* b a ==> accp r b"
krauss@26748
   572
  apply (blast dest: accp_downwards_aux)
krauss@26748
   573
  done
krauss@26748
   574
krauss@26748
   575
theorem accp_wfPI: "\<forall>x. accp r x ==> wfP r"
krauss@26748
   576
  apply (rule wfPUNIVI)
krauss@26748
   577
  apply (induct_tac P x rule: accp_induct)
krauss@26748
   578
   apply blast
krauss@26748
   579
  apply blast
krauss@26748
   580
  done
krauss@26748
   581
krauss@26748
   582
theorem accp_wfPD: "wfP r ==> accp r x"
krauss@26748
   583
  apply (erule wfP_induct_rule)
krauss@26748
   584
  apply (rule accp.accI)
krauss@26748
   585
  apply blast
krauss@26748
   586
  done
krauss@26748
   587
krauss@26748
   588
theorem wfP_accp_iff: "wfP r = (\<forall>x. accp r x)"
krauss@26748
   589
  apply (blast intro: accp_wfPI dest: accp_wfPD)
krauss@26748
   590
  done
krauss@26748
   591
krauss@26748
   592
krauss@26748
   593
text {* Smaller relations have bigger accessible parts: *}
krauss@26748
   594
krauss@26748
   595
lemma accp_subset:
krauss@26748
   596
  assumes sub: "R1 \<le> R2"
krauss@26748
   597
  shows "accp R2 \<le> accp R1"
berghofe@26803
   598
proof (rule predicate1I)
krauss@26748
   599
  fix x assume "accp R2 x"
krauss@26748
   600
  then show "accp R1 x"
krauss@26748
   601
  proof (induct x)
krauss@26748
   602
    fix x
krauss@26748
   603
    assume ih: "\<And>y. R2 y x \<Longrightarrow> accp R1 y"
krauss@26748
   604
    with sub show "accp R1 x"
krauss@26748
   605
      by (blast intro: accp.accI)
krauss@26748
   606
  qed
krauss@26748
   607
qed
krauss@26748
   608
krauss@26748
   609
krauss@26748
   610
text {* This is a generalized induction theorem that works on
krauss@26748
   611
  subsets of the accessible part. *}
krauss@26748
   612
krauss@26748
   613
lemma accp_subset_induct:
krauss@26748
   614
  assumes subset: "D \<le> accp R"
krauss@26748
   615
    and dcl: "\<And>x z. \<lbrakk>D x; R z x\<rbrakk> \<Longrightarrow> D z"
krauss@26748
   616
    and "D x"
krauss@26748
   617
    and istep: "\<And>x. \<lbrakk>D x; (\<And>z. R z x \<Longrightarrow> P z)\<rbrakk> \<Longrightarrow> P x"
krauss@26748
   618
  shows "P x"
krauss@26748
   619
proof -
krauss@26748
   620
  from subset and `D x`
krauss@26748
   621
  have "accp R x" ..
krauss@26748
   622
  then show "P x" using `D x`
krauss@26748
   623
  proof (induct x)
krauss@26748
   624
    fix x
krauss@26748
   625
    assume "D x"
krauss@26748
   626
      and "\<And>y. R y x \<Longrightarrow> D y \<Longrightarrow> P y"
krauss@26748
   627
    with dcl and istep show "P x" by blast
krauss@26748
   628
  qed
krauss@26748
   629
qed
krauss@26748
   630
krauss@26748
   631
krauss@26748
   632
text {* Set versions of the above theorems *}
krauss@26748
   633
krauss@26748
   634
lemmas acc_induct = accp_induct [to_set]
krauss@26748
   635
krauss@26748
   636
lemmas acc_induct_rule = acc_induct [rule_format, induct set: acc]
krauss@26748
   637
krauss@26748
   638
lemmas acc_downward = accp_downward [to_set]
krauss@26748
   639
krauss@26748
   640
lemmas not_acc_down = not_accp_down [to_set]
krauss@26748
   641
krauss@26748
   642
lemmas acc_downwards_aux = accp_downwards_aux [to_set]
krauss@26748
   643
krauss@26748
   644
lemmas acc_downwards = accp_downwards [to_set]
krauss@26748
   645
krauss@26748
   646
lemmas acc_wfI = accp_wfPI [to_set]
krauss@26748
   647
krauss@26748
   648
lemmas acc_wfD = accp_wfPD [to_set]
krauss@26748
   649
krauss@26748
   650
lemmas wf_acc_iff = wfP_accp_iff [to_set]
krauss@26748
   651
berghofe@26803
   652
lemmas acc_subset = accp_subset [to_set pred_subset_eq]
krauss@26748
   653
berghofe@26803
   654
lemmas acc_subset_induct = accp_subset_induct [to_set pred_subset_eq]
krauss@26748
   655
krauss@26748
   656
krauss@26748
   657
subsection {* Tools for building wellfounded relations *}
krauss@26748
   658
krauss@26748
   659
text {* Inverse Image *}
krauss@26748
   660
krauss@26748
   661
lemma wf_inv_image [simp,intro!]: "wf(r) ==> wf(inv_image r (f::'a=>'b))"
krauss@26748
   662
apply (simp (no_asm_use) add: inv_image_def wf_eq_minimal)
krauss@26748
   663
apply clarify
krauss@26748
   664
apply (subgoal_tac "EX (w::'b) . w : {w. EX (x::'a) . x: Q & (f x = w) }")
krauss@26748
   665
prefer 2 apply (blast del: allE)
krauss@26748
   666
apply (erule allE)
krauss@26748
   667
apply (erule (1) notE impE)
krauss@26748
   668
apply blast
krauss@26748
   669
done
krauss@26748
   670
haftmann@31775
   671
text {* Measure Datatypes into @{typ nat} *}
krauss@26748
   672
krauss@26748
   673
definition measure :: "('a => nat) => ('a * 'a)set"
krauss@26748
   674
where "measure == inv_image less_than"
krauss@26748
   675
krauss@26748
   676
lemma in_measure[simp]: "((x,y) : measure f) = (f x < f y)"
krauss@26748
   677
  by (simp add:measure_def)
krauss@26748
   678
krauss@26748
   679
lemma wf_measure [iff]: "wf (measure f)"
krauss@26748
   680
apply (unfold measure_def)
krauss@26748
   681
apply (rule wf_less_than [THEN wf_inv_image])
krauss@26748
   682
done
krauss@26748
   683
krauss@26748
   684
text{* Lexicographic combinations *}
krauss@26748
   685
krauss@26748
   686
definition
krauss@26748
   687
 lex_prod  :: "[('a*'a)set, ('b*'b)set] => (('a*'b)*('a*'b))set"
krauss@26748
   688
               (infixr "<*lex*>" 80)
krauss@26748
   689
where
krauss@26748
   690
    "ra <*lex*> rb == {((a,b),(a',b')). (a,a') : ra | a=a' & (b,b') : rb}"
krauss@26748
   691
krauss@26748
   692
lemma wf_lex_prod [intro!]: "[| wf(ra); wf(rb) |] ==> wf(ra <*lex*> rb)"
krauss@26748
   693
apply (unfold wf_def lex_prod_def) 
krauss@26748
   694
apply (rule allI, rule impI)
krauss@26748
   695
apply (simp (no_asm_use) only: split_paired_All)
krauss@26748
   696
apply (drule spec, erule mp) 
krauss@26748
   697
apply (rule allI, rule impI)
krauss@26748
   698
apply (drule spec, erule mp, blast) 
krauss@26748
   699
done
krauss@26748
   700
krauss@26748
   701
lemma in_lex_prod[simp]: 
krauss@26748
   702
  "(((a,b),(a',b')): r <*lex*> s) = ((a,a'): r \<or> (a = a' \<and> (b, b') : s))"
krauss@26748
   703
  by (auto simp:lex_prod_def)
krauss@26748
   704
krauss@26748
   705
text{* @{term "op <*lex*>"} preserves transitivity *}
krauss@26748
   706
krauss@26748
   707
lemma trans_lex_prod [intro!]: 
krauss@26748
   708
    "[| trans R1; trans R2 |] ==> trans (R1 <*lex*> R2)"
krauss@26748
   709
by (unfold trans_def lex_prod_def, blast) 
krauss@26748
   710
haftmann@31775
   711
text {* lexicographic combinations with measure Datatypes *}
krauss@26748
   712
krauss@26748
   713
definition 
krauss@26748
   714
  mlex_prod :: "('a \<Rightarrow> nat) \<Rightarrow> ('a \<times> 'a) set \<Rightarrow> ('a \<times> 'a) set" (infixr "<*mlex*>" 80)
krauss@26748
   715
where
krauss@26748
   716
  "f <*mlex*> R = inv_image (less_than <*lex*> R) (%x. (f x, x))"
krauss@26748
   717
krauss@26748
   718
lemma wf_mlex: "wf R \<Longrightarrow> wf (f <*mlex*> R)"
krauss@26748
   719
unfolding mlex_prod_def
krauss@26748
   720
by auto
krauss@26748
   721
krauss@26748
   722
lemma mlex_less: "f x < f y \<Longrightarrow> (x, y) \<in> f <*mlex*> R"
krauss@26748
   723
unfolding mlex_prod_def by simp
krauss@26748
   724
krauss@26748
   725
lemma mlex_leq: "f x \<le> f y \<Longrightarrow> (x, y) \<in> R \<Longrightarrow> (x, y) \<in> f <*mlex*> R"
krauss@26748
   726
unfolding mlex_prod_def by auto
krauss@26748
   727
krauss@26748
   728
text {* proper subset relation on finite sets *}
krauss@26748
   729
krauss@26748
   730
definition finite_psubset  :: "('a set * 'a set) set"
krauss@26748
   731
where "finite_psubset == {(A,B). A < B & finite B}"
krauss@26748
   732
krauss@28260
   733
lemma wf_finite_psubset[simp]: "wf(finite_psubset)"
krauss@26748
   734
apply (unfold finite_psubset_def)
krauss@26748
   735
apply (rule wf_measure [THEN wf_subset])
krauss@26748
   736
apply (simp add: measure_def inv_image_def less_than_def less_eq)
krauss@26748
   737
apply (fast elim!: psubset_card_mono)
krauss@26748
   738
done
krauss@26748
   739
krauss@26748
   740
lemma trans_finite_psubset: "trans finite_psubset"
berghofe@26803
   741
by (simp add: finite_psubset_def less_le trans_def, blast)
krauss@26748
   742
krauss@28260
   743
lemma in_finite_psubset[simp]: "(A, B) \<in> finite_psubset = (A < B & finite B)"
krauss@28260
   744
unfolding finite_psubset_def by auto
krauss@26748
   745
krauss@28735
   746
text {* max- and min-extension of order to finite sets *}
krauss@28735
   747
krauss@28735
   748
inductive_set max_ext :: "('a \<times> 'a) set \<Rightarrow> ('a set \<times> 'a set) set" 
krauss@28735
   749
for R :: "('a \<times> 'a) set"
krauss@28735
   750
where
krauss@28735
   751
  max_extI[intro]: "finite X \<Longrightarrow> finite Y \<Longrightarrow> Y \<noteq> {} \<Longrightarrow> (\<And>x. x \<in> X \<Longrightarrow> \<exists>y\<in>Y. (x, y) \<in> R) \<Longrightarrow> (X, Y) \<in> max_ext R"
krauss@28735
   752
krauss@28735
   753
lemma max_ext_wf:
krauss@28735
   754
  assumes wf: "wf r"
krauss@28735
   755
  shows "wf (max_ext r)"
krauss@28735
   756
proof (rule acc_wfI, intro allI)
krauss@28735
   757
  fix M show "M \<in> acc (max_ext r)" (is "_ \<in> ?W")
krauss@28735
   758
  proof cases
krauss@28735
   759
    assume "finite M"
krauss@28735
   760
    thus ?thesis
krauss@28735
   761
    proof (induct M)
krauss@28735
   762
      show "{} \<in> ?W"
krauss@28735
   763
        by (rule accI) (auto elim: max_ext.cases)
krauss@28735
   764
    next
krauss@28735
   765
      fix M a assume "M \<in> ?W" "finite M"
krauss@28735
   766
      with wf show "insert a M \<in> ?W"
krauss@28735
   767
      proof (induct arbitrary: M)
krauss@28735
   768
        fix M a
krauss@28735
   769
        assume "M \<in> ?W"  and  [intro]: "finite M"
krauss@28735
   770
        assume hyp: "\<And>b M. (b, a) \<in> r \<Longrightarrow> M \<in> ?W \<Longrightarrow> finite M \<Longrightarrow> insert b M \<in> ?W"
krauss@28735
   771
        {
krauss@28735
   772
          fix N M :: "'a set"
krauss@28735
   773
          assume "finite N" "finite M"
krauss@28735
   774
          then
krauss@28735
   775
          have "\<lbrakk>M \<in> ?W ; (\<And>y. y \<in> N \<Longrightarrow> (y, a) \<in> r)\<rbrakk> \<Longrightarrow>  N \<union> M \<in> ?W"
krauss@28735
   776
            by (induct N arbitrary: M) (auto simp: hyp)
krauss@28735
   777
        }
krauss@28735
   778
        note add_less = this
krauss@28735
   779
        
krauss@28735
   780
        show "insert a M \<in> ?W"
krauss@28735
   781
        proof (rule accI)
krauss@28735
   782
          fix N assume Nless: "(N, insert a M) \<in> max_ext r"
krauss@28735
   783
          hence asm1: "\<And>x. x \<in> N \<Longrightarrow> (x, a) \<in> r \<or> (\<exists>y \<in> M. (x, y) \<in> r)"
krauss@28735
   784
            by (auto elim!: max_ext.cases)
krauss@28735
   785
krauss@28735
   786
          let ?N1 = "{ n \<in> N. (n, a) \<in> r }"
krauss@28735
   787
          let ?N2 = "{ n \<in> N. (n, a) \<notin> r }"
krauss@28735
   788
          have N: "?N1 \<union> ?N2 = N" by (rule set_ext) auto
krauss@28735
   789
          from Nless have "finite N" by (auto elim: max_ext.cases)
krauss@28735
   790
          then have finites: "finite ?N1" "finite ?N2" by auto
krauss@28735
   791
          
krauss@28735
   792
          have "?N2 \<in> ?W"
krauss@28735
   793
          proof cases
krauss@28735
   794
            assume [simp]: "M = {}"
krauss@28735
   795
            have Mw: "{} \<in> ?W" by (rule accI) (auto elim: max_ext.cases)
krauss@28735
   796
krauss@28735
   797
            from asm1 have "?N2 = {}" by auto
krauss@28735
   798
            with Mw show "?N2 \<in> ?W" by (simp only:)
krauss@28735
   799
          next
krauss@28735
   800
            assume "M \<noteq> {}"
krauss@28735
   801
            have N2: "(?N2, M) \<in> max_ext r" 
krauss@28735
   802
              by (rule max_extI[OF _ _ `M \<noteq> {}`]) (insert asm1, auto intro: finites)
krauss@28735
   803
            
krauss@28735
   804
            with `M \<in> ?W` show "?N2 \<in> ?W" by (rule acc_downward)
krauss@28735
   805
          qed
krauss@28735
   806
          with finites have "?N1 \<union> ?N2 \<in> ?W" 
krauss@28735
   807
            by (rule add_less) simp
krauss@28735
   808
          then show "N \<in> ?W" by (simp only: N)
krauss@28735
   809
        qed
krauss@28735
   810
      qed
krauss@28735
   811
    qed
krauss@28735
   812
  next
krauss@28735
   813
    assume [simp]: "\<not> finite M"
krauss@28735
   814
    show ?thesis
krauss@28735
   815
      by (rule accI) (auto elim: max_ext.cases)
krauss@28735
   816
  qed
krauss@28735
   817
qed
krauss@28735
   818
krauss@29125
   819
lemma max_ext_additive: 
krauss@29125
   820
 "(A, B) \<in> max_ext R \<Longrightarrow> (C, D) \<in> max_ext R \<Longrightarrow>
krauss@29125
   821
  (A \<union> C, B \<union> D) \<in> max_ext R"
krauss@29125
   822
by (force elim!: max_ext.cases)
krauss@29125
   823
krauss@28735
   824
krauss@28735
   825
definition
krauss@28735
   826
  min_ext :: "('a \<times> 'a) set \<Rightarrow> ('a set \<times> 'a set) set" 
krauss@28735
   827
where
krauss@28735
   828
  [code del]: "min_ext r = {(X, Y) | X Y. X \<noteq> {} \<and> (\<forall>y \<in> Y. (\<exists>x \<in> X. (x, y) \<in> r))}"
krauss@28735
   829
krauss@28735
   830
lemma min_ext_wf:
krauss@28735
   831
  assumes "wf r"
krauss@28735
   832
  shows "wf (min_ext r)"
krauss@28735
   833
proof (rule wfI_min)
krauss@28735
   834
  fix Q :: "'a set set"
krauss@28735
   835
  fix x
krauss@28735
   836
  assume nonempty: "x \<in> Q"
krauss@28735
   837
  show "\<exists>m \<in> Q. (\<forall> n. (n, m) \<in> min_ext r \<longrightarrow> n \<notin> Q)"
krauss@28735
   838
  proof cases
krauss@28735
   839
    assume "Q = {{}}" thus ?thesis by (simp add: min_ext_def)
krauss@28735
   840
  next
krauss@28735
   841
    assume "Q \<noteq> {{}}"
krauss@28735
   842
    with nonempty
krauss@28735
   843
    obtain e x where "x \<in> Q" "e \<in> x" by force
krauss@28735
   844
    then have eU: "e \<in> \<Union>Q" by auto
krauss@28735
   845
    with `wf r` 
krauss@28735
   846
    obtain z where z: "z \<in> \<Union>Q" "\<And>y. (y, z) \<in> r \<Longrightarrow> y \<notin> \<Union>Q" 
krauss@28735
   847
      by (erule wfE_min)
krauss@28735
   848
    from z obtain m where "m \<in> Q" "z \<in> m" by auto
krauss@28735
   849
    from `m \<in> Q`
krauss@28735
   850
    show ?thesis
krauss@28735
   851
    proof (rule, intro bexI allI impI)
krauss@28735
   852
      fix n
krauss@28735
   853
      assume smaller: "(n, m) \<in> min_ext r"
krauss@28735
   854
      with `z \<in> m` obtain y where y: "y \<in> n" "(y, z) \<in> r" by (auto simp: min_ext_def)
krauss@28735
   855
      then show "n \<notin> Q" using z(2) by auto
krauss@28735
   856
    qed      
krauss@28735
   857
  qed
krauss@28735
   858
qed
krauss@26748
   859
krauss@26748
   860
krauss@26748
   861
subsection{*Weakly decreasing sequences (w.r.t. some well-founded order) 
krauss@26748
   862
   stabilize.*}
krauss@26748
   863
krauss@26748
   864
text{*This material does not appear to be used any longer.*}
krauss@26748
   865
krauss@28845
   866
lemma sequence_trans: "[| ALL i. (f (Suc i), f i) : r^* |] ==> (f (i+k), f i) : r^*"
krauss@28845
   867
by (induct k) (auto intro: rtrancl_trans)
krauss@26748
   868
krauss@28845
   869
lemma wf_weak_decr_stable: 
krauss@28845
   870
  assumes as: "ALL i. (f (Suc i), f i) : r^*" "wf (r^+)"
krauss@28845
   871
  shows "EX i. ALL k. f (i+k) = f i"
krauss@28845
   872
proof -
krauss@28845
   873
  have lem: "!!x. [| ALL i. (f (Suc i), f i) : r^*; wf (r^+) |]  
krauss@26748
   874
      ==> ALL m. f m = x --> (EX i. ALL k. f (m+i+k) = f (m+i))"
krauss@28845
   875
  apply (erule wf_induct, clarify)
krauss@28845
   876
  apply (case_tac "EX j. (f (m+j), f m) : r^+")
krauss@28845
   877
   apply clarify
krauss@28845
   878
   apply (subgoal_tac "EX i. ALL k. f ((m+j) +i+k) = f ( (m+j) +i) ")
krauss@28845
   879
    apply clarify
krauss@28845
   880
    apply (rule_tac x = "j+i" in exI)
krauss@28845
   881
    apply (simp add: add_ac, blast)
krauss@28845
   882
  apply (rule_tac x = 0 in exI, clarsimp)
krauss@28845
   883
  apply (drule_tac i = m and k = k in sequence_trans)
krauss@28845
   884
  apply (blast elim: rtranclE dest: rtrancl_into_trancl1)
krauss@28845
   885
  done
krauss@26748
   886
krauss@28845
   887
  from lem[OF as, THEN spec, of 0, simplified] 
krauss@28845
   888
  show ?thesis by auto
krauss@28845
   889
qed
krauss@26748
   890
krauss@26748
   891
(* special case of the theorem above: <= *)
krauss@26748
   892
lemma weak_decr_stable:
krauss@26748
   893
     "ALL i. f (Suc i) <= ((f i)::nat) ==> EX i. ALL k. f (i+k) = f i"
krauss@26748
   894
apply (rule_tac r = pred_nat in wf_weak_decr_stable)
krauss@26748
   895
apply (simp add: pred_nat_trancl_eq_le)
krauss@26748
   896
apply (intro wf_trancl wf_pred_nat)
krauss@26748
   897
done
krauss@26748
   898
krauss@26748
   899
krauss@26748
   900
subsection {* size of a datatype value *}
krauss@26748
   901
haftmann@31775
   902
use "Tools/Function/size.ML"
krauss@26748
   903
krauss@26748
   904
setup Size.setup
krauss@26748
   905
haftmann@28562
   906
lemma size_bool [code]:
haftmann@27823
   907
  "size (b\<Colon>bool) = 0" by (cases b) auto
haftmann@27823
   908
haftmann@28562
   909
lemma nat_size [simp, code]: "size (n\<Colon>nat) = n"
krauss@26748
   910
  by (induct n) simp_all
krauss@26748
   911
haftmann@27823
   912
declare "prod.size" [noatp]
krauss@26748
   913
haftmann@30430
   914
lemma [code]:
haftmann@30430
   915
  "size (P :: 'a Predicate.pred) = 0" by (cases P) simp
haftmann@30430
   916
haftmann@30430
   917
lemma [code]:
haftmann@30430
   918
  "pred_size f P = 0" by (cases P) simp
haftmann@30430
   919
krauss@26748
   920
end