Theorem karatsuba 12961

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 9650. (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 9389	. . . . 5 |- A e. CC
3		10nn0 9603	. . . . . . 7 |- ; 1 0 e. NN0
4	3	nn0cni 9389	. . . . . 6 |- ; 1 0 e. CC
5		karatsuba.m	. . . . . 6 |- M e. NN0
6		expcl 10787	. . . . . 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 8159	. . . 4 |- ( A x. (; 1 0 ^ M ) ) e. CC
9		karatsuba.b	. . . . 5 |- B e. NN0
10	9	nn0cni 9389	. . . 4 |- B e. CC
11		karatsuba.c	. . . . . 6 |- C e. NN0
12	11	nn0cni 9389	. . . . 5 |- C e. CC
13	12, 7	mulcli 8159	. . . 4 |- ( C x. (; 1 0 ^ M ) ) e. CC
14		karatsuba.d	. . . . 5 |- D e. NN0
15	14	nn0cni 9389	. . . 4 |- D e. CC
16	8, 10, 13, 15	muladdi 8563	. . 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 8159	. . . 4 |- ( ( A x. (; 1 0 ^ M ) ) x. ( C x. (; 1 0 ^ M ) ) ) e. CC
18	15, 10	mulcli 8159	. . . 4 |- ( D x. B ) e. CC
19	8, 15	mulcli 8159	. . . . 5 |- ( ( A x. (; 1 0 ^ M ) ) x. D ) e. CC
20	13, 10	mulcli 8159	. . . . 5 |- ( ( C x. (; 1 0 ^ M ) ) x. B ) e. CC
21	19, 20	addcli 8158	. . . 4 |- ( ( ( A x. (; 1 0 ^ M ) ) x. D ) + ( ( C x. (; 1 0 ^ M ) ) x. B ) ) e. CC
22	17, 18, 21	add32i 8318	. . 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 8159	. . . . . 6 |- ( ( A x. (; 1 0 ^ M ) ) x. C ) e. CC
24		karatsuba.s	. . . . . . 7 |- S e. NN0
25	24	nn0cni 9389	. . . . . 6 |- S e. CC
26	23, 25, 7	adddiri 8165	. . . . 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 8301	. . . . . . . . 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 6017	. . . . . . . . 9 |- ( ( A x. C ) x. (; 1 0 ^ M ) ) = ( R x. (; 1 0 ^ M ) )
30	27, 29	eqtri 2250	. . . . . . . 8 |- ( ( A x. (; 1 0 ^ M ) ) x. C ) = ( R x. (; 1 0 ^ M ) )
31	30	oveq1i 6017	. . . . . . 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 2250	. . . . . 6 |- ( ( ( A x. (; 1 0 ^ M ) ) x. C ) + S ) = W
34	33	oveq1i 6017	. . . . 5 |- ( ( ( ( A x. (; 1 0 ^ M ) ) x. C ) + S ) x. (; 1 0 ^ M ) ) = ( W x. (; 1 0 ^ M ) )
35	8, 12, 7	mulassi 8163	. . . . . 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 8159	. . . . . . . . . . . 12 |- ( A x. C ) e. CC
37	36, 18, 25	add32i 8318	. . . . . . . . . . 11 |- ( ( ( A x. C ) + ( D x. B ) ) + S ) = ( ( ( A x. C ) + S ) + ( D x. B ) )
38	28	oveq1i 6017	. . . . . . . . . . . 12 |- ( ( A x. C ) + S ) = ( R + S )
39		karatsuba.t	. . . . . . . . . . . . 13 |- ( B x. D ) = T
40	10, 15, 39	mulcomli 8161	. . . . . . . . . . . 12 |- ( D x. B ) = T
41	38, 40	oveq12i 6019	. . . . . . . . . . 11 |- ( ( ( A x. C ) + S ) + ( D x. B ) ) = ( ( R + S ) + T )
42	37, 41	eqtri 2250	. . . . . . . . . 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 8563	. . . . . . . . . 10 |- ( ( A + B ) x. ( C + D ) ) = ( ( ( A x. C ) + ( D x. B ) ) + ( ( A x. D ) + ( C x. B ) ) )
45	42, 43, 44	3eqtr2i 2256	. . . . . . . . 9 |- ( ( ( A x. C ) + ( D x. B ) ) + S ) = ( ( ( A x. C ) + ( D x. B ) ) + ( ( A x. D ) + ( C x. B ) ) )
46	36, 18	addcli 8158	. . . . . . . . . 10 |- ( ( A x. C ) + ( D x. B ) ) e. CC
47	2, 15	mulcli 8159	. . . . . . . . . . 11 |- ( A x. D ) e. CC
48	12, 10	mulcli 8159	. . . . . . . . . . 11 |- ( C x. B ) e. CC
49	47, 48	addcli 8158	. . . . . . . . . 10 |- ( ( A x. D ) + ( C x. B ) ) e. CC
50	46, 25, 49	addcani 8336	. . . . . . . . 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 6017	. . . . . . 7 |- ( S x. (; 1 0 ^ M ) ) = ( ( ( A x. D ) + ( C x. B ) ) x. (; 1 0 ^ M ) )
53	47, 48, 7	adddiri 8165	. . . . . . 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 8301	. . . . . . . 8 |- ( ( A x. D ) x. (; 1 0 ^ M ) ) = ( ( A x. (; 1 0 ^ M ) ) x. D )
55	12, 10, 7	mul32i 8301	. . . . . . . 8 |- ( ( C x. B ) x. (; 1 0 ^ M ) ) = ( ( C x. (; 1 0 ^ M ) ) x. B )
56	54, 55	oveq12i 6019	. . . . . . 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 2254	. . . . . 6 |- ( S x. (; 1 0 ^ M ) ) = ( ( ( A x. (; 1 0 ^ M ) ) x. D ) + ( ( C x. (; 1 0 ^ M ) ) x. B ) )
58	35, 57	oveq12i 6019	. . . . 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 2259	. . . 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 6019	. . 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 2254	. 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 6019	. 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 2258	1 |- ( X x. Y ) = Z

Syntax hints:   = wceq 1395   e. wcel 2200  (class class class)co 6007   CCcc 8005   0cc0 8007   1c1 8008   + caddc 8010   x. cmul 8012   NN0cn0 9377  ;cdc 9586   ^cexp 10768

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 617  ax-in2 618  ax-io 714  ax-5 1493  ax-7 1494  ax-gen 1495  ax-ie1 1539  ax-ie2 1540  ax-8 1550  ax-10 1551  ax-11 1552  ax-i12 1553  ax-bndl 1555  ax-4 1556  ax-17 1572  ax-i9 1576  ax-ial 1580  ax-i5r 1581  ax-13 2202  ax-14 2203  ax-ext 2211  ax-coll 4199  ax-sep 4202  ax-nul 4210  ax-pow 4258  ax-pr 4293  ax-un 4524  ax-setind 4629  ax-iinf 4680  ax-cnex 8098  ax-resscn 8099  ax-1cn 8100  ax-1re 8101  ax-icn 8102  ax-addcl 8103  ax-addrcl 8104  ax-mulcl 8105  ax-mulrcl 8106  ax-addcom 8107  ax-mulcom 8108  ax-addass 8109  ax-mulass 8110  ax-distr 8111  ax-i2m1 8112  ax-0lt1 8113  ax-1rid 8114  ax-0id 8115  ax-rnegex 8116  ax-precex 8117  ax-cnre 8118  ax-pre-ltirr 8119  ax-pre-ltwlin 8120  ax-pre-lttrn 8121  ax-pre-apti 8122  ax-pre-ltadd 8123  ax-pre-mulgt0 8124  ax-pre-mulext 8125

This theorem depends on definitions:  df-bi 117  df-dc 840  df-3or 1003  df-3an 1004  df-tru 1398  df-fal 1401  df-nf 1507  df-sb 1809  df-eu 2080  df-mo 2081  df-clab 2216  df-cleq 2222  df-clel 2225  df-nfc 2361  df-ne 2401  df-nel 2496  df-ral 2513  df-rex 2514  df-reu 2515  df-rmo 2516  df-rab 2517  df-v 2801  df-sbc 3029  df-csb 3125  df-dif 3199  df-un 3201  df-in 3203  df-ss 3210  df-nul 3492  df-if 3603  df-pw 3651  df-sn 3672  df-pr 3673  df-op 3675  df-uni 3889  df-int 3924  df-iun 3967  df-br 4084  df-opab 4146  df-mpt 4147  df-tr 4183  df-id 4384  df-po 4387  df-iso 4388  df-iord 4457  df-on 4459  df-ilim 4460  df-suc 4462  df-iom 4683  df-xp 4725  df-rel 4726  df-cnv 4727  df-co 4728  df-dm 4729  df-rn 4730  df-res 4731  df-ima 4732  df-iota 5278  df-fun 5320  df-fn 5321  df-f 5322  df-f1 5323  df-fo 5324  df-f1o 5325  df-fv 5326  df-riota 5960  df-ov 6010  df-oprab 6011  df-mpo 6012  df-1st 6292  df-2nd 6293  df-recs 6457  df-frec 6543  df-pnf 8191  df-mnf 8192  df-xr 8193  df-ltxr 8194  df-le 8195  df-sub 8327  df-neg 8328  df-reap 8730  df-ap 8737  df-div 8828  df-inn 9119  df-2 9177  df-3 9178  df-4 9179  df-5 9180  df-6 9181  df-7 9182  df-8 9183  df-9 9184  df-n0 9378  df-z 9455  df-dec 9587  df-uz 9731  df-seqfrec 10678  df-exp 10769

This theorem is referenced by: (None)