秀爾演算法（英語：Shor's algorithm）是一個于1994年發現的，以數學家彼得·秀爾命名，針對整數分解題目的的量子演算法（在量子計算機上面運作的演算法）。不正式地說，它解決的題目是：給定一個整數 ${\displaystyle N}$，找出它的質因數。在一個量子計算機上面，要分解整數${\displaystyle N}$，秀爾演算法的運作需要多項式時間（時間是 ${\displaystyle \log N}$的某個多項式這麼長，${\displaystyle \log N}$ 在這裡的意義是輸入的檔案長度）。准确来说，该演算法花費 ${\displaystyle O((\log N)^{3})}$ 的時間，展示出質因數分解問題可以使用量子計算機以多項式時間解出，因此在複雜度類BQP裡面。這比傳統上已知的最快的因數分解演算法普通數域篩選法所花費的次指數時間——大約 ${\displaystyle O(e^{1.9(\log N)^{\frac {1}{3}}(\log \log N)^{\frac {2}{3}}})}$還要快了一個指數。

秀爾演算法非常重要，它意味着：假如有可用的量子計算機，我們可以破解基于公開密鑰加密的算法（比如RSA加密演算法）。RSA演算法的基礎在於假設了我們不能有效率的分解一個已知的整數。就目前所知，這假設對傳統（即非量子）的電腦為真；沒有已知的傳統演算法可以在多項式時間內解決這個問題。但秀爾演算法展示了因數分解問題在量子計算機上可以很有效率的解決，所以一個足夠大的量子計算機可以破解RSA。這對於建立量子計算機和研究新的量子計算機演算法是一個强大的動力。

2001年，IBM的一個小組展示了秀爾演算法的實例，使用NMR實驗的量子計算機及7個量子位元，將15分解成3×5。^[1]不过，对于IBM的實驗是否是量子計算的真實展示，有一些疑慮出現，因為沒有发现纏結現象。^[2]IBM的實驗之後，也有其他的團隊以光學量子位元實驗秀爾演算法，並強調可觀察到其纏結現象。^[3][4]

試著解決的問題是：給定一個合成數${\displaystyle N}$，找到整數${\displaystyle p}$在${\displaystyle 1}$和${\displaystyle N}$之間且不包含${\displaystyle 1}$和${\displaystyle N}$，並且${\displaystyle N}$整除於${\displaystyle p}$。

秀爾演算法包含兩個部份：

1. 一個以傳統電腦運作的簡化演算法，將因數分解簡化成搜尋階問題。
2. 一個量子演算法，解決搜尋階問題。

1. 選擇任意數字${\displaystyle a}$ < ${\displaystyle N}$
2. 計算gcd（ ${\displaystyle a}$ , ${\displaystyle N}$ ）。這裡可以使用輾轉相除法來計算。
3. 若gcd( ${\displaystyle a}$ , ${\displaystyle N}$ ) ≠ 1，則我們有了一個 ${\displaystyle N}$ 非平凡的因數，因此這部份結束了。
4. 否則，利用下面的週期尋找子程序（Period-finding subroutine，下面會列出）來找出下面這個函數的週期 ${\displaystyle r}$ ：

   ${\displaystyle f(x)=a^{x}\ {\mbox{mod}}\ N}$,

   換句話說，找出${\displaystyle a}$在${\displaystyle \mathbb {Z} _{N}}$裡面的階${\displaystyle r}$，或者最小的正整數 ${\displaystyle r}$ 令 ${\displaystyle f(x+r)=f(x)}$。
5. 若 ${\displaystyle r}$ 是奇數，回到第一步。
6. 若 ${\displaystyle a}$ ^{${\displaystyle r}$ /2} ≡ -1 (mod ${\displaystyle N}$ ),回到第一步。
7. gcd(a ^{${\displaystyle r}$ /2}+1, ${\displaystyle N}$ )與gcd(a ^{${\displaystyle r}$ /2}-1, ${\displaystyle N}$ )至少有一個是 ${\displaystyle N}$ 非平凡的因數。分解完成。

這個演算法使用的量子線路是為了在給定一個固定常數 ${\displaystyle N}$ 以及一個任意常數 ${\displaystyle a}$ 之下，找出${\displaystyle f(x)=a^{x}\mod {N}}$ 所設定的。給定 ${\displaystyle N}$ ,找出 ${\displaystyle Q}$ = 2 ^{${\displaystyle q}$} 且合乎${\displaystyle N^{2}\leq Q<2N^{2}}$（這同時表示${\displaystyle Q/r>N}$）。輸入和輸出量子位元暫存器需要儲存從0到 ${\displaystyle Q}$ -1所有值的疊加態，因此分別需要 ${\displaystyle q}$ 個量子位元。這裡使用看起來比所需的數量還要更多一倍的量子位元，保證了即使週期 ${\displaystyle r}$ 的大小逼近 ${\displaystyle N}$ /2，也至少有 ${\displaystyle N}$ 個不同的 ${\displaystyle x}$ 會產生相同的 ${\displaystyle f(x)}$ 。

