數學中,正整數的階乘(英語:factorial)是所有小於等於該數的正整數,記為,例如5的階乘表示為,其值為120

並定義,1的階乘和0的階乘都為1,其中0的階乘表示一個空積[2]

實數範圍內的階乘函數,負整數除外[註 1]

1808年,基斯頓·卡曼引進這個表示法:,符號表示連續乘積,亦即。階乘亦可以遞歸方式定義:。除了自然數之外,階乘亦可定義於整個實數(負整數除外),其與伽瑪函數的關係為:

階乘應用在許多數學領域中,最常應用在組合數學代數學和數學分析中。在組合數學中,階乘代表的意義為個相異物件任意排列的數量,例如前述例子,其代表了5個相異物件共有120種排列法。在正整數的情形下,的階乘又可以稱為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]

定義

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

 

用連乘積符號可表示為:

 

從上述公式中,可以推導出遞歸關係式

 

但遞歸定義須給出起點,因此需要定義零的階乘。 除此之外,遞歸關係在階乘函數中各個值皆成立,例如:

 

0的階乘

為了將遞歸關係式擴展 ,因此需要定義0的階乘:

 

可以得到:

 

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

  •  的情況, 定義為「沒有任何數字相乘的結果」,所以更廣泛之慣例的例子是以不存在任何因數的乘法單位元素來當作其解。(參閱空積
  • 對於零個物品只有一種排列方式,因為沒有任何東西可以置換,唯一的重新排列就是什麼都不做。
  • 它使組合數學中的許多恆等式對所有適用的值皆有效,例如從空集合中選擇0個元素的方法數,可由二項式系數給出:
 .
而從空集合中選擇0個元素的方法數為一種,即沒有任何東西可以取,唯一的取法就是什麼都不做。定義 可以滿足:
 .
更一般地,在 個相異元素的集合中取出 個相異元素的方法數,可由二項式系數給出:
 .
其方法數只有一種,即全部取出。定義 可以滿足:
 
  • 此定義允許將許多公式更嚴謹地表達為冪級數,例如指數函數:
 

性質

 質因子分解 ,如 [9]

計算

階乘與斯特靈公式
 (藍色)、 (橘色),數字越大 會越趨近 。但 在負值則會因為出現虛數而無法使用。

計算 時,若 不太,普通的科學計算機都可以計算,能夠處理不超過 古高爾)數值的計算機可以計算至 ,而雙精度浮點數的計算機則可計算至 

 時,可用斯特林公式估計:  
更精確的估計是:  
其中  

部分函數值

