Theorem eucalg 11991

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 9500	. . . . . . . 8 |- NN0 = ( ZZ>= ` 0 )
2		eucalg.2	. . . . . . . 8 |- R = seq 0 ( ( E o. 1st ) , ( NN0 X. { A } ) )
3		0zd 9203	. . . . . . . 8 |- ( ( M e. NN0 /\ N e. NN0 ) -> 0 e. ZZ )
4		eucalg.3	. . . . . . . . 9 |- A = <. M , N >.
5		opelxpi 4636	. . . . . . . . 9 |- ( ( M e. NN0 /\ N e. NN0 ) -> <. M , N >. e. ( NN0 X. NN0 ) )
6	4, 5	eqeltrid 2253	. . . . . . . 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 11987	. . . . . . . . 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 11977	. . . . . . 7 |- ( ( M e. NN0 /\ N e. NN0 ) -> R : NN0 --> ( NN0 X. NN0 ) )
11		ffvelrn 5618	. . . . . . 7 |- ( ( R : NN0 --> ( NN0 X. NN0 ) /\ N e. NN0 ) -> ( R ` N ) e. ( NN0 X. NN0 ) )
12	10, 11	sylancom 417	. . . . . 6 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( R ` N ) e. ( NN0 X. NN0 ) )
13		1st2nd2 6143	. . . . . 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 5490	. . . 4 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( gcd ` ( R ` N ) ) = ( gcd ` <. ( 1st ` ( R ` N ) ) , ( 2nd ` ( R ` N ) ) >. ) )
16		df-ov 5845	. . . 4 |- ( ( 1st ` ( R ` N ) ) gcd ( 2nd ` ( R ` N ) ) ) = ( gcd ` <. ( 1st ` ( R ` N ) ) , ( 2nd ` ( R ` N ) ) >. )
17	15, 16	eqtr4di 2217	. . 3 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( gcd ` ( R ` N ) ) = ( ( 1st ` ( R ` N ) ) gcd ( 2nd ` ( R ` N ) ) ) )
18	4	fveq2i 5489	. . . . . . . 8 |- ( 2nd ` A ) = ( 2nd ` <. M , N >. )
19		op2ndg 6119	. . . . . . . 8 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( 2nd ` <. M , N >. ) = N )
20	18, 19	syl5eq 2211	. . . . . . 7 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( 2nd ` A ) = N )
21	20	fveq2d 5490	. . . . . 6 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( R ` ( 2nd ` A ) ) = ( R ` N ) )
22	21	fveq2d 5490	. . . . 5 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( 2nd ` ( R ` ( 2nd ` A ) ) ) = ( 2nd ` ( R ` N ) ) )
23		xp2nd 6134	. . . . . . . . 9 |- ( A e. ( NN0 X. NN0 ) -> ( 2nd ` A ) e. NN0 )
24	23	nn0zd 9311	. . . . . . . 8 |- ( A e. ( NN0 X. NN0 ) -> ( 2nd ` A ) e. ZZ )
25		uzid 9480	. . . . . . . 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 2165	. . . . . . . 8 |- ( 2nd ` A ) = ( 2nd ` A )
28	7, 2, 27	eucalgcvga 11990	. . . . . . 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 2200	. . . 4 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( 2nd ` ( R ` N ) ) = 0 )
32	31	oveq2d 5858	. . 3 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( ( 1st ` ( R ` N ) ) gcd ( 2nd ` ( R ` N ) ) ) = ( ( 1st ` ( R ` N ) ) gcd 0 ) )
33		xp1st 6133	. . . 4 |- ( ( R ` N ) e. ( NN0 X. NN0 ) -> ( 1st ` ( R ` N ) ) e. NN0 )
34		nn0gcdid0 11914	. . . 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 2203	. 2 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( 1st ` ( R ` N ) ) = ( gcd ` ( R ` N ) ) )
37	7	eucalginv 11988	. . . . . 6 |- ( z e. ( NN0 X. NN0 ) -> ( gcd ` ( E ` z ) ) = ( gcd ` z ) )
38	8	ffvelrni 5619	. . . . . . 7 |- ( z e. ( NN0 X. NN0 ) -> ( E ` z ) e. ( NN0 X. NN0 ) )
39		fvres 5510	. . . . . . 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 5510	. . . . . 6 |- ( z e. ( NN0 X. NN0 ) -> ( ( gcd |` ( NN0 X. NN0 ) ) ` z ) = ( gcd ` z ) )
42	37, 40, 41	3eqtr4d 2208	. . . . 5 |- ( z e. ( NN0 X. NN0 ) -> ( ( gcd |` ( NN0 X. NN0 ) ) ` ( E ` z ) ) = ( ( gcd |` ( NN0 X. NN0 ) ) ` z ) )
43	2, 8, 42	alginv 11979	. . . 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 417	. . 3 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( ( gcd |` ( NN0 X. NN0 ) ) ` ( R ` N ) ) = ( ( gcd |` ( NN0 X. NN0 ) ) ` ( R ` 0 ) ) )
45		fvres 5510	. . . 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 9129	. . . . 5 |- 0 e. NN0
48		ffvelrn 5618	. . . . 5 |- ( ( R : NN0 --> ( NN0 X. NN0 ) /\ 0 e. NN0 ) -> ( R ` 0 ) e. ( NN0 X. NN0 ) )
49	10, 47, 48	sylancl 410	. . . 4 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( R ` 0 ) e. ( NN0 X. NN0 ) )
50		fvres 5510	. . . 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 2206	. 2 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( gcd ` ( R ` N ) ) = ( gcd ` ( R ` 0 ) ) )
53	1, 2, 3, 6, 9	ialgr0 11976	. . . . 5 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( R ` 0 ) = A )
54	53, 4	eqtrdi 2215	. . . 4 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( R ` 0 ) = <. M , N >. )
55	54	fveq2d 5490	. . 3 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( gcd ` ( R ` 0 ) ) = ( gcd ` <. M , N >. ) )
56		df-ov 5845	. . 3 |- ( M gcd N ) = ( gcd ` <. M , N >. )
57	55, 56	eqtr4di 2217	. 2 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( gcd ` ( R ` 0 ) ) = ( M gcd N ) )
58	36, 52, 57	3eqtrd 2202	1 |- ( ( M e. NN0 /\ N e. NN0 ) -> ( 1st ` ( R ` N ) ) = ( M gcd N ) )

