Theorem karatsuba 13211

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 9843. (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 9577	. . . . 5 ⊢ 𝐴 ∈ ℂ
3		10nn0 9796	. . . . . . 7 ⊢ ;10 ∈ ℕ0
4	3	nn0cni 9577	. . . . . 6 ⊢ ;10 ∈ ℂ
5		karatsuba.m	. . . . . 6 ⊢ 𝑀 ∈ ℕ0
6		expcl 10996	. . . . . 6 ⊢ ((;10 ∈ ℂ ∧ 𝑀 ∈ ℕ0) → (;10↑𝑀) ∈ ℂ)
7	4, 5, 6	mp2an 430	. . . . 5 ⊢ (;10↑𝑀) ∈ ℂ
8	2, 7	mulcli 8331	. . . 4 ⊢ (𝐴 · (;10↑𝑀)) ∈ ℂ
9		karatsuba.b	. . . . 5 ⊢ 𝐵 ∈ ℕ0
10	9	nn0cni 9577	. . . 4 ⊢ 𝐵 ∈ ℂ
11		karatsuba.c	. . . . . 6 ⊢ 𝐶 ∈ ℕ0
12	11	nn0cni 9577	. . . . 5 ⊢ 𝐶 ∈ ℂ
13	12, 7	mulcli 8331	. . . 4 ⊢ (𝐶 · (;10↑𝑀)) ∈ ℂ
14		karatsuba.d	. . . . 5 ⊢ 𝐷 ∈ ℕ0
15	14	nn0cni 9577	. . . 4 ⊢ 𝐷 ∈ ℂ
16	8, 10, 13, 15	muladdi 8736	. . 3 ⊢ (((𝐴 · (;10↑𝑀)) + 𝐵) · ((𝐶 · (;10↑𝑀)) + 𝐷)) = ((((𝐴 · (;10↑𝑀)) · (𝐶 · (;10↑𝑀))) + (𝐷 · 𝐵)) + (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵)))
17	8, 13	mulcli 8331	. . . 4 ⊢ ((𝐴 · (;10↑𝑀)) · (𝐶 · (;10↑𝑀))) ∈ ℂ
18	15, 10	mulcli 8331	. . . 4 ⊢ (𝐷 · 𝐵) ∈ ℂ
19	8, 15	mulcli 8331	. . . . 5 ⊢ ((𝐴 · (;10↑𝑀)) · 𝐷) ∈ ℂ
20	13, 10	mulcli 8331	. . . . 5 ⊢ ((𝐶 · (;10↑𝑀)) · 𝐵) ∈ ℂ
21	19, 20	addcli 8330	. . . 4 ⊢ (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵)) ∈ ℂ
22	17, 18, 21	add32i 8490	. . 3 ⊢ ((((𝐴 · (;10↑𝑀)) · (𝐶 · (;10↑𝑀))) + (𝐷 · 𝐵)) + (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵))) = ((((𝐴 · (;10↑𝑀)) · (𝐶 · (;10↑𝑀))) + (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵))) + (𝐷 · 𝐵))
23	8, 12	mulcli 8331	. . . . . 6 ⊢ ((𝐴 · (;10↑𝑀)) · 𝐶) ∈ ℂ
24		karatsuba.s	. . . . . . 7 ⊢ 𝑆 ∈ ℕ0
25	24	nn0cni 9577	. . . . . 6 ⊢ 𝑆 ∈ ℂ
26	23, 25, 7	adddiri 8337	. . . . 5 ⊢ ((((𝐴 · (;10↑𝑀)) · 𝐶) + 𝑆) · (;10↑𝑀)) = ((((𝐴 · (;10↑𝑀)) · 𝐶) · (;10↑𝑀)) + (𝑆 · (;10↑𝑀)))
27	2, 7, 12	mul32i 8473	. . . . . . . . 9 ⊢ ((𝐴 · (;10↑𝑀)) · 𝐶) = ((𝐴 · 𝐶) · (;10↑𝑀))
28		karatsuba.r	. . . . . . . . . 10 ⊢ (𝐴 · 𝐶) = 𝑅
29	28	oveq1i 6095	. . . . . . . . 9 ⊢ ((𝐴 · 𝐶) · (;10↑𝑀)) = (𝑅 · (;10↑𝑀))
30	27, 29	eqtri 2259	. . . . . . . 8 ⊢ ((𝐴 · (;10↑𝑀)) · 𝐶) = (𝑅 · (;10↑𝑀))
31	30	oveq1i 6095	. . . . . . 7 ⊢ (((𝐴 · (;10↑𝑀)) · 𝐶) + 𝑆) = ((𝑅 · (;10↑𝑀)) + 𝑆)
32		karatsuba.w	. . . . . . 7 ⊢ ((𝑅 · (;10↑𝑀)) + 𝑆) = 𝑊
33	31, 32	eqtri 2259	. . . . . 6 ⊢ (((𝐴 · (;10↑𝑀)) · 𝐶) + 𝑆) = 𝑊
34	33	oveq1i 6095	. . . . 5 ⊢ ((((𝐴 · (;10↑𝑀)) · 𝐶) + 𝑆) · (;10↑𝑀)) = (𝑊 · (;10↑𝑀))
35	8, 12, 7	mulassi 8335	. . . . . 6 ⊢ (((𝐴 · (;10↑𝑀)) · 𝐶) · (;10↑𝑀)) = ((𝐴 · (;10↑𝑀)) · (𝐶 · (;10↑𝑀)))
36	2, 12	mulcli 8331	. . . . . . . . . . . 12 ⊢ (𝐴 · 𝐶) ∈ ℂ
37	36, 18, 25	add32i 8490	. . . . . . . . . . 11 ⊢ (((𝐴 · 𝐶) + (𝐷 · 𝐵)) + 𝑆) = (((𝐴 · 𝐶) + 𝑆) + (𝐷 · 𝐵))
38	28	oveq1i 6095	. . . . . . . . . . . 12 ⊢ ((𝐴 · 𝐶) + 𝑆) = (𝑅 + 𝑆)
39		karatsuba.t	. . . . . . . . . . . . 13 ⊢ (𝐵 · 𝐷) = 𝑇
40	10, 15, 39	mulcomli 8333	. . . . . . . . . . . 12 ⊢ (𝐷 · 𝐵) = 𝑇
41	38, 40	oveq12i 6097	. . . . . . . . . . 11 ⊢ (((𝐴 · 𝐶) + 𝑆) + (𝐷 · 𝐵)) = ((𝑅 + 𝑆) + 𝑇)
42	37, 41	eqtri 2259	. . . . . . . . . 10 ⊢ (((𝐴 · 𝐶) + (𝐷 · 𝐵)) + 𝑆) = ((𝑅 + 𝑆) + 𝑇)
43		karatsuba.e	. . . . . . . . . 10 ⊢ ((𝐴 + 𝐵) · (𝐶 + 𝐷)) = ((𝑅 + 𝑆) + 𝑇)
44	2, 10, 12, 15	muladdi 8736	. . . . . . . . . 10 ⊢ ((𝐴 + 𝐵) · (𝐶 + 𝐷)) = (((𝐴 · 𝐶) + (𝐷 · 𝐵)) + ((𝐴 · 𝐷) + (𝐶 · 𝐵)))
45	42, 43, 44	3eqtr2i 2265	. . . . . . . . 9 ⊢ (((𝐴 · 𝐶) + (𝐷 · 𝐵)) + 𝑆) = (((𝐴 · 𝐶) + (𝐷 · 𝐵)) + ((𝐴 · 𝐷) + (𝐶 · 𝐵)))
46	36, 18	addcli 8330	. . . . . . . . . 10 ⊢ ((𝐴 · 𝐶) + (𝐷 · 𝐵)) ∈ ℂ
47	2, 15	mulcli 8331	. . . . . . . . . . 11 ⊢ (𝐴 · 𝐷) ∈ ℂ
48	12, 10	mulcli 8331	. . . . . . . . . . 11 ⊢ (𝐶 · 𝐵) ∈ ℂ
49	47, 48	addcli 8330	. . . . . . . . . 10 ⊢ ((𝐴 · 𝐷) + (𝐶 · 𝐵)) ∈ ℂ
50	46, 25, 49	addcani 8508	. . . . . . . . 9 ⊢ ((((𝐴 · 𝐶) + (𝐷 · 𝐵)) + 𝑆) = (((𝐴 · 𝐶) + (𝐷 · 𝐵)) + ((𝐴 · 𝐷) + (𝐶 · 𝐵))) ↔ 𝑆 = ((𝐴 · 𝐷) + (𝐶 · 𝐵)))
51	45, 50	mpbi 145	. . . . . . . 8 ⊢ 𝑆 = ((𝐴 · 𝐷) + (𝐶 · 𝐵))
52	51	oveq1i 6095	. . . . . . 7 ⊢ (𝑆 · (;10↑𝑀)) = (((𝐴 · 𝐷) + (𝐶 · 𝐵)) · (;10↑𝑀))
53	47, 48, 7	adddiri 8337	. . . . . . 7 ⊢ (((𝐴 · 𝐷) + (𝐶 · 𝐵)) · (;10↑𝑀)) = (((𝐴 · 𝐷) · (;10↑𝑀)) + ((𝐶 · 𝐵) · (;10↑𝑀)))
54	2, 15, 7	mul32i 8473	. . . . . . . 8 ⊢ ((𝐴 · 𝐷) · (;10↑𝑀)) = ((𝐴 · (;10↑𝑀)) · 𝐷)
55	12, 10, 7	mul32i 8473	. . . . . . . 8 ⊢ ((𝐶 · 𝐵) · (;10↑𝑀)) = ((𝐶 · (;10↑𝑀)) · 𝐵)
56	54, 55	oveq12i 6097	. . . . . . 7 ⊢ (((𝐴 · 𝐷) · (;10↑𝑀)) + ((𝐶 · 𝐵) · (;10↑𝑀))) = (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵))
57	52, 53, 56	3eqtri 2263	. . . . . 6 ⊢ (𝑆 · (;10↑𝑀)) = (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵))
58	35, 57	oveq12i 6097	. . . . 5 ⊢ ((((𝐴 · (;10↑𝑀)) · 𝐶) · (;10↑𝑀)) + (𝑆 · (;10↑𝑀))) = (((𝐴 · (;10↑𝑀)) · (𝐶 · (;10↑𝑀))) + (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵)))
59	26, 34, 58	3eqtr3ri 2268	. . . 4 ⊢ (((𝐴 · (;10↑𝑀)) · (𝐶 · (;10↑𝑀))) + (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵))) = (𝑊 · (;10↑𝑀))
60	59, 40	oveq12i 6097	. . 3 ⊢ ((((𝐴 · (;10↑𝑀)) · (𝐶 · (;10↑𝑀))) + (((𝐴 · (;10↑𝑀)) · 𝐷) + ((𝐶 · (;10↑𝑀)) · 𝐵))) + (𝐷 · 𝐵)) = ((𝑊 · (;10↑𝑀)) + 𝑇)
61	16, 22, 60	3eqtri 2263	. 2 ⊢ (((𝐴 · (;10↑𝑀)) + 𝐵) · ((𝐶 · (;10↑𝑀)) + 𝐷)) = ((𝑊 · (;10↑𝑀)) + 𝑇)
62		karatsuba.x	. . 3 ⊢ ((𝐴 · (;10↑𝑀)) + 𝐵) = 𝑋
63		karatsuba.y	. . 3 ⊢ ((𝐶 · (;10↑𝑀)) + 𝐷) = 𝑌
64	62, 63	oveq12i 6097	. 2 ⊢ (((𝐴 · (;10↑𝑀)) + 𝐵) · ((𝐶 · (;10↑𝑀)) + 𝐷)) = (𝑋 · 𝑌)
65		karatsuba.z	. 2 ⊢ ((𝑊 · (;10↑𝑀)) + 𝑇) = 𝑍
66	61, 64, 65	3eqtr3i 2267	1 ⊢ (𝑋 · 𝑌) = 𝑍

