Theorem karatsuba 13158

Description: The Karatsuba multiplication algorithm. If X and Y are decomposed into two groups of digits of length M (only the lower group is known to be this size but the algorithm is most efficient when the partition is chosen near the middle of the digit string), then X Y can be written in three groups of digits, where each group needs only one multiplication. Thus, we can halve both inputs with only three multiplications on the smaller operands, yielding an asymptotic improvement of n^(log₂ 3) instead of n^2 for the "naive" algorithm decmul1c 9795. (Contributed by Mario Carneiro, 16-Jul-2015.) (Revised by AV, 9-Sep-2021.)

Hypotheses

Ref	Expression
karatsuba.a	|- A e. NN0
karatsuba.b	|- B e. NN0
karatsuba.c	|- C e. NN0
karatsuba.d	|- D e. NN0
karatsuba.s	|- S e. NN0
karatsuba.m	|- M e. NN0
karatsuba.r	|- ( A x. C ) = R
karatsuba.t	|- ( B x. D ) = T
karatsuba.e	|- ( ( A + B ) x. ( C + D ) ) = ( ( R + S ) + T )
karatsuba.x	|- ( ( A x. (; 1 0 ^ M ) ) + B ) = X
karatsuba.y	|- ( ( C x. (; 1 0 ^ M ) ) + D ) = Y
karatsuba.w	|- ( ( R x. (; 1 0 ^ M ) ) + S ) = W
karatsuba.z	|- ( ( W x. (; 1 0 ^ M ) ) + T ) = Z

Assertion

Ref	Expression
karatsuba	|- ( X x. Y ) = Z

Proof of Theorem karatsuba

