Theorem karatsuba 17144

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 12782. (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 12517	. . . . 5 ⊢ 𝐴 ∈ ℂ
3		10nn0 12734	. . . . . . 7 ⊢ ;10 ∈ ℕ0
4	3	nn0cni 12517	. . . . . 6 ⊢ ;10 ∈ ℂ
5		karatsuba.m	. . . . . 6 ⊢ 𝑀 ∈ ℕ0
6		expcl 14117	. . . . . 6 ⊢ ((;10 ∈ ℂ ∧ 𝑀 ∈ ℕ0) → (;10↑𝑀) ∈ ℂ)
7	4, 5, 6	mp2an 704	. . . . 5 ⊢ (;10↑𝑀) ∈ ℂ
8	2, 7	mulcli 11217	. . . 4 ⊢ (𝐴 · (;10↑𝑀)) ∈ ℂ
9		karatsuba.b	. . . . 5 ⊢ 𝐵 ∈ ℕ0
10	9	nn0cni 12517	. . . 4 ⊢ 𝐵 ∈ ℂ
11		karatsuba.c	. . . . . 6 ⊢ 𝐶 ∈ ℕ0
12	11	nn0cni 12517	. . . . 5 ⊢ 𝐶 ∈ ℂ
13	12, 7	mulcli 11217	. . . 4 ⊢ (𝐶 · (;10↑𝑀)) ∈ ℂ
14		karatsuba.d	. . . . 5 ⊢ 𝐷 ∈ ℕ0
15	14	nn0cni 12517	. . . 4 ⊢ 𝐷 ∈ ℂ
16	8, 10, 13, 15	muladdi 11666	. . 3 ⊢ (((𝐴 · (;10↑𝑀)) + 𝐵) · ((𝐶 · (;10↑𝑀)) + 𝐷)) = ((((𝐴 · (;10↑𝑀)) · (𝐶 · (;10↑𝑀))) + (𝐷 · 𝐵)) + (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵)))
17	8, 13	mulcli 11217	. . . 4 ⊢ ((𝐴 · (;10↑𝑀)) · (𝐶 · (;10↑𝑀))) ∈ ℂ
18	15, 10	mulcli 11217	. . . 4 ⊢ (𝐷 · 𝐵) ∈ ℂ
19	8, 15	mulcli 11217	. . . . 5 ⊢ ((𝐴 · (;10↑𝑀)) · 𝐷) ∈ ℂ
20	13, 10	mulcli 11217	. . . . 5 ⊢ ((𝐶 · (;10↑𝑀)) · 𝐵) ∈ ℂ
21	19, 20	addcli 11216	. . . 4 ⊢ (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵)) ∈ ℂ
22	17, 18, 21	add32i 11435	. . 3 ⊢ ((((𝐴 · (;10↑𝑀)) · (𝐶 · (;10↑𝑀))) + (𝐷 · 𝐵)) + (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵))) = ((((𝐴 · (;10↑𝑀)) · (𝐶 · (;10↑𝑀))) + (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵))) + (𝐷 · 𝐵))
23	8, 12	mulcli 11217	. . . . . 6 ⊢ ((𝐴 · (;10↑𝑀)) · 𝐶) ∈ ℂ
24		karatsuba.s	. . . . . . 7 ⊢ 𝑆 ∈ ℕ0
25	24	nn0cni 12517	. . . . . 6 ⊢ 𝑆 ∈ ℂ
26	23, 25, 7	adddiri 11223	. . . . 5 ⊢ ((((𝐴 · (;10↑𝑀)) · 𝐶) + 𝑆) · (;10↑𝑀)) = ((((𝐴 · (;10↑𝑀)) · 𝐶) · (;10↑𝑀)) + (𝑆 · (;10↑𝑀)))
27	2, 7, 12	mul32i 11407	. . . . . . . . 9 ⊢ ((𝐴 · (;10↑𝑀)) · 𝐶) = ((𝐴 · 𝐶) · (;10↑𝑀))
28		karatsuba.r	. . . . . . . . . 10 ⊢ (𝐴 · 𝐶) = 𝑅
29	28	oveq1i 7422	. . . . . . . . 9 ⊢ ((𝐴 · 𝐶) · (;10↑𝑀)) = (𝑅 · (;10↑𝑀))
30	27, 29	eqtri 2786	. . . . . . . 8 ⊢ ((𝐴 · (;10↑𝑀)) · 𝐶) = (𝑅 · (;10↑𝑀))
31	30	oveq1i 7422	. . . . . . 7 ⊢ (((𝐴 · (;10↑𝑀)) · 𝐶) + 𝑆) = ((𝑅 · (;10↑𝑀)) + 𝑆)
32		karatsuba.w	. . . . . . 7 ⊢ ((𝑅 · (;10↑𝑀)) + 𝑆) = 𝑊
33	31, 32	eqtri 2786	. . . . . 6 ⊢ (((𝐴 · (;10↑𝑀)) · 𝐶) + 𝑆) = 𝑊
34	33	oveq1i 7422	. . . . 5 ⊢ ((((𝐴 · (;10↑𝑀)) · 𝐶) + 𝑆) · (;10↑𝑀)) = (𝑊 · (;10↑𝑀))
35	8, 12, 7	mulassi 11221	. . . . . 6 ⊢ (((𝐴 · (;10↑𝑀)) · 𝐶) · (;10↑𝑀)) = ((𝐴 · (;10↑𝑀)) · (𝐶 · (;10↑𝑀)))
36	2, 12	mulcli 11217	. . . . . . . . . . . 12 ⊢ (𝐴 · 𝐶) ∈ ℂ
37	36, 18, 25	add32i 11435	. . . . . . . . . . 11 ⊢ (((𝐴 · 𝐶) + (𝐷 · 𝐵)) + 𝑆) = (((𝐴 · 𝐶) + 𝑆) + (𝐷 · 𝐵))
38	28	oveq1i 7422	. . . . . . . . . . . 12 ⊢ ((𝐴 · 𝐶) + 𝑆) = (𝑅 + 𝑆)
39		karatsuba.t	. . . . . . . . . . . . 13 ⊢ (𝐵 · 𝐷) = 𝑇
40	10, 15, 39	mulcomli 11219	. . . . . . . . . . . 12 ⊢ (𝐷 · 𝐵) = 𝑇
41	38, 40	oveq12i 7424	. . . . . . . . . . 11 ⊢ (((𝐴 · 𝐶) + 𝑆) + (𝐷 · 𝐵)) = ((𝑅 + 𝑆) + 𝑇)
42	37, 41	eqtri 2786	. . . . . . . . . 10 ⊢ (((𝐴 · 𝐶) + (𝐷 · 𝐵)) + 𝑆) = ((𝑅 + 𝑆) + 𝑇)
43		karatsuba.e	. . . . . . . . . 10 ⊢ ((𝐴 + 𝐵) · (𝐶 + 𝐷)) = ((𝑅 + 𝑆) + 𝑇)
44	2, 10, 12, 15	muladdi 11666	. . . . . . . . . 10 ⊢ ((𝐴 + 𝐵) · (𝐶 + 𝐷)) = (((𝐴 · 𝐶) + (𝐷 · 𝐵)) + ((𝐴 · 𝐷) + (𝐶 · 𝐵)))
45	42, 43, 44	3eqtr2i 2792	. . . . . . . . 9 ⊢ (((𝐴 · 𝐶) + (𝐷 · 𝐵)) + 𝑆) = (((𝐴 · 𝐶) + (𝐷 · 𝐵)) + ((𝐴 · 𝐷) + (𝐶 · 𝐵)))
46	36, 18	addcli 11216	. . . . . . . . . 10 ⊢ ((𝐴 · 𝐶) + (𝐷 · 𝐵)) ∈ ℂ
47	2, 15	mulcli 11217	. . . . . . . . . . 11 ⊢ (𝐴 · 𝐷) ∈ ℂ
48	12, 10	mulcli 11217	. . . . . . . . . . 11 ⊢ (𝐶 · 𝐵) ∈ ℂ
49	47, 48	addcli 11216	. . . . . . . . . 10 ⊢ ((𝐴 · 𝐷) + (𝐶 · 𝐵)) ∈ ℂ
50	46, 25, 49	addcani 11404	. . . . . . . . 9 ⊢ ((((𝐴 · 𝐶) + (𝐷 · 𝐵)) + 𝑆) = (((𝐴 · 𝐶) + (𝐷 · 𝐵)) + ((𝐴 · 𝐷) + (𝐶 · 𝐵))) ↔ 𝑆 = ((𝐴 · 𝐷) + (𝐶 · 𝐵)))
51	45, 50	mpbi 233	. . . . . . . 8 ⊢ 𝑆 = ((𝐴 · 𝐷) + (𝐶 · 𝐵))
52	51	oveq1i 7422	. . . . . . 7 ⊢ (𝑆 · (;10↑𝑀)) = (((𝐴 · 𝐷) + (𝐶 · 𝐵)) · (;10↑𝑀))
53	47, 48, 7	adddiri 11223	. . . . . . 7 ⊢ (((𝐴 · 𝐷) + (𝐶 · 𝐵)) · (;10↑𝑀)) = (((𝐴 · 𝐷) · (;10↑𝑀)) + ((𝐶 · 𝐵) · (;10↑𝑀)))
54	2, 15, 7	mul32i 11407	. . . . . . . 8 ⊢ ((𝐴 · 𝐷) · (;10↑𝑀)) = ((𝐴 · (;10↑𝑀)) · 𝐷)
55	12, 10, 7	mul32i 11407	. . . . . . . 8 ⊢ ((𝐶 · 𝐵) · (;10↑𝑀)) = ((𝐶 · (;10↑𝑀)) · 𝐵)
56	54, 55	oveq12i 7424	. . . . . . 7 ⊢ (((𝐴 · 𝐷) · (;10↑𝑀)) + ((𝐶 · 𝐵) · (;10↑𝑀))) = (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵))
57	52, 53, 56	3eqtri 2790	. . . . . 6 ⊢ (𝑆 · (;10↑𝑀)) = (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵))
58	35, 57	oveq12i 7424	. . . . 5 ⊢ ((((𝐴 · (;10↑𝑀)) · 𝐶) · (;10↑𝑀)) + (𝑆 · (;10↑𝑀))) = (((𝐴 · (;10↑𝑀)) · (𝐶 · (;10↑𝑀))) + (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵)))
59	26, 34, 58	3eqtr3ri 2795	. . . 4 ⊢ (((𝐴 · (;10↑𝑀)) · (𝐶 · (;10↑𝑀))) + (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵))) = (𝑊 · (;10↑𝑀))
60	59, 40	oveq12i 7424	. . 3 ⊢ ((((𝐴 · (;10↑𝑀)) · (𝐶 · (;10↑𝑀))) + (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵))) + (𝐷 · 𝐵)) = ((𝑊 · (;10↑𝑀)) + 𝑇)
61	16, 22, 60	3eqtri 2790	. 2 ⊢ (((𝐴 · (;10↑𝑀)) + 𝐵) · ((𝐶 · (;10↑𝑀)) + 𝐷)) = ((𝑊 · (;10↑𝑀)) + 𝑇)
62		karatsuba.x	. . 3 ⊢ ((𝐴 · (;10↑𝑀)) + 𝐵) = 𝑋
63		karatsuba.y	. . 3 ⊢ ((𝐶 · (;10↑𝑀)) + 𝐷) = 𝑌
64	62, 63	oveq12i 7424	. 2 ⊢ (((𝐴 · (;10↑𝑀)) + 𝐵) · ((𝐶 · (;10↑𝑀)) + 𝐷)) = (𝑋 · 𝑌)
65		karatsuba.z	. 2 ⊢ ((𝑊 · (;10↑𝑀)) + 𝑇) = 𝑍
66	61, 64, 65	3eqtr3i 2794	1 ⊢ (𝑋 · 𝑌) = 𝑍