This proof depends on syntax axioms:   = wceq 1402   ∈ wcel 2209  (class class class)co 6085  ℂcc 8177  0cc0 8179  1c1 8180   + caddc 8182   · cmul 8184  ℕ₀cn0 9565  ;cdc 9779  ↑cexp 10977

This proof depends on axioms:  ax-mp 5  ax-1 6  ax-2 7  ax-ia1 106  ax-ia2 107  ax-ia3 108  ax-in1 623  ax-in2 624  ax-io 721  ax-5 1500  ax-7 1501  ax-gen 1502  ax-ie1 1546  ax-ie2 1547  ax-8 1557  ax-10 1558  ax-11 1559  ax-i12 1560  ax-bndl 1562  ax-4 1563  ax-17 1579  ax-i9 1583  ax-ial 1587  ax-i5r 1588  ax-14 2212  ax-ext 2220  ax-coll 4246  ax-sep 4249  ax-nul 4259  ax-pow 4311  ax-pr 4346  ax-un 4578  ax-setind 4684  ax-iinf 4735  ax-cnex 8270  ax-resscn 8271  ax-1cn 8272  ax-1re 8273  ax-icn 8274  ax-addcl 8275  ax-addrcl 8276  ax-mulcl 8277  ax-mulrcl 8278  ax-addcom 8279  ax-mulcom 8280  ax-addass 8281  ax-mulass 8282  ax-distr 8283  ax-i2m1 8284  ax-0lt1 8285  ax-1rid 8286  ax-0id 8287  ax-rnegex 8288  ax-precex 8289  ax-cnre 8290  ax-pre-ltirr 8291  ax-pre-ltwlin 8292  ax-pre-lttrn 8293  ax-pre-apti 8294  ax-pre-ltadd 8295  ax-pre-mulgt0 8296  ax-pre-mulext 8297

