Theorem eucalg 12862

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 9973	. . . . . . . 8 |- NN0 = ( ZZ>= ` 0 )
2		eucalg.2	. . . . . . . 8 |- R = seq 0 ( ( E o. 1st ) , ( NN0 X. { A } ) )
3		0zd 9667	. . . . . . . 8 |- ( ( M e. NN0 /\ N e. NN0 ) -> 0 e. ZZ )
4		eucalg.3	. . . . . . . . 9 |- A = <. M , N >.
5		opelxpi 4808	. . . . . . . . 9 |- ( ( M e. NN0 /\ N e. NN0 ) -> <. M , N >. e. ( NN0 X. NN0 ) )
6	4, 5	eqeltrid 2325	. . . . . . . 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 12858	. . . . . . . . 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 12848	. . . . . . 7 |- ( ( M e. NN0 /\ N e. NN0 ) -> R : NN0 --> ( NN0 X. NN0 ) )
11		ffvelcdm 5844	. . . . . . 7 |- ( ( R : NN0 --> ( NN0 X. NN0 ) /\ N e. NN0 ) -> ( R ` N ) e. ( NN0 X. NN0 ) )
12	10, 11	sylancom 424	. . . . . 6 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( R ` N ) e. ( NN0 X. NN0 ) )
13		1st2nd2 6413	. . . . . 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 5702	. . . 4 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( gcd ` ( R ` N ) ) = ( gcd ` <. ( 1st ` ( R ` N ) ) , ( 2nd ` ( R ` N ) ) >. ) )
16		df-ov 6091	. . . 4 |- ( ( 1st ` ( R ` N ) ) gcd ( 2nd ` ( R ` N ) ) ) = ( gcd ` <. ( 1st ` ( R ` N ) ) , ( 2nd ` ( R ` N ) ) >. )
17	15, 16	eqtr4di 2289	. . 3 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( gcd ` ( R ` N ) ) = ( ( 1st ` ( R ` N ) ) gcd ( 2nd ` ( R ` N ) ) ) )
18	4	fveq2i 5701	. . . . . . . 8 |- ( 2nd ` A ) = ( 2nd ` <. M , N >. )
19		op2ndg 6389	. . . . . . . 8 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( 2nd ` <. M , N >. ) = N )
20	18, 19	eqtrid 2283	. . . . . . 7 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( 2nd ` A ) = N )
21	20	fveq2d 5702	. . . . . 6 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( R ` ( 2nd ` A ) ) = ( R ` N ) )
22	21	fveq2d 5702	. . . . 5 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( 2nd ` ( R ` ( 2nd ` A ) ) ) = ( 2nd ` ( R ` N ) ) )
23		xp2nd 6404	. . . . . . . . 9 |- ( A e. ( NN0 X. NN0 ) -> ( 2nd ` A ) e. NN0 )
24	23	nn0zd 9777	. . . . . . . 8 |- ( A e. ( NN0 X. NN0 ) -> ( 2nd ` A ) e. ZZ )
25		uzid 9952	. . . . . . . 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 2238	. . . . . . . 8 |- ( 2nd ` A ) = ( 2nd ` A )
28	7, 2, 27	eucalgcvga 12861	. . . . . . 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 2273	. . . 4 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( 2nd ` ( R ` N ) ) = 0 )
32	31	oveq2d 6104	. . 3 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( ( 1st ` ( R ` N ) ) gcd ( 2nd ` ( R ` N ) ) ) = ( ( 1st ` ( R ` N ) ) gcd 0 ) )
33		xp1st 6403	. . . 4 |- ( ( R ` N ) e. ( NN0 X. NN0 ) -> ( 1st ` ( R ` N ) ) e. NN0 )
34		nn0gcdid0 12783	. . . 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 2276	. 2 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( 1st ` ( R ` N ) ) = ( gcd ` ( R ` N ) ) )
37	7	eucalginv 12859	. . . . . 6 |- ( z e. ( NN0 X. NN0 ) -> ( gcd ` ( E ` z ) ) = ( gcd ` z ) )
38	8	ffvelcdmi 5845	. . . . . . 7 |- ( z e. ( NN0 X. NN0 ) -> ( E ` z ) e. ( NN0 X. NN0 ) )
39		fvres 5722	. . . . . . 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 5722	. . . . . 6 |- ( z e. ( NN0 X. NN0 ) -> ( ( gcd |` ( NN0 X. NN0 ) ) ` z ) = ( gcd ` z ) )
42	37, 40, 41	3eqtr4d 2281	. . . . 5 |- ( z e. ( NN0 X. NN0 ) -> ( ( gcd |` ( NN0 X. NN0 ) ) ` ( E ` z ) ) = ( ( gcd |` ( NN0 X. NN0 ) ) ` z ) )
43	2, 8, 42	alginv 12850	. . . 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 424	. . 3 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( ( gcd |` ( NN0 X. NN0 ) ) ` ( R ` N ) ) = ( ( gcd |` ( NN0 X. NN0 ) ) ` ( R ` 0 ) ) )
45		fvres 5722	. . . 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 9589	. . . . 5 |- 0 e. NN0
48		ffvelcdm 5844	. . . . 5 |- ( ( R : NN0 --> ( NN0 X. NN0 ) /\ 0 e. NN0 ) -> ( R ` 0 ) e. ( NN0 X. NN0 ) )
49	10, 47, 48	sylancl 417	. . . 4 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( R ` 0 ) e. ( NN0 X. NN0 ) )
50		fvres 5722	. . . 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 2279	. 2 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( gcd ` ( R ` N ) ) = ( gcd ` ( R ` 0 ) ) )
53	1, 2, 3, 6, 9	ialgr0 12847	. . . . 5 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( R ` 0 ) = A )
54	53, 4	eqtrdi 2287	. . . 4 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( R ` 0 ) = <. M , N >. )
55	54	fveq2d 5702	. . 3 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( gcd ` ( R ` 0 ) ) = ( gcd ` <. M , N >. ) )
56		df-ov 6091	. . 3 |- ( M gcd N ) = ( gcd ` <. M , N >. )
57	55, 56	eqtr4di 2289	. 2 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( gcd ` ( R ` 0 ) ) = ( M gcd N ) )
58	36, 52, 57	3eqtrd 2275	1 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( 1st ` ( R ` N ) ) = ( M gcd N ) )

