階乘

在數學中，正整数的階乘（英語：Factorial）是所有小於等於該數的正整數的積，計為$n!$，例如5的階乘表示為${\displaystyle 5!}$，其值為120：

${\displaystyle 5!={{{{{5}\times {4}}\times {3}}\times {2}}\times {1}}=120}$

並定義，1的階乘${\displaystyle 1!}$和0的階乘${\displaystyle 0!}$都為1，其中0的階乘表示一個空積^[2]。

实数范围内的階乘函数，负整数除外^{[註 1]}

1808年，基斯頓·卡曼引進這個表示法：$n!=\prod _{k=1}^{n}k\quad \forall n\geq 1$，符號$\Pi$表示連續乘積，亦即${\displaystyle n!=1\times 2\times 3\times \cdots \times n}$。階乘亦可以遞迴方式定義：$0!=1$，${\displaystyle n!=(n-1)!\times n}$。除了自然數之外，階乘亦可定義于整個實數（負整數除外），其与伽瑪函數的关系为：

$z!=\Gamma (z+1)=\int _{0}^{\infty }t^{z}e^{-t}\,dt$

階乘應用在許多數學領域中，最常應用在組合數學、代數學和数学分析中。在組合數學中，階乘代表的意義為$n$個相異物件任意排列的數量，例如前述例子，${\displaystyle 5!=120}$其代表了5個相異物件共有120種排列法。在正整數的情形下，$n$的階乘又可以稱為n的排列數。

歷史

早在12世紀，印度學者就已有使用階乘的概念來計算排列數的紀錄^[3]。1677年時，法比安·斯特德曼使用Change ringing（英语：Change ringing）來解釋階乘的概念^[5]。在描述遞迴方法之後，斯特德將階乘描述為：「現在這些方法的本質是這樣的：一個數字的變化數包含了所有比他小的數字（包括本身）的所有變化數……因為一個數字的完全變化數是將較小數字的變化數視為一個整體，並透過將所有數字的完整變化聯合起來。」，其原文如下：

Now the nature of these methods is such, that the changes on one number comprehends [includes] the changes on all lesser numbers ... insomuch that a compleat Peal of changes on one number seemeth to be formed by uniting of the compleat Peals on all lesser numbers into one entire body.^[6]

而符號n!是由法國數學家克里斯蒂安·克蘭普在1808年使用^[8]。

定義

階乘可透過連乘積來定義：

${\displaystyle n!=1\cdot 2\cdot 3\cdots (n-2)\cdot (n-1)\cdot n,}$ 

用連乘積符號可表示為：

${\displaystyle n!=\prod _{i=1}^{n}i.\quad \forall n\geq 1}$ 

從上述公式中，可以推導出遞迴關係式：

${\displaystyle n!=n\cdot (n-1)!}$ 

但遞迴定義須給出base case，因此需要定義零的階乘。 除此之外，遞迴關係在階乘函數中各個值皆成立，例如：

${\displaystyle {\begin{aligned}5!&=5\cdot 4!\\6!&=6\cdot 5!\\50!&=50\cdot 49!\end{aligned}}}$ 

0的階乘

為了將遞迴關係式擴展到$n=0$ ，因此需要定義0的階乘：

${\displaystyle 0!=1}$ 

可以得到：

${\displaystyle 1!=1\cdot 0!=1}$ 

有幾個獨立的理由認為這個定義是和諧的。 其中包括：

• 在$n=0$ 的情況，$n!$ 定義為「沒有任何數字相乘的結果」，所以更廣泛之慣例的例子是以不存在任何因數的乘法單位元素來當作其解。（參閱空積）
• 對於零個物品只有一種排列方式，因為沒有任何東西可以置換，唯一的重新排列就是什麼都不做。
• 它使組合數學中的許多恆等式對所有適用的值皆有效，例如從空集合中選擇0個元素的方法數，可由二項式係數給出：

${\displaystyle {\binom {0}{0}}=1}$ .

而從空集合中選擇0個元素的方法數為一種，即沒有任何東西可以取，唯一的取法就是什麼都不做。定義${\displaystyle 0!=1}$ 可以滿足：

${\displaystyle {\binom {0}{0}}={\frac {0!}{0!0!}}=1}$ .

更一般地，在$n$ 個相異元素的集合中取出$n$ 個相異元素的方法數，可由二項式係數給出：

${\displaystyle {\binom {n}{n}}=1}$ .

其方法數只有一種，即全部取出。定義${\displaystyle 0!=1}$ 可以滿足：

