(* This code is copyrighted by its authors; it is distributed under  *)
(* the terms of the LGPL license (see LICENSE and description files) *)
(* Contribution to the Coq Library   V6.3 (July 1999)                    *)
(****************************************************************************)
(*                                                                          *)
(*                                                                          *)
(*                   Solange Coupet-Grimal & Line Jakubiec                  *)
(*                                                                          *)
(*                                                                          *)
(*              Laboratoire d'Informatique de Marseille                     *)
(*               CMI-Technopole de Chateau-Gombert                          *)
(*                   39, Rue F. Joliot Curie                                *)
(*                   13453 MARSEILLE Cedex 13                               *)
(*           e-mail:{Solange. Coupet,Line. Jakubiec}@lim.univ-mrs.fr          *)
(*                                                                          *)
(*                                                                          *)
(*                                Coq V5.10                                 *)
(*                              May 30th 1996                               *)
(*                                                                          *)
(****************************************************************************)
(*                            Lib_Div_Even_Odd.v                            *)
(****************************************************************************)
Require Export Lib_Mult.
(****************************** even/odd *********************************)
Inductive even (n : nat) : Prop :=
even_intro : forall q : nat, n = 2 * q -> even n.
Inductive odd (n : nat) : Prop :=
odd_intro : forall q : nat, n = 2 * q + 1 -> odd n.
(************************* quotients and divisors ***********************)
simple induction q.
intro.
absurd (0 = n).
apply lt_O_neq; auto with arith.
rewrite (mult_n_O m); auto with arith.
auto with arith.
Qed.
intros.
rewrite (plus_n_O q).
rewrite H0.
elim plus_mult.
apply plus_lt_compat_l.
apply (no_zero_div n 2 H q H0).
Qed.
Hint Resolve lt_quotient2_n.
intros.
rewrite H0.
elim (mult_O_le b a); auto with arith.
intro.
absurd (n = 0).
unfold not in |- *; intro.
elim (lt_irrefl n).
pattern n at 1 in |- *.
rewrite H2; auto with arith.
rewrite H0.
rewrite H1; auto with arith.
Qed.
(************************************************************************)
simple induction n.
left; apply (even_intro 0 0); auto with arith.
intros y hrec.
elim hrec.
intro eveny.
right; elim eveny.
intros q hy; apply (odd_intro (S y) q).
elim hy.
elim plus_n_Sm.
auto with arith.
intro oddy; left.
elim oddy.
intros q hy; apply (even_intro (S y) (S q)).
rewrite hy.
rewrite plus_n_Sm.
elim (mult_comm (S q) 2).
elim (mult_comm q 2).
elim plus_comm.
auto with arith.
Qed.
unfold not in |- *; intros.
elim H0; elim H.
intros k1 evena k2 odda.
apply (le_Sn_n 1).
apply (less_div 1 (k1 - k2) 2); auto with arith.
rewrite mult_comm.
replace (2 * (k1 - k2)) with (2 * k1 - 2 * k2).
elim evena.
rewrite odda.
auto with arith.
elim mult_comm.
replace (2 * k1) with (k1 * 2).
2: apply mult_comm.
replace (2 * k2) with (k2 * 2).
2: apply mult_comm.
(* Rewrite (mult_minus_distr k1 k2 (S(S O))). *)
elim (mult_comm 2 k1); elim (mult_comm 2 k2).
apply sym_equal; auto with arith.
Qed.
unfold not in |- *; intros.
apply (even_odd a H0).
auto with arith.
Qed.
(********************** S,pred,plus,odd and even ************************)
intros.
elim H0; elim H.
intros.
rewrite H2; rewrite H1.
elim (mult_plus_distr_left q q0 2).
apply (even_intro (2 * (q + q0)) (q + q0)).
auto with arith.
Qed.
intros.
elim H; intros.
rewrite H0.
rewrite (plus_n_Sm (2 * q) 1).
apply (plus_even_even (2 * q) 2).
apply (even_intro (2 * q) q); auto with arith.
apply (even_intro 2 1).
auto with arith.
Qed.
intros.
elim H; intros.
rewrite H0.
elim (plus_n_Sm (2 * q) 0).
elim (pred_Sn (2 * q + 0)).
elim (plus_n_O (2 * q)).
apply (even_intro (2 * q) q); auto with arith.
Qed.
intros.
elim H0; elim H; intros.
rewrite H2; rewrite H1.
rewrite (plus_assoc (2 * q) (2 * q0) 1).
elim (mult_plus_distr_left q q0 2).
apply (odd_intro (2 * (q + q0) + 1) (q + q0)).
auto with arith.
Qed.
Hint Resolve plus_even_odd.
intros.
rewrite plus_comm.
auto with arith.
Qed.
intros.
elim H; intros.
rewrite H0.
apply (odd_intro (S (2 * q)) q).
elim (plus_n_Sm (2 * q) 0).
auto with arith.
Qed.
intros.
elim H0; elim H; intros.
rewrite H2; rewrite H1; intros.
rewrite (plus_comm (2 * q0) 1).
rewrite (plus_assoc (2 * q + 1) 1 (2 * q0)).
rewrite (plus_assoc_reverse (2 * q) 1 1).
apply (plus_even_even (2 * q + (1 + 1)) (2 * q0)).
apply (plus_even_even (2 * q) (1 + 1)).
apply (even_intro (2 * q) q); auto with arith.
apply (even_intro (1 + 1) 1); auto with arith.
apply (even_intro (2 * q0) q0); auto with arith.
Qed.
(************************** mult, even and odd ***************************)
intros.
elim H.
intros.
apply (even_intro (a * b) (q * b)).
rewrite H0.
auto with arith.
Qed.
intros.
elim H0; intros.
elim H; intros.
clear H0; clear H.
rewrite H2; rewrite H1.
rewrite (mult_plus_distr_r (2 * q0) 1 (2 * q + 1)).
apply (plus_even_odd (2 * q0 * (2 * q + 1)) (1 * (2 * q + 1))).
apply (mult_even (2 * q0) (2 * q + 1)).
apply (even_intro (2 * q0) q0); auto with arith.
apply (odd_intro (1 * (2 * q + 1)) q); auto with arith.
Qed.
(************************** div, even and odd ****************************)
Definition div (d a : nat) := exists k : nat, a = d * k.
(*************************************************************************)
intros.
elim H.
intros.
cut (even x).
intro.
elim H2.
intros q evenx.
unfold div in |- *.
exists q.
apply (mult_reg_l_bis n (d * q) 2).
auto with arith.
rewrite H1.
rewrite evenx.
elim (mult_assoc_reverse 2 d q).
rewrite (mult_comm 2 d).
auto with arith.
elim (even_or_odd x).
auto with arith.
intros.
absurd (even (d * x)).
apply odd_even.
apply mult_odd_odd; auto with arith.
cut (d * x = 2 * n); auto with arith.
intro.
apply (even_intro (d * x) n).
auto with arith.
Qed.

Lemma div_odd_odd : forall n : nat, odd n -> forall d : nat, div d n -> odd d.