This proof depends on syntax axioms:   -> wi 4   /\ wa 104   = wceq 1402   e. wcel 2209   ifcif 3638   {csn 3709   <.cop 3712   X. cxp 4774   |` cres 4778   o. ccom 4780   -->wf 5376   ` cfv 5380  (class class class)co 6088   e. cmpo 6090   1stc1st 6376   2ndc2nd 6377   0cc0 8185   NN0cn0 9574   ZZcz 9655   ZZ>=cuz 9937   mod cmo 10780   seqcseq 10905   gcd cgcd 12755

This proof depends on axioms:  ax-mp 5  ax-1 6  ax-2 7  ax-ia1 106  ax-ia2 107  ax-ia3 108  ax-in1 623  ax-in2 624  ax-io 721  ax-5 1500  ax-7 1501  ax-gen 1502  ax-ie1 1546  ax-ie2 1547  ax-8 1557  ax-10 1558  ax-11 1559  ax-i12 1560  ax-bndl 1562  ax-4 1563  ax-17 1579  ax-i9 1583  ax-ial 1587  ax-i5r 1588  ax-14 2212  ax-ext 2220  ax-coll 4248  ax-sep 4251  ax-nul 4261  ax-pow 4313  ax-pr 4348  ax-un 4580  ax-setind 4686  ax-iinf 4737  ax-cnex 8276  ax-resscn 8277  ax-1cn 8278  ax-1re 8279  ax-icn 8280  ax-addcl 8281  ax-addrcl 8282  ax-mulcl 8283  ax-mulrcl 8284  ax-addcom 8285  ax-mulcom 8286  ax-addass 8287  ax-mulass 8288  ax-distr 8289  ax-i2m1 8290  ax-0lt1 8291  ax-1rid 8292  ax-0id 8293  ax-rnegex 8294  ax-precex 8295  ax-cnre 8296  ax-pre-ltirr 8297  ax-pre-ltwlin 8298  ax-pre-lttrn 8299  ax-pre-apti 8300  ax-pre-ltadd 8301  ax-pre-mulgt0 8302  ax-pre-mulext 8303  ax-arch 8304  ax-caucvg 8305

This proof depends on definitions:  df-bi 117  df-stab 843  df-dc 847  df-3or 1010  df-3an 1011  df-tru 1405  df-fal 1408  df-nf 1514  df-sb 1816  df-eu 2089  df-mo 2090  df-clab 2225  df-cleq 2231  df-clel 2234  df-nfc 2381  df-ne 2421  df-nel 2516  df-ral 2533  df-rex 2534  df-reu 2535  df-rmo 2536  df-rab 2537  df-v 2823  df-sbc 3052  df-csb 3148  df-dif 3222  df-un 3224  df-in 3226  df-ss 3233  df-nul 3521  df-if 3639  df-pw 3690  df-sn 3715  df-pr 3716  df-op 3718  df-uni 3936  df-int 3971  df-iun 4014  df-br 4132  df-opab 4194  df-mpt 4195  df-tr 4232  df-id 4440  df-po 4443  df-iso 4444  df-iord 4513  df-on 4515  df-ilim 4516  df-suc 4518  df-iom 4740  df-xp 4782  df-rel 4783  df-cnv 4784  df-co 4785  df-dm 4786  df-rn 4787  df-res 4788  df-ima 4789  df-iota 5340  df-fun 5382  df-fn 5383  df-f 5384  df-f1 5385  df-fo 5386  df-f1o 5387  df-fv 5388  df-riota 6041  df-ov 6091  df-oprab 6092  df-mpo 6093  df-1st 6378  df-2nd 6379  df-recs 6580  df-frec 6666  df-sup 7330  df-pnf 8368  df-mnf 8369  df-xr 8370  df-ltxr 8371  df-le 8372  df-sub 8506  df-neg 8507  df-reap 8911  df-ap 8918  df-div 9011  df-inn 9314  df-2 9372  df-3 9373  df-4 9374  df-n0 9575  df-z 9656  df-uz 9938  df-q 10036  df-rp 10072  df-fz 10429  df-fzo 10567  df-fl 10722  df-mod 10781  df-seqfrec 10906  df-exp 10997  df-cj 11629  df-re 11630  df-im 11631  df-rsqrt 11786  df-abs 11787  df-dvds 12580  df-gcd 12756

This theorem is used by: (None)