Theorem eucalg 12013

Description: Euclid's Algorithm computes the greatest common divisor of two nonnegative integers by repeatedly replacing the larger of them with its remainder modulo the smaller until the remainder is 0. Theorem 1.15 in [ApostolNT] p. 20.

Upon halting, the 1st member of the final state ( R ` N ) is equal to the gcd of the values comprising the input state <. M , N >.. This is Metamath 100 proof #69 (greatest common divisor algorithm). (Contributed by Paul Chapman, 31-Mar-2011.) (Proof shortened by Mario Carneiro, 29-May-2014.)

Hypotheses

Ref	Expression
eucalgval.1	|- E = ( x e. NN0 , y e. NN0 |-> if ( y = 0 , <. x , y >. , <. y , ( x mod y ) >. ) )
eucalg.2	|- R = seq 0 ( ( E o. 1st ) , ( NN0 X. { A } ) )
eucalg.3	|- A = <. M , N >.

Assertion

Ref	Expression
eucalg	|- ( ( M e. NN0 /\ N e. NN0 ) -> ( 1st ` ( R ` N ) ) = ( M gcd N ) )

Distinct variable groups:   x, y, M   x, N, y   x, A, y   x, R

Allowed substitution hints:   R( y)   E( x, y)

Proof of Theorem eucalg

Dummy variable z is distinct from all other variables.

Step	Hyp	Ref	Expression
1		nn0uz 9521	. . . . . . . 8 |- NN0 = ( ZZ>= ` 0 )
2		eucalg.2	. . . . . . . 8 |- R = seq 0 ( ( E o. 1st ) , ( NN0 X. { A } ) )
3		0zd 9224	. . . . . . . 8 |- ( ( M e. NN0 /\ N e. NN0 ) -> 0 e. ZZ )
4		eucalg.3	. . . . . . . . 9 |- A = <. M , N >.
5		opelxpi 4643	. . . . . . . . 9 |- ( ( M e. NN0 /\ N e. NN0 ) -> <. M , N >. e. ( NN0 X. NN0 ) )
6	4, 5	eqeltrid 2257	. . . . . . . 8 |- ( ( M e. NN0 /\ N e. NN0 ) -> A e. ( NN0 X. NN0 ) )
7		eucalgval.1	. . . . . . . . . 10 |- E = ( x e. NN0 , y e. NN0 |-> if ( y = 0 , <. x , y >. , <. y , ( x mod y ) >. ) )
8	7	eucalgf 12009	. . . . . . . . 9 |- E : ( NN0 X. NN0 ) --> ( NN0 X. NN0 )
9	8	a1i 9	. . . . . . . 8 |- ( ( M e. NN0 /\ N e. NN0 ) -> E : ( NN0 X. NN0 ) --> ( NN0 X. NN0 ) )
10	1, 2, 3, 6, 9	algrf 11999	. . . . . . 7 |- ( ( M e. NN0 /\ N e. NN0 ) -> R : NN0 --> ( NN0 X. NN0 ) )
11		ffvelrn 5629	. . . . . . 7 |- ( ( R : NN0 --> ( NN0 X. NN0 ) /\ N e. NN0 ) -> ( R ` N ) e. ( NN0 X. NN0 ) )
12	10, 11	sylancom 418	. . . . . 6 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( R ` N ) e. ( NN0 X. NN0 ) )
13		1st2nd2 6154	. . . . . 6 |- ( ( R ` N ) e. ( NN0 X. NN0 ) -> ( R ` N ) = <. ( 1st ` ( R ` N ) ) , ( 2nd ` ( R ` N ) ) >. )
14	12, 13	syl 14	. . . . 5 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( R ` N ) = <. ( 1st ` ( R ` N ) ) , ( 2nd ` ( R ` N ) ) >. )
15	14	fveq2d 5500	. . . 4 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( gcd ` ( R ` N ) ) = ( gcd ` <. ( 1st ` ( R ` N ) ) , ( 2nd ` ( R ` N ) ) >. ) )
16		df-ov 5856	. . . 4 |- ( ( 1st ` ( R ` N ) ) gcd ( 2nd ` ( R ` N ) ) ) = ( gcd ` <. ( 1st ` ( R ` N ) ) , ( 2nd ` ( R ` N ) ) >. )
17	15, 16	eqtr4di 2221	. . 3 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( gcd ` ( R ` N ) ) = ( ( 1st ` ( R ` N ) ) gcd ( 2nd ` ( R ` N ) ) ) )
18	4	fveq2i 5499	. . . . . . . 8 |- ( 2nd ` A ) = ( 2nd ` <. M , N >. )
19		op2ndg 6130	. . . . . . . 8 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( 2nd ` <. M , N >. ) = N )
20	18, 19	eqtrid 2215	. . . . . . 7 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( 2nd ` A ) = N )
21	20	fveq2d 5500	. . . . . 6 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( R ` ( 2nd ` A ) ) = ( R ` N ) )
22	21	fveq2d 5500	. . . . 5 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( 2nd ` ( R ` ( 2nd ` A ) ) ) = ( 2nd ` ( R ` N ) ) )
23		xp2nd 6145	. . . . . . . . 9 |- ( A e. ( NN0 X. NN0 ) -> ( 2nd ` A ) e. NN0 )
24	23	nn0zd 9332	. . . . . . . 8 |- ( A e. ( NN0 X. NN0 ) -> ( 2nd ` A ) e. ZZ )
25		uzid 9501	. . . . . . . 8 |- ( ( 2nd ` A ) e. ZZ -> ( 2nd ` A ) e. ( ZZ>= ` ( 2nd ` A ) ) )
26	24, 25	syl 14	. . . . . . 7 |- ( A e. ( NN0 X. NN0 ) -> ( 2nd ` A ) e. ( ZZ>= ` ( 2nd ` A ) ) )
27		eqid 2170	. . . . . . . 8 |- ( 2nd ` A ) = ( 2nd ` A )
28	7, 2, 27	eucalgcvga 12012	. . . . . . 7 |- ( A e. ( NN0 X. NN0 ) -> ( ( 2nd ` A ) e. ( ZZ>= ` ( 2nd ` A ) ) -> ( 2nd ` ( R ` ( 2nd ` A ) ) ) = 0 ) )
29	26, 28	mpd 13	. . . . . 6 |- ( A e. ( NN0 X. NN0 ) -> ( 2nd ` ( R ` ( 2nd ` A ) ) ) = 0 )
30	6, 29	syl 14	. . . . 5 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( 2nd ` ( R ` ( 2nd ` A ) ) ) = 0 )
31	22, 30	eqtr3d 2205	. . . 4 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( 2nd ` ( R ` N ) ) = 0 )
32	31	oveq2d 5869	. . 3 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( ( 1st ` ( R ` N ) ) gcd ( 2nd ` ( R ` N ) ) ) = ( ( 1st ` ( R ` N ) ) gcd 0 ) )
33		xp1st 6144	. . . 4 |- ( ( R ` N ) e. ( NN0 X. NN0 ) -> ( 1st ` ( R ` N ) ) e. NN0 )
34		nn0gcdid0 11936	. . . 4 |- ( ( 1st ` ( R ` N ) ) e. NN0 -> ( ( 1st ` ( R ` N ) ) gcd 0 ) = ( 1st ` ( R ` N ) ) )
35	12, 33, 34	3syl 17	. . 3 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( ( 1st ` ( R ` N ) ) gcd 0 ) = ( 1st ` ( R ` N ) ) )
36	17, 32, 35	3eqtrrd 2208	. 2 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( 1st ` ( R ` N ) ) = ( gcd ` ( R ` N ) ) )
37	7	eucalginv 12010	. . . . . 6 |- ( z e. ( NN0 X. NN0 ) -> ( gcd ` ( E ` z ) ) = ( gcd ` z ) )
38	8	ffvelrni 5630	. . . . . . 7 |- ( z e. ( NN0 X. NN0 ) -> ( E ` z ) e. ( NN0 X. NN0 ) )
39		fvres 5520	. . . . . . 7 |- ( ( E ` z ) e. ( NN0 X. NN0 ) -> ( ( gcd |` ( NN0 X. NN0 ) ) ` ( E ` z ) ) = ( gcd ` ( E ` z ) ) )
40	38, 39	syl 14	. . . . . 6 |- ( z e. ( NN0 X. NN0 ) -> ( ( gcd |` ( NN0 X. NN0 ) ) ` ( E ` z ) ) = ( gcd ` ( E ` z ) ) )
41		fvres 5520	. . . . . 6 |- ( z e. ( NN0 X. NN0 ) -> ( ( gcd |` ( NN0 X. NN0 ) ) ` z ) = ( gcd ` z ) )
42	37, 40, 41	3eqtr4d 2213	. . . . 5 |- ( z e. ( NN0 X. NN0 ) -> ( ( gcd |` ( NN0 X. NN0 ) ) ` ( E ` z ) ) = ( ( gcd |` ( NN0 X. NN0 ) ) ` z ) )
43	2, 8, 42	alginv 12001	. . . 4 |- ( ( A e. ( NN0 X. NN0 ) /\ N e. NN0 ) -> ( ( gcd |` ( NN0 X. NN0 ) ) ` ( R ` N ) ) = ( ( gcd |` ( NN0 X. NN0 ) ) ` ( R ` 0 ) ) )
44	6, 43	sylancom 418	. . 3 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( ( gcd |` ( NN0 X. NN0 ) ) ` ( R ` N ) ) = ( ( gcd |` ( NN0 X. NN0 ) ) ` ( R ` 0 ) ) )
45		fvres 5520	. . . 4 |- ( ( R ` N ) e. ( NN0 X. NN0 ) -> ( ( gcd |` ( NN0 X. NN0 ) ) ` ( R ` N ) ) = ( gcd ` ( R ` N ) ) )
46	12, 45	syl 14	. . 3 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( ( gcd |` ( NN0 X. NN0 ) ) ` ( R ` N ) ) = ( gcd ` ( R ` N ) ) )
47		0nn0 9150	. . . . 5 |- 0 e. NN0
48		ffvelrn 5629	. . . . 5 |- ( ( R : NN0 --> ( NN0 X. NN0 ) /\ 0 e. NN0 ) -> ( R ` 0 ) e. ( NN0 X. NN0 ) )
49	10, 47, 48	sylancl 411	. . . 4 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( R ` 0 ) e. ( NN0 X. NN0 ) )
50		fvres 5520	. . . 4 |- ( ( R ` 0 ) e. ( NN0 X. NN0 ) -> ( ( gcd |` ( NN0 X. NN0 ) ) ` ( R ` 0 ) ) = ( gcd ` ( R ` 0 ) ) )
51	49, 50	syl 14	. . 3 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( ( gcd |` ( NN0 X. NN0 ) ) ` ( R ` 0 ) ) = ( gcd ` ( R ` 0 ) ) )
52	44, 46, 51	3eqtr3d 2211	. 2 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( gcd ` ( R ` N ) ) = ( gcd ` ( R ` 0 ) ) )
53	1, 2, 3, 6, 9	ialgr0 11998	. . . . 5 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( R ` 0 ) = A )
54	53, 4	eqtrdi 2219	. . . 4 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( R ` 0 ) = <. M , N >. )
55	54	fveq2d 5500	. . 3 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( gcd ` ( R ` 0 ) ) = ( gcd ` <. M , N >. ) )
56		df-ov 5856	. . 3 |- ( M gcd N ) = ( gcd ` <. M , N >. )
57	55, 56	eqtr4di 2221	. 2 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( gcd ` ( R ` 0 ) ) = ( M gcd N ) )
58	36, 52, 57	3eqtrd 2207	1 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( 1st ` ( R ` N ) ) = ( M gcd N ) )

