Theorem mul4sq 12832

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 12831. (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 12830	. 2 |- ( A e. S <-> E. a e. ZZ[_i] E. b e. ZZ[_i] A = ( ( ( abs ` a ) ^ 2 ) + ( ( abs ` b ) ^ 2 ) ) )
3	1	4sqlem4 12830	. 2 |- ( B e. S <-> E. c e. ZZ[_i] E. d e. ZZ[_i] B = ( ( ( abs ` c ) ^ 2 ) + ( ( abs ` d ) ^ 2 ) ) )
4		reeanv 2678	. . 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 2678	. . . . 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 12819	. . . . . . . . . . . . 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 529	. . . . . . . . . . . . 13 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> b e. ZZ[_i] )
10		gzabssqcl 12819	. . . . . . . . . . . . 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 9387	. . . . . . . . . . 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 9385	. . . . . . . . . 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 8888	. . . . . . . . 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 528	. . . . . . . . . . . . 13 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> c e. ZZ[_i] )
16		gzabssqcl 12819	. . . . . . . . . . . . 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 531	. . . . . . . . . . . . 13 |- ( ( ( a e. ZZ[_i] /\ c e. ZZ[_i] ) /\ ( b e. ZZ[_i] /\ d e. ZZ[_i] ) ) -> d e. ZZ[_i] )
19		gzabssqcl 12819	. . . . . . . . . . . . 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 9387	. . . . . . . . . . 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 9385	. . . . . . . . . 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 8888	. . . . . . . . 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 5985	. . . . . . . 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 2207	. . . . . . . . 9 |- ( ( ( abs ` a ) ^ 2 ) + ( ( abs ` b ) ^ 2 ) ) = ( ( ( abs ` a ) ^ 2 ) + ( ( abs ` b ) ^ 2 ) )
26		eqid 2207	. . . . . . . . 9 |- ( ( ( abs ` c ) ^ 2 ) + ( ( abs ` d ) ^ 2 ) ) = ( ( ( abs ` c ) ^ 2 ) + ( ( abs ` d ) ^ 2 ) )
27		1nn 9082	. . . . . . . . . 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 12818	. . . . . . . . . . . . 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 12810	. . . . . . . . . . . 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 8888	. . . . . . . . . 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 2284	. . . . . . . . 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 12818	. . . . . . . . . . . . 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 12810	. . . . . . . . . . . 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 8888	. . . . . . . . . 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 2284	. . . . . . . . 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 2284	. . . . . . . . 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 12831	. . . . . . . 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 2285	. . . . . . 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 5976	. . . . . . . 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 2276	. . . . . . 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 2633	. . . . 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 2631	. . 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 1373   e. wcel 2178   {cab 2193   E.wrex 2487   ` cfv 5290  (class class class)co 5967   CCcc 7958   1c1 7961   + caddc 7963   x. cmul 7965   - cmin 8278   / cdiv 8780   NNcn 9071   2c2 9122   NN0cn0 9330   ZZcz 9407   ^cexp 10720   abscabs 11423   ZZ[_i]cgz 12807

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 615  ax-in2 616  ax-io 711  ax-5 1471  ax-7 1472  ax-gen 1473  ax-ie1 1517  ax-ie2 1518  ax-8 1528  ax-10 1529  ax-11 1530  ax-i12 1531  ax-bndl 1533  ax-4 1534  ax-17 1550  ax-i9 1554  ax-ial 1558  ax-i5r 1559  ax-13 2180  ax-14 2181  ax-ext 2189  ax-coll 4175  ax-sep 4178  ax-nul 4186  ax-pow 4234  ax-pr 4269  ax-un 4498  ax-setind 4603  ax-iinf 4654  ax-cnex 8051  ax-resscn 8052  ax-1cn 8053  ax-1re 8054  ax-icn 8055  ax-addcl 8056  ax-addrcl 8057  ax-mulcl 8058  ax-mulrcl 8059  ax-addcom 8060  ax-mulcom 8061  ax-addass 8062  ax-mulass 8063  ax-distr 8064  ax-i2m1 8065  ax-0lt1 8066  ax-1rid 8067  ax-0id 8068  ax-rnegex 8069  ax-precex 8070  ax-cnre 8071  ax-pre-ltirr 8072  ax-pre-ltwlin 8073  ax-pre-lttrn 8074  ax-pre-apti 8075  ax-pre-ltadd 8076  ax-pre-mulgt0 8077  ax-pre-mulext 8078  ax-arch 8079  ax-caucvg 8080

This theorem depends on definitions:  df-bi 117  df-dc 837  df-3or 982  df-3an 983  df-tru 1376  df-fal 1379  df-nf 1485  df-sb 1787  df-eu 2058  df-mo 2059  df-clab 2194  df-cleq 2200  df-clel 2203  df-nfc 2339  df-ne 2379  df-nel 2474  df-ral 2491  df-rex 2492  df-reu 2493  df-rmo 2494  df-rab 2495  df-v 2778  df-sbc 3006  df-csb 3102  df-dif 3176  df-un 3178  df-in 3180  df-ss 3187  df-nul 3469  df-if 3580  df-pw 3628  df-sn 3649  df-pr 3650  df-op 3652  df-uni 3865  df-int 3900  df-iun 3943  df-br 4060  df-opab 4122  df-mpt 4123  df-tr 4159  df-id 4358  df-po 4361  df-iso 4362  df-iord 4431  df-on 4433  df-ilim 4434  df-suc 4436  df-iom 4657  df-xp 4699  df-rel 4700  df-cnv 4701  df-co 4702  df-dm 4703  df-rn 4704  df-res 4705  df-ima 4706  df-iota 5251  df-fun 5292  df-fn 5293  df-f 5294  df-f1 5295  df-fo 5296  df-f1o 5297  df-fv 5298  df-riota 5922  df-ov 5970  df-oprab 5971  df-mpo 5972  df-1st 6249  df-2nd 6250  df-recs 6414  df-frec 6500  df-pnf 8144  df-mnf 8145  df-xr 8146  df-ltxr 8147  df-le 8148  df-sub 8280  df-neg 8281  df-reap 8683  df-ap 8690  df-div 8781  df-inn 9072  df-2 9130  df-3 9131  df-4 9132  df-n0 9331  df-z 9408  df-uz 9684  df-rp 9811  df-seqfrec 10630  df-exp 10721  df-cj 11268  df-re 11269  df-im 11270  df-rsqrt 11424  df-abs 11425  df-gz 12808

This theorem is referenced by:  4sqlem19  12847