程序如下：

1. 將暫存器初始化成

   ${\displaystyle Q^{-1/2}\sum _{x=0}^{Q-1}\left|x\right\rangle \left|0\right\rangle }$

   ${\displaystyle x}$ 從0到 ${\displaystyle Q}$ − 1。所以這一個初始態是 ${\displaystyle Q}$ 個狀態的疊加。
2. 建立量子函式版本的 ${\displaystyle f}$ ( ${\displaystyle x}$ )，並且應用於上面的疊加態，得到

   ${\displaystyle Q^{-1/2}\sum _{x}\left|x\right\rangle \left|f(x)\right\rangle }$.

   這裡仍舊是 ${\displaystyle Q}$ 個狀態的疊加。
3. 對輸入暫存器進行量子傅立葉變換。這個變換（操作於二的冪次—— ${\displaystyle Q=2^{q}}$ 個疊加態上面） 使用一個 ${\displaystyle Q}$次單位根例如${\displaystyle \omega =e^{2\pi i/Q}}$將任意給定態${\displaystyle \left|x\right\rangle }$的振幅平均分佈在所有 ${\displaystyle Q}$ 個${\displaystyle \left|y\right\rangle }$態上。另一個方法是對於每個不同的 ${\displaystyle x}$ ：

   ${\displaystyle U_{QFT}\left|x\right\rangle =Q^{-1/2}\sum _{y}\omega ^{xy}\left|y\right\rangle }$。

   由此得到最終狀態：

   ${\displaystyle Q^{-1}\sum _{x}\sum _{y}\omega ^{xy}\left|y\right\rangle \left|f(x)\right\rangle }$.

   這是一個遠多過 ${\displaystyle Q}$ 個狀態的疊加態，但是遠低過 ${\displaystyle Q}$ ²個。雖然在和中有 ${\displaystyle Q}$ ²項，但只要 ${\displaystyle x_{0}}$ 和 ${\displaystyle x}$ 的值相同，態${\displaystyle \left|y\right\rangle \left|f(x_{0})\right\rangle }$就可被提出來。令

   ${\displaystyle \omega =e^{2\pi i/Q}}$為 ${\displaystyle Q^{th}}$的一個單位根，

   ${\displaystyle r}$ 為 ${\displaystyle f}$ 的週期，

   ${\displaystyle x}$ ₀為一個產生相同 ${\displaystyle f}$ ( ${\displaystyle x}$ )的 ${\displaystyle x}$ 的集裡面的最小元素（我們已經有 ${\displaystyle x}$ ₀ < ${\displaystyle r}$ ），以及

   b在0到${\displaystyle \lfloor (Q-x_{0}-1)/r\rfloor }$之間使得${\displaystyle x_{0}+rb<Q}$。那麼${\displaystyle \omega ^{ry}}$則是復平面的一個單位向量（${\displaystyle \omega }$是一個單位根， ${\displaystyle r}$ 和y是整數），而${\displaystyle Q^{-1}\left|y\right\rangle \left|f(x_{0})\right\rangle }$在最終狀態下的係數則為

   ${\displaystyle \sum _{x:\,f(x)=f(x_{0})}\omega ^{xy}=\sum _{b}\omega ^{(x_{0}+rb)y}=\omega ^{x_{0}y}\sum _{b}\omega ^{rby}}$。這一求和的每一項代表一個獲得相同結果的不同路徑，而量子干涉發生。在單位向量${\displaystyle \omega ^{ryb}}$幾乎與復平面指向同一方向（要求${\displaystyle \omega ^{ry}}$指向正實數軸）時，干涉將是相長的。

4. 進行測量。我們由輸入寄存器取得結果y，由輸出寄存器取得${\displaystyle f(x_{0})}$。而既然 ${\displaystyle f}$ 是週期，對某對y和${\displaystyle f(x_{0})}$進行測量的概率則由

   ${\displaystyle \left|Q^{-1}\sum _{x:\,f(x)=f(x_{0})}\omega ^{xy}\right|^{2}=Q^{-2}\left|\sum _{b}\omega ^{(x_{0}+rb)y}\right|^{2}}$

   給出。分析顯示這個概率越高，單位向量${\displaystyle \omega ^{ry}}$就越接近正實數軸，或者${\displaystyle {\frac {ry}{Q}}}$越接近一個整數。除非 ${\displaystyle r}$ 是2的乘方，否則它不會是${\displaystyle Q}$的因子。