${\displaystyle {\binom {n}{n}}={\frac {n!}{n!0!}}=1}$ 

• 此定義允許將許多公式更嚴謹地表達為冪級數，例如指數函數：

${\displaystyle e^{x}=\sum _{n=0}^{\infty }{\frac {x^{n}}{n!}}.}$ 

性質

$n!$ 可质因子分解为$\prod _{p\leq n}p^{\sum _{r=1}^{n}[{\frac {n}{p^{r}}}]}$ ，如${\displaystyle 6!=2^{4}\times 3^{2}\times 5^{1}}$ 。^[9]

計算

階乘與斯特靈公式

$n!$ （藍色）、${\displaystyle {\sqrt {2\pi n}}\left({\frac {n}{e}}\right)^{n}}$ （橘色），數字越大${\displaystyle {\sqrt {2\pi n}}\left({\frac {n}{e}}\right)^{n},}$ 會越趨近$n!$ 。但${\displaystyle {\sqrt {2\pi n}}\left({\frac {n}{e}}\right)^{n}}$ 在負值則會因為出現虛數而無法使用。

計算${\displaystyle n!}$ 時，若${\displaystyle n}$ 不太大，普通的科學計算機都可以計算，能夠處理不超過${\displaystyle 10^{100}}$ （古高爾）數值的計算機可以計算至${\displaystyle 69!}$ ，而雙精度浮點數的計算機則可計算至${\displaystyle 170!}$ 。

當${\displaystyle n}$ 很大時，可用斯特林公式估計： ${\displaystyle n!\approx {\sqrt {2\pi n}}\;\left({\frac {n}{e}}\right)^{n}}$ 
更精确的估计是： ${\displaystyle n!={\sqrt {2\pi n}}\;\left({\frac {n}{e}}\right)^{n}e^{\lambda _{n}}}$ 
其中 ${\displaystyle {\frac {1}{12n+1}}<\lambda _{n}<{\frac {1}{12n}}}$ 

部分函數值

部分的階乘值 （OEIS數列A000142）

n	n!
0	1
1	1
2	2
3	6
4	24
5	120
6	720
7	7003504000000000000♠5040
8	7004403200000000000♠40320
9	7005362880000000000♠362880
10	7006362880000000000♠3628800
11	7007399168000000000♠39916800
12	7008479001600000000♠479001600
13	7009622702080000000♠6227020800
14	7010871782912000000♠87178291200
15	7012130767436800000♠1307674368000
16	7013209227898880000♠20922789888000
17	7014355687428096000♠355687428096000
18	7015640237370572800♠6402373705728000
19	7017121645100408832♠121645100408832000
20	7018243290200817664♠2432902008176640000
25	7025155112100433310♠1.5511210043331×10^25
50	7064304140932017130♠3.0414093201713×10^64
70	7100119785716699700♠1.197857166997×10^100
100	7157933262154400000♠9.332621544×10^157
450	9000000000000000000♠1.733368733×10^1000
7003100000000000000♠1000	9000000000000000000♠4.023872601×10^2567
7003324900000000000♠3249	9000000000000000000♠6.412337688×10^10000
7004100000000000000♠10000	9000000000000000000♠2.846259681×10^35659
7004252060000000000♠25206	9000000000000000000♠1.205703438×10^100000
7005100000000000000♠100000	9000000000000000000♠2.824229408×10^456573
7005205023000000000♠205023	9000000000000000000♠2.503898932×10^1000004
7006100000000000000♠1000000	9000000000000000000♠8.263931688×10^5565708
7100100000000000000♠10^100	10^10101.9981097754820

非正整數的階乘

階乘原始的定義是在整數，為離散，然而在部分領域如機率論要探討到連續或其他需求（如組合數當取出的數量大於原有的數量會出現負階乘）時，則需要將階乘從正整數推廣到實數，甚至是複數。

${\displaystyle \Gamma }$ 函数和${\displaystyle \Pi }$ 函数

 

伽馬函數將階乘函數為非整數插值。主要線索是階乘函數的遞歸關係在連續的伽馬函數中也存在。

除了非負整數之外，還可以為非整數值定義階乘函數，但這需要使用更高級的數值分析方法。

可以透過插值的方式將階乘兩整數之間填入數值，但其插入的數值必須也要滿足階乘的遞迴定義。一個良好的插值結果是${\displaystyle \Gamma }$ 函数，其為所有非負整數和複數給出了定義，而當${\displaystyle z}$ 的實部為正時，可以透過下列瑕積分來計算${\displaystyle \Gamma }$ 函数值：

