Theorem karatsuba 17052

Description: The Karatsuba multiplication algorithm. If 𝑋 and 𝑌 are decomposed into two groups of digits of length 𝑀 (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 𝑋𝑌 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 12707. (Contributed by Mario Carneiro, 16-Jul-2015.) (Revised by AV, 9-Sep-2021.)

Hypotheses

Ref	Expression
karatsuba.a	⊢ 𝐴 ∈ ℕ0
karatsuba.b	⊢ 𝐵 ∈ ℕ0
karatsuba.c	⊢ 𝐶 ∈ ℕ0
karatsuba.d	⊢ 𝐷 ∈ ℕ0
karatsuba.s	⊢ 𝑆 ∈ ℕ0
karatsuba.m	⊢ 𝑀 ∈ ℕ0
karatsuba.r	⊢ (𝐴 · 𝐶) = 𝑅
karatsuba.t	⊢ (𝐵 · 𝐷) = 𝑇
karatsuba.e	⊢ ((𝐴 + 𝐵) · (𝐶 + 𝐷)) = ((𝑅 + 𝑆) + 𝑇)
karatsuba.x	⊢ ((𝐴 · (;10↑𝑀)) + 𝐵) = 𝑋
karatsuba.y	⊢ ((𝐶 · (;10↑𝑀)) + 𝐷) = 𝑌
karatsuba.w	⊢ ((𝑅 · (;10↑𝑀)) + 𝑆) = 𝑊
karatsuba.z	⊢ ((𝑊 · (;10↑𝑀)) + 𝑇) = 𝑍

Assertion

Ref	Expression
karatsuba	⊢ (𝑋 · 𝑌) = 𝑍

Proof of Theorem karatsuba

Step	Hyp	Ref	Expression
1		karatsuba.a	. . . . . 6 ⊢ 𝐴 ∈ ℕ0
2	1	nn0cni 12447	. . . . 5 ⊢ 𝐴 ∈ ℂ
3		10nn0 12660	. . . . . . 7 ⊢ ;10 ∈ ℕ0
4	3	nn0cni 12447	. . . . . 6 ⊢ ;10 ∈ ℂ
5		karatsuba.m	. . . . . 6 ⊢ 𝑀 ∈ ℕ0
6		expcl 14039	. . . . . 6 ⊢ ((;10 ∈ ℂ ∧ 𝑀 ∈ ℕ0) → (;10↑𝑀) ∈ ℂ)
7	4, 5, 6	mp2an 698	. . . . 5 ⊢ (;10↑𝑀) ∈ ℂ
8	2, 7	mulcli 11150	. . . 4 ⊢ (𝐴 · (;10↑𝑀)) ∈ ℂ
9		karatsuba.b	. . . . 5 ⊢ 𝐵 ∈ ℕ0
10	9	nn0cni 12447	. . . 4 ⊢ 𝐵 ∈ ℂ
11		karatsuba.c	. . . . . 6 ⊢ 𝐶 ∈ ℕ0
12	11	nn0cni 12447	. . . . 5 ⊢ 𝐶 ∈ ℂ
13	12, 7	mulcli 11150	. . . 4 ⊢ (𝐶 · (;10↑𝑀)) ∈ ℂ
14		karatsuba.d	. . . . 5 ⊢ 𝐷 ∈ ℕ0
15	14	nn0cni 12447	. . . 4 ⊢ 𝐷 ∈ ℂ
16	8, 10, 13, 15	muladdi 11599	. . 3 ⊢ (((𝐴 · (;10↑𝑀)) + 𝐵) · ((𝐶 · (;10↑𝑀)) + 𝐷)) = ((((𝐴 · (;10↑𝑀)) · (𝐶 · (;10↑𝑀))) + (𝐷 · 𝐵)) + (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵)))
17	8, 13	mulcli 11150	. . . 4 ⊢ ((𝐴 · (;10↑𝑀)) · (𝐶 · (;10↑𝑀))) ∈ ℂ
18	15, 10	mulcli 11150	. . . 4 ⊢ (𝐷 · 𝐵) ∈ ℂ
19	8, 15	mulcli 11150	. . . . 5 ⊢ ((𝐴 · (;10↑𝑀)) · 𝐷) ∈ ℂ
20	13, 10	mulcli 11150	. . . . 5 ⊢ ((𝐶 · (;10↑𝑀)) · 𝐵) ∈ ℂ
21	19, 20	addcli 11149	. . . 4 ⊢ (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵)) ∈ ℂ
22	17, 18, 21	add32i 11368	. . 3 ⊢ ((((𝐴 · (;10↑𝑀)) · (𝐶 · (;10↑𝑀))) + (𝐷 · 𝐵)) + (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵))) = ((((𝐴 · (;10↑𝑀)) · (𝐶 · (;10↑𝑀))) + (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵))) + (𝐷 · 𝐵))
23	8, 12	mulcli 11150	. . . . . 6 ⊢ ((𝐴 · (;10↑𝑀)) · 𝐶) ∈ ℂ
24		karatsuba.s	. . . . . . 7 ⊢ 𝑆 ∈ ℕ0
25	24	nn0cni 12447	. . . . . 6 ⊢ 𝑆 ∈ ℂ
26	23, 25, 7	adddiri 11156	. . . . 5 ⊢ ((((𝐴 · (;10↑𝑀)) · 𝐶) + 𝑆) · (;10↑𝑀)) = ((((𝐴 · (;10↑𝑀)) · 𝐶) · (;10↑𝑀)) + (𝑆 · (;10↑𝑀)))
27	2, 7, 12	mul32i 11340	. . . . . . . . 9 ⊢ ((𝐴 · (;10↑𝑀)) · 𝐶) = ((𝐴 · 𝐶) · (;10↑𝑀))
28		karatsuba.r	. . . . . . . . . 10 ⊢ (𝐴 · 𝐶) = 𝑅
29	28	oveq1i 7373	. . . . . . . . 9 ⊢ ((𝐴 · 𝐶) · (;10↑𝑀)) = (𝑅 · (;10↑𝑀))
30	27, 29	eqtri 2763	. . . . . . . 8 ⊢ ((𝐴 · (;10↑𝑀)) · 𝐶) = (𝑅 · (;10↑𝑀))
31	30	oveq1i 7373	. . . . . . 7 ⊢ (((𝐴 · (;10↑𝑀)) · 𝐶) + 𝑆) = ((𝑅 · (;10↑𝑀)) + 𝑆)
32		karatsuba.w	. . . . . . 7 ⊢ ((𝑅 · (;10↑𝑀)) + 𝑆) = 𝑊
33	31, 32	eqtri 2763	. . . . . 6 ⊢ (((𝐴 · (;10↑𝑀)) · 𝐶) + 𝑆) = 𝑊
34	33	oveq1i 7373	. . . . 5 ⊢ ((((𝐴 · (;10↑𝑀)) · 𝐶) + 𝑆) · (;10↑𝑀)) = (𝑊 · (;10↑𝑀))
35	8, 12, 7	mulassi 11154	. . . . . 6 ⊢ (((𝐴 · (;10↑𝑀)) · 𝐶) · (;10↑𝑀)) = ((𝐴 · (;10↑𝑀)) · (𝐶 · (;10↑𝑀)))
36	2, 12	mulcli 11150	. . . . . . . . . . . 12 ⊢ (𝐴 · 𝐶) ∈ ℂ
37	36, 18, 25	add32i 11368	. . . . . . . . . . 11 ⊢ (((𝐴 · 𝐶) + (𝐷 · 𝐵)) + 𝑆) = (((𝐴 · 𝐶) + 𝑆) + (𝐷 · 𝐵))
38	28	oveq1i 7373	. . . . . . . . . . . 12 ⊢ ((𝐴 · 𝐶) + 𝑆) = (𝑅 + 𝑆)
39		karatsuba.t	. . . . . . . . . . . . 13 ⊢ (𝐵 · 𝐷) = 𝑇
40	10, 15, 39	mulcomli 11152	. . . . . . . . . . . 12 ⊢ (𝐷 · 𝐵) = 𝑇
41	38, 40	oveq12i 7375	. . . . . . . . . . 11 ⊢ (((𝐴 · 𝐶) + 𝑆) + (𝐷 · 𝐵)) = ((𝑅 + 𝑆) + 𝑇)
42	37, 41	eqtri 2763	. . . . . . . . . 10 ⊢ (((𝐴 · 𝐶) + (𝐷 · 𝐵)) + 𝑆) = ((𝑅 + 𝑆) + 𝑇)
43		karatsuba.e	. . . . . . . . . 10 ⊢ ((𝐴 + 𝐵) · (𝐶 + 𝐷)) = ((𝑅 + 𝑆) + 𝑇)
44	2, 10, 12, 15	muladdi 11599	. . . . . . . . . 10 ⊢ ((𝐴 + 𝐵) · (𝐶 + 𝐷)) = (((𝐴 · 𝐶) + (𝐷 · 𝐵)) + ((𝐴 · 𝐷) + (𝐶 · 𝐵)))
45	42, 43, 44	3eqtr2i 2769	. . . . . . . . 9 ⊢ (((𝐴 · 𝐶) + (𝐷 · 𝐵)) + 𝑆) = (((𝐴 · 𝐶) + (𝐷 · 𝐵)) + ((𝐴 · 𝐷) + (𝐶 · 𝐵)))
46	36, 18	addcli 11149	. . . . . . . . . 10 ⊢ ((𝐴 · 𝐶) + (𝐷 · 𝐵)) ∈ ℂ
47	2, 15	mulcli 11150	. . . . . . . . . . 11 ⊢ (𝐴 · 𝐷) ∈ ℂ
48	12, 10	mulcli 11150	. . . . . . . . . . 11 ⊢ (𝐶 · 𝐵) ∈ ℂ
49	47, 48	addcli 11149	. . . . . . . . . 10 ⊢ ((𝐴 · 𝐷) + (𝐶 · 𝐵)) ∈ ℂ
50	46, 25, 49	addcani 11337	. . . . . . . . 9 ⊢ ((((𝐴 · 𝐶) + (𝐷 · 𝐵)) + 𝑆) = (((𝐴 · 𝐶) + (𝐷 · 𝐵)) + ((𝐴 · 𝐷) + (𝐶 · 𝐵))) ↔ 𝑆 = ((𝐴 · 𝐷) + (𝐶 · 𝐵)))
51	45, 50	mpbi 231	. . . . . . . 8 ⊢ 𝑆 = ((𝐴 · 𝐷) + (𝐶 · 𝐵))
52	51	oveq1i 7373	. . . . . . 7 ⊢ (𝑆 · (;10↑𝑀)) = (((𝐴 · 𝐷) + (𝐶 · 𝐵)) · (;10↑𝑀))
53	47, 48, 7	adddiri 11156	. . . . . . 7 ⊢ (((𝐴 · 𝐷) + (𝐶 · 𝐵)) · (;10↑𝑀)) = (((𝐴 · 𝐷) · (;10↑𝑀)) + ((𝐶 · 𝐵) · (;10↑𝑀)))
54	2, 15, 7	mul32i 11340	. . . . . . . 8 ⊢ ((𝐴 · 𝐷) · (;10↑𝑀)) = ((𝐴 · (;10↑𝑀)) · 𝐷)
55	12, 10, 7	mul32i 11340	. . . . . . . 8 ⊢ ((𝐶 · 𝐵) · (;10↑𝑀)) = ((𝐶 · (;10↑𝑀)) · 𝐵)
56	54, 55	oveq12i 7375	. . . . . . 7 ⊢ (((𝐴 · 𝐷) · (;10↑𝑀)) + ((𝐶 · 𝐵) · (;10↑𝑀))) = (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵))
57	52, 53, 56	3eqtri 2767	. . . . . 6 ⊢ (𝑆 · (;10↑𝑀)) = (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵))
58	35, 57	oveq12i 7375	. . . . 5 ⊢ ((((𝐴 · (;10↑𝑀)) · 𝐶) · (;10↑𝑀)) + (𝑆 · (;10↑𝑀))) = (((𝐴 · (;10↑𝑀)) · (𝐶 · (;10↑𝑀))) + (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵)))
59	26, 34, 58	3eqtr3ri 2772	. . . 4 ⊢ (((𝐴 · (;10↑𝑀)) · (𝐶 · (;10↑𝑀))) + (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵))) = (𝑊 · (;10↑𝑀))
60	59, 40	oveq12i 7375	. . 3 ⊢ ((((𝐴 · (;10↑𝑀)) · (𝐶 · (;10↑𝑀))) + (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵))) + (𝐷 · 𝐵)) = ((𝑊 · (;10↑𝑀)) + 𝑇)
61	16, 22, 60	3eqtri 2767	. 2 ⊢ (((𝐴 · (;10↑𝑀)) + 𝐵) · ((𝐶 · (;10↑𝑀)) + 𝐷)) = ((𝑊 · (;10↑𝑀)) + 𝑇)
62		karatsuba.x	. . 3 ⊢ ((𝐴 · (;10↑𝑀)) + 𝐵) = 𝑋
63		karatsuba.y	. . 3 ⊢ ((𝐶 · (;10↑𝑀)) + 𝐷) = 𝑌
64	62, 63	oveq12i 7375	. 2 ⊢ (((𝐴 · (;10↑𝑀)) + 𝐵) · ((𝐶 · (;10↑𝑀)) + 𝐷)) = (𝑋 · 𝑌)
65		karatsuba.z	. 2 ⊢ ((𝑊 · (;10↑𝑀)) + 𝑇) = 𝑍
66	61, 64, 65	3eqtr3i 2771	1 ⊢ (𝑋 · 𝑌) = 𝑍