Step	Hyp	Ref	Expression
1		karatsuba.a	. . . . . 6 |- A e. NN0
2	1	nn0cni 9529	. . . . 5 |- A e. CC
3		10nn0 9748	. . . . . . 7 |- ; 1 0 e. NN0
4	3	nn0cni 9529	. . . . . 6 |- ; 1 0 e. CC
5		karatsuba.m	. . . . . 6 |- M e. NN0
6		expcl 10947	. . . . . 6 |- ( (; 1 0 e. CC /\ M e. NN0 ) -> (; 1 0 ^ M ) e. CC )
7	4, 5, 6	mp2an 426	. . . . 5 |- (; 1 0 ^ M ) e. CC
8	2, 7	mulcli 8296	. . . 4 |- ( A x. (; 1 0 ^ M ) ) e. CC
9		karatsuba.b	. . . . 5 |- B e. NN0
10	9	nn0cni 9529	. . . 4 |- B e. CC
11		karatsuba.c	. . . . . 6 |- C e. NN0
12	11	nn0cni 9529	. . . . 5 |- C e. CC
13	12, 7	mulcli 8296	. . . 4 |- ( C x. (; 1 0 ^ M ) ) e. CC
14		karatsuba.d	. . . . 5 |- D e. NN0
15	14	nn0cni 9529	. . . 4 |- D e. CC
16	8, 10, 13, 15	muladdi 8701	. . 3 |- ( ( ( A x. (; 1 0 ^ M ) ) + B ) x. ( ( C x. (; 1 0 ^ M ) ) + D ) ) = ( ( ( ( A x. (; 1 0 ^ M ) ) x. ( C x. (; 1 0 ^ M ) ) ) + ( D x. B ) ) + ( ( ( A x. (; 1 0 ^ M ) ) x. D ) + ( ( C x. (; 1 0 ^ M ) ) x. B ) ) )
17	8, 13	mulcli 8296	. . . 4 |- ( ( A x. (; 1 0 ^ M ) ) x. ( C x. (; 1 0 ^ M ) ) ) e. CC
18	15, 10	mulcli 8296	. . . 4 |- ( D x. B ) e. CC
19	8, 15	mulcli 8296	. . . . 5 |- ( ( A x. (; 1 0 ^ M ) ) x. D ) e. CC
20	13, 10	mulcli 8296	. . . . 5 |- ( ( C x. (; 1 0 ^ M ) ) x. B ) e. CC
21	19, 20	addcli 8295	. . . 4 |- ( ( ( A x. (; 1 0 ^ M ) ) x. D ) + ( ( C x. (; 1 0 ^ M ) ) x. B ) ) e. CC
22	17, 18, 21	add32i 8455	. . 3 |- ( ( ( ( A x. (; 1 0 ^ M ) ) x. ( C x. (; 1 0 ^ M ) ) ) + ( D x. B ) ) + ( ( ( A x. (; 1 0 ^ M ) ) x. D ) + ( ( C x. (; 1 0 ^ M ) ) x. B ) ) ) = ( ( ( ( A x. (; 1 0 ^ M ) ) x. ( C x. (; 1 0 ^ M ) ) ) + ( ( ( A x. (; 1 0 ^ M ) ) x. D ) + ( ( C x. (; 1 0 ^ M ) ) x. B ) ) ) + ( D x. B ) )
23	8, 12	mulcli 8296	. . . . . 6 |- ( ( A x. (; 1 0 ^ M ) ) x. C ) e. CC
24		karatsuba.s	. . . . . . 7 |- S e. NN0
25	24	nn0cni 9529	. . . . . 6 |- S e. CC
26	23, 25, 7	adddiri 8302	. . . . 5 |- ( ( ( ( A x. (; 1 0 ^ M ) ) x. C ) + S ) x. (; 1 0 ^ M ) ) = ( ( ( ( A x. (; 1 0 ^ M ) ) x. C ) x. (; 1 0 ^ M ) ) + ( S x. (; 1 0 ^ M ) ) )
27	2, 7, 12	mul32i 8438	. . . . . . . . 9 |- ( ( A x. (; 1 0 ^ M ) ) x. C ) = ( ( A x. C ) x. (; 1 0 ^ M ) )
28		karatsuba.r	. . . . . . . . . 10 |- ( A x. C ) = R
29	28	oveq1i 6069	. . . . . . . . 9 |- ( ( A x. C ) x. (; 1 0 ^ M ) ) = ( R x. (; 1 0 ^ M ) )
30	27, 29	eqtri 2255	. . . . . . . 8 |- ( ( A x. (; 1 0 ^ M ) ) x. C ) = ( R x. (; 1 0 ^ M ) )
31	30	oveq1i 6069	. . . . . . 7 |- ( ( ( A x. (; 1 0 ^ M ) ) x. C ) + S ) = ( ( R x. (; 1 0 ^ M ) ) + S )
32		karatsuba.w	. . . . . . 7 |- ( ( R x. (; 1 0 ^ M ) ) + S ) = W
33	31, 32	eqtri 2255	. . . . . 6 |- ( ( ( A x. (; 1 0 ^ M ) ) x. C ) + S ) = W
34	33	oveq1i 6069	. . . . 5 |- ( ( ( ( A x. (; 1 0 ^ M ) ) x. C ) + S ) x. (; 1 0 ^ M ) ) = ( W x. (; 1 0 ^ M ) )
35	8, 12, 7	mulassi 8300	. . . . . 6 |- ( ( ( A x. (; 1 0 ^ M ) ) x. C ) x. (; 1 0 ^ M ) ) = ( ( A x. (; 1 0 ^ M ) ) x. ( C x. (; 1 0 ^ M ) ) )
36	2, 12	mulcli 8296	. . . . . . . . . . . 12 |- ( A x. C ) e. CC
37	36, 18, 25	add32i 8455	. . . . . . . . . . 11 |- ( ( ( A x. C ) + ( D x. B ) ) + S ) = ( ( ( A x. C ) + S ) + ( D x. B ) )
38	28	oveq1i 6069	. . . . . . . . . . . 12 |- ( ( A x. C ) + S ) = ( R + S )
39		karatsuba.t	. . . . . . . . . . . . 13 |- ( B x. D ) = T
40	10, 15, 39	mulcomli 8298	. . . . . . . . . . . 12 |- ( D x. B ) = T
41	38, 40	oveq12i 6071	. . . . . . . . . . 11 |- ( ( ( A x. C ) + S ) + ( D x. B ) ) = ( ( R + S ) + T )
42	37, 41	eqtri 2255	. . . . . . . . . 10 |- ( ( ( A x. C ) + ( D x. B ) ) + S ) = ( ( R + S ) + T )
43		karatsuba.e	. . . . . . . . . 10 |- ( ( A + B ) x. ( C + D ) ) = ( ( R + S ) + T )
44	2, 10, 12, 15	muladdi 8701	. . . . . . . . . 10 |- ( ( A + B ) x. ( C + D ) ) = ( ( ( A x. C ) + ( D x. B ) ) + ( ( A x. D ) + ( C x. B ) ) )
45	42, 43, 44	3eqtr2i 2261	. . . . . . . . 9 |- ( ( ( A x. C ) + ( D x. B ) ) + S ) = ( ( ( A x. C ) + ( D x. B ) ) + ( ( A x. D ) + ( C x. B ) ) )
46	36, 18	addcli 8295	. . . . . . . . . 10 |- ( ( A x. C ) + ( D x. B ) ) e. CC
47	2, 15	mulcli 8296	. . . . . . . . . . 11 |- ( A x. D ) e. CC
48	12, 10	mulcli 8296	. . . . . . . . . . 11 |- ( C x. B ) e. CC
49	47, 48	addcli 8295	. . . . . . . . . 10 |- ( ( A x. D ) + ( C x. B ) ) e. CC
50	46, 25, 49	addcani 8473	. . . . . . . . 9 |- ( ( ( ( A x. C ) + ( D x. B ) ) + S ) = ( ( ( A x. C ) + ( D x. B ) ) + ( ( A x. D ) + ( C x. B ) ) ) <-> S = ( ( A x. D ) + ( C x. B ) ) )
51	45, 50	mpbi 145	. . . . . . . 8 |- S = ( ( A x. D ) + ( C x. B ) )
52	51	oveq1i 6069	. . . . . . 7 |- ( S x. (; 1 0 ^ M ) ) = ( ( ( A x. D ) + ( C x. B ) ) x. (; 1 0 ^ M ) )
53	47, 48, 7	adddiri 8302	. . . . . . 7 |- ( ( ( A x. D ) + ( C x. B ) ) x. (; 1 0 ^ M ) ) = ( ( ( A x. D ) x. (; 1 0 ^ M ) ) + ( ( C x. B ) x. (; 1 0 ^ M ) ) )
54	2, 15, 7	mul32i 8438	. . . . . . . 8 |- ( ( A x. D ) x. (; 1 0 ^ M ) ) = ( ( A x. (; 1 0 ^ M ) ) x. D )
55	12, 10, 7	mul32i 8438	. . . . . . . 8 |- ( ( C x. B ) x. (; 1 0 ^ M ) ) = ( ( C x. (; 1 0 ^ M ) ) x. B )
56	54, 55	oveq12i 6071	. . . . . . 7 |- ( ( ( A x. D ) x. (; 1 0 ^ M ) ) + ( ( C x. B ) x. (; 1 0 ^ M ) ) ) = ( ( ( A x. (; 1 0 ^ M ) ) x. D ) + ( ( C x. (; 1 0 ^ M ) ) x. B ) )
57	52, 53, 56	3eqtri 2259	. . . . . 6 |- ( S x. (; 1 0 ^ M ) ) = ( ( ( A x. (; 1 0 ^ M ) ) x. D ) + ( ( C x. (; 1 0 ^ M ) ) x. B ) )
58	35, 57	oveq12i 6071	. . . . 5 |- ( ( ( ( A x. (; 1 0 ^ M ) ) x. C ) x. (; 1 0 ^ M ) ) + ( S x. (; 1 0 ^ M ) ) ) = ( ( ( A x. (; 1 0 ^ M ) ) x. ( C x. (; 1 0 ^ M ) ) ) + ( ( ( A x. (; 1 0 ^ M ) ) x. D ) + ( ( C x. (; 1 0 ^ M ) ) x. B ) ) )
59	26, 34, 58	3eqtr3ri 2264	. . . 4 |- ( ( ( A x. (; 1 0 ^ M ) ) x. ( C x. (; 1 0 ^ M ) ) ) + ( ( ( A x. (; 1 0 ^ M ) ) x. D ) + ( ( C x. (; 1 0 ^ M ) ) x. B ) ) ) = ( W x. (; 1 0 ^ M ) )
60	59, 40	oveq12i 6071	. . 3 |- ( ( ( ( A x. (; 1 0 ^ M ) ) x. ( C x. (; 1 0 ^ M ) ) ) + ( ( ( A x. (; 1 0 ^ M ) ) x. D ) + ( ( C x. (; 1 0 ^ M ) ) x. B ) ) ) + ( D x. B ) ) = ( ( W x. (; 1 0 ^ M ) ) + T )
61	16, 22, 60	3eqtri 2259	. 2 |- ( ( ( A x. (; 1 0 ^ M ) ) + B ) x. ( ( C x. (; 1 0 ^ M ) ) + D ) ) = ( ( W x. (; 1 0 ^ M ) ) + T )
62		karatsuba.x	. . 3 |- ( ( A x. (; 1 0 ^ M ) ) + B ) = X
63		karatsuba.y	. . 3 |- ( ( C x. (; 1 0 ^ M ) ) + D ) = Y
64	62, 63	oveq12i 6071	. 2 |- ( ( ( A x. (; 1 0 ^ M ) ) + B ) x. ( ( C x. (; 1 0 ^ M ) ) + D ) ) = ( X x. Y )
65		karatsuba.z	. 2 |- ( ( W x. (; 1 0 ^ M ) ) + T ) = Z
66	61, 64, 65	3eqtr3i 2263	1 |- ( X x. Y ) = Z

