Theorem karatsuba 17290

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 12911. (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 12645	. . . . 5 ⊢ 𝐴 ∈ ℂ
3		10nn0 12863	. . . . . . 7 ⊢ ;10 ∈ ℕ0
4	3	nn0cni 12645	. . . . . 6 ⊢ ;10 ∈ ℂ
5		karatsuba.m	. . . . . 6 ⊢ 𝑀 ∈ ℕ0
6		expcl 14249	. . . . . 6 ⊢ ((;10 ∈ ℂ ∧ 𝑀 ∈ ℕ0) → (;10↑𝑀) ∈ ℂ)
7	4, 5, 6	mp2an 705	. . . . 5 ⊢ (;10↑𝑀) ∈ ℂ
8	2, 7	mulcli 11343	. . . 4 ⊢ (𝐴 · (;10↑𝑀)) ∈ ℂ
9		karatsuba.b	. . . . 5 ⊢ 𝐵 ∈ ℕ0
10	9	nn0cni 12645	. . . 4 ⊢ 𝐵 ∈ ℂ
11		karatsuba.c	. . . . . 6 ⊢ 𝐶 ∈ ℕ0
12	11	nn0cni 12645	. . . . 5 ⊢ 𝐶 ∈ ℂ
13	12, 7	mulcli 11343	. . . 4 ⊢ (𝐶 · (;10↑𝑀)) ∈ ℂ
14		karatsuba.d	. . . . 5 ⊢ 𝐷 ∈ ℕ0
15	14	nn0cni 12645	. . . 4 ⊢ 𝐷 ∈ ℂ
16	8, 10, 13, 15	muladdi 11794	. . 3 ⊢ (((𝐴 · (;10↑𝑀)) + 𝐵) · ((𝐶 · (;10↑𝑀)) + 𝐷)) = ((((𝐴 · (;10↑𝑀)) · (𝐶 · (;10↑𝑀))) + (𝐷 · 𝐵)) + (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵)))
17	8, 13	mulcli 11343	. . . 4 ⊢ ((𝐴 · (;10↑𝑀)) · (𝐶 · (;10↑𝑀))) ∈ ℂ
18	15, 10	mulcli 11343	. . . 4 ⊢ (𝐷 · 𝐵) ∈ ℂ
19	8, 15	mulcli 11343	. . . . 5 ⊢ ((𝐴 · (;10↑𝑀)) · 𝐷) ∈ ℂ
20	13, 10	mulcli 11343	. . . . 5 ⊢ ((𝐶 · (;10↑𝑀)) · 𝐵) ∈ ℂ
21	19, 20	addcli 11342	. . . 4 ⊢ (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵)) ∈ ℂ
22	17, 18, 21	add32i 11561	. . 3 ⊢ ((((𝐴 · (;10↑𝑀)) · (𝐶 · (;10↑𝑀))) + (𝐷 · 𝐵)) + (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵))) = ((((𝐴 · (;10↑𝑀)) · (𝐶 · (;10↑𝑀))) + (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵))) + (𝐷 · 𝐵))
23	8, 12	mulcli 11343	. . . . . 6 ⊢ ((𝐴 · (;10↑𝑀)) · 𝐶) ∈ ℂ
24		karatsuba.s	. . . . . . 7 ⊢ 𝑆 ∈ ℕ0
25	24	nn0cni 12645	. . . . . 6 ⊢ 𝑆 ∈ ℂ
26	23, 25, 7	adddiri 11349	. . . . 5 ⊢ ((((𝐴 · (;10↑𝑀)) · 𝐶) + 𝑆) · (;10↑𝑀)) = ((((𝐴 · (;10↑𝑀)) · 𝐶) · (;10↑𝑀)) + (𝑆 · (;10↑𝑀)))
27	2, 7, 12	mul32i 11533	. . . . . . . . 9 ⊢ ((𝐴 · (;10↑𝑀)) · 𝐶) = ((𝐴 · 𝐶) · (;10↑𝑀))
28		karatsuba.r	. . . . . . . . . 10 ⊢ (𝐴 · 𝐶) = 𝑅
29	28	oveq1i 7453	. . . . . . . . 9 ⊢ ((𝐴 · 𝐶) · (;10↑𝑀)) = (𝑅 · (;10↑𝑀))
30	27, 29	eqtri 2786	. . . . . . . 8 ⊢ ((𝐴 · (;10↑𝑀)) · 𝐶) = (𝑅 · (;10↑𝑀))
31	30	oveq1i 7453	. . . . . . 7 ⊢ (((𝐴 · (;10↑𝑀)) · 𝐶) + 𝑆) = ((𝑅 · (;10↑𝑀)) + 𝑆)
32		karatsuba.w	. . . . . . 7 ⊢ ((𝑅 · (;10↑𝑀)) + 𝑆) = 𝑊
33	31, 32	eqtri 2786	. . . . . 6 ⊢ (((𝐴 · (;10↑𝑀)) · 𝐶) + 𝑆) = 𝑊
34	33	oveq1i 7453	. . . . 5 ⊢ ((((𝐴 · (;10↑𝑀)) · 𝐶) + 𝑆) · (;10↑𝑀)) = (𝑊 · (;10↑𝑀))
35	8, 12, 7	mulassi 11347	. . . . . 6 ⊢ (((𝐴 · (;10↑𝑀)) · 𝐶) · (;10↑𝑀)) = ((𝐴 · (;10↑𝑀)) · (𝐶 · (;10↑𝑀)))
36	2, 12	mulcli 11343	. . . . . . . . . . . 12 ⊢ (𝐴 · 𝐶) ∈ ℂ
37	36, 18, 25	add32i 11561	. . . . . . . . . . 11 ⊢ (((𝐴 · 𝐶) + (𝐷 · 𝐵)) + 𝑆) = (((𝐴 · 𝐶) + 𝑆) + (𝐷 · 𝐵))
38	28	oveq1i 7453	. . . . . . . . . . . 12 ⊢ ((𝐴 · 𝐶) + 𝑆) = (𝑅 + 𝑆)
39		karatsuba.t	. . . . . . . . . . . . 13 ⊢ (𝐵 · 𝐷) = 𝑇
40	10, 15, 39	mulcomli 11345	. . . . . . . . . . . 12 ⊢ (𝐷 · 𝐵) = 𝑇
41	38, 40	oveq12i 7455	. . . . . . . . . . 11 ⊢ (((𝐴 · 𝐶) + 𝑆) + (𝐷 · 𝐵)) = ((𝑅 + 𝑆) + 𝑇)
42	37, 41	eqtri 2786	. . . . . . . . . 10 ⊢ (((𝐴 · 𝐶) + (𝐷 · 𝐵)) + 𝑆) = ((𝑅 + 𝑆) + 𝑇)
43		karatsuba.e	. . . . . . . . . 10 ⊢ ((𝐴 + 𝐵) · (𝐶 + 𝐷)) = ((𝑅 + 𝑆) + 𝑇)
44	2, 10, 12, 15	muladdi 11794	. . . . . . . . . 10 ⊢ ((𝐴 + 𝐵) · (𝐶 + 𝐷)) = (((𝐴 · 𝐶) + (𝐷 · 𝐵)) + ((𝐴 · 𝐷) + (𝐶 · 𝐵)))
45	42, 43, 44	3eqtr2i 2792	. . . . . . . . 9 ⊢ (((𝐴 · 𝐶) + (𝐷 · 𝐵)) + 𝑆) = (((𝐴 · 𝐶) + (𝐷 · 𝐵)) + ((𝐴 · 𝐷) + (𝐶 · 𝐵)))
46	36, 18	addcli 11342	. . . . . . . . . 10 ⊢ ((𝐴 · 𝐶) + (𝐷 · 𝐵)) ∈ ℂ
47	2, 15	mulcli 11343	. . . . . . . . . . 11 ⊢ (𝐴 · 𝐷) ∈ ℂ
48	12, 10	mulcli 11343	. . . . . . . . . . 11 ⊢ (𝐶 · 𝐵) ∈ ℂ
49	47, 48	addcli 11342	. . . . . . . . . 10 ⊢ ((𝐴 · 𝐷) + (𝐶 · 𝐵)) ∈ ℂ
50	46, 25, 49	addcani 11530	. . . . . . . . 9 ⊢ ((((𝐴 · 𝐶) + (𝐷 · 𝐵)) + 𝑆) = (((𝐴 · 𝐶) + (𝐷 · 𝐵)) + ((𝐴 · 𝐷) + (𝐶 · 𝐵))) ↔ 𝑆 = ((𝐴 · 𝐷) + (𝐶 · 𝐵)))
51	45, 50	mpbi 233	. . . . . . . 8 ⊢ 𝑆 = ((𝐴 · 𝐷) + (𝐶 · 𝐵))
52	51	oveq1i 7453	. . . . . . 7 ⊢ (𝑆 · (;10↑𝑀)) = (((𝐴 · 𝐷) + (𝐶 · 𝐵)) · (;10↑𝑀))
53	47, 48, 7	adddiri 11349	. . . . . . 7 ⊢ (((𝐴 · 𝐷) + (𝐶 · 𝐵)) · (;10↑𝑀)) = (((𝐴 · 𝐷) · (;10↑𝑀)) + ((𝐶 · 𝐵) · (;10↑𝑀)))
54	2, 15, 7	mul32i 11533	. . . . . . . 8 ⊢ ((𝐴 · 𝐷) · (;10↑𝑀)) = ((𝐴 · (;10↑𝑀)) · 𝐷)
55	12, 10, 7	mul32i 11533	. . . . . . . 8 ⊢ ((𝐶 · 𝐵) · (;10↑𝑀)) = ((𝐶 · (;10↑𝑀)) · 𝐵)
56	54, 55	oveq12i 7455	. . . . . . 7 ⊢ (((𝐴 · 𝐷) · (;10↑𝑀)) + ((𝐶 · 𝐵) · (;10↑𝑀))) = (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵))
57	52, 53, 56	3eqtri 2790	. . . . . 6 ⊢ (𝑆 · (;10↑𝑀)) = (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵))
58	35, 57	oveq12i 7455	. . . . 5 ⊢ ((((𝐴 · (;10↑𝑀)) · 𝐶) · (;10↑𝑀)) + (𝑆 · (;10↑𝑀))) = (((𝐴 · (;10↑𝑀)) · (𝐶 · (;10↑𝑀))) + (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵)))
59	26, 34, 58	3eqtr3ri 2795	. . . 4 ⊢ (((𝐴 · (;10↑𝑀)) · (𝐶 · (;10↑𝑀))) + (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵))) = (𝑊 · (;10↑𝑀))
60	59, 40	oveq12i 7455	. . 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 7455	. 2 ⊢ (((𝐴 · (;10↑𝑀)) + 𝐵) · ((𝐶 · (;10↑𝑀)) + 𝐷)) = (𝑋 · 𝑌)
65		karatsuba.z	. 2 ⊢ ((𝑊 · (;10↑𝑀)) + 𝑇) = 𝑍
66	61, 64, 65	3eqtr3i 2794	1 ⊢ (𝑋 · 𝑌) = 𝑍