Syntax hints:   -> wi 4   /\ wa 103   = wceq 1348   e. wcel 2141   ifcif 3526   {csn 3583   <.cop 3586   X. cxp 4609   |` cres 4613   o. ccom 4615   -->wf 5194   ` cfv 5198  (class class class)co 5853   e. cmpo 5855   1stc1st 6117   2ndc2nd 6118   0cc0 7774   NN0cn0 9135   ZZcz 9212   ZZ>=cuz 9487   mod cmo 10278   seqcseq 10401   gcd cgcd 11897

This theorem was proved from axioms:  ax-mp 5  ax-1 6  ax-2 7  ax-ia1 105  ax-ia2 106  ax-ia3 107  ax-in1 609  ax-in2 610  ax-io 704  ax-5 1440  ax-7 1441  ax-gen 1442  ax-ie1 1486  ax-ie2 1487  ax-8 1497  ax-10 1498  ax-11 1499  ax-i12 1500  ax-bndl 1502  ax-4 1503  ax-17 1519  ax-i9 1523  ax-ial 1527  ax-i5r 1528  ax-13 2143  ax-14 2144  ax-ext 2152  ax-coll 4104  ax-sep 4107  ax-nul 4115  ax-pow 4160  ax-pr 4194  ax-un 4418  ax-setind 4521  ax-iinf 4572  ax-cnex 7865  ax-resscn 7866  ax-1cn 7867  ax-1re 7868  ax-icn 7869  ax-addcl 7870  ax-addrcl 7871  ax-mulcl 7872  ax-mulrcl 7873  ax-addcom 7874  ax-mulcom 7875  ax-addass 7876  ax-mulass 7877  ax-distr 7878  ax-i2m1 7879  ax-0lt1 7880  ax-1rid 7881  ax-0id 7882  ax-rnegex 7883  ax-precex 7884  ax-cnre 7885  ax-pre-ltirr 7886  ax-pre-ltwlin 7887  ax-pre-lttrn 7888  ax-pre-apti 7889  ax-pre-ltadd 7890  ax-pre-mulgt0 7891  ax-pre-mulext 7892  ax-arch 7893  ax-caucvg 7894

This theorem depends on definitions:  df-bi 116  df-stab 826  df-dc 830  df-3or 974  df-3an 975  df-tru 1351  df-fal 1354  df-nf 1454  df-sb 1756  df-eu 2022  df-mo 2023  df-clab 2157  df-cleq 2163  df-clel 2166  df-nfc 2301  df-ne 2341  df-nel 2436  df-ral 2453  df-rex 2454  df-reu 2455  df-rmo 2456  df-rab 2457  df-v 2732  df-sbc 2956  df-csb 3050  df-dif 3123  df-un 3125  df-in 3127  df-ss 3134  df-nul 3415  df-if 3527  df-pw 3568  df-sn 3589  df-pr 3590  df-op 3592  df-uni 3797  df-int 3832  df-iun 3875  df-br 3990  df-opab 4051  df-mpt 4052  df-tr 4088  df-id 4278  df-po 4281  df-iso 4282  df-iord 4351  df-on 4353  df-ilim 4354  df-suc 4356  df-iom 4575  df-xp 4617  df-rel 4618  df-cnv 4619  df-co 4620  df-dm 4621  df-rn 4622  df-res 4623  df-ima 4624  df-iota 5160  df-fun 5200  df-fn 5201  df-f 5202  df-f1 5203  df-fo 5204  df-f1o 5205  df-fv 5206  df-riota 5809  df-ov 5856  df-oprab 5857  df-mpo 5858  df-1st 6119  df-2nd 6120  df-recs 6284  df-frec 6370  df-sup 6961  df-pnf 7956  df-mnf 7957  df-xr 7958  df-ltxr 7959  df-le 7960  df-sub 8092  df-neg 8093  df-reap 8494  df-ap 8501  df-div 8590  df-inn 8879  df-2 8937  df-3 8938  df-4 8939  df-n0 9136  df-z 9213  df-uz 9488  df-q 9579  df-rp 9611  df-fz 9966  df-fzo 10099  df-fl 10226  df-mod 10279  df-seqfrec 10402  df-exp 10476  df-cj 10806  df-re 10807  df-im 10808  df-rsqrt 10962  df-abs 10963  df-dvds 11750  df-gcd 11898

This theorem is referenced by: (None)