Syntax hints:   = wceq 1547   ∈ wcel 2119  (class class class)co 7363  ℂcc 11034  0cc0 11036  1c1 11037   + caddc 11039   · cmul 11041  ℕ₀cn0 12435  ;cdc 12642  ↑cexp 14021

This theorem was proved from axioms:  ax-mp 5  ax-1 6  ax-2 7  ax-3 8  ax-gen 1802  ax-4 1816  ax-5 1917  ax-6 1974  ax-7 2015  ax-8 2121  ax-9 2129  ax-10 2152  ax-11 2168  ax-12 2189  ax-ext 2712  ax-sep 5225  ax-nul 5235  ax-pow 5301  ax-pr 5369  ax-un 7685  ax-cnex 11092  ax-resscn 11093  ax-1cn 11094  ax-icn 11095  ax-addcl 11096  ax-addrcl 11097  ax-mulcl 11098  ax-mulrcl 11099  ax-mulcom 11100  ax-addass 11101  ax-mulass 11102  ax-distr 11103  ax-i2m1 11104  ax-1ne0 11105  ax-1rid 11106  ax-rnegex 11107  ax-rrecex 11108  ax-cnre 11109  ax-pre-lttri 11110  ax-pre-lttrn 11111  ax-pre-ltadd 11112  ax-pre-mulgt0 11113