${\displaystyle \Gamma (z)=\int _{0}^{\infty }t^{z-1}e^{-t}\,dt\,.}$ 

它與階乘的關係是對於任何自然數n滿足：

${\displaystyle n!=\Gamma (n+1)\,.}$ 

另外，我们也可利用此式以计算任意大于-1的实数的阶乘：

${\displaystyle x!=\lim _{N\to \infty }N^{x}\prod _{k=1}^{N}{\frac {k}{x+k}}=\int _{0}^{1}(-\ln {(t)})^{x}\,dt\,.}$ 

複數的階乘

 

複數階乘之模與輻角的等值線

可以透過${\displaystyle \Gamma }$ 函數來計算複數的階乘。右圖顯示了複數階乘之模與輻角的等值線

令${\displaystyle f}$ 為：

${\displaystyle f=\rho e^{i\varphi }=(x+{\rm {i}}y)!=\Gamma (x+iy+1)}$ 

右圖顯示了幾個模（絕對值）${\displaystyle \rho }$ 與輻角${\displaystyle \varphi }$ 的幾個等級，圖表的繪製範圍為${\displaystyle -3\leq x\leq 3}$ , ${\displaystyle -2\leq y\leq 2}$ 個單位長。較粗的鉛直線為輻角值為${\displaystyle \varphi =\pm \pi }$ 的等值線。

細線表示模或輻角相等之函數值的位置。在每個負整數的位置為奇點，無法定義其模和輻角，並且在離奇點越近的地方，等值線的密度就越密集。

在|z| < 1時，可使用泰勒級數來計算：

${\displaystyle z!=\sum _{n=0}^{\infty }g_{n}z^{n}\,.}$ 

其泰勒級數的前幾項係數為：

n	gn	近似值
0	1	1
1	${\displaystyle -\gamma }$	3000422784335100000♠−0.5772156649
2	${\displaystyle {\frac {\pi ^{2}}{12}}+{\frac {\gamma ^{2}}{2}}}$	6999989055995500000♠0.9890559955
3	${\displaystyle -{\frac {\zeta (3)}{3}}-{\frac {\pi ^{2}}{12}}-{\frac {\gamma ^{3}}{6}}}$	3000092520924000000♠−0.9074790760

其中，γ為歐拉-馬斯刻若尼常數

ζ(z)為黎曼ζ函數

部分計算機代數的系統存在可以直接產生這些展開式係數的語法，例如SageMath。

此種方式甚至可以將階乘推廣至四元數甚至其他數學結構。

z	z!
實數
1、2、3、4、5	1、2、6、24、120 （OEIS數列A000142）
${\displaystyle {\frac {1}{2}}}$	${\displaystyle {\frac {\sqrt {\pi }}{2}}\approx \,}$ ${\displaystyle 0.88622692545276}$  （OEIS數列A019704）
複數
${\displaystyle i}$	${\displaystyle 0.49801566811836-0.15494982830181i}$  （OEIS數列A212877）、（OEIS數列A212878）
${\displaystyle 2i}$	${\displaystyle 0.15190400267024+0.019804880162337i}$
${\displaystyle 1+i}$	${\displaystyle 0.65296549642017+0.34306583981655i}$
四元數
${\displaystyle j}$	${\displaystyle 0.49801566811836-0.15494982830181j}$
${\displaystyle k}$	${\displaystyle 0.49801566811836-0.15494982830181k}$
${\displaystyle 1+i+j}$	${\displaystyle 0.31694069797431-0.045151191260681i-0.045151191260681j}$

 

階乘的色相環複變函數圖形。顏色越深代表絕對值越接近零；顏色越接近白色代表絕對值趨於無窮。其中紅色為正實數、青藍色為負實數。

較大的階乘值可透過双伽玛函数積分的連續分數來近似，這個方法由T. J. Stieltjes於1894提出。

將階乘寫為${\displaystyle z!=e^{P(z)}}$ ，其中${\displaystyle P(z)}$ 為：

${\displaystyle P(z)=p(z)+{\frac {\ln 2\pi }{2}}-z+\left(z+{\tfrac {1}{2}}\right)\ln(z)\,,}$ 

Stieltjes給出了其連分數值：

${\displaystyle p(z)={\cfrac {a_{0}}{z+{\cfrac {a_{1}}{z+{\cfrac {a_{2}}{z+{\cfrac {a_{3}}{z+\ddots }}}}}}}}}$ 

前幾項係數${\displaystyle a_{n}}$ 為^[10]：