This proof depends on syntax axioms:   = wceq 1570   ∈ wcel 2147  (class class class)co 7443  ℂcc 11225  0cc0 11227  1c1 11228   + caddc 11230   · cmul 11232  ℕ₀cn0 12633  ;cdc 12841  ↑cexp 14231

This proof depends on axioms:  ax-mp 5  ax-1 6  ax-2 7  ax-3 8  ax-gen 1828  ax-4 1842  ax-5 1945  ax-6 2002  ax-7 2043  ax-8 2149  ax-9 2157  ax-10 2180  ax-11 2197  ax-12 2216  ax-ext 2735  ax-sep 5252  ax-nul 5263  ax-pow 5331  ax-pr 5395  ax-un 7775  ax-cnex 11283  ax-resscn 11284  ax-1cn 11285  ax-icn 11286  ax-addcl 11287  ax-addrcl 11288  ax-mulcl 11289  ax-mulrcl 11290  ax-mulcom 11291  ax-addass 11292  ax-mulass 11293  ax-distr 11294  ax-i2m1 11295  ax-1ne0 11296  ax-1rid 11297  ax-rnegex 11298  ax-rrecex 11299  ax-cnre 11300  ax-pre-lttri 11301  ax-pre-lttrn 11302  ax-pre-ltadd 11303  ax-pre-mulgt0 11304

