Theorem mul4sq 13117

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 13116. (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 13115	. 2 |- ( A e. S <-> E. a e. ZZ[_i] E. b e. ZZ[_i] A = ( ( ( abs ` a ) ^ 2 ) + ( ( abs ` b ) ^ 2 ) ) )
3	1	4sqlem4 13115	. 2 |- ( B e. S <-> E. c e. ZZ[_i] E. d e. ZZ[_i] B = ( ( ( abs ` c ) ^ 2 ) + ( ( abs ` d ) ^ 2 ) ) )
4		reeanv 2715	. . 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 2715	. . . . 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 527	. . . . . . . . . . . . 13 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> a e. ZZ[_i] )
7		gzabssqcl 13104	. . . . . . . . . . . . 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 531	. . . . . . . . . . . . 13 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> b e. ZZ[_i] )
10		gzabssqcl 13104	. . . . . . . . . . . . 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 9574	. . . . . . . . . . 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 9572	. . . . . . . . . 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 9071	. . . . . . . . 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 529	. . . . . . . . . . . . 13 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> c e. ZZ[_i] )
16		gzabssqcl 13104	. . . . . . . . . . . . 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 533	. . . . . . . . . . . . 13 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> d e. ZZ[_i] )
19		gzabssqcl 13104	. . . . . . . . . . . . 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 9574	. . . . . . . . . . 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 9572	. . . . . . . . . 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 9071	. . . . . . . . 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 6076	. . . . . . . 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 2234	. . . . . . . . 9 |- ( ( ( abs ` a ) ^ 2 ) + ( ( abs ` b ) ^ 2 ) ) = ( ( ( abs ` a ) ^ 2 ) + ( ( abs ` b ) ^ 2 ) )
26		eqid 2234	. . . . . . . . 9 |- ( ( ( abs ` c ) ^ 2 ) + ( ( abs ` d ) ^ 2 ) ) = ( ( ( abs ` c ) ^ 2 ) + ( ( abs ` d ) ^ 2 ) )
27		1nn 9265	. . . . . . . . . 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 13103	. . . . . . . . . . . . 13 |- ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) -> ( a - c ) e. ZZ[_i] )
30	29	adantr 276	. . . . . . . . . . . 12 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> ( a - c ) e. ZZ[_i] )
31		gzcn 13095	. . . . . . . . . . . 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 9071	. . . . . . . . . 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 2311	. . . . . . . . 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 13103	. . . . . . . . . . . . 13 |- ( ( b e. ZZ[_i] /\ d e. ZZ[_i] ) -> ( b - d ) e. ZZ[_i] )
36	35	adantl 277	. . . . . . . . . . . 12 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> ( b - d ) e. ZZ[_i] )
37		gzcn 13095	. . . . . . . . . . . 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 9071	. . . . . . . . . 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 2311	. . . . . . . . 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 2311	. . . . . . . . 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 13116	. . . . . . . 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 2312	. . . . . . 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 6067	. . . . . . . 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 2303	. . . . . . 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 157	. . . . . 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 2670	. . . . 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	biimtrrid 153	. . . 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 2668	. . 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 135	. 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 291	1 |- ( ( A e. S /\ B e. S ) -> ( A x. B ) e. S )

Syntax hints:   -> wi 4   /\ wa 104   = wceq 1398   e. wcel 2205   {cab 2220   E.wrex 2523   ` cfv 5357  (class class class)co 6058   CCcc 8141   1c1 8144   + caddc 8146   x. cmul 8148   - cmin 8460   / cdiv 8963   NNcn 9254   2c2 9305   NN0cn0 9513   ZZcz 9594   ^cexp 10924   abscabs 11707   ZZ[_i]cgz 13092

This theorem was proved from axioms:  ax-mp 5  ax-1 6  ax-2 7  ax-ia1 106  ax-ia2 107  ax-ia3 108  ax-in1 619  ax-in2 620  ax-io 717  ax-5 1496  ax-7 1497  ax-gen 1498  ax-ie1 1542  ax-ie2 1543  ax-8 1553  ax-10 1554  ax-11 1555  ax-i12 1556  ax-bndl 1558  ax-4 1559  ax-17 1575  ax-i9 1579  ax-ial 1583  ax-i5r 1584  ax-13 2207  ax-14 2208  ax-ext 2216  ax-coll 4230  ax-sep 4233  ax-nul 4241  ax-pow 4292  ax-pr 4327  ax-un 4559  ax-setind 4664  ax-iinf 4715  ax-cnex 8234  ax-resscn 8235  ax-1cn 8236  ax-1re 8237  ax-icn 8238  ax-addcl 8239  ax-addrcl 8240  ax-mulcl 8241  ax-mulrcl 8242  ax-addcom 8243  ax-mulcom 8244  ax-addass 8245  ax-mulass 8246  ax-distr 8247  ax-i2m1 8248  ax-0lt1 8249  ax-1rid 8250  ax-0id 8251  ax-rnegex 8252  ax-precex 8253  ax-cnre 8254  ax-pre-ltirr 8255  ax-pre-ltwlin 8256  ax-pre-lttrn 8257  ax-pre-apti 8258  ax-pre-ltadd 8259  ax-pre-mulgt0 8260  ax-pre-mulext 8261  ax-arch 8262  ax-caucvg 8263

This theorem depends on definitions:  df-bi 117  df-dc 843  df-3or 1006  df-3an 1007  df-tru 1401  df-fal 1404  df-nf 1510  df-sb 1812  df-eu 2085  df-mo 2086  df-clab 2221  df-cleq 2227  df-clel 2230  df-nfc 2375  df-ne 2415  df-nel 2510  df-ral 2527  df-rex 2528  df-reu 2529  df-rmo 2530  df-rab 2531  df-v 2817  df-sbc 3046  df-csb 3142  df-dif 3216  df-un 3218  df-in 3220  df-ss 3227  df-nul 3513  df-if 3625  df-pw 3676  df-sn 3700  df-pr 3701  df-op 3703  df-uni 3920  df-int 3955  df-iun 3998  df-br 4115  df-opab 4177  df-mpt 4178  df-tr 4214  df-id 4419  df-po 4422  df-iso 4423  df-iord 4492  df-on 4494  df-ilim 4495  df-suc 4497  df-iom 4718  df-xp 4760  df-rel 4761  df-cnv 4762  df-co 4763  df-dm 4764  df-rn 4765  df-res 4766  df-ima 4767  df-iota 5317  df-fun 5359  df-fn 5360  df-f 5361  df-f1 5362  df-fo 5363  df-f1o 5364  df-fv 5365  df-riota 6011  df-ov 6061  df-oprab 6062  df-mpo 6063  df-1st 6347  df-2nd 6348  df-recs 6549  df-frec 6635  df-pnf 8326  df-mnf 8327  df-xr 8328  df-ltxr 8329  df-le 8330  df-sub 8462  df-neg 8463  df-reap 8866  df-ap 8873  df-div 8964  df-inn 9255  df-2 9313  df-3 9314  df-4 9315  df-n0 9514  df-z 9595  df-uz 9872  df-rp 10005  df-seqfrec 10834  df-exp 10925  df-cj 11552  df-re 11553  df-im 11554  df-rsqrt 11708  df-abs 11709  df-gz 13093

This theorem is referenced by:  4sqlem19  13132