Theorem mul4sq 12346

Description: Euler's four-square identity: The product of two sums of four squares is also a sum of four squares. This is usually quoted as an explicit formula involving eight real variables; we save some time by working with complex numbers (gaussian integers) instead, so that we only have to work with four variables, and also hiding the actual formula for the product in the proof of mul4sqlem 12345. (For the curious, the explicit formula that is used is ( | a | ^ 2 + | b | ^ 2 ) ( | c | ^ 2 + | d | ^ 2 ) = | a * x. c + b x. d * | ^ 2 + | a * x. d - b x. c * | ^ 2.) (Contributed by Mario Carneiro, 14-Jul-2014.)

Hypothesis

Ref	Expression
4sq.1	|- S = { n | E. x e. ZZ E. y e. ZZ E. z e. ZZ E. w e. ZZ n = ( ( ( x ^ 2 ) + ( y ^ 2 ) ) + ( ( z ^ 2 ) + ( w ^ 2 ) ) ) }

Assertion

Ref	Expression
mul4sq	|- ( ( A e. S /\ B e. S ) -> ( A x. B ) e. S )

Distinct variable groups:   w, n, x, y, z   B, n   A, n   S, n

Allowed substitution hints:   A( x, y, z, w)   B( x, y, z, w)   S( x, y, z, w)

Proof of Theorem mul4sq

Dummy variables a b c d are mutually distinct and distinct from all other variables.