This theorem depends on definitions:  df-bi 208  df-an 397  df-or 854  df-3or 1093  df-3an 1094  df-tru 1550  df-fal 1560  df-ex 1787  df-nf 1791  df-sb 2074  df-mo 2543  df-eu 2573  df-clab 2719  df-cleq 2732  df-clel 2815  df-nfc 2889  df-ne 2936  df-nel 3040  df-ral 3055  df-rex 3065  df-reu 3346  df-rab 3393  df-v 3434  df-sbc 3731  df-csb 3839  df-dif 3893  df-un 3895  df-in 3897  df-ss 3907  df-pss 3910  df-nul 4269  df-if 4462  df-pw 4538  df-sn 4563  df-pr 4565  df-op 4569  df-uni 4846  df-iun 4930  df-br 5080  df-opab 5142  df-mpt 5161  df-tr 5187  df-id 5520  df-eprel 5525  df-po 5533  df-so 5534  df-fr 5578  df-we 5580  df-xp 5631  df-rel 5632  df-cnv 5633  df-co 5634  df-dm 5635  df-rn 5636  df-res 5637  df-ima 5638  df-pred 6259  df-ord 6320  df-on 6321  df-lim 6322  df-suc 6323  df-iota 6448  df-fun 6494  df-fn 6495  df-f 6496  df-f1 6497  df-fo 6498  df-f1o 6499  df-fv 6500  df-riota 7320  df-ov 7366  df-oprab 7367  df-mpo 7368  df-om 7814  df-2nd 7939  df-frecs 8228  df-wrecs 8259  df-recs 8308  df-rdg 8346  df-er 8640  df-en 8891  df-dom 8892  df-sdom 8893  df-pnf 11179  df-mnf 11180  df-xr 11181  df-ltxr 11182  df-le 11183  df-sub 11377  df-neg 11378  df-nn 12173  df-2 12242  df-3 12243  df-4 12244  df-5 12245  df-6 12246  df-7 12247  df-8 12248  df-9 12249  df-n0 12436  df-z 12523  df-dec 12643  df-uz 12787  df-seq 13962  df-exp 14022

This theorem is referenced by:  dpmul4  32999