Syntax hints:   = wceq 1398   e. wcel 2205  (class class class)co 6059   CCcc 8142   0cc0 8144   1c1 8145   + caddc 8147   x. cmul 8149   NN0cn0 9517  ;cdc 9731   ^cexp 10928

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 4231  ax-sep 4234  ax-nul 4242  ax-pow 4293  ax-pr 4328  ax-un 4560  ax-setind 4665  ax-iinf 4716  ax-cnex 8235  ax-resscn 8236  ax-1cn 8237  ax-1re 8238  ax-icn 8239  ax-addcl 8240  ax-addrcl 8241  ax-mulcl 8242  ax-mulrcl 8243  ax-addcom 8244  ax-mulcom 8245  ax-addass 8246  ax-mulass 8247  ax-distr 8248  ax-i2m1 8249  ax-0lt1 8250  ax-1rid 8251  ax-0id 8252  ax-rnegex 8253  ax-precex 8254  ax-cnre 8255  ax-pre-ltirr 8256  ax-pre-ltwlin 8257  ax-pre-lttrn 8258  ax-pre-apti 8259  ax-pre-ltadd 8260  ax-pre-mulgt0 8261  ax-pre-mulext 8262

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 3626  df-pw 3677  df-sn 3701  df-pr 3702  df-op 3704  df-uni 3921  df-int 3956  df-iun 3999  df-br 4116  df-opab 4178  df-mpt 4179  df-tr 4215  df-id 4420  df-po 4423  df-iso 4424  df-iord 4493  df-on 4495  df-ilim 4496  df-suc 4498  df-iom 4719  df-xp 4761  df-rel 4762  df-cnv 4763  df-co 4764  df-dm 4765  df-rn 4766  df-res 4767  df-ima 4768  df-iota 5318  df-fun 5360  df-fn 5361  df-f 5362  df-f1 5363  df-fo 5364  df-f1o 5365  df-fv 5366  df-riota 6012  df-ov 6062  df-oprab 6063  df-mpo 6064  df-1st 6348  df-2nd 6349  df-recs 6550  df-frec 6636  df-pnf 8327  df-mnf 8328  df-xr 8329  df-ltxr 8330  df-le 8331  df-sub 8464  df-neg 8465  df-reap 8868  df-ap 8875  df-div 8968  df-inn 9259  df-2 9317  df-3 9318  df-4 9319  df-5 9320  df-6 9321  df-7 9322  df-8 9323  df-9 9324  df-n0 9518  df-z 9599  df-dec 9732  df-uz 9876  df-seqfrec 10838  df-exp 10929

This theorem is referenced by: (None)