n	an
0	${\displaystyle {\frac {1}{12}}}$
1	${\displaystyle {\frac {1}{30}}}$
2	${\displaystyle {\frac {53}{210}}}$
3	${\displaystyle {\frac {195}{371}}}$
4	${\displaystyle {\frac {22999}{22737}}}$
5	${\displaystyle {\frac {29944523}{19733142}}}$
6	${\displaystyle {\frac {109535241009}{48264275462}}}$

負整數的階乘

負整數的階乘可透過階乘的遞迴定義${\displaystyle n!=n\times (n-1)!}$ 逆推而得：

${\displaystyle (n-1)!={\frac {n!}{n}}.}$ 

但由於在此定義下計算負一的階乘會出現除以零（即${\displaystyle (0-1)!={\frac {0!}{0}}}$ ），因此無法直接給出負整數的階乘。

其他數學結構的階乘

透過伽瑪函數或其展開式亦可以將階乘擴展到其他能定義加法和乘法等基本運算的數學結構，如矩陣^[11]。

矩陣的階乘具有如下性質：

${\displaystyle A!=\Gamma (A+I)=A\Gamma (A)=A(A-I)!}$ 。

並且${\displaystyle \Gamma (I)=I}$ ，其中，${\displaystyle I}$ 是單位矩陣、${\displaystyle A}$ 是一個方陣，同時${\displaystyle A!}$ 是一個非奇異矩陣^[12]。

換句話說，即矩陣${\displaystyle A}$ 為單位矩陣的純量${\displaystyle n}$ 倍，其階乘為${\displaystyle A!=(nI)!=n!I}$ ，例如${\displaystyle {\bigl (}{\begin{smallmatrix}n&0\\0&n\end{smallmatrix}}{\bigr )}!=n!I={\bigl (}{\begin{smallmatrix}n!&0\\0&n!\end{smallmatrix}}{\bigr )}}$ 

對於一個可對角化矩陣${\displaystyle {\bigl (}{\begin{smallmatrix}a&b\\c&d\end{smallmatrix}}{\bigr )}}$ 其階乘為：

${\displaystyle \left.{\begin{pmatrix}a&b\\c&d\end{pmatrix}}\right.!=\Gamma \left({\bigl (}{\begin{smallmatrix}a+1&b\\c&d+1\end{smallmatrix}}{\bigr )}\right)={\frac {1}{2\Omega }}{\begin{pmatrix}\Gamma (\lambda _{1})\left(d-a+\Omega \right)+\Gamma (\lambda _{2})\left(a-d+\Omega \right)&-2b\left(\Gamma (\lambda _{1})-\Gamma (\lambda _{2})\right)\\-2c\left(\Gamma (\lambda _{1})-\Gamma (\lambda _{2})\right)&\Gamma (\lambda _{1})\left(a-d+\Omega \right)+\Gamma (\lambda _{2})\left(d-a+\Omega \right)\end{pmatrix}}}$ ^[12]

其中，${\displaystyle \lambda _{1}}$ 和${\displaystyle \lambda _{2}}$ 是${\displaystyle {\bigl (}{\begin{smallmatrix}a+1&b\\c&d+1\end{smallmatrix}}{\bigr )}}$ 的特徵值，分別為${\displaystyle \lambda _{1}=1+{\begin{smallmatrix}{\frac {\left(a+d-\Omega \right)}{2}}\end{smallmatrix}}}$ 和${\displaystyle \lambda _{2}=1+{\begin{smallmatrix}{\frac {\left(a+d+\Omega \right)}{2}}\end{smallmatrix}}}$ ，其中，${\displaystyle \Omega ={\begin{smallmatrix}{\sqrt {(a-d)^{2}+4bc}}\end{smallmatrix}}}$ ^[12]

變化

定义扩展

 

伽瑪函數

階乘的定義可推廣到複數，其与伽瑪函數的关系为：

${\displaystyle z!=\Gamma (z+1)=\int _{0}^{\infty }t^{z}e^{-t}\,\mathrm {d} t.\!}$ 。

伽瑪函數滿足${\displaystyle \Gamma (n+1)=(n)\Gamma (n)}$ ，

另一種定义扩展是阿達馬伽瑪函數，但由於其不在所有實數上皆能滿足階乘的遞迴定義，只有在正整數上滿足階乘的遞迴定義${\displaystyle n!=n\times (n-1)!}$ 因此比較少被拿出來討論。

${\displaystyle H(x+1)=x\,H(x)+{\frac {1}{\Gamma (1-x)}}}$ 

