| author | wenzelm |
| Tue, 05 Nov 2024 23:01:09 +0100 | |
| changeset 81351 | 95cb584cb777 |
| parent 80934 | 8e72f55295fd |
| child 81595 | ed264056f5dc |
| permissions | -rw-r--r-- |
| 58101 | 1 |
(* Author: Tobias Nipkow, TU Muenchen *) |
2 |
||
| 60758 | 3 |
section \<open>Sum and product over lists\<close> |
| 58101 | 4 |
|
5 |
theory Groups_List |
|
6 |
imports List |
|
7 |
begin |
|
8 |
||
| 58320 | 9 |
locale monoid_list = monoid |
10 |
begin |
|
|
67489
f1ba59ddd9a6
drop redundant cong rules
Lars Hupel <lars.hupel@mytum.de>
parents:
67399
diff
changeset
|
11 |
|
| 58320 | 12 |
definition F :: "'a list \<Rightarrow> 'a" |
13 |
where |
|
|
63290
9ac558ab0906
boldify syntax in abstract algebraic structures, to avoid clashes with concrete syntax in corresponding type classes
haftmann
parents:
63101
diff
changeset
|
14 |
eq_foldr [code]: "F xs = foldr f xs \<^bold>1" |
|
67489
f1ba59ddd9a6
drop redundant cong rules
Lars Hupel <lars.hupel@mytum.de>
parents:
67399
diff
changeset
|
15 |
|
| 58320 | 16 |
lemma Nil [simp]: |
|
63290
9ac558ab0906
boldify syntax in abstract algebraic structures, to avoid clashes with concrete syntax in corresponding type classes
haftmann
parents:
63101
diff
changeset
|
17 |
"F [] = \<^bold>1" |
| 58320 | 18 |
by (simp add: eq_foldr) |
|
67489
f1ba59ddd9a6
drop redundant cong rules
Lars Hupel <lars.hupel@mytum.de>
parents:
67399
diff
changeset
|
19 |
|
| 58320 | 20 |
lemma Cons [simp]: |
|
63290
9ac558ab0906
boldify syntax in abstract algebraic structures, to avoid clashes with concrete syntax in corresponding type classes
haftmann
parents:
63101
diff
changeset
|
21 |
"F (x # xs) = x \<^bold>* F xs" |
| 58320 | 22 |
by (simp add: eq_foldr) |
|
67489
f1ba59ddd9a6
drop redundant cong rules
Lars Hupel <lars.hupel@mytum.de>
parents:
67399
diff
changeset
|
23 |
|
| 58320 | 24 |
lemma append [simp]: |
|
63290
9ac558ab0906
boldify syntax in abstract algebraic structures, to avoid clashes with concrete syntax in corresponding type classes
haftmann
parents:
63101
diff
changeset
|
25 |
"F (xs @ ys) = F xs \<^bold>* F ys" |
| 58320 | 26 |
by (induct xs) (simp_all add: assoc) |
|
67489
f1ba59ddd9a6
drop redundant cong rules
Lars Hupel <lars.hupel@mytum.de>
parents:
67399
diff
changeset
|
27 |
|
| 58320 | 28 |
end |
| 58101 | 29 |
|
| 58320 | 30 |
locale comm_monoid_list = comm_monoid + monoid_list |
31 |
begin |
|
|
67489
f1ba59ddd9a6
drop redundant cong rules
Lars Hupel <lars.hupel@mytum.de>
parents:
67399
diff
changeset
|
32 |
|
| 58320 | 33 |
lemma rev [simp]: |
34 |
"F (rev xs) = F xs" |
|
35 |
by (simp add: eq_foldr foldr_fold fold_rev fun_eq_iff assoc left_commute) |
|
|
67489
f1ba59ddd9a6
drop redundant cong rules
Lars Hupel <lars.hupel@mytum.de>
parents:
67399
diff
changeset
|
36 |
|
| 58320 | 37 |
end |
|
67489
f1ba59ddd9a6
drop redundant cong rules
Lars Hupel <lars.hupel@mytum.de>
parents:
67399
diff
changeset
|
38 |
|
| 58320 | 39 |
locale comm_monoid_list_set = list: comm_monoid_list + set: comm_monoid_set |
40 |
begin |
|
| 58101 | 41 |
|
| 58320 | 42 |
lemma distinct_set_conv_list: |
43 |
"distinct xs \<Longrightarrow> set.F g (set xs) = list.F (map g xs)" |
|
44 |
by (induct xs) simp_all |
|
| 58101 | 45 |
|
| 58320 | 46 |
lemma set_conv_list [code]: |
47 |
"set.F g (set xs) = list.F (map g (remdups xs))" |
|
48 |
by (simp add: distinct_set_conv_list [symmetric]) |
|
49 |
||
|
80061
4c1347e172b1
moved over material from AFP; most importantly on algebraic numbers and algebraically closed fields
Manuel Eberl <eberlm@in.tum.de>
parents:
79017
diff
changeset
|
50 |
lemma list_conv_set_nth: |
|
4c1347e172b1
moved over material from AFP; most importantly on algebraic numbers and algebraically closed fields
Manuel Eberl <eberlm@in.tum.de>
parents:
79017
diff
changeset
|
51 |
"list.F xs = set.F (\<lambda>i. xs ! i) {0..<length xs}"
|
|
4c1347e172b1
moved over material from AFP; most importantly on algebraic numbers and algebraically closed fields
Manuel Eberl <eberlm@in.tum.de>
parents:
79017
diff
changeset
|
52 |
proof - |
|
4c1347e172b1
moved over material from AFP; most importantly on algebraic numbers and algebraically closed fields
Manuel Eberl <eberlm@in.tum.de>
parents:
79017
diff
changeset
|
53 |
have "xs = map (\<lambda>i. xs ! i) [0..<length xs]" |
|
4c1347e172b1
moved over material from AFP; most importantly on algebraic numbers and algebraically closed fields
Manuel Eberl <eberlm@in.tum.de>
parents:
79017
diff
changeset
|
54 |
by (simp add: map_nth) |
|
4c1347e172b1
moved over material from AFP; most importantly on algebraic numbers and algebraically closed fields
Manuel Eberl <eberlm@in.tum.de>
parents:
79017
diff
changeset
|
55 |
also have "list.F \<dots> = set.F (\<lambda>i. xs ! i) {0..<length xs}"
|
|
4c1347e172b1
moved over material from AFP; most importantly on algebraic numbers and algebraically closed fields
Manuel Eberl <eberlm@in.tum.de>
parents:
79017
diff
changeset
|
56 |
by (subst distinct_set_conv_list [symmetric]) auto |
|
4c1347e172b1
moved over material from AFP; most importantly on algebraic numbers and algebraically closed fields
Manuel Eberl <eberlm@in.tum.de>
parents:
79017
diff
changeset
|
57 |
finally show ?thesis . |
|
4c1347e172b1
moved over material from AFP; most importantly on algebraic numbers and algebraically closed fields
Manuel Eberl <eberlm@in.tum.de>
parents:
79017
diff
changeset
|
58 |
qed |
|
4c1347e172b1
moved over material from AFP; most importantly on algebraic numbers and algebraically closed fields
Manuel Eberl <eberlm@in.tum.de>
parents:
79017
diff
changeset
|
59 |
|
| 58320 | 60 |
end |
61 |
||
62 |
||
| 60758 | 63 |
subsection \<open>List summation\<close> |
| 58320 | 64 |
|
65 |
context monoid_add |
|
66 |
begin |
|
67 |
||
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
68 |
sublocale sum_list: monoid_list plus 0 |
| 61776 | 69 |
defines |
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
70 |
sum_list = sum_list.F .. |
|
67489
f1ba59ddd9a6
drop redundant cong rules
Lars Hupel <lars.hupel@mytum.de>
parents:
67399
diff
changeset
|
71 |
|
| 58320 | 72 |
end |
73 |
||
74 |
context comm_monoid_add |
|
75 |
begin |
|
76 |
||
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
77 |
sublocale sum_list: comm_monoid_list plus 0 |
|
61566
c3d6e570ccef
Keyword 'rewrites' identifies rewrite morphisms.
ballarin
parents:
61378
diff
changeset
|
78 |
rewrites |
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
79 |
"monoid_list.F plus 0 = sum_list" |
| 58320 | 80 |
proof - |
81 |
show "comm_monoid_list plus 0" .. |
|
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
82 |
then interpret sum_list: comm_monoid_list plus 0 . |
|
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
83 |
from sum_list_def show "monoid_list.F plus 0 = sum_list" by simp |
| 58101 | 84 |
qed |
85 |
||
| 64267 | 86 |
sublocale sum: comm_monoid_list_set plus 0 |
|
61566
c3d6e570ccef
Keyword 'rewrites' identifies rewrite morphisms.
ballarin
parents:
61378
diff
changeset
|
87 |
rewrites |
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
88 |
"monoid_list.F plus 0 = sum_list" |
| 64267 | 89 |
and "comm_monoid_set.F plus 0 = sum" |
| 58320 | 90 |
proof - |
91 |
show "comm_monoid_list_set plus 0" .. |
|
| 64267 | 92 |
then interpret sum: comm_monoid_list_set plus 0 . |
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
93 |
from sum_list_def show "monoid_list.F plus 0 = sum_list" by simp |
| 64267 | 94 |
from sum_def show "comm_monoid_set.F plus 0 = sum" by (auto intro: sym) |
| 58320 | 95 |
qed |
96 |
||
97 |
end |
|
98 |
||
| 60758 | 99 |
text \<open>Some syntactic sugar for summing a function over a list:\<close> |
|
61955
e96292f32c3c
former "xsymbols" syntax is used by default, and ASCII replacement syntax with print mode "ASCII";
wenzelm
parents:
61799
diff
changeset
|
100 |
syntax (ASCII) |
| 80934 | 101 |
"_sum_list" :: "pttrn => 'a list => 'b => 'b" (\<open>(\<open>indent=3 notation=\<open>binder SUM\<close>\<close>SUM _<-_. _)\<close> [0, 51, 10] 10) |
| 58101 | 102 |
syntax |
| 80934 | 103 |
"_sum_list" :: "pttrn => 'a list => 'b => 'b" (\<open>(\<open>indent=3 notation=\<open>binder \<Sum>\<close>\<close>\<Sum>_\<leftarrow>_. _)\<close> [0, 51, 10] 10) |
| 80760 | 104 |
syntax_consts |
105 |
"_sum_list" == sum_list |
|
| 61799 | 106 |
translations \<comment> \<open>Beware of argument permutation!\<close> |
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
107 |
"\<Sum>x\<leftarrow>xs. b" == "CONST sum_list (CONST map (\<lambda>x. b) xs)" |
| 58101 | 108 |
|
| 70928 | 109 |
context |
110 |
includes lifting_syntax |
|
111 |
begin |
|
112 |
||
113 |
lemma sum_list_transfer [transfer_rule]: |
|
114 |
"(list_all2 A ===> A) sum_list sum_list" |
|
115 |
if [transfer_rule]: "A 0 0" "(A ===> A ===> A) (+) (+)" |
|
116 |
unfolding sum_list.eq_foldr [abs_def] |
|
117 |
by transfer_prover |
|
118 |
||
119 |
end |
|
120 |
||
| 60758 | 121 |
text \<open>TODO duplicates\<close> |
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
122 |
lemmas sum_list_simps = sum_list.Nil sum_list.Cons |
|
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
123 |
lemmas sum_list_append = sum_list.append |
|
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
124 |
lemmas sum_list_rev = sum_list.rev |
| 58320 | 125 |
|
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
126 |
lemma (in monoid_add) fold_plus_sum_list_rev: |
|
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
127 |
"fold plus xs = plus (sum_list (rev xs))" |
| 58320 | 128 |
proof |
129 |
fix x |
|
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
130 |
have "fold plus xs x = sum_list (rev xs @ [x])" |
|
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
131 |
by (simp add: foldr_conv_fold sum_list.eq_foldr) |
|
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
132 |
also have "\<dots> = sum_list (rev xs) + x" |
| 58320 | 133 |
by simp |
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
134 |
finally show "fold plus xs x = sum_list (rev xs) + x" |
| 58320 | 135 |
. |
136 |
qed |
|
137 |
||
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
138 |
lemma (in comm_monoid_add) sum_list_map_remove1: |
|
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
139 |
"x \<in> set xs \<Longrightarrow> sum_list (map f xs) = f x + sum_list (map f (remove1 x xs))" |
| 58101 | 140 |
by (induct xs) (auto simp add: ac_simps) |
141 |
||
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
142 |
lemma (in monoid_add) size_list_conv_sum_list: |
|
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
143 |
"size_list f xs = sum_list (map f xs) + size xs" |
| 58101 | 144 |
by (induct xs) auto |
145 |
||
146 |
lemma (in monoid_add) length_concat: |
|
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
147 |
"length (concat xss) = sum_list (map length xss)" |
| 58101 | 148 |
by (induct xss) simp_all |
149 |
||
150 |
lemma (in monoid_add) length_product_lists: |
|
|
69064
5840724b1d71
Prefix form of infix with * on either side no longer needs special treatment
nipkow
parents:
67489
diff
changeset
|
151 |
"length (product_lists xss) = foldr (*) (map length xss) 1" |
| 58101 | 152 |
proof (induct xss) |
153 |
case (Cons xs xss) then show ?case by (induct xs) (auto simp: length_concat o_def) |
|
154 |
qed simp |
|
155 |
||
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
156 |
lemma (in monoid_add) sum_list_map_filter: |
| 58101 | 157 |
assumes "\<And>x. x \<in> set xs \<Longrightarrow> \<not> P x \<Longrightarrow> f x = 0" |
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
158 |
shows "sum_list (map f (filter P xs)) = sum_list (map f xs)" |
| 58101 | 159 |
using assms by (induct xs) auto |
160 |
||
| 69231 | 161 |
lemma sum_list_filter_le_nat: |
162 |
fixes f :: "'a \<Rightarrow> nat" |
|
163 |
shows "sum_list (map f (filter P xs)) \<le> sum_list (map f xs)" |
|
164 |
by(induction xs; simp) |
|
165 |
||
| 64267 | 166 |
lemma (in comm_monoid_add) distinct_sum_list_conv_Sum: |
167 |
"distinct xs \<Longrightarrow> sum_list xs = Sum (set xs)" |
|
| 58101 | 168 |
by (induct xs) simp_all |
169 |
||
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
170 |
lemma sum_list_upt[simp]: |
|
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
171 |
"m \<le> n \<Longrightarrow> sum_list [m..<n] = \<Sum> {m..<n}"
|
| 64267 | 172 |
by(simp add: distinct_sum_list_conv_Sum) |
| 58995 | 173 |
|
| 66311 | 174 |
context ordered_comm_monoid_add |
175 |
begin |
|
176 |
||
177 |
lemma sum_list_nonneg: "(\<And>x. x \<in> set xs \<Longrightarrow> 0 \<le> x) \<Longrightarrow> 0 \<le> sum_list xs" |
|
178 |
by (induction xs) auto |
|
179 |
||
180 |
lemma sum_list_nonpos: "(\<And>x. x \<in> set xs \<Longrightarrow> x \<le> 0) \<Longrightarrow> sum_list xs \<le> 0" |
|
181 |
by (induction xs) (auto simp: add_nonpos_nonpos) |
|
| 58101 | 182 |
|
| 66311 | 183 |
lemma sum_list_nonneg_eq_0_iff: |
184 |
"(\<And>x. x \<in> set xs \<Longrightarrow> 0 \<le> x) \<Longrightarrow> sum_list xs = 0 \<longleftrightarrow> (\<forall>x\<in> set xs. x = 0)" |
|
185 |
by (induction xs) (simp_all add: add_nonneg_eq_0_iff sum_list_nonneg) |
|
186 |
||
187 |
end |
|
188 |
||
189 |
context canonically_ordered_monoid_add |
|
190 |
begin |
|
| 58101 | 191 |
|
| 66311 | 192 |
lemma sum_list_eq_0_iff [simp]: |
193 |
"sum_list ns = 0 \<longleftrightarrow> (\<forall>n \<in> set ns. n = 0)" |
|
194 |
by (simp add: sum_list_nonneg_eq_0_iff) |
|
195 |
||
196 |
lemma member_le_sum_list: |
|
197 |
"x \<in> set xs \<Longrightarrow> x \<le> sum_list xs" |
|
198 |
by (induction xs) (auto simp: add_increasing add_increasing2) |
|
| 58101 | 199 |
|
| 66311 | 200 |
lemma elem_le_sum_list: |
201 |
"k < size ns \<Longrightarrow> ns ! k \<le> sum_list (ns)" |
|
202 |
by (rule member_le_sum_list) simp |
|
203 |
||
204 |
end |
|
205 |
||
206 |
lemma (in ordered_cancel_comm_monoid_diff) sum_list_update: |
|
207 |
"k < size xs \<Longrightarrow> sum_list (xs[k := x]) = sum_list xs + x - xs ! k" |
|
208 |
apply(induction xs arbitrary:k) |
|
209 |
apply (auto simp: add_ac split: nat.split) |
|
210 |
apply(drule elem_le_sum_list) |
|
211 |
by (simp add: local.add_diff_assoc local.add_increasing) |
|
| 58101 | 212 |
|
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
213 |
lemma (in monoid_add) sum_list_triv: |
| 58101 | 214 |
"(\<Sum>x\<leftarrow>xs. r) = of_nat (length xs) * r" |
215 |
by (induct xs) (simp_all add: distrib_right) |
|
216 |
||
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
217 |
lemma (in monoid_add) sum_list_0 [simp]: |
| 58101 | 218 |
"(\<Sum>x\<leftarrow>xs. 0) = 0" |
219 |
by (induct xs) (simp_all add: distrib_right) |
|
220 |
||
| 61799 | 221 |
text\<open>For non-Abelian groups \<open>xs\<close> needs to be reversed on one side:\<close> |
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
222 |
lemma (in ab_group_add) uminus_sum_list_map: |
|
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
223 |
"- sum_list (map f xs) = sum_list (map (uminus \<circ> f) xs)" |
| 58101 | 224 |
by (induct xs) simp_all |
225 |
||
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
226 |
lemma (in comm_monoid_add) sum_list_addf: |
|
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
227 |
"(\<Sum>x\<leftarrow>xs. f x + g x) = sum_list (map f xs) + sum_list (map g xs)" |
| 58101 | 228 |
by (induct xs) (simp_all add: algebra_simps) |
229 |
||
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
230 |
lemma (in ab_group_add) sum_list_subtractf: |
|
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
231 |
"(\<Sum>x\<leftarrow>xs. f x - g x) = sum_list (map f xs) - sum_list (map g xs)" |
| 58101 | 232 |
by (induct xs) (simp_all add: algebra_simps) |
233 |
||
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
234 |
lemma (in semiring_0) sum_list_const_mult: |
| 58101 | 235 |
"(\<Sum>x\<leftarrow>xs. c * f x) = c * (\<Sum>x\<leftarrow>xs. f x)" |
236 |
by (induct xs) (simp_all add: algebra_simps) |
|
237 |
||
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
238 |
lemma (in semiring_0) sum_list_mult_const: |
| 58101 | 239 |
"(\<Sum>x\<leftarrow>xs. f x * c) = (\<Sum>x\<leftarrow>xs. f x) * c" |
240 |
by (induct xs) (simp_all add: algebra_simps) |
|
241 |
||
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
242 |
lemma (in ordered_ab_group_add_abs) sum_list_abs: |
|
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
243 |
"\<bar>sum_list xs\<bar> \<le> sum_list (map abs xs)" |
| 58101 | 244 |
by (induct xs) (simp_all add: order_trans [OF abs_triangle_ineq]) |
245 |
||
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
246 |
lemma sum_list_mono: |
| 58101 | 247 |
fixes f g :: "'a \<Rightarrow> 'b::{monoid_add, ordered_ab_semigroup_add}"
|
248 |
shows "(\<And>x. x \<in> set xs \<Longrightarrow> f x \<le> g x) \<Longrightarrow> (\<Sum>x\<leftarrow>xs. f x) \<le> (\<Sum>x\<leftarrow>xs. g x)" |
|
| 69231 | 249 |
by (induct xs) (simp, simp add: add_mono) |
250 |
||
251 |
lemma sum_list_strict_mono: |
|
252 |
fixes f g :: "'a \<Rightarrow> 'b::{monoid_add, strict_ordered_ab_semigroup_add}"
|
|
253 |
shows "\<lbrakk> xs \<noteq> []; \<And>x. x \<in> set xs \<Longrightarrow> f x < g x \<rbrakk> |
|
254 |
\<Longrightarrow> sum_list (map f xs) < sum_list (map g xs)" |
|
255 |
proof (induction xs) |
|
256 |
case Nil thus ?case by simp |
|
257 |
next |
|
258 |
case C: (Cons _ xs) |
|
259 |
show ?case |
|
260 |
proof (cases xs) |
|
261 |
case Nil thus ?thesis using C.prems by simp |
|
262 |
next |
|
263 |
case Cons thus ?thesis using C by(simp add: add_strict_mono) |
|
264 |
qed |
|
265 |
qed |
|
| 58101 | 266 |
|
| 75693 | 267 |
text \<open>A much more general version of this monotonicity lemma |
268 |
can be formulated with multisets and the multiset order\<close> |
|
269 |
||
270 |
lemma sum_list_mono2: fixes xs :: "'a ::ordered_comm_monoid_add list" |
|
271 |
shows "\<lbrakk> length xs = length ys; \<And>i. i < length xs \<longrightarrow> xs!i \<le> ys!i \<rbrakk> |
|
272 |
\<Longrightarrow> sum_list xs \<le> sum_list ys" |
|
273 |
apply(induction xs ys rule: list_induct2) |
|
274 |
by(auto simp: nth_Cons' less_Suc_eq_0_disj imp_ex add_mono) |
|
275 |
||
| 64267 | 276 |
lemma (in monoid_add) sum_list_distinct_conv_sum_set: |
277 |
"distinct xs \<Longrightarrow> sum_list (map f xs) = sum f (set xs)" |
|
| 58101 | 278 |
by (induct xs) simp_all |
279 |
||
| 64267 | 280 |
lemma (in monoid_add) interv_sum_list_conv_sum_set_nat: |
281 |
"sum_list (map f [m..<n]) = sum f (set [m..<n])" |
|
282 |
by (simp add: sum_list_distinct_conv_sum_set) |
|
| 58101 | 283 |
|
| 64267 | 284 |
lemma (in monoid_add) interv_sum_list_conv_sum_set_int: |
285 |
"sum_list (map f [k..l]) = sum f (set [k..l])" |
|
286 |
by (simp add: sum_list_distinct_conv_sum_set) |
|
| 58101 | 287 |
|
| 69593 | 288 |
text \<open>General equivalence between \<^const>\<open>sum_list\<close> and \<^const>\<open>sum\<close>\<close> |
| 64267 | 289 |
lemma (in monoid_add) sum_list_sum_nth: |
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
290 |
"sum_list xs = (\<Sum> i = 0 ..< length xs. xs ! i)" |
| 67399 | 291 |
using interv_sum_list_conv_sum_set_nat [of "(!) xs" 0 "length xs"] by (simp add: map_nth) |
| 58101 | 292 |
|
| 64267 | 293 |
lemma sum_list_map_eq_sum_count: |
294 |
"sum_list (map f xs) = sum (\<lambda>x. count_list xs x * f x) (set xs)" |
|
| 59728 | 295 |
proof(induction xs) |
296 |
case (Cons x xs) |
|
297 |
show ?case (is "?l = ?r") |
|
298 |
proof cases |
|
299 |
assume "x \<in> set xs" |
|
| 60541 | 300 |
have "?l = f x + (\<Sum>x\<in>set xs. count_list xs x * f x)" by (simp add: Cons.IH) |
| 60758 | 301 |
also have "set xs = insert x (set xs - {x})" using \<open>x \<in> set xs\<close>by blast
|
| 60541 | 302 |
also have "f x + (\<Sum>x\<in>insert x (set xs - {x}). count_list xs x * f x) = ?r"
|
| 64267 | 303 |
by (simp add: sum.insert_remove eq_commute) |
| 59728 | 304 |
finally show ?thesis . |
305 |
next |
|
306 |
assume "x \<notin> set xs" |
|
307 |
hence "\<And>xa. xa \<in> set xs \<Longrightarrow> x \<noteq> xa" by blast |
|
| 60758 | 308 |
thus ?thesis by (simp add: Cons.IH \<open>x \<notin> set xs\<close>) |
| 59728 | 309 |
qed |
310 |
qed simp |
|
311 |
||
| 64267 | 312 |
lemma sum_list_map_eq_sum_count2: |
| 59728 | 313 |
assumes "set xs \<subseteq> X" "finite X" |
| 64267 | 314 |
shows "sum_list (map f xs) = sum (\<lambda>x. count_list xs x * f x) X" |
| 59728 | 315 |
proof- |
| 60541 | 316 |
let ?F = "\<lambda>x. count_list xs x * f x" |
| 64267 | 317 |
have "sum ?F X = sum ?F (set xs \<union> (X - set xs))" |
| 59728 | 318 |
using Un_absorb1[OF assms(1)] by(simp) |
| 64267 | 319 |
also have "\<dots> = sum ?F (set xs)" |
| 59728 | 320 |
using assms(2) |
| 64267 | 321 |
by(simp add: sum.union_disjoint[OF _ _ Diff_disjoint] del: Un_Diff_cancel) |
322 |
finally show ?thesis by(simp add:sum_list_map_eq_sum_count) |
|
| 59728 | 323 |
qed |
324 |
||
| 72545 | 325 |
lemma sum_list_replicate: "sum_list (replicate n c) = of_nat n * c" |
326 |
by(induction n)(auto simp add: distrib_right) |
|
327 |
||
328 |
||
|
67489
f1ba59ddd9a6
drop redundant cong rules
Lars Hupel <lars.hupel@mytum.de>
parents:
67399
diff
changeset
|
329 |
lemma sum_list_nonneg: |
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
330 |
"(\<And>x. x \<in> set xs \<Longrightarrow> (x :: 'a :: ordered_comm_monoid_add) \<ge> 0) \<Longrightarrow> sum_list xs \<ge> 0" |
|
63099
af0e964aad7b
Moved material from AFP/Randomised_Social_Choice to distribution
eberlm
parents:
61955
diff
changeset
|
331 |
by (induction xs) simp_all |
|
af0e964aad7b
Moved material from AFP/Randomised_Social_Choice to distribution
eberlm
parents:
61955
diff
changeset
|
332 |
|
| 69231 | 333 |
lemma sum_list_Suc: |
334 |
"sum_list (map (\<lambda>x. Suc(f x)) xs) = sum_list (map f xs) + length xs" |
|
335 |
by(induction xs; simp) |
|
336 |
||
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
337 |
lemma (in monoid_add) sum_list_map_filter': |
|
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
338 |
"sum_list (map f (filter P xs)) = sum_list (map (\<lambda>x. if P x then f x else 0) xs)" |
|
63099
af0e964aad7b
Moved material from AFP/Randomised_Social_Choice to distribution
eberlm
parents:
61955
diff
changeset
|
339 |
by (induction xs) simp_all |
|
af0e964aad7b
Moved material from AFP/Randomised_Social_Choice to distribution
eberlm
parents:
61955
diff
changeset
|
340 |
|
|
67489
f1ba59ddd9a6
drop redundant cong rules
Lars Hupel <lars.hupel@mytum.de>
parents:
67399
diff
changeset
|
341 |
text \<open>Summation of a strictly ascending sequence with length \<open>n\<close> |
|
66434
5d7e770c7d5d
added sorted_wrt to List; added Data_Structures/Binomial_Heap.thy
nipkow
parents:
66311
diff
changeset
|
342 |
can be upper-bounded by summation over \<open>{0..<n}\<close>.\<close>
|
|
5d7e770c7d5d
added sorted_wrt to List; added Data_Structures/Binomial_Heap.thy
nipkow
parents:
66311
diff
changeset
|
343 |
|
|
5d7e770c7d5d
added sorted_wrt to List; added Data_Structures/Binomial_Heap.thy
nipkow
parents:
66311
diff
changeset
|
344 |
lemma sorted_wrt_less_sum_mono_lowerbound: |
|
5d7e770c7d5d
added sorted_wrt to List; added Data_Structures/Binomial_Heap.thy
nipkow
parents:
66311
diff
changeset
|
345 |
fixes f :: "nat \<Rightarrow> ('b::ordered_comm_monoid_add)"
|
|
5d7e770c7d5d
added sorted_wrt to List; added Data_Structures/Binomial_Heap.thy
nipkow
parents:
66311
diff
changeset
|
346 |
assumes mono: "\<And>x y. x\<le>y \<Longrightarrow> f x \<le> f y" |
| 67399 | 347 |
shows "sorted_wrt (<) ns \<Longrightarrow> |
|
66434
5d7e770c7d5d
added sorted_wrt to List; added Data_Structures/Binomial_Heap.thy
nipkow
parents:
66311
diff
changeset
|
348 |
(\<Sum>i\<in>{0..<length ns}. f i) \<le> (\<Sum>i\<leftarrow>ns. f i)"
|
|
5d7e770c7d5d
added sorted_wrt to List; added Data_Structures/Binomial_Heap.thy
nipkow
parents:
66311
diff
changeset
|
349 |
proof (induction ns rule: rev_induct) |
|
5d7e770c7d5d
added sorted_wrt to List; added Data_Structures/Binomial_Heap.thy
nipkow
parents:
66311
diff
changeset
|
350 |
case Nil |
|
5d7e770c7d5d
added sorted_wrt to List; added Data_Structures/Binomial_Heap.thy
nipkow
parents:
66311
diff
changeset
|
351 |
then show ?case by simp |
|
5d7e770c7d5d
added sorted_wrt to List; added Data_Structures/Binomial_Heap.thy
nipkow
parents:
66311
diff
changeset
|
352 |
next |
|
5d7e770c7d5d
added sorted_wrt to List; added Data_Structures/Binomial_Heap.thy
nipkow
parents:
66311
diff
changeset
|
353 |
case (snoc n ns) |
|
67489
f1ba59ddd9a6
drop redundant cong rules
Lars Hupel <lars.hupel@mytum.de>
parents:
67399
diff
changeset
|
354 |
have "sum f {0..<length (ns @ [n])}
|
|
f1ba59ddd9a6
drop redundant cong rules
Lars Hupel <lars.hupel@mytum.de>
parents:
67399
diff
changeset
|
355 |
= sum f {0..<length ns} + f (length ns)"
|
|
66434
5d7e770c7d5d
added sorted_wrt to List; added Data_Structures/Binomial_Heap.thy
nipkow
parents:
66311
diff
changeset
|
356 |
by simp |
|
5d7e770c7d5d
added sorted_wrt to List; added Data_Structures/Binomial_Heap.thy
nipkow
parents:
66311
diff
changeset
|
357 |
also have "sum f {0..<length ns} \<le> sum_list (map f ns)"
|
|
5d7e770c7d5d
added sorted_wrt to List; added Data_Structures/Binomial_Heap.thy
nipkow
parents:
66311
diff
changeset
|
358 |
using snoc by (auto simp: sorted_wrt_append) |
|
5d7e770c7d5d
added sorted_wrt to List; added Data_Structures/Binomial_Heap.thy
nipkow
parents:
66311
diff
changeset
|
359 |
also have "length ns \<le> n" |
|
67489
f1ba59ddd9a6
drop redundant cong rules
Lars Hupel <lars.hupel@mytum.de>
parents:
67399
diff
changeset
|
360 |
using sorted_wrt_less_idx[OF snoc.prems(1), of "length ns"] by auto |
|
66434
5d7e770c7d5d
added sorted_wrt to List; added Data_Structures/Binomial_Heap.thy
nipkow
parents:
66311
diff
changeset
|
361 |
finally have "sum f {0..<length (ns @ [n])} \<le> sum_list (map f ns) + f n"
|
|
5d7e770c7d5d
added sorted_wrt to List; added Data_Structures/Binomial_Heap.thy
nipkow
parents:
66311
diff
changeset
|
362 |
using mono add_mono by blast |
|
5d7e770c7d5d
added sorted_wrt to List; added Data_Structures/Binomial_Heap.thy
nipkow
parents:
66311
diff
changeset
|
363 |
thus ?case by simp |
|
67489
f1ba59ddd9a6
drop redundant cong rules
Lars Hupel <lars.hupel@mytum.de>
parents:
67399
diff
changeset
|
364 |
qed |
|
66434
5d7e770c7d5d
added sorted_wrt to List; added Data_Structures/Binomial_Heap.thy
nipkow
parents:
66311
diff
changeset
|
365 |
|
| 58101 | 366 |
|
| 72024 | 367 |
subsection \<open>Horner sums\<close> |
368 |
||
369 |
context comm_semiring_0 |
|
370 |
begin |
|
371 |
||
372 |
definition horner_sum :: \<open>('b \<Rightarrow> 'a) \<Rightarrow> 'a \<Rightarrow> 'b list \<Rightarrow> 'a\<close>
|
|
373 |
where horner_sum_foldr: \<open>horner_sum f a xs = foldr (\<lambda>x b. f x + a * b) xs 0\<close> |
|
374 |
||
375 |
lemma horner_sum_simps [simp]: |
|
376 |
\<open>horner_sum f a [] = 0\<close> |
|
377 |
\<open>horner_sum f a (x # xs) = f x + a * horner_sum f a xs\<close> |
|
378 |
by (simp_all add: horner_sum_foldr) |
|
379 |
||
380 |
lemma horner_sum_eq_sum_funpow: |
|
381 |
\<open>horner_sum f a xs = (\<Sum>n = 0..<length xs. ((*) a ^^ n) (f (xs ! n)))\<close> |
|
382 |
proof (induction xs) |
|
383 |
case Nil |
|
384 |
then show ?case |
|
385 |
by simp |
|
386 |
next |
|
387 |
case (Cons x xs) |
|
388 |
then show ?case |
|
389 |
by (simp add: sum.atLeast0_lessThan_Suc_shift sum_distrib_left del: sum.op_ivl_Suc) |
|
390 |
qed |
|
391 |
||
392 |
end |
|
393 |
||
394 |
context |
|
395 |
includes lifting_syntax |
|
396 |
begin |
|
397 |
||
398 |
lemma horner_sum_transfer [transfer_rule]: |
|
399 |
\<open>((B ===> A) ===> A ===> list_all2 B ===> A) horner_sum horner_sum\<close> |
|
400 |
if [transfer_rule]: \<open>A 0 0\<close> |
|
401 |
and [transfer_rule]: \<open>(A ===> A ===> A) (+) (+)\<close> |
|
402 |
and [transfer_rule]: \<open>(A ===> A ===> A) (*) (*)\<close> |
|
403 |
by (unfold horner_sum_foldr) transfer_prover |
|
404 |
||
405 |
end |
|
406 |
||
407 |
context comm_semiring_1 |
|
408 |
begin |
|
409 |
||
410 |
lemma horner_sum_eq_sum: |
|
411 |
\<open>horner_sum f a xs = (\<Sum>n = 0..<length xs. f (xs ! n) * a ^ n)\<close> |
|
412 |
proof - |
|
413 |
have \<open>(*) a ^^ n = (*) (a ^ n)\<close> for n |
|
414 |
by (induction n) (simp_all add: ac_simps) |
|
415 |
then show ?thesis |
|
416 |
by (simp add: horner_sum_eq_sum_funpow ac_simps) |
|
417 |
qed |
|
418 |
||
| 72619 | 419 |
lemma horner_sum_append: |
420 |
\<open>horner_sum f a (xs @ ys) = horner_sum f a xs + a ^ length xs * horner_sum f a ys\<close> |
|
421 |
using sum.atLeastLessThan_shift_bounds [of _ 0 \<open>length xs\<close> \<open>length ys\<close>] |
|
422 |
atLeastLessThan_add_Un [of 0 \<open>length xs\<close> \<open>length ys\<close>] |
|
423 |
by (simp add: horner_sum_eq_sum sum_distrib_left sum.union_disjoint ac_simps nth_append power_add) |
|
424 |
||
| 72024 | 425 |
end |
426 |
||
| 75662 | 427 |
context linordered_semidom |
428 |
begin |
|
429 |
||
430 |
lemma horner_sum_nonnegative: |
|
431 |
\<open>0 \<le> horner_sum of_bool 2 bs\<close> |
|
432 |
by (induction bs) simp_all |
|
433 |
||
434 |
end |
|
435 |
||
|
78935
5e788ff7a489
explicit type class for discrete linordered semidoms
haftmann
parents:
75693
diff
changeset
|
436 |
context discrete_linordered_semidom |
| 75662 | 437 |
begin |
438 |
||
439 |
lemma horner_sum_bound: |
|
440 |
\<open>horner_sum of_bool 2 bs < 2 ^ length bs\<close> |
|
441 |
proof (induction bs) |
|
442 |
case Nil |
|
443 |
then show ?case |
|
444 |
by simp |
|
445 |
next |
|
446 |
case (Cons b bs) |
|
447 |
moreover define a where \<open>a = 2 ^ length bs - horner_sum of_bool 2 bs\<close> |
|
448 |
ultimately have *: \<open>2 ^ length bs = horner_sum of_bool 2 bs + a\<close> |
|
449 |
by simp |
|
|
78935
5e788ff7a489
explicit type class for discrete linordered semidoms
haftmann
parents:
75693
diff
changeset
|
450 |
have \<open>0 < a\<close> |
|
5e788ff7a489
explicit type class for discrete linordered semidoms
haftmann
parents:
75693
diff
changeset
|
451 |
using Cons * by simp |
|
5e788ff7a489
explicit type class for discrete linordered semidoms
haftmann
parents:
75693
diff
changeset
|
452 |
moreover have \<open>1 \<le> a\<close> |
|
5e788ff7a489
explicit type class for discrete linordered semidoms
haftmann
parents:
75693
diff
changeset
|
453 |
using \<open>0 < a\<close> by (simp add: less_eq_iff_succ_less) |
|
5e788ff7a489
explicit type class for discrete linordered semidoms
haftmann
parents:
75693
diff
changeset
|
454 |
ultimately have \<open>0 + 1 < a + a\<close> |
|
5e788ff7a489
explicit type class for discrete linordered semidoms
haftmann
parents:
75693
diff
changeset
|
455 |
by (rule add_less_le_mono) |
|
5e788ff7a489
explicit type class for discrete linordered semidoms
haftmann
parents:
75693
diff
changeset
|
456 |
then have \<open>1 < a * 2\<close> |
|
5e788ff7a489
explicit type class for discrete linordered semidoms
haftmann
parents:
75693
diff
changeset
|
457 |
by (simp add: mult_2_right) |
| 75662 | 458 |
with Cons show ?case |
|
78935
5e788ff7a489
explicit type class for discrete linordered semidoms
haftmann
parents:
75693
diff
changeset
|
459 |
by (simp add: * algebra_simps) |
| 75662 | 460 |
qed |
461 |
||
| 79017 | 462 |
lemma horner_sum_of_bool_2_less: |
463 |
\<open>(horner_sum of_bool 2 bs) < 2 ^ length bs\<close> |
|
464 |
by (fact horner_sum_bound) |
|
465 |
||
| 75662 | 466 |
end |
467 |
||
468 |
lemma nat_horner_sum [simp]: |
|
469 |
\<open>nat (horner_sum of_bool 2 bs) = horner_sum of_bool 2 bs\<close> |
|
470 |
by (induction bs) (auto simp add: nat_add_distrib horner_sum_nonnegative) |
|
471 |
||
|
78935
5e788ff7a489
explicit type class for discrete linordered semidoms
haftmann
parents:
75693
diff
changeset
|
472 |
context discrete_linordered_semidom |
| 75662 | 473 |
begin |
474 |
||
475 |
lemma horner_sum_less_eq_iff_lexordp_eq: |
|
476 |
\<open>horner_sum of_bool 2 bs \<le> horner_sum of_bool 2 cs \<longleftrightarrow> lexordp_eq (rev bs) (rev cs)\<close> |
|
477 |
if \<open>length bs = length cs\<close> |
|
478 |
proof - |
|
479 |
have \<open>horner_sum of_bool 2 (rev bs) \<le> horner_sum of_bool 2 (rev cs) \<longleftrightarrow> lexordp_eq bs cs\<close> |
|
480 |
if \<open>length bs = length cs\<close> for bs cs |
|
481 |
using that proof (induction bs cs rule: list_induct2) |
|
482 |
case Nil |
|
483 |
then show ?case |
|
484 |
by simp |
|
485 |
next |
|
486 |
case (Cons b bs c cs) |
|
487 |
with horner_sum_nonnegative [of \<open>rev bs\<close>] horner_sum_nonnegative [of \<open>rev cs\<close>] |
|
488 |
horner_sum_bound [of \<open>rev bs\<close>] horner_sum_bound [of \<open>rev cs\<close>] |
|
489 |
show ?case |
|
490 |
by (auto simp add: horner_sum_append not_le Cons intro: add_strict_increasing2 add_increasing) |
|
491 |
qed |
|
492 |
from that this [of \<open>rev bs\<close> \<open>rev cs\<close>] show ?thesis |
|
493 |
by simp |
|
494 |
qed |
|
495 |
||
496 |
lemma horner_sum_less_iff_lexordp: |
|
497 |
\<open>horner_sum of_bool 2 bs < horner_sum of_bool 2 cs \<longleftrightarrow> ord_class.lexordp (rev bs) (rev cs)\<close> |
|
498 |
if \<open>length bs = length cs\<close> |
|
499 |
proof - |
|
500 |
have \<open>horner_sum of_bool 2 (rev bs) < horner_sum of_bool 2 (rev cs) \<longleftrightarrow> ord_class.lexordp bs cs\<close> |
|
501 |
if \<open>length bs = length cs\<close> for bs cs |
|
502 |
using that proof (induction bs cs rule: list_induct2) |
|
503 |
case Nil |
|
504 |
then show ?case |
|
505 |
by simp |
|
506 |
next |
|
507 |
case (Cons b bs c cs) |
|
508 |
with horner_sum_nonnegative [of \<open>rev bs\<close>] horner_sum_nonnegative [of \<open>rev cs\<close>] |
|
509 |
horner_sum_bound [of \<open>rev bs\<close>] horner_sum_bound [of \<open>rev cs\<close>] |
|
510 |
show ?case |
|
511 |
by (auto simp add: horner_sum_append not_less Cons intro: add_strict_increasing2 add_increasing) |
|
512 |
qed |
|
513 |
from that this [of \<open>rev bs\<close> \<open>rev cs\<close>] show ?thesis |
|
514 |
by simp |
|
515 |
qed |
|
516 |
||
517 |
end |
|
518 |
||
| 72024 | 519 |
|
| 69593 | 520 |
subsection \<open>Further facts about \<^const>\<open>List.n_lists\<close>\<close> |
| 58101 | 521 |
|
522 |
lemma length_n_lists: "length (List.n_lists n xs) = length xs ^ n" |
|
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
523 |
by (induct n) (auto simp add: comp_def length_concat sum_list_triv) |
| 58101 | 524 |
|
525 |
lemma distinct_n_lists: |
|
526 |
assumes "distinct xs" |
|
527 |
shows "distinct (List.n_lists n xs)" |
|
528 |
proof (rule card_distinct) |
|
529 |
from assms have card_length: "card (set xs) = length xs" by (rule distinct_card) |
|
530 |
have "card (set (List.n_lists n xs)) = card (set xs) ^ n" |
|
531 |
proof (induct n) |
|
532 |
case 0 then show ?case by simp |
|
533 |
next |
|
534 |
case (Suc n) |
|
535 |
moreover have "card (\<Union>ys\<in>set (List.n_lists n xs). (\<lambda>y. y # ys) ` set xs) |
|
536 |
= (\<Sum>ys\<in>set (List.n_lists n xs). card ((\<lambda>y. y # ys) ` set xs))" |
|
537 |
by (rule card_UN_disjoint) auto |
|
538 |
moreover have "\<And>ys. card ((\<lambda>y. y # ys) ` set xs) = card (set xs)" |
|
539 |
by (rule card_image) (simp add: inj_on_def) |
|
540 |
ultimately show ?case by auto |
|
541 |
qed |
|
542 |
also have "\<dots> = length xs ^ n" by (simp add: card_length) |
|
543 |
finally show "card (set (List.n_lists n xs)) = length (List.n_lists n xs)" |
|
544 |
by (simp add: length_n_lists) |
|
545 |
qed |
|
546 |
||
547 |
||
| 60758 | 548 |
subsection \<open>Tools setup\<close> |
| 58101 | 549 |
|
| 64267 | 550 |
lemmas sum_code = sum.set_conv_list |
| 58320 | 551 |
|
| 64267 | 552 |
lemma sum_set_upto_conv_sum_list_int [code_unfold]: |
553 |
"sum f (set [i..j::int]) = sum_list (map f [i..j])" |
|
554 |
by (simp add: interv_sum_list_conv_sum_set_int) |
|
| 58101 | 555 |
|
| 64267 | 556 |
lemma sum_set_upt_conv_sum_list_nat [code_unfold]: |
557 |
"sum f (set [m..<n]) = sum_list (map f [m..<n])" |
|
558 |
by (simp add: interv_sum_list_conv_sum_set_nat) |
|
| 58101 | 559 |
|
| 58368 | 560 |
|
| 60758 | 561 |
subsection \<open>List product\<close> |
| 58368 | 562 |
|
563 |
context monoid_mult |
|
564 |
begin |
|
565 |
||
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
566 |
sublocale prod_list: monoid_list times 1 |
| 61776 | 567 |
defines |
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
568 |
prod_list = prod_list.F .. |
| 58368 | 569 |
|
| 58320 | 570 |
end |
| 58368 | 571 |
|
572 |
context comm_monoid_mult |
|
573 |
begin |
|
574 |
||
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
575 |
sublocale prod_list: comm_monoid_list times 1 |
|
61566
c3d6e570ccef
Keyword 'rewrites' identifies rewrite morphisms.
ballarin
parents:
61378
diff
changeset
|
576 |
rewrites |
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
577 |
"monoid_list.F times 1 = prod_list" |
| 58368 | 578 |
proof - |
579 |
show "comm_monoid_list times 1" .. |
|
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
580 |
then interpret prod_list: comm_monoid_list times 1 . |
|
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
581 |
from prod_list_def show "monoid_list.F times 1 = prod_list" by simp |
| 58368 | 582 |
qed |
583 |
||
| 64272 | 584 |
sublocale prod: comm_monoid_list_set times 1 |
|
61566
c3d6e570ccef
Keyword 'rewrites' identifies rewrite morphisms.
ballarin
parents:
61378
diff
changeset
|
585 |
rewrites |
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
586 |
"monoid_list.F times 1 = prod_list" |
| 64272 | 587 |
and "comm_monoid_set.F times 1 = prod" |
| 58368 | 588 |
proof - |
589 |
show "comm_monoid_list_set times 1" .. |
|
| 64272 | 590 |
then interpret prod: comm_monoid_list_set times 1 . |
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
591 |
from prod_list_def show "monoid_list.F times 1 = prod_list" by simp |
| 64272 | 592 |
from prod_def show "comm_monoid_set.F times 1 = prod" by (auto intro: sym) |
| 58368 | 593 |
qed |
594 |
||
595 |
end |
|
596 |
||
| 60758 | 597 |
text \<open>Some syntactic sugar:\<close> |
| 58368 | 598 |
|
|
61955
e96292f32c3c
former "xsymbols" syntax is used by default, and ASCII replacement syntax with print mode "ASCII";
wenzelm
parents:
61799
diff
changeset
|
599 |
syntax (ASCII) |
| 80934 | 600 |
"_prod_list" :: "pttrn => 'a list => 'b => 'b" (\<open>(\<open>indent=3 notation=\<open>binder PROD\<close>\<close>PROD _<-_. _)\<close> [0, 51, 10] 10) |
|
61955
e96292f32c3c
former "xsymbols" syntax is used by default, and ASCII replacement syntax with print mode "ASCII";
wenzelm
parents:
61799
diff
changeset
|
601 |
syntax |
| 80934 | 602 |
"_prod_list" :: "pttrn => 'a list => 'b => 'b" (\<open>(\<open>indent=3 notation=\<open>binder \<Prod>\<close>\<close>\<Prod>_\<leftarrow>_. _)\<close> [0, 51, 10] 10) |
| 80760 | 603 |
syntax_consts |
604 |
"_prod_list" \<rightleftharpoons> prod_list |
|
| 61799 | 605 |
translations \<comment> \<open>Beware of argument permutation!\<close> |
|
63882
018998c00003
renamed listsum -> sum_list, listprod ~> prod_list
nipkow
parents:
63343
diff
changeset
|
606 |
"\<Prod>x\<leftarrow>xs. b" \<rightleftharpoons> "CONST prod_list (CONST map (\<lambda>x. b) xs)" |
| 58368 | 607 |
|
| 70928 | 608 |
context |
609 |
includes lifting_syntax |
|
610 |
begin |
|
611 |
||
612 |
lemma prod_list_transfer [transfer_rule]: |
|
613 |
"(list_all2 A ===> A) prod_list prod_list" |
|
614 |
if [transfer_rule]: "A 1 1" "(A ===> A ===> A) (*) (*)" |
|
615 |
unfolding prod_list.eq_foldr [abs_def] |
|
616 |
by transfer_prover |
|
617 |
||
| 58368 | 618 |
end |
| 70928 | 619 |
|
620 |
lemma prod_list_zero_iff: |
|
621 |
"prod_list xs = 0 \<longleftrightarrow> (0 :: 'a :: {semiring_no_zero_divisors, semiring_1}) \<in> set xs"
|
|
622 |
by (induction xs) simp_all |
|
623 |
||
624 |
end |