Section contents.
Variables (ee: E) (ol: list ee) (i j: Index ee ol).
Variable n: nat.
Variable v: Vector.t (Index ee ol) (S n).
Variable iltj: (i < j)%nat.
Let flt x0 cr := filter (fun f: Index ee ol => unsum_bool (cmp_cmp (Ecmp (UE ee ol) f (vec.nth v x0)) cr)) (vec.to_list (vec.remove v x0)).
Variable IH: forall (x0: natBelow (S n)) (cr: comparison) (b: nat),
IndexSeq b (flt x0 cr) ->
monoid_expec (U.ijcount i j) (qs (U.cmp (e:=ee) (ol:=ol)) U.pick (flt x0 cr)) <= 2 / INR (S (j - i)).
Variables (b: nat) (is: IndexSeq b (vec.to_list v)).
apply IndexSeq_NoDup with b.
assumption.
Qed.
Variable iin: In i (vec.to_list v).
Variable jin: In j (vec.to_list v).
Variable pi: natBelow (S n).
Proof with auto.
unfold flt.
intros.
intro.
destruct (proj1_conj (filter_In _ _ _) H). clear H.
simpl in H0.
simpl in H2.
destruct (Ecmp ee (subscript k) (subscript (vec.nth v pi))); destruct dr; auto; simpl in H2; try discriminate...
Qed.
Hint Immediate ndi_flt.
Hint Immediate U.hom_ijcount.
Hint Immediate vec.remove_perm.
Proof with auto.
intros H.
pose proof ndi as H0.
rewrite (Permutation.Permutation_sym (vec.List_Permutation (vec.remove_perm pi v))) in H0.
inversion_clear H0...
destruct (In_filter_inv _ _ _ H)...
Qed.
Proof with auto.
intros.
inversion_clear H.
apply NoDup_map'...
intros.
unfold U.unordered_nat_pair.
intro.
apply H3.
unfold Index in *.
apply natBelow_unique.
destruct (le_lt_dec x0 x); destruct (le_lt_dec y x); inversion H4; reflexivity.
Qed.
Proof with auto with real.
intros.
rewrite partition_0...
Focus 2.
intro.
rewrite H0 in H.
apply (lt_asym _ _ H iltj).
rewrite Rplus_0_l.
rewrite monoid_expec_plus... Focus 2. intros. repeat rewrite monoid_expec_plus...
rewrite monoid_expec_plus...
rewrite monoid_expec_ret...
rewrite Rplus_0_r.
rewrite sound_cmp_expec_0...
Focus 2.
left.
apply not_In_flt...
intro.
apply (lt_asym _ _ H).
symmetry in H0.
apply (IndicesCorrect _ _ H0).
rewrite Rplus_0_l.
case_eq (Ecmp ee (subscript i) (subscript (vec.nth v pi))); intro.
(* i in the Eq part *)
unfold Rdiv.
rewrite sound_cmp_expec_0...
left.
apply not_In_flt...
simpl.
simpl in H0.
rewrite H0. intro. discriminate.
(* in the Lt part, impossible *)
elimtype False.
apply (lt_irrefl i).
apply lt_trans with (vec.nth v pi)...
apply (IndicesCorrect i (vec.nth v pi))...
(* i in the Gt part. recursive sorting of upper part: *)
assert (IndexSeq b (vec.nth v pi :: (vec.to_list (vec.remove v pi)))).
apply IndexSeq_perm with (vec.to_list v)...
cset (vec.List_Permutation (vec.remove_perm pi v))...
apply (IH _ Gt (@InvIndexSeq_filterGt' ee ol (vec.nth v pi) (vec.to_list (vec.remove v pi)) b H1))...
Qed.
Proof with auto with real.
intros.
rewrite partition_0...
Focus 2.
intro.
cset (natBelow_unique _ _ H0).
subst i.
apply (lt_asym _ _ H iltj).
rewrite Rplus_0_l.
rewrite monoid_expec_plus...
rewrite monoid_expec_plus...
rewrite monoid_expec_ret...
rewrite Rplus_0_r.
rewrite (@sound_cmp_expec_0 ee ol i j (flt pi Gt))...
Focus 2.
right.
apply not_In_flt.
intro.
apply (lt_asym _ _ H).
cset (Ecmp_sym ee (subscript (vec.nth v pi)) (subscript j)).
simpl in H0.
simpl in H1.
rewrite <- H0 in H1.
simpl in H1.
apply (IndicesCorrect _ _ H1).
rewrite Rplus_0_r.
case_eq (Ecmp ee (subscript j) (subscript (vec.nth v pi))); intro.
(* j in the Eq part *)
unfold Rdiv.
rewrite sound_cmp_expec_0...
right.
apply not_In_flt...
simpl.
simpl in H0. simpl.
rewrite H0. intro. discriminate.
Focus 2.
(* in the Gt part, impossible *)
elimtype False.
cset (Ecmp_sym ee (subscript (vec.nth v pi)) (subscript j)).
rewrite H0 in H1.
simpl in H1.
apply (lt_asym _ _ H).
apply (IndicesCorrect _ _ H1)...
(* in the Lt part *)
apply IH with b...
unfold flt.
assert (length (vec.nth v pi :: (vec.to_list (vec.remove v pi))) = S n).
simpl @length. rewrite vec.length...
assert (IndexSeq b (vec.nth v pi :: (vec.to_list (vec.remove v pi)))).
apply IndexSeq_perm with (vec.to_list v)...
apply (vec.List_Permutation (vec.remove_perm pi v)).
cset (@IndexSeq_filterLt ee ol (vec.nth v pi) (S n) b (vec.nth v pi :: (vec.to_list (vec.remove v pi))) H1 H2).
simpl filter in H3.
rewrite Ecmp_refl in H3...
intros.
repeat rewrite monoid_expec_plus...
Qed.

Lemma case_C: (i < vec.nth v pi)%nat -> (vec.nth v pi < j)%nat ->
INR (U.ijcount i j (map (fun i0: Index ee ol => U.unordered_nat_pair i0 (vec.nth v pi)) (vec.to_list (vec.remove v pi)))) +
monoid_expec (U.ijcount i j)
(foo <- @U.qs ee ol (flt pi Lt);
bar <- @U.qs ee ol (flt pi Gt);
ret (foo ++ (vec.nth v pi :: flt pi Eq) ++ bar))
= 0.