其後面的項${\displaystyle {\frac {1}{\Gamma (1-x)}}}$ 只有在正整數的情形為零。也因為其有加上一項，也因此，此擴展在描述負階乘時不會有除以零的情況，而使阿達馬伽瑪函數是一個處處連續、無奇點的函數。

遞進/遞降階乘

• 遞進階乘：${\displaystyle (x)_{n}=x^{\overline {n}}=x(x+1)...(x+n-1)}$ 
• 遞降階乘：${\displaystyle x^{\underline {n}}=x(x-1)...(x-n+1)}$ 
• ${\displaystyle x^{\overline {n}}=(-1)^{n}(-x)^{\underline {n}}}$ 

雙階乘

正整數的雙階乘表示小於等於該數的所有具相同奇偶性的正整數的乘積，即：

${\displaystyle {\begin{cases}(2n-1)!!=1\times 3\times 5\times \cdots \times (2n-1)\\(2n)!!=2\times 4\times 6\times \cdots \times (2n)\end{cases}},n\in \mathbb {N} }$ 

廣義的雙階乘

無視上述定義的${\displaystyle n!!}$ 因為即使值的${\displaystyle N}$ ，雙階乘為奇數可擴展到最實數和複數${\displaystyle z}$ 的注意到，當${\displaystyle z}$ 是一個正的奇數則：

${\displaystyle z!!=z(z-2)\cdots (3)=2^{\frac {z-1}{2}}\left({\frac {z}{2}}\right)\left({\frac {z-2}{2}}\right)\cdots \left({\frac {3}{2}}\right)=2^{\frac {z-1}{2}}{\frac {\Gamma \left({\frac {z}{2}}+1\right)}{\Gamma \left({\frac {1}{2}}+1\right)}}={\sqrt {\frac {2^{z+1}}{\pi }}}\Gamma \left({\frac {z}{2}}+1\right)\,.}$ 

獲得的表達接受一個以上公式${\displaystyle (2n+1)!!}$ 和${\displaystyle (2n-1)!!}$ 並表示在條件發生的階乘函數的${\displaystyle \gamma }$ 既可以看出（使用乘法定理）等同於一個給定在這裡。

${\displaystyle z!!}$ 定義為所有複數除負偶數。

比較上式與${\displaystyle (2n)!!}$ 的原始定義，廣義的雙階乘在${\displaystyle (2n)!!}$ 的計算上須包含${\displaystyle 0!!}$ ，即

${\displaystyle (2n)!!=2n\times (2n-2)\times (2n-4)\times \cdots \times 4\times 2\times 0!!}$ 

其中 ${\displaystyle 0!!={\sqrt {\frac {2}{\pi }}}}$ 

使用它的定義，半徑為${\displaystyle R}$ 的n維超球其體積可表示為：

${\displaystyle V_{n}={\frac {2(2\pi )^{\frac {n-1}{2}}}{n!!}}R^{n}.}$  n=1,3,5,...

${\displaystyle V_{n}={\frac {(\pi )^{\frac {n}{2}}}{{\frac {n}{2}}!}}R^{n}.}$  n=2,4,6,...

多重阶乘

${\displaystyle n!^{(k)}}$ 被称为${\displaystyle n}$ 的${\displaystyle k}$ 重阶乘，定义为：