This proof depends on definitions:  df-bi 117  df-dc 847  df-3or 1010  df-3an 1011  df-tru 1405  df-fal 1408  df-nf 1514  df-sb 1816  df-eu 2089  df-mo 2090  df-clab 2225  df-cleq 2231  df-clel 2234  df-nfc 2381  df-ne 2421  df-nel 2516  df-ral 2533  df-rex 2534  df-reu 2535  df-rmo 2536  df-rab 2537  df-v 2823  df-sbc 3052  df-csb 3148  df-dif 3222  df-un 3224  df-in 3226  df-ss 3233  df-nul 3521  df-if 3639  df-pw 3690  df-sn 3715  df-pr 3716  df-op 3718  df-uni 3936  df-int 3971  df-iun 4014  df-br 4131  df-opab 4193  df-mpt 4194  df-tr 4230  df-id 4438  df-po 4441  df-iso 4442  df-iord 4511  df-on 4513  df-ilim 4514  df-suc 4516  df-iom 4738  df-xp 4780  df-rel 4781  df-cnv 4782  df-co 4783  df-dm 4784  df-rn 4785  df-res 4786  df-ima 4787  df-iota 5337  df-fun 5379  df-fn 5380  df-f 5381  df-f1 5382  df-fo 5383  df-f1o 5384  df-fv 5385  df-riota 6038  df-ov 6088  df-oprab 6089  df-mpo 6090  df-1st 6374  df-2nd 6375  df-recs 6576  df-frec 6662  df-pnf 8362  df-mnf 8363  df-xr 8364  df-ltxr 8365  df-le 8366  df-sub 8499  df-neg 8500  df-reap 8904  df-ap 8911  df-div 9004  df-inn 9306  df-2 9364  df-3 9365  df-4 9366  df-5 9367  df-6 9368  df-7 9369  df-8 9370  df-9 9371  df-n0 9566  df-z 9647  df-dec 9780  df-uz 9924  df-seqfrec 10887  df-exp 10978

This theorem is used by: (None)