Syntax hints:   = wceq 1570   ∈ wcel 2143  (class class class)co 7412  ℂcc 11099  0cc0 11101  1c1 11102   + caddc 11104   · cmul 11106  ℕ₀cn0 12505  ;cdc 12712  ↑cexp 14099

This theorem was proved from axioms:  ax-mp 5  ax-1 6  ax-2 7  ax-3 8  ax-gen 1825  ax-4 1839  ax-5 1940  ax-6 1997  ax-7 2038  ax-8 2145  ax-9 2153  ax-10 2176  ax-11 2192  ax-12 2213  ax-ext 2735  ax-sep 5258  ax-nul 5270  ax-pow 5338  ax-pr 5406  ax-un 7734  ax-cnex 11157  ax-resscn 11158  ax-1cn 11159  ax-icn 11160  ax-addcl 11161  ax-addrcl 11162  ax-mulcl 11163  ax-mulrcl 11164  ax-mulcom 11165  ax-addass 11166  ax-mulass 11167  ax-distr 11168  ax-i2m1 11169  ax-1ne0 11170  ax-1rid 11171  ax-rnegex 11172  ax-rrecex 11173  ax-cnre 11174  ax-pre-lttri 11175  ax-pre-lttrn 11176  ax-pre-ltadd 11177  ax-pre-mulgt0 11178

