使用cuFFT实现大整数乘法

序言

---

在某些场合，我们可能需要使用远超内置整型范围的整数进行运算，比如公钥加密等。如果使用最原始的竖式计算，那么时间复杂度是$T(n^2)$，其中n是相乘的两个整数的位数。使用Karatsuba算法优化，时间复杂度可以降至$T(n^{\log_23})\approx T(n^{1.585})$。而如果使用快速傅里叶变换（FFT），则可以优化到$T(n\log n)$。

原理

---

整数乘法与多项式

一个整数，我们可以把它表示成一个多项式。比如长度为n的十进制整数$a_0a_1a_2a_3...a_{n-1}$写成多项式：

$$A(x)=a_0x^{n-1}+a_1x^{n-2}+...+a_{n-2}x+a_{n-1}$$
整数2345，写成多项式就是：
$$A(x)=2x^3+3x^2+4x+5$$
这样一来，整数乘法就变成了多项式的乘法。比如2345*123，就变成了：

×+=2x5+2x5+4x4+