Step	Hyp	Ref	Expression
1		4sq.1	. . 3 |- S = { n | E. x e. ZZ E. y e. ZZ E. z e. ZZ E. w e. ZZ n = ( ( ( x ^ 2 ) + ( y ^ 2 ) ) + ( ( z ^ 2 ) + ( w ^ 2 ) ) ) }
2	1	4sqlem4 12344	. 2 |- ( A e. S <-> E. a e. ZZ[_i] E. b e. ZZ[_i] A = ( ( ( abs ` a ) ^ 2 ) + ( ( abs ` b ) ^ 2 ) ) )
3	1	4sqlem4 12344	. 2 |- ( B e. S <-> E. c e. ZZ[_i] E. d e. ZZ[_i] B = ( ( ( abs ` c ) ^ 2 ) + ( ( abs ` d ) ^ 2 ) ) )
4		reeanv 2639	. . 3 |- ( E. a e. ZZ[_i] E. c e. ZZ[_i] ( E. b e. ZZ[_i] A = ( ( ( abs ` a ) ^ 2 ) + ( ( abs ` b ) ^ 2 ) ) /\ E. d e. ZZ[_i] B = ( ( ( abs ` c ) ^ 2 ) + ( ( abs ` d ) ^ 2 ) ) ) <-> ( E. a e. ZZ[_i] E. b e. ZZ[_i] A = ( ( ( abs ` a ) ^ 2 ) + ( ( abs ` b ) ^ 2 ) ) /\ E. c e. ZZ[_i] E. d e. ZZ[_i] B = ( ( ( abs ` c ) ^ 2 ) + ( ( abs ` d ) ^ 2 ) ) ) )
5		reeanv 2639	. . . . 5 |- ( E. b e. ZZ[_i] E. d e. ZZ[_i] ( A = ( ( ( abs ` a ) ^ 2 ) + ( ( abs ` b ) ^ 2 ) ) /\ B = ( ( ( abs ` c ) ^ 2 ) + ( ( abs ` d ) ^ 2 ) ) ) <-> ( E. b e. ZZ[_i] A = ( ( ( abs ` a ) ^ 2 ) + ( ( abs ` b ) ^ 2 ) ) /\ E. d e. ZZ[_i] B = ( ( ( abs ` c ) ^ 2 ) + ( ( abs ` d ) ^ 2 ) ) ) )
6		simpll 524	. . . . . . . . . . . . 13 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> a e. ZZ[_i] )
7		gzabssqcl 12333	. . . . . . . . . . . . 13 |- ( a e. ZZ[_i] -> ( ( abs ` a ) ^ 2 ) e. NN0 )
8	6, 7	syl 14	. . . . . . . . . . . 12 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> ( ( abs ` a ) ^ 2 ) e. NN0 )
9		simprl 526	. . . . . . . . . . . . 13 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> b e. ZZ[_i] )
10		gzabssqcl 12333	. . . . . . . . . . . . 13 |- ( b e. ZZ[_i] -> ( ( abs ` b ) ^ 2 ) e. NN0 )
11	9, 10	syl 14	. . . . . . . . . . . 12 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> ( ( abs ` b ) ^ 2 ) e. NN0 )
12	8, 11	nn0addcld 9192	. . . . . . . . . . 11 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> ( ( ( abs ` a ) ^ 2 ) + ( ( abs ` b ) ^ 2 ) ) e. NN0 )
13	12	nn0cnd 9190	. . . . . . . . . 10 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> ( ( ( abs ` a ) ^ 2 ) + ( ( abs ` b ) ^ 2 ) ) e. CC )
14	13	div1d 8697	. . . . . . . . 9 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> ( ( ( ( abs ` a ) ^ 2 ) + ( ( abs ` b ) ^ 2 ) ) / 1 ) = ( ( ( abs ` a ) ^ 2 ) + ( ( abs ` b ) ^ 2 ) ) )
15		simplr 525	. . . . . . . . . . . . 13 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> c e. ZZ[_i] )
16		gzabssqcl 12333	. . . . . . . . . . . . 13 |- ( c e. ZZ[_i] -> ( ( abs ` c ) ^ 2 ) e. NN0 )
17	15, 16	syl 14	. . . . . . . . . . . 12 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> ( ( abs ` c ) ^ 2 ) e. NN0 )
18		simprr 527	. . . . . . . . . . . . 13 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> d e. ZZ[_i] )
19		gzabssqcl 12333	. . . . . . . . . . . . 13 |- ( d e. ZZ[_i] -> ( ( abs ` d ) ^ 2 ) e. NN0 )
20	18, 19	syl 14	. . . . . . . . . . . 12 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> ( ( abs ` d ) ^ 2 ) e. NN0 )
21	17, 20	nn0addcld 9192	. . . . . . . . . . 11 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> ( ( ( abs ` c ) ^ 2 ) + ( ( abs ` d ) ^ 2 ) ) e. NN0 )
22	21	nn0cnd 9190	. . . . . . . . . 10 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> ( ( ( abs ` c ) ^ 2 ) + ( ( abs ` d ) ^ 2 ) ) e. CC )
23	22	div1d 8697	. . . . . . . . 9 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> ( ( ( ( abs ` c ) ^ 2 ) + ( ( abs ` d ) ^ 2 ) ) / 1 ) = ( ( ( abs ` c ) ^ 2 ) + ( ( abs ` d ) ^ 2 ) ) )
24	14, 23	oveq12d 5871	. . . . . . . 8 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> ( ( ( ( ( abs ` a ) ^ 2 ) + ( ( abs ` b ) ^ 2 ) ) / 1 ) x. ( ( ( ( abs ` c ) ^ 2 ) + ( ( abs ` d ) ^ 2 ) ) / 1 ) ) = ( ( ( ( abs ` a ) ^ 2 ) + ( ( abs ` b ) ^ 2 ) ) x. ( ( ( abs ` c ) ^ 2 ) + ( ( abs ` d ) ^ 2 ) ) ) )
25		eqid 2170	. . . . . . . . 9 |- ( ( ( abs ` a ) ^ 2 ) + ( ( abs ` b ) ^ 2 ) ) = ( ( ( abs ` a ) ^ 2 ) + ( ( abs ` b ) ^ 2 ) )
26		eqid 2170	. . . . . . . . 9 |- ( ( ( abs ` c ) ^ 2 ) + ( ( abs ` d ) ^ 2 ) ) = ( ( ( abs ` c ) ^ 2 ) + ( ( abs ` d ) ^ 2 ) )
27		1nn 8889	. . . . . . . . . 10 |- 1 e. NN
28	27	a1i 9	. . . . . . . . 9 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> 1 e. NN )
29		gzsubcl 12332	. . . . . . . . . . . . 13 |- ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) -> ( a - c ) e. ZZ[_i] )
30	29	adantr 274	. . . . . . . . . . . 12 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> ( a - c ) e. ZZ[_i] )
31		gzcn 12324	. . . . . . . . . . . 12 |- ( ( a - c ) e. ZZ[_i] -> ( a - c ) e. CC )
32	30, 31	syl 14	. . . . . . . . . . 11 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> ( a - c ) e. CC )
33	32	div1d 8697	. . . . . . . . . 10 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> ( ( a - c ) / 1 ) = ( a - c ) )
34	33, 30	eqeltrd 2247	. . . . . . . . 9 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> ( ( a - c ) / 1 ) e. ZZ[_i] )
35		gzsubcl 12332	. . . . . . . . . . . . 13 |- ( ( b e. ZZ[_i] /\ d e. ZZ[_i] ) -> ( b - d ) e. ZZ[_i] )
36	35	adantl 275	. . . . . . . . . . . 12 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> ( b - d ) e. ZZ[_i] )
37		gzcn 12324	. . . . . . . . . . . 12 |- ( ( b - d ) e. ZZ[_i] -> ( b - d ) e. CC )
38	36, 37	syl 14	. . . . . . . . . . 11 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> ( b - d ) e. CC )
39	38	div1d 8697	. . . . . . . . . 10 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> ( ( b - d ) / 1 ) = ( b - d ) )
40	39, 36	eqeltrd 2247	. . . . . . . . 9 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> ( ( b - d ) / 1 ) e. ZZ[_i] )
41	14, 12	eqeltrd 2247	. . . . . . . . 9 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> ( ( ( ( abs ` a ) ^ 2 ) + ( ( abs ` b ) ^ 2 ) ) / 1 ) e. NN0 )
42	1, 6, 9, 15, 18, 25, 26, 28, 34, 40, 41	mul4sqlem 12345	. . . . . . . 8 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> ( ( ( ( ( abs ` a ) ^ 2 ) + ( ( abs ` b ) ^ 2 ) ) / 1 ) x. ( ( ( ( abs ` c ) ^ 2 ) + ( ( abs ` d ) ^ 2 ) ) / 1 ) ) e. S )
43	24, 42	eqeltrrd 2248	. . . . . . 7 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> ( ( ( ( abs ` a ) ^ 2 ) + ( ( abs ` b ) ^ 2 ) ) x. ( ( ( abs ` c ) ^ 2 ) + ( ( abs ` d ) ^ 2 ) ) ) e. S )
44		oveq12 5862	. . . . . . . 8 |- ( ( A = ( ( ( abs ` a ) ^ 2 ) + ( ( abs ` b ) ^ 2 ) ) /\ B = ( ( ( abs ` c ) ^ 2 ) + ( ( abs ` d ) ^ 2 ) ) ) -> ( A x. B ) = ( ( ( ( abs ` a ) ^ 2 ) + ( ( abs ` b ) ^ 2 ) ) x. ( ( ( abs ` c ) ^ 2 ) + ( ( abs ` d ) ^ 2 ) ) ) )
45	44	eleq1d 2239	. . . . . . 7 |- ( ( A = ( ( ( abs ` a ) ^ 2 ) + ( ( abs ` b ) ^ 2 ) ) /\ B = ( ( ( abs ` c ) ^ 2 ) + ( ( abs ` d ) ^ 2 ) ) ) -> ( ( A x. B ) e. S <-> ( ( ( ( abs ` a ) ^ 2 ) + ( ( abs ` b ) ^ 2 ) ) x. ( ( ( abs ` c ) ^ 2 ) + ( ( abs ` d ) ^ 2 ) ) ) e. S ) )
46	43, 45	syl5ibrcom 156	. . . . . 6 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> ( ( A = ( ( ( abs ` a ) ^ 2 ) + ( ( abs ` b ) ^ 2 ) ) /\ B = ( ( ( abs ` c ) ^ 2 ) + ( ( abs ` d ) ^ 2 ) ) ) -> ( A x. B ) e. S ) )
47	46	rexlimdvva 2595	. . . . 5 |- ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) -> ( E. b e. ZZ[_i] E. d e. ZZ[_i] ( A = ( ( ( abs ` a ) ^ 2 ) + ( ( abs ` b ) ^ 2 ) ) /\ B = ( ( ( abs ` c ) ^ 2 ) + ( ( abs ` d ) ^ 2 ) ) ) -> ( A x. B ) e. S ) )
48	5, 47	syl5bir 152	. . . 4 |- ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) -> ( ( E. b e. ZZ[_i] A = ( ( ( abs ` a ) ^ 2 ) + ( ( abs ` b ) ^ 2 ) ) /\ E. d e. ZZ[_i] B = ( ( ( abs ` c ) ^ 2 ) + ( ( abs ` d ) ^ 2 ) ) ) -> ( A x. B ) e. S ) )
49	48	rexlimivv 2593	. . 3 |- ( E. a e. ZZ[_i] E. c e. ZZ[_i] ( E. b e. ZZ[_i] A = ( ( ( abs ` a ) ^ 2 ) + ( ( abs ` b ) ^ 2 ) ) /\ E. d e. ZZ[_i] B = ( ( ( abs ` c ) ^ 2 ) + ( ( abs ` d ) ^ 2 ) ) ) -> ( A x. B ) e. S )
50	4, 49	sylbir 134	. 2 |- ( ( E. a e. ZZ[_i] E. b e. ZZ[_i] A = ( ( ( abs ` a ) ^ 2 ) + ( ( abs ` b ) ^ 2 ) ) /\ E. c e. ZZ[_i] E. d e. ZZ[_i] B = ( ( ( abs ` c ) ^ 2 ) + ( ( abs ` d ) ^ 2 ) ) ) -> ( A x. B ) e. S )
51	2, 3, 50	syl2anb 289	1 |- ( ( A e. S /\ B e. S ) -> ( A x. B ) e. S )

Syntax hints:   -> wi 4   /\ wa 103   = wceq 1348   e. wcel 2141   {cab 2156   E.wrex 2449   ` cfv 5198  (class class class)co 5853   CCcc 7772   1c1 7775   + caddc 7777   x. cmul 7779   - cmin 8090   / cdiv 8589   NNcn 8878   2c2 8929   NN0cn0 9135   ZZcz 9212   ^cexp 10475   abscabs 10961   ZZ[_i]cgz 12321

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-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-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-rp 9611  df-seqfrec 10402  df-exp 10476  df-cj 10806  df-re 10807  df-im 10808  df-rsqrt 10962  df-abs 10963  df-gz 12322

This theorem is referenced by: (None)