This proof depends on definitions:  df-bi 210  df-an 402  df-or 862  df-3or 1104  df-3an 1105  df-tru 1573  df-fal 1583  df-ex 1813  df-nf 1817  df-sb 2102  df-mo 2567  df-eu 2597  df-clab 2742  df-cleq 2755  df-clel 2839  df-nfc 2913  df-ne 2960  df-nel 3066  df-ral 3081  df-rex 3091  df-reu 3370  df-rab 3417  df-v 3456  df-sbc 3743  df-csb 3852  df-ss 3901  df-pss 4017  df-dif 4050  df-in 4094  df-un 4154  df-nul 4282  df-if 4486  df-pw 4562  df-sn 4589  df-pr 4591  df-op 4839  df-uni 4871  df-iun 4956  df-br 5107  df-opab 5171  df-mpt 5190  df-tr 5216  df-id 5550  df-eprel 5555  df-po 5563  df-so 5564  df-fr 5608  df-we 5610  df-xp 5661  df-rel 5662  df-cnv 5663  df-co 5664  df-dm 5665  df-rn 5666  df-res 5667  df-ima 5668  df-pred 6322  df-ord 6383  df-on 6384  df-lim 6385  df-suc 6386  df-iota 6512  df-fun 6560  df-fn 6561  df-f 6562  df-f1 6563  df-fo 6564  df-f1o 6565  df-fv 6566  df-riota 7400  df-ov 7446  df-oprab 7447  df-mpo 7448  df-om 7902  df-2nd 8026  df-frecs 8327  df-wrecs 8358  df-recs 8407  df-rdg 8446  df-er 8745  df-en 9002  df-dom 9003  df-sdom 9004  df-pnf 11372  df-mnf 11373  df-xr 11374  df-ltxr 11375  df-le 11376  df-sub 11570  df-neg 11571  df-nn 12363  df-2 12432  df-3 12433  df-4 12434  df-5 12435  df-6 12436  df-7 12437  df-8 12438  df-9 12439  df-n0 12634  df-z 12721  df-dec 12842  df-uz 12993  df-seq 14172  df-exp 14232

This theorem is used by:  dpmul4  33550