5. 對${\displaystyle {\frac {y}{Q}}}$進行連分數展開來計算其近似值，並生成滿足下列兩個條件的${\displaystyle {\frac {c}{r'}}}$：

   ${\displaystyle A:r'<N}$

   ${\displaystyle B:|{\frac {y}{Q}}-{\frac {c}{r'}}|<{\frac {1}{2Q}}}$

   藉著滿足這一些條件， ${\displaystyle r}$ ′有很高的機率會是我們要找的週期 ${\displaystyle r}$ 。
6. 檢查 ${\displaystyle f}$ ( ${\displaystyle x}$ ) = ${\displaystyle f}$ ( ${\displaystyle x}$ + ${\displaystyle r}$ ′) ${\displaystyle \Leftrightarrow }$ ${\displaystyle a^{r}\equiv 1{\pmod {N}}}$。如果成功了，我們就完成了。
7. 否則，以接近y左右的數值作為 ${\displaystyle r}$ 的候選，或者說多取幾個 ${\displaystyle r'}$.如果任何候選成功了，我們就完成了。
8. 否則，回到第一步驟（也就是全部重新做一次）。

此演算法包含兩個部份。演算法的第一部份是將因數分解問題轉成尋找一個函式的週期，而且這部份可以用傳統方式實作。第二部份則是使用量子傅立葉變換來搜尋這個函式的週期，而且這一部份是量子加速這整個演算法的理由。

小於N且互質於N的整數組成一個有限大，且對乘法同餘N的群。在步驟3之後，我們有一個屬於這個群的整數a。既然這個群是有限的，a必有一個有限大的階r,也就是最小的正整數令

${\displaystyle a^{r}\equiv 1\ {\mbox{mod}}\ N.\,}$

因此可知，N是a ^r − 1的因數。先假設我們有能力獲得r，而且r是偶數。則

${\displaystyle a^{r}-1=(a^{r/2}-1)(a^{r/2}+1)\equiv 0\ {\mbox{mod}}\ N}$

${\displaystyle \Rightarrow N\ |(a^{r/2}-1)(a^{r/2}+1).\,}$

r是令a ^r ≡ 1最小的正整數，所以（a ^{r / 2} − 1）必定不能整除於N。若（a ^{r / 2} + 1）也不整除於N的話，則N必定與（a ^{r / 2} − 1）或者（a ^{r / 2} + 1）有一個非顯然的公因數。

證明：為了簡化，我們分別用u和v表示（a ^{r / 2} − 1）和（a ^{r / 2} + 1）。N | uv，故kN = uv（k是某個整數）。假設gcd(v, N) = 1：則必定存在某個整數m和n令mv + nN = 1 (這是最大公因數的一個性質)。兩邊同乘以u，我們得到mkN + nuN = u, 所以 N | u，从而產生矛盾（前文提到（a ^{r / 2} − 1）必定不能整除於N）。故假设不成立，即gcd(v, N) ≠ 1。用類似的方法，可以得知若（a ^{r / 2} + 1）不能整除於N, gcd(u, N) ≠ 1。故得證u和v不同時與N互質。

這給予我們一個N的因數分解。若N是兩個質數的乘積，則這是唯一可能的分解。

秀爾的週期尋找演算法非常倚賴量子計算機同時處於許多狀態的能力。物理學家稱呼這個特性為狀態的「疊加」。在計算函數f的週期時，我們會同時估計這個函數的所有點。

不過，量子物理不允許我們直接取得這些資訊。測量會令觀測結果塌陷到其中一個可能的結果，並摧毀所有其他可能。如果不是因為不可複製原理，我們可以先測量f(x)而非測量x，再製造幾個結果態的拷貝（將會是多個有著同樣的f(x)的態的疊加）。對這些態上的x的測量將會給出能給出相同f(x)的不同的x，由此推導出週期來。因為我們不能複製完全相同的量子狀態，這個方法行不通。因此在這裡我們需要小心的轉變疊加態，使其成為可以以高機率回傳正確答案的狀態。這可使用量子傅立葉變換來達成。

因此秀爾在這裡必須解決三個「實作」的問題，而且速度必須夠快，也就是可這些問題可以用量子閘個數為${\displaystyle \log N}$的多項式來實作。

1. 製造狀態的疊加。 這可以藉著對所有輸入的量子位元使用哈達瑪閘（Hadamard gate）來達成。另一個方法是使用量子傅立葉變換（如下述）。
2. 以量子變換實作函數f。 要達成這件事情，秀爾使用了反覆平方以完成模的取幂變換。值得注意的是，這一個步驟比起量子傅立葉變換還難以實作，需要更多輔助的量子位元以及明顯更多的閘來完成。
3. 進行量子傅立葉變換。 利用受控旋轉閘（controlled rotation gate）和哈達瑪閘，秀爾設計了一個量子傅立葉變換的線路，只使用了${\displaystyle q(q-1)/2=O((\log Q)^{2})}$個閘（Q = 2^q）。^[5]

在這一些變換之後，測量會給出週期r的近似值。為簡明起見，假設存在y令yr/Q是整數，則測量到y的機率是1（理論上）。要推導出這個，我們首先注意到對任何整數b，

${\displaystyle e^{-2\pi ibyr/Q}=1}$。因此Q/r的平方能告訴我們測量y的概率，因為b大致上取Q/r個值，因此概率為${\displaystyle 1/r^{2}}$。有r個y使得yr/Q為整數，${\displaystyle f(x_{0})}$也有r種可能，所以機率總和為1。

Note:另一種解釋秀爾演算法的方式是將之視為是喬裝過後的量子相位估計演算法。

• 泡利矩陣
• 量子退火