部分的階乘值(OEIS數列A000142
n n!
0 1
1 1
2 2
3 6
4 24
5 120
6 720
7 5040
8 40320
9 362880
10 3628800
11 39916800
12 479001600
13 6227020800
14 87178291200
15 1307674368000
16 20922789888000
17 355687428096000
18 6402373705728000
19 121645100408832000
20 2432902008176640000
25 1.5511210043331×1025
50 3.0414093201713×1064
70 1.197857166997×10100
100 9.332621544×10157
450 1.733368733×101000
1000 4.023872601×102567
3249 6.412337688×1010000
10000 2.846259681×1035659
25206 1.205703438×10100000
100000 2.824229408×10456573
205023 2.503898932×101000004
1000000 8.263931688×105565708
10100 1010101.9981097754820

非正整數的階乘

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

Γ函數和Π函數

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

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

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

 

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

 

另外,我們也可利用此式以計算任意大於-1的實數的階乘:

 

複數的階乘

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

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

 為:

 

右圖顯示了幾個模(絕對值) 與輻角 的幾個等級,圖表的繪製範圍為 ,  個單位長。較粗的鉛直線為輻角值為 的等值線。

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

|z| < 1時,可使用泰勒級數來計算:

 

其中,γ歐拉-馬斯刻若尼常數ζ(z)黎曼ζ函數。部分計算機代數的系統存在可以直接產生這些展開式系數的語法。

z z!
實數
1、2、3、4、5 1、2、6、24、120 (OEIS數列A000142
    OEIS數列A019704
複數
   OEIS數列A212877)、(OEIS數列A212878
   
   
四元數
   
   
   
 
階乘的色相環複變函數圖形。顏色越深代表絕對值越接近零;顏色越接近白色代表絕對值趨於無窮。其中紅色為正實數、青藍色為負實數。

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

將階乘寫為 ,其中 為:

 

Stieltjes給出了其連分數值:

 

前幾項系數 [10]

n an
0  
1  
2  
3  
4  
5  
6  

負整數的階乘

負整數的階乘可透過階乘的遞歸定義 逆推而得:

 

但由於在此定義下計算負一的階乘會出現除以零(即 ),因此無法直接給出負整數的階乘。

其他數學結構的階乘

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

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

 

並且 ,其中, 是單位矩陣、 是一個方陣,同時 是一個非奇異矩陣[12]

換句話說,即矩陣 為單位矩陣的純量 倍,其階乘為 ,例如 

對於一個可對角化矩陣 其階乘為:

 [12]

其中,   特徵值,分別為  ,其中, [12]

變化

定義擴展

 
伽瑪函數

階乘的定義可推廣到複數,其與伽瑪函數的關係為:

 

伽瑪函數滿足 

另一種定義擴展是阿達馬伽瑪函數,但由於其不在所有實數上皆能滿足階乘的遞歸定義,只有在正整數上滿足階乘的遞歸定義 因此比較少被拿出來討論。

 

其後面的項 只有在正整數的情形為零。也因為其有加上一項,也因此,此擴展在描述負階乘時不會有除以零的情況,而使阿達馬伽瑪函數是一個處處連續、無奇異點的函數。

遞進/遞降階乘

  • 遞進階乘: 
  • 遞降階乘: 
  •  

雙階乘

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

 

廣義的雙階乘

無視上述定義的 因為即使值的 ,雙階乘為奇數可擴展到最實數和複數 的注意到,當 是一個正的奇數則:

 

獲得的表達接受一個以上公式  並表示在條件發生的階乘函數的 既可以看出(使用乘法定理)等同於一個給定在這裏。

 定義為所有複數除負偶數。

比較上式與 的原始定義,廣義的雙階乘在 的計算上須包含 ,即

 

其中  

使用它的定義,半徑為 n維超球其體積可表示為:

  n=1,3,5,...
  n=2,4,6,...

多重階乘

 被稱為  重階乘,定義為:

 

廣義的多重階乘

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

 

四次階乘

所謂的四次階乘(又稱四重階乘) 不是  ,而是  ,前幾個四次階乘

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

它也等於

 

過階乘

hyperfactorial(有時譯作過階乘)寫作 ,其定義為:

 

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

1, 4, 108, 27648, 86400000, ... (OEIS數列A002109

超階乘

1995年,尼爾·斯洛恩西蒙·普勞夫定義了超階乘(superfactorial)為首 個階乘的積。即 。一般來說

 

前幾項的超階乘為:

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

另一種定義

柯利弗德·皮寇弗在他的書Key to Infinity定義了另一個超階乘,寫作  為!和S重疊在一起): (4),表示hyper4,使用高德納箭號表示法 。這個數列:

 
 
 ,讀作6個6重冪。
  =  ,一直寫24個24,讀作24個24重冪。

質數階乘

質數階乘是所有小於或等於該數且大於或等於2的質數的積,自然數 質數階乘,寫作 

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

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

 

其中,  質數計數函數,小於或等於某個實數 的質數的個數的函數 

自然數階冪

階冪也稱疊冪或者重冪記作 (感嘆號!寫在自然數的右上角),它的定義是將自然數1至 的數由大到小作冪指數重疊排列,數學定義如下:

 

其中 ,前幾項的重冪數為:

1 , 2 , 9 , 262144 , ... (OEIS數列A049384

第5個重冪數是一個有183231位阿拉伯數字組成的超大自然數[13][14],其值約為 

另外一種定義則是每個階冪都先取一次階乘:

 
前幾個階乘階冪為:
1, 2, 36, 48708493958471199415506599153950129703565945470976, ... (OEIS數列A073581
第5個階乘階冪值已大於 [15][16],其值約為 

二次階冪:

 
前幾個二次階冪為:
1, 2, 81...
第4個階乘階冪值已大於 ,其值約為 

相應地, 次階冪定義如下:

 

其中  ,且 

倒數階乘

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

 

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

 

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

 

這個值又稱為弗朗桑-羅賓遜常數英語Fransén–Robinson_constant[18]

反階乘

 
反階乘的複變函數圖形

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

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

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

 

因此,反階乘也可以寫成如下的漸近分析形式:[19]

 

其中 朗伯W函數。這個公式是利用史特靈公式求逆得到的,因此也可以展開為漸近級數。

符號史

  • 瑞士數學家歐拉(Euler, L.)於1751年用大寫字母 表示 階乘 
  • 意大利數學家魯菲尼(Ruffini, P.)在1799年出版的方程著述中,用小寫字母 表示 階乘。
  • 德國數學家高斯(Gauss, C.F)於1818年則用 表示n階乘。
  • 用符號 表示 階乘的方法起源於英國,尚不能確定其創始人,1827年,由雅來特(Jarrett)的建議得以流行,現代有時亦用此階乘符號。
  • 現在通用的階乘符號 是法國數學家克拉姆(Kramp, C.)於1808年最先提出來的,後經德國數學家、物理學家格奧爾格·歐姆(Ohm, M.)等人的倡議而流行起來,直用到現在。

參見

註釋

  1. ^ 例如:   
  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, 具體數學, 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. ^ 4.0 4.1 4.2 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. 2015, 3: 63––68. 
  12. ^ 12.0 12.1 12.2 Cardoso, Joao R and Sadeghi, Amir. Computation of matrix gamma function. arXiv preprint arXiv:1806.10554. {2018. 
  13. ^ print len(str(6**(5**(4**(3**2)))))頁面存檔備份,存於互聯網檔案館[連結內容執行超時]
  14. ^ Wolfram, Stephen. " ". from Wolfram Alpha: Computational Knowledge Engine, Wolfram Research. [2018-11-19] (英語). 
  15. ^ Wolfram, Stephen. " ". 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).