${\displaystyle n!^{(k)}=\left\{{\begin{matrix}1,\qquad \qquad \ &&{\mbox{if }}0\leq n<k;\\n(n-k)!^{(k)},&&{\mbox{if }}n\geq k.\quad \ \ \,\end{matrix}}\right.}$ 

廣義的多重階乘

能將多重階乘推廣到複數（甚至是四元數）

${\displaystyle z!^{(k)}=z(z-k)\cdots (k+1)=k^{\frac {z-1}{k}}\left({\frac {z}{k}}\right)\left({\frac {z-k}{k}}\right)\cdots \left({\frac {k+1}{k}}\right)=k^{\frac {z-1}{k}}{\frac {\Gamma \left({\frac {z}{k}}+1\right)}{\Gamma \left({\frac {1}{k}}+1\right)}}\,.}$ 

四次階乘

所謂的四次阶乘（又称四重阶乘） 不是 ${\displaystyle n!^{4}}$ ，而是 ${\displaystyle {\frac {(2n)!}{n!}}}$ ，前幾個四次階乘為

1, 2, 12, 120, 1680, 30240, 665280, ....

它也等於

${\displaystyle {\begin{aligned}2^{n}{\frac {(2n)!}{n!2^{n}}}&=2^{n}{\frac {(2\cdot 4\cdots 2n)[1\cdot 3\cdots (2n-1)]}{2\cdot 4\cdots 2n}}\\[8pt]&=(1\cdot 2)\cdot (3\cdot 2)\cdots [(2n-1)\cdot 2]=(4n-2)!^{(4)}.\end{aligned}}}$ 

過階乘

hyperfactorial（有時譯作過階乘）寫作${\displaystyle H(n)}$ ，其定義為：

${\displaystyle H(n)=\prod _{k=1}^{n}k^{k}=1^{1}\cdot 2^{2}\cdot 3^{3}\cdots (n-1)^{n-1}\cdot n^{n}}$ 

hyper階乘和階乘差不多，但產生更大的數。hyper階乘的增長速度卻並非跟一般階乘在大小上相差很遠。 前幾項的hyper階乘為：

1, 4, 108, 27648, 86400000, ... （OEIS數列A002109）

超階乘

1995年，尼爾·斯洛恩和西蒙·普勞夫定義了超階乘（superfactorial）為首${\displaystyle n}$ 個階乘的積。即${\displaystyle \mathrm {sf} (n)=1!\times 2!\times 3!\times \cdots \times n!}$ 。一般來說

${\displaystyle \mathrm {sf} (n)=\prod _{k=1}^{n}k!=\prod _{k=1}^{n}k^{n-k+1}=1^{n}\cdot 2^{n-1}\cdot 3^{n-2}\cdots (n-1)^{2}\cdot n^{1}.}$ 

前幾項的超階乘為：

1, 2, 12, 288, 34560, 24883200, ... （OEIS數列A000178）

另一種定義

柯利弗德·皮寇弗在他的書Key to Infinity定義了另一個超階乘，寫作${\displaystyle n\mathrm {S} \!\!\!\!\!\;\,{!}}$ （ ${\displaystyle \mathrm {S} \!\!\!\!\!\;\,{!}}$ 為!和S重疊在一起）：${\displaystyle n\mathrm {S} \!\!\!\!\!\;\,{!}=n^{(4)}n}$ ^(4),表示hyper4，使用高德納箭號表示法即${\displaystyle n\mathrm {S} \!\!\!\!\!\;\,{!}=(n!)\uparrow \uparrow (n!)}$ 。這個數列：

${\displaystyle 1\mathrm {S} \!\!\!\!\!\;\,{!}=1}$ 

${\displaystyle 2\mathrm {S} \!\!\!\!\!\;\,{!}=2^{2}=4}$ 

${\displaystyle 3\mathrm {S} \!\!\!\!\!\;\,{!}=6\uparrow \uparrow 6=6^{6^{6^{6^{6^{6}}}}}}$ ，读作6个6重幂。

${\displaystyle 4\mathrm {S} \!\!\!\!\!\;\,{!}=(4!)\uparrow \uparrow (4!)=24\uparrow \uparrow 24}$  = ${\displaystyle {\begin{matrix}{24_{}^{24^{{}^{.\,^{.\,^{.\,^{24}}}}}}}\\\end{matrix}}}$ ，一直写24个24，读作24个24重幂。

質數階乘

主条目：質數階乘

質數階乘是所有小於或等於該數且大於或等於2的質數的積，自然數${\displaystyle n}$ 的質數階乘，寫作${\displaystyle n\#}$ 。

目前質數階乘只能用遞迴方式定義，因為尚未找到一個能用基本函數表示所有質數的函數或一條包含所有質數的曲線

一般情況下質數階乘定義為：

${\displaystyle n\#=\prod _{i=1}^{\pi (n)}p_{i}=p_{\pi (n)}\#}$ 

其中, ${\displaystyle \pi (n)}$ 是質數計數函數，小於或等於某個實數${\displaystyle n}$ 的質數的個數的函數${\displaystyle \leq n}$ 。

自然数阶幂

阶幂也称叠幂或者重幂记作${\displaystyle n^{!}}$ （感叹号!写在自然数的右上角），它的定义是将自然数1至${\displaystyle n}$ 的数由大到小作幂指数重叠排列，数学定义如下：

${\displaystyle n^{!}=n^{{(n-1)}^{!}}=n_{}^{(n-1)^{{(n-2)}^{.\,^{.\,^{.\,^{{3}^{{2}^{1}}}}}}}}}$ 

其中${\displaystyle n\geq 1}$ ，前几项的重幂数为：

1 , 2 , 9 , 262144 , ... （OEIS數列A049384）

第5个重幂数是一个有183231位阿拉伯数字组成的超大自然数^[13][14]，其值約為${\displaystyle 6.20606987866\times 10^{183230}}$ 

另外一種定義則是每個阶幂都先取一次階乘：

${\displaystyle n!^{{(n-1)!}^{!}}=n!_{}^{(n-1)!^{{(n-2)!}^{.\,^{.\,^{.\,^{{3!}^{{2!}^{1!}}}}}}}}}$ 

前幾個阶乘阶幂為：

1, 2, 36, 48708493958471199415506599153950129703565945470976, ... （OEIS數列A073581）

第5个阶乘阶幂值已大於${\displaystyle 10^{10^{50}}}$ ^[15][16]，其值約為${\displaystyle 4.3056\times 10^{1.01274\times 10^{50}}\approx 10^{10^{50.00549705084703}}}$ 

二次阶幂：

${\displaystyle n^{!!}=n^{{!}^{2}}={n^{{!}{(n-1)^{{!}{{(n-2)}^{{!}{{.}^{{.}^{{.}^{3^{{!}{2^{{!}{1^{!}}}}}}}}}}}}}}}}$ 

前幾個二次阶幂為：

1, 2, 81...

第4个阶乘阶幂值已大於${\displaystyle 10^{438}}$ ，其值約為${\displaystyle 7.975\times 10^{438}}$ 。

相应地，${\displaystyle m}$ 次阶幂定义如下：

${\displaystyle n^{{!}^{m}}=n^{{!}^{(m-1)}{(n-1)}^{!^{m}}}={n^{{!^{(m-1)}}{(n-1)^{{!^{(m-1)}}{{(n-2)}^{{!^{(m-1)}}{{.}^{{.}^{{.}^{3^{{!^{(m-1)}}{2^{{!^{(m-1)}}{1^{!^{(m-1)}}}}}}}}}}}}}}}}}$ 

其中${\displaystyle n}$ ，${\displaystyle m\geq 1}$ ，且${\displaystyle n,m\in Z}$ 。

倒數階乘

更多信息：倒數伽瑪函數 § 階乘倒數

倒數階乘是指所有小於及等於該數的正整數之倒數的積，其值與階乘的倒數相同：

${\displaystyle \prod _{k=1}^{n}{\frac {1}{k}}={\frac {1}{n!}}\quad \forall n\geq 1}$ 

其無窮級數收斂在e^[17]：

${\displaystyle \sum _{n=0}^{\infty }\prod _{k=1}^{n}{\frac {1}{k}}=e}$ 

考量階乘可以表示為連續的伽瑪函數，則有

${\displaystyle \int _{-1}^{\infty }{\frac {dx}{x!}}\,=\int _{0}^{\infty }{\frac {dx}{\Gamma (x)}}\,\approx 2.80777024,}$ 

這個值又稱為弗朗桑－羅賓遜常數（英语：Fransén–Robinson_constant）。^[18]

反階乘

 

反階乘的複變函數圖形

更多信息：反伽瑪函數 § 反階乘

反階乘是階乘的反函數，用於求解指定的數是哪個數的階乘。例如120的反階乘為5，因為5的階乘為120。反階乘可以透過泰勒級數或反伽瑪函數來評估與計算。

反階乘可以用了推算某個數大約是多少的階乘。

由於階乘與伽瑪函數之間的關聯，反階乘也可以透過反伽瑪函數近似公式來估計：

${\displaystyle \mathrm {ArcFactorial} \left(z\right)\approx -1+\alpha +{\sqrt {\frac {2\left(x-\Gamma \left(\alpha \right)\right)}{\Psi \left(1,\ \alpha \right)\Gamma \left(\alpha \right)}}}.}$ 

因此，反階乘也可以寫成如下的渐近分析形式：^[19]

${\displaystyle \mathrm {ArcFactorial} \left(x\right)\sim {\frac {\ln \left({\frac {x}{\sqrt {2\pi }}}\right)}{W_{0}\left(e^{-1}\ln \left({\frac {x}{\sqrt {2\pi }}}\right)\right)}}-{\frac {1}{2}}}$ 

其中${\displaystyle W_{0}(x)}$ 是朗伯W函数。這個公式是利用史特靈公式求逆得到的，因此也可以展開為漸近級數。

符号史

• 瑞士数学家欧拉（Euler, L.）于1751年用大写字母${\displaystyle M}$ 表示${\displaystyle m}$ 阶乘${\displaystyle M=1\cdot 2\cdot 3\cdot \cdots \cdot m}$ 。
• 意大利数学家鲁菲尼（Ruffini, P.）在1799年出版的方程著述中，用小写字母${\displaystyle \pi }$ 表示${\displaystyle m}$ 阶乘。
• 德国数学家高斯（Gauss, C.F）于1818年则用${\displaystyle \Pi (n)}$ 表示n阶乘。
• 用符号${\displaystyle {\underline {\mid n}}}$ 表示${\displaystyle n}$ 阶乘的方法起源于英国，尚不能确定其创始人，1827年，由雅来特（Jarrett）的建议得以流行，现代有时亦用此阶乘符号。
• 现在通用的阶乘符号${\displaystyle n!}$ 是法国数学家克拉姆（Kramp, C.）于1808年最先提出来的，后经德国数学家、物理学家格奥尔格·欧姆（Ohm, M.）等人的倡议而流行起来，直用到现在。

參見

• 伽瑪函數
• 斯特靈公式
• 階乘倒數
• 排列组合
• 威尔逊定理

註釋

1. 例如：${\displaystyle 1!=0!=1\,}$ ，$(-0.5)!={\sqrt {\pi }}$ ，$0.5!=0.5{\sqrt {\pi }}.$ 
2. The publisher of Stedman 1677^[4] is given as "W.S." who may have been William Smith, possibly acting as agent for the Society of College Youths, to which society the "Dedicatory" is addressed.

參考文獻

1. ^埃里克·韦斯坦因. Factorial. MathWorld. 

1. Graham, Ronald L.; Knuth, Donald E.; Patashnik, Oren, Concrete Mathematics, Reading, MA: Addison-Wesley, 1988, ISBN 0-201-14236-8 
2. Graham，Knuth & Patashnik 1988^[1], p.111
3. Biggs, Norman L. The roots of combinatorics. Historia Mathematica. May 1979, 6 (2): 109–136. ISSN 0315-0860. doi:10.1016/0315-0860(79)90074-0 –通过ScienceDirect. 
4. Stedman, Fabian, Campanalogia, London, 1677 ^{[註 2]}
5. Stedman 1677^[4], pp. 6–9.
6. Stedman 1677^[4], p. 8.
7. Higgins, Peter, Number Story: From Counting to Cryptography, New York: Copernicus, 2008, ISBN 978-1-84800-000-1 
8. Higgins 2008,^[7] p. 12
9. 潘承洞. 《数论基础》. 现代数学基础 （丛书）. 高等教育出版社. 2012年12月. ISBN 9787040364729 （中文（中国大陆））. 
10. 5.10. Digital Library of Mathematical Functions. [2010-10-17]. （原始内容于2010-05-29）. 
11. Abul-Dahab, MA and Bakhet, AK. A certain generalized gamma matrix functions and their properties. J Anal Number Theory: 63––68. 2015.  参数|journal=与模板{{cite web}}不匹配（建议改用{{cite journal}}或|website=） (帮助); |volume=被忽略 (帮助); 缺少或|url=为空 (帮助)
12. Cardoso, Joao R and Sadeghi, Amir. Computation of matrix gamma function. arXiv preprint arXiv:1806.10554. {2018.  参数|journal=与模板{{cite web}}不匹配（建议改用{{cite journal}}或|website=） (帮助); 请检查|date=中的日期值 (帮助); 缺少或|url=为空 (帮助)
13. print len(str(6**(5**(4**(3**2))))) （页面存档备份，存于互联网档案馆）^{[連結內容執行超時]}
14. Wolfram, Stephen. "${\displaystyle 5^{262144}}$ ". from Wolfram Alpha: Computational Knowledge Engine, Wolfram Research. [2018-11-19] （英语）. 
15. Wolfram, Stephen. "${\displaystyle 120^{24^{6^{2^{1}}}}}$ ". from Wolfram Alpha: Computational Knowledge Engine, Wolfram Research. [2018-11-19] （英语）. 
16. Sloane, N.J.A. (编). Sequence A073581 (Factorials successively exponentiated). The On-Line Encyclopedia of Integer Sequences. OEIS Foundation. a(5) > 10^(10^50). 
17. Iwanami Sūgaku Jiten Fourth, Tokyo: Iwanami Shoten, 2007, ISBN 978-4-00-080309-0, MR 2383190 （日语）  142.D
18. Finch, S. R. "Fransén-Robinson Constant." §4.6 in Mathematical Constants. Cambridge, England: Cambridge University Press, pp. 262-264, 2003.
19. Amenyou, Folitse Komla; Jeffrey, David. "Properties and Computation of the inverse of the Gamma Function" (学位论文): 28. 2018 [2023-10-03]. （原始内容于2022-05-09）.