Theorem karatsuba 13064

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 9718. (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 9457	. . . . 5 |- A e. CC
3		10nn0 9671	. . . . . . 7 |- ; 1 0 e. NN0
4	3	nn0cni 9457	. . . . . 6 |- ; 1 0 e. CC
5		karatsuba.m	. . . . . 6 |- M e. NN0
6		expcl 10863	. . . . . 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 8227	. . . 4 |- ( A x. (; 1 0 ^ M ) ) e. CC
9		karatsuba.b	. . . . 5 |- B e. NN0
10	9	nn0cni 9457	. . . 4 |- B e. CC
11		karatsuba.c	. . . . . 6 |- C e. NN0
12	11	nn0cni 9457	. . . . 5 |- C e. CC
13	12, 7	mulcli 8227	. . . 4 |- ( C x. (; 1 0 ^ M ) ) e. CC
14		karatsuba.d	. . . . 5 |- D e. NN0
15	14	nn0cni 9457	. . . 4 |- D e. CC
16	8, 10, 13, 15	muladdi 8631	. . 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 8227	. . . 4 |- ( ( A x. (; 1 0 ^ M ) ) x. ( C x. (; 1 0 ^ M ) ) ) e. CC
18	15, 10	mulcli 8227	. . . 4 |- ( D x. B ) e. CC
19	8, 15	mulcli 8227	. . . . 5 |- ( ( A x. (; 1 0 ^ M ) ) x. D ) e. CC
20	13, 10	mulcli 8227	. . . . 5 |- ( ( C x. (; 1 0 ^ M ) ) x. B ) e. CC
21	19, 20	addcli 8226	. . . 4 |- ( ( ( A x. (; 1 0 ^ M ) ) x. D ) + ( ( C x. (; 1 0 ^ M ) ) x. B ) ) e. CC
22	17, 18, 21	add32i 8386	. . 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 8227	. . . . . 6 |- ( ( A x. (; 1 0 ^ M ) ) x. C ) e. CC
24		karatsuba.s	. . . . . . 7 |- S e. NN0
25	24	nn0cni 9457	. . . . . 6 |- S e. CC
26	23, 25, 7	adddiri 8233	. . . . 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 8369	. . . . . . . . 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 6038	. . . . . . . . 9 |- ( ( A x. C ) x. (; 1 0 ^ M ) ) = ( R x. (; 1 0 ^ M ) )
30	27, 29	eqtri 2252	. . . . . . . 8 |- ( ( A x. (; 1 0 ^ M ) ) x. C ) = ( R x. (; 1 0 ^ M ) )
31	30	oveq1i 6038	. . . . . . 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 2252	. . . . . 6 |- ( ( ( A x. (; 1 0 ^ M ) ) x. C ) + S ) = W
34	33	oveq1i 6038	. . . . 5 |- ( ( ( ( A x. (; 1 0 ^ M ) ) x. C ) + S ) x. (; 1 0 ^ M ) ) = ( W x. (; 1 0 ^ M ) )
35	8, 12, 7	mulassi 8231	. . . . . 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 8227	. . . . . . . . . . . 12 |- ( A x. C ) e. CC
37	36, 18, 25	add32i 8386	. . . . . . . . . . 11 |- ( ( ( A x. C ) + ( D x. B ) ) + S ) = ( ( ( A x. C ) + S ) + ( D x. B ) )
38	28	oveq1i 6038	. . . . . . . . . . . 12 |- ( ( A x. C ) + S ) = ( R + S )
39		karatsuba.t	. . . . . . . . . . . . 13 |- ( B x. D ) = T
40	10, 15, 39	mulcomli 8229	. . . . . . . . . . . 12 |- ( D x. B ) = T
41	38, 40	oveq12i 6040	. . . . . . . . . . 11 |- ( ( ( A x. C ) + S ) + ( D x. B ) ) = ( ( R + S ) + T )
42	37, 41	eqtri 2252	. . . . . . . . . 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 8631	. . . . . . . . . 10 |- ( ( A + B ) x. ( C + D ) ) = ( ( ( A x. C ) + ( D x. B ) ) + ( ( A x. D ) + ( C x. B ) ) )
45	42, 43, 44	3eqtr2i 2258	. . . . . . . . 9 |- ( ( ( A x. C ) + ( D x. B ) ) + S ) = ( ( ( A x. C ) + ( D x. B ) ) + ( ( A x. D ) + ( C x. B ) ) )
46	36, 18	addcli 8226	. . . . . . . . . 10 |- ( ( A x. C ) + ( D x. B ) ) e. CC
47	2, 15	mulcli 8227	. . . . . . . . . . 11 |- ( A x. D ) e. CC
48	12, 10	mulcli 8227	. . . . . . . . . . 11 |- ( C x. B ) e. CC
49	47, 48	addcli 8226	. . . . . . . . . 10 |- ( ( A x. D ) + ( C x. B ) ) e. CC
50	46, 25, 49	addcani 8404	. . . . . . . . 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 6038	. . . . . . 7 |- ( S x. (; 1 0 ^ M ) ) = ( ( ( A x. D ) + ( C x. B ) ) x. (; 1 0 ^ M ) )
53	47, 48, 7	adddiri 8233	. . . . . . 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 8369	. . . . . . . 8 |- ( ( A x. D ) x. (; 1 0 ^ M ) ) = ( ( A x. (; 1 0 ^ M ) ) x. D )
55	12, 10, 7	mul32i 8369	. . . . . . . 8 |- ( ( C x. B ) x. (; 1 0 ^ M ) ) = ( ( C x. (; 1 0 ^ M ) ) x. B )
56	54, 55	oveq12i 6040	. . . . . . 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 2256	. . . . . 6 |- ( S x. (; 1 0 ^ M ) ) = ( ( ( A x. (; 1 0 ^ M ) ) x. D ) + ( ( C x. (; 1 0 ^ M ) ) x. B ) )
58	35, 57	oveq12i 6040	. . . . 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 2261	. . . 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 6040	. . 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 2256	. 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 6040	. 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 2260	1 |- ( X x. Y ) = Z

Syntax hints:   = wceq 1398   e. wcel 2202  (class class class)co 6028   CCcc 8073   0cc0 8075   1c1 8076   + caddc 8078   x. cmul 8080   NN0cn0 9445  ;cdc 9654   ^cexp 10844

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 2204  ax-14 2205  ax-ext 2213  ax-coll 4209  ax-sep 4212  ax-nul 4220  ax-pow 4270  ax-pr 4305  ax-un 4536  ax-setind 4641  ax-iinf 4692  ax-cnex 8166  ax-resscn 8167  ax-1cn 8168  ax-1re 8169  ax-icn 8170  ax-addcl 8171  ax-addrcl 8172  ax-mulcl 8173  ax-mulrcl 8174  ax-addcom 8175  ax-mulcom 8176  ax-addass 8177  ax-mulass 8178  ax-distr 8179  ax-i2m1 8180  ax-0lt1 8181  ax-1rid 8182  ax-0id 8183  ax-rnegex 8184  ax-precex 8185  ax-cnre 8186  ax-pre-ltirr 8187  ax-pre-ltwlin 8188  ax-pre-lttrn 8189  ax-pre-apti 8190  ax-pre-ltadd 8191  ax-pre-mulgt0 8192  ax-pre-mulext 8193

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 1811  df-eu 2082  df-mo 2083  df-clab 2218  df-cleq 2224  df-clel 2227  df-nfc 2364  df-ne 2404  df-nel 2499  df-ral 2516  df-rex 2517  df-reu 2518  df-rmo 2519  df-rab 2520  df-v 2805  df-sbc 3033  df-csb 3129  df-dif 3203  df-un 3205  df-in 3207  df-ss 3214  df-nul 3497  df-if 3608  df-pw 3658  df-sn 3679  df-pr 3680  df-op 3682  df-uni 3899  df-int 3934  df-iun 3977  df-br 4094  df-opab 4156  df-mpt 4157  df-tr 4193  df-id 4396  df-po 4399  df-iso 4400  df-iord 4469  df-on 4471  df-ilim 4472  df-suc 4474  df-iom 4695  df-xp 4737  df-rel 4738  df-cnv 4739  df-co 4740  df-dm 4741  df-rn 4742  df-res 4743  df-ima 4744  df-iota 5293  df-fun 5335  df-fn 5336  df-f 5337  df-f1 5338  df-fo 5339  df-f1o 5340  df-fv 5341  df-riota 5981  df-ov 6031  df-oprab 6032  df-mpo 6033  df-1st 6312  df-2nd 6313  df-recs 6514  df-frec 6600  df-pnf 8259  df-mnf 8260  df-xr 8261  df-ltxr 8262  df-le 8263  df-sub 8395  df-neg 8396  df-reap 8798  df-ap 8805  df-div 8896  df-inn 9187  df-2 9245  df-3 9246  df-4 9247  df-5 9248  df-6 9249  df-7 9250  df-8 9251  df-9 9252  df-n0 9446  df-z 9523  df-dec 9655  df-uz 9799  df-seqfrec 10754  df-exp 10845

This theorem is referenced by: (None)