This theorem depends on definitions:  df-bi 210  df-an 401  df-or 861  df-3or 1104  df-3an 1105  df-tru 1573  df-fal 1583  df-ex 1810  df-nf 1814  df-sb 2097  df-mo 2567  df-eu 2597  df-clab 2742  df-cleq 2755  df-clel 2838  df-nfc 2912  df-ne 2959  df-nel 3065  df-ral 3080  df-rex 3090  df-reu 3370  df-rab 3417  df-v 3457  df-sbc 3746  df-csb 3855  df-dif 3909  df-un 3911  df-in 3913  df-ss 3923  df-pss 3926  df-nul 4288  df-if 4489  df-pw 4565  df-sn 4591  df-pr 4593  df-op 4597  df-uni 4874  df-iun 4959  df-br 5111  df-opab 5175  df-mpt 5194  df-tr 5220  df-id 5558  df-eprel 5563  df-po 5571  df-so 5572  df-fr 5616  df-we 5618  df-xp 5669  df-rel 5670  df-cnv 5671  df-co 5672  df-dm 5673  df-rn 5674  df-res 5675  df-ima 5676  df-pred 6304  df-ord 6365  df-on 6366  df-lim 6367  df-suc 6368  df-iota 6494  df-fun 6540  df-fn 6541  df-f 6542  df-f1 6543  df-fo 6544  df-f1o 6545  df-fv 6546  df-riota 7369  df-ov 7415  df-oprab 7416  df-mpo 7417  df-om 7864  df-2nd 7988  df-frecs 8279  df-wrecs 8310  df-recs 8359  df-rdg 8398  df-er 8695  df-en 8945  df-dom 8946  df-sdom 8947  df-pnf 11246  df-mnf 11247  df-xr 11248  df-ltxr 11249  df-le 11250  df-sub 11444  df-neg 11445  df-nn 12235  df-2 12304  df-3 12305  df-4 12306  df-5 12307  df-6 12308  df-7 12309  df-8 12310  df-9 12311  df-n0 12506  df-z 12593  df-dec 12713  df-uz 12864  df-seq 14040  df-exp 14100

This theorem is referenced by:  dpmul4  33211