Syntax hints:   -> wi 4   /\ wa 103   = wceq 1343   e. wcel 2136   ifcif 3520   {csn 3576   <.cop 3579   X. cxp 4602   |` cres 4606   o. ccom 4608   -->wf 5184   ` cfv 5188  (class class class)co 5842   e. cmpo 5844   1stc1st 6106   2ndc2nd 6107   0cc0 7753   NN0cn0 9114   ZZcz 9191   ZZ>=cuz 9466   mod cmo 10257   seqcseq 10380   gcd cgcd 11875

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 604  ax-in2 605  ax-io 699  ax-5 1435  ax-7 1436  ax-gen 1437  ax-ie1 1481  ax-ie2 1482  ax-8 1492  ax-10 1493  ax-11 1494  ax-i12 1495  ax-bndl 1497  ax-4 1498  ax-17 1514  ax-i9 1518  ax-ial 1522  ax-i5r 1523  ax-13 2138  ax-14 2139  ax-ext 2147  ax-coll 4097  ax-sep 4100  ax-nul 4108  ax-pow 4153  ax-pr 4187  ax-un 4411  ax-setind 4514  ax-iinf 4565  ax-cnex 7844  ax-resscn 7845  ax-1cn 7846  ax-1re 7847  ax-icn 7848  ax-addcl 7849  ax-addrcl 7850  ax-mulcl 7851  ax-mulrcl 7852  ax-addcom 7853  ax-mulcom 7854  ax-addass 7855  ax-mulass 7856  ax-distr 7857  ax-i2m1 7858  ax-0lt1 7859  ax-1rid 7860  ax-0id 7861  ax-rnegex 7862  ax-precex 7863  ax-cnre 7864  ax-pre-ltirr 7865  ax-pre-ltwlin 7866  ax-pre-lttrn 7867  ax-pre-apti 7868  ax-pre-ltadd 7869  ax-pre-mulgt0 7870  ax-pre-mulext 7871  ax-arch 7872  ax-caucvg 7873

This theorem depends on definitions:  df-bi 116  df-stab 821  df-dc 825  df-3or 969  df-3an 970  df-tru 1346  df-fal 1349  df-nf 1449  df-sb 1751  df-eu 2017  df-mo 2018  df-clab 2152  df-cleq 2158  df-clel 2161  df-nfc 2297  df-ne 2337  df-nel 2432  df-ral 2449  df-rex 2450  df-reu 2451  df-rmo 2452  df-rab 2453  df-v 2728  df-sbc 2952  df-csb 3046  df-dif 3118  df-un 3120  df-in 3122  df-ss 3129  df-nul 3410  df-if 3521  df-pw 3561  df-sn 3582  df-pr 3583  df-op 3585  df-uni 3790  df-int 3825  df-iun 3868  df-br 3983  df-opab 4044  df-mpt 4045  df-tr 4081  df-id 4271  df-po 4274  df-iso 4275  df-iord 4344  df-on 4346  df-ilim 4347  df-suc 4349  df-iom 4568  df-xp 4610  df-rel 4611  df-cnv 4612  df-co 4613  df-dm 4614  df-rn 4615  df-res 4616  df-ima 4617  df-iota 5153  df-fun 5190  df-fn 5191  df-f 5192  df-f1 5193  df-fo 5194  df-f1o 5195  df-fv 5196  df-riota 5798  df-ov 5845  df-oprab 5846  df-mpo 5847  df-1st 6108  df-2nd 6109  df-recs 6273  df-frec 6359  df-sup 6949  df-pnf 7935  df-mnf 7936  df-xr 7937  df-ltxr 7938  df-le 7939  df-sub 8071  df-neg 8072  df-reap 8473  df-ap 8480  df-div 8569  df-inn 8858  df-2 8916  df-3 8917  df-4 8918  df-n0 9115  df-z 9192  df-uz 9467  df-q 9558  df-rp 9590  df-fz 9945  df-fzo 10078  df-fl 10205  df-mod 10258  df-seqfrec 10381  df-exp 10455  df-cj 10784  df-re 10785  df-im 10786  df-rsqrt 10940  df-abs 10941  df-dvds 11728  df-gcd 11876

This theorem is referenced by: (None)