A-16.1 Background and Motivation

Given two $(n-1)$-degree polynomials:

$A(X)=\underset{i=0}{\overset{n-1}{\sum <!--\; FUNCTION\; APPLICATION\; -->}}{a}_i{X}^i$, ... $B(X)=\underset{i=0}{\overset{n-1}{\sum <!--\; FUNCTION\; APPLICATION\; -->}}{b}_i{X}^i$

, the polynomial multiplication $C(X)=A(X)\cdot B(X)$ is computed as follows:

$C(X)=\underset{i=0}{\overset{2n-1}{\sum <!--\; FUNCTION\; APPLICATION\; -->}}{c}_i{X}^i$, where $c_i=\underset{k=0}{\overset{i}{\sum <!--\; FUNCTION\; APPLICATION\; -->}}{a}_k{b}_{i-k}$

This operation of computing $\overrightarrow{c}=(c_0,c_1,\cdots ,c_{2n-1})$ is also called the convolution of $\overrightarrow{a}$ and $\overrightarrow{b}$ denoted as $\overrightarrow{c}=\overrightarrow{a}\otimes \overrightarrow{b}$. The time complexity of this operation (i.e., the total number of multiplications between two numbers) is $O(n^2)$.

Another way of multiplying two polynomials is based on point-value representation. The point-value representation of an $(n-1)$-degree (or lesser degree) polynomial $A(X)$ is a set of $n$ coordinates $\{(x_0,y_0),(x_1,y_1),\cdots (x_{n-1},y_{n-1})\}$, where each $x_i$ is a distinct $X$ coordinate (whereas each $y_i$ is not necessarily a distinct $Y$ coordinate). Given a point-value representation of an $(n-1)$-degree (or lesser degree) polynomial, we can use polynomial interpolation (§A-15) to derive the polynomial. Let’s denote the point-value representation of $(n-1)$-degree (or lesser degree) polynomial $A(X)$ and $B(X)$ as follows:

$A(X)$ : $((x_0,y_0^{⟨a⟩}),(x_1,y_1^{⟨a⟩}),\cdots (x_{n-1},y_{n-1}^{⟨a⟩}))$

$B(X)$ : $((x_0,y_0^{⟨b⟩}),(x_1,y_1^{⟨b⟩}),\cdots (x_{n-1},y_{n-1}^{⟨b⟩}))$

Then, the point-value representation of polynomial $C(X)=A(X)\cdot B(X)$ can be computed as a Hadamard product (Definition A-10.1 in §A-10.1) of the $y$ values of point-value representation of $A(X)$ and $B(X)$ as follows:

$C(X)$ : $((x_0,y_0^{⟨c⟩}),(x_1,y_1^{⟨c⟩}),\cdots (x_{n-1},y_{n-1}^{⟨c⟩}))$, where $y_i^{⟨c⟩}=y_i^{⟨a⟩}\cdot y_i^{⟨b⟩}$

However, we cannot derive polynomial $C(X)$ based on these $n$ coordinates, because the degree of $C(X)$ is $2n-2$ (or lesser than $2n-2$). But if we regard all polynomials (including $A(X),B(X)$, and $C(X)$) to be in the polynomial ring $ℝ[X]∕X^n+1$ (or $\mathbb{Z}{[X]}_p∕X^n+1$), then we can reduce the $(2n-2)$-degree polynomial $C(X)$ to a congruent $(n-1)$-degree (or lesser degree) polynomial in the ring. Then, the $n-1$ coordinates of $C(X)$ are sufficient to derive $C(X)$.

However, the time complexity of this new method is still $O(n^2)$. The Hadamard product between two polynomials’ point-value representation takes $O(n)$, but evaluating a polynomial at $n$ distinct $x$ values takes $O(n^2)$ (because each polynomial has $n$ terms and we have to compute each term for $n$ distinct $x$ values), and the polynomial interpolation of deriving $C(X)$ also takes $O(n^2)$.

To solve this efficiency problem, this section will introduce an efficient technique of polynomial evaluation which can evaluate a polynomial at $n$ distinct roots of unity in $O(n\log <!--\; FUNCTION\; APPLICATION\; -->n)$. This technique is classified into 2 types: Fast Fourier Transform (FFT) and Number-theoretic Transform(NTT). These two types are technically almost the same with the only difference that FFT assumes a polynomial ring over complex numbers (§A-7), whereas NTT assumes a polynomial ring over the integer ring (§A-9). Polynomial multiplication based on FFT (or NTT) comprises 3 steps: (1) forward FFT (or NTT); (2) point-value multiplication; and (3) inverse FFT (or NTT).