朱世傑一 一范德蒙公式的發展簡介
羅見今
摘要: 朱世傑−范德蒙公式 (the Chu-Vandermonde formula) 是現代組合計數 理論的一個基本公式。 本文闡明了該公式的組合意義, 分別論述了朱世傑 (1303) 和 A. T. Vandermonde (1772) 的歷史貢獻, 追蹤錢寶琮 1923 年首次解讀朱世傑的 成果和傳播海外的歷程, 並對朱−范公式的現代發展作出概要介紹。
關鍵詞: 組合數學、 組合計數、 朱世傑−范德蒙公式、 數學史。
現在國內外組合分析家所發表的一些論文和著作中, 當涉及對兩組合的乘積求和或 “卷積 型” 粧等式時, 時常提到 “the Chu-Vandermonde formula” (下文簡稱 “朱—范公式”)。 這 裡的 “Chu” 即朱世傑 (字漢卿, 號松庭), 我國宋元數學四大家之一, 他在 1303 年出版了 hh四 元玉鑒ii, 對中世紀世界數學卓有貢獻。 Vandermonde 就是學高等代數的人熟知的法國數學家 范德蒙 (A. T. Vandermonde, 1735∼1796, 或譯為旺德蒙德)。 二百多年來, 朱−范公式已有 Rothe (1793), Gauss, Hagen (1891), Gould (1956), Handa 和 Mohanty (1969) 等人的 多種擴充, 應用甚廣。 因此, 有必要對朱世傑、 Vandermonde 的工作和後來的進展作一概述。
1. hh四元玉鑒ii 中的組合粧等式
朱世傑在 hh四元玉鑒ii[1] (1303) “茭草形段”、 “如象招數” 和 “果垛疊藏” 中用傳統的垛 積招差術獲得了一批重要成果, 其中有朱世傑粧等式:
n
X
k=1
k + p − 1 p
=n + p p + 1
(1) 它是賈憲三角形中每斜行前 n 項和的計數公式, 表明了組合的一項基本性質, 作者認為與帕斯 卡粧等式具有同樣的重要性。 式 (1) 右邊的組合數 (在賈憲和朱世傑時僅為計算二項式係數的 演算法, 不含從多少元素中取出多少元素不同取法總數的組合意義) 在中國數學史上叫做三角 垛系列, 朱世傑給出 p = 2, 3, 4, 5, 6 時的名目依次為: 茭草垛、 三角垛、 撒星形垛、 撒星更落 一形垛、 三角撒星更落一形垛。
66
在通常的組合數學書和組合粧等式表[2]中常見的等式是
m
X
k=0
k p
=m + 1 p + 1
(2)
及 q
X
k=0
p + k k
=p + q + 1 q
(3) 其實兩者均為朱世傑粧等式的等價公式, 表示賈憲三角形中的同一性質。
證明: 由 (3) 知
q
X
k=0
p + k p
=p + q + 1 p + 1
, 令 m = p + q, q = m − p, 代入得
m−p
X
k=0
p + k p
=m + 1 p + 1
, 亦
m
X
k=0
k p
=m + 1 p + 1
, 即式 (3) 變換成 (2)。
再證 (1) (2) 兩式等價。 令 m = n + p − 1,
n
X
k=1
k + p − 1 p
=
m−p+1
X
k=1
k + p − 1 p
=
m
X
k=p
k p
=
m
X
k=0
k p
=m + 1 p + 1
式 (1) 即變換成 (2) 式; 反之, (2) (3) 式亦可變換成 (1) 式。 證完。
據清代學者羅士琳的解釋[3], 可將 “茭草形段” 第四題給出的關係表示為
n
X
k=1
k + 1 2
n + 2 − k 2
=
n
X
k=1
k + 3 4
=n + 4 5
(4) 在原著中 “撒星更落一形垛” 可用式右的組合表示, 是式 (1) 中 p = 5 時的情況。 羅士琳將 ak = k(k + 1)/2 叫做 “三角積”, bk = (n − k + 2)(n − k + 1)/2 叫做 “逆列三角積”。 顯然 Σakbk 構成卷積。
朱世傑在 “果垛疊藏” 第六題中依據同樣的關係和方法求得
n
X
k=1
k + 2 3
n + 2 − k 2
=
n
X
k=1
k + 4 5
=n + 5 6
(5) 朱世傑命名的 “撒星更落一形果子積”, 其組合形式應如式 (5) 右邊所示, 是式 (1) 中 p = 6 時 的情況。 羅士琳將 ak = (k + 2)(k + 1)k/6 叫做 “三角積積數”, bk(n + 2 − k)(n + 1 − k)/2 叫做 “逆列三角積”, 顯然 Σakbk 構成卷積。 不難看出, (4)、 (5) 兩式屬於更廣泛一類公式的兩 個特例, 本文將其推廣, 獲得
n
X
k=1
k + p − 1 p
n + q − k q
=
n
X
k=1
k + p + q − 1 p + q
=n + p + q p + q + 1
(6)
當 p = 2, q = 2 時便是式 (4), 當 p = 3, q = 2 時便是式 (5), 思路較為清晰。 有理由 認為朱世傑已掌握這類公式, 推廣它是題中應有之義。 我們將式 (4)、 (5)、 (6) 稱為朱世傑組合 卷積公式。 推測它的立術之原, 不大可能是從代數變換而得, 朱世傑一定是從大量操演垛積實物 獲得了其間的數量關係, 因此可以說, 垛積成為中國傳統數學早期一種有效的數學模型。
2. Vandermonde 公式的組合意義
A. T. Vandermonde 是法國科學院院士 (1771), 對高等代數的發展有重要貢獻。 Van- dermonde 公式在三種中譯本的 hh范氏大代數ii 中皆有介紹, 並出現於近年出版的組合數學專 著、 教材的例題和習題中[4,5]。 然而, 大概是缺乏歷史感之故, 有的書卻未注明范德蒙公式。
在研究怎樣對兩組合的乘積求和、 同時考慮遊動變數可能出現的各種位置時, 范德蒙 1772 年向法國科學院提交了一篇論文[6], 其中的一個組合公式具有一般性, 在同類公式中最為重要, 通常被稱為 “Vandermonde 綜合式”, 原式為:
X
k
r k
s n − k
=r + s n
(7a) 可以這樣來理解這個公式: 式右是從 r 個男人和 s 個女人中來選擇 n 個人的方法數, 而 式左中的每一項是從男人中選出 k 個同時從女人中選出 n − k 個人的方法數。 設 α, β 為任 意實數, 則 Vandermonde 公式的一般形式為
n
X
k=0
α k
β n − k
=α + β n
(7b) 特別當 α, β 為正整數時, 它的組合意義是: 從 α 中取 k (k = 0, 1, . . . , n) 個, 同時從 β 中取 n − k (n − k = n, n − 1, . . . , 1, 0) 個, 其不同取法的總數, 與從 α + β 中取 n 個不同 取法的總數相等。
這是一個組合卷積 (convolution, 以前曾譯為 “摺積”) 公式。 設兩函數 f 與 g 構成卷 積, 記作 f ∗ g, 則其意思即指
f (0)g(n) + f (1)g(n − 1) + · · · + f (n)g(0) = (f ∗ g)(n)
亦即 (f ∗ g)(n) =
n
X
k=0
f (k)g(n − k) (8)
此處展開式中共含 n + 1 項。
可以看出, 上述垛積術中正列 “三角積” 與 “逆列三角積” 逐項對應相乘求和, 恰是卷積 的形式, 朱世傑的垛積公式 (4)、 (5)、 (6) 即屬於卷積型的。
演算法和程式設計技術的先驅者、 美國斯坦福大學教授 Donald E. Knuth (克努特, 中 文名高德納) 對中國數學感興趣, 他曾關注賈憲、 楊輝、 朱世傑以及明安圖[7]的計數成就。 在他 的名著 The Art of Computer Programming hh電腦程式設計藝術ii (1998) 中, 講到 (6)、
(7) 兩式之間的關係, 指明朱世傑的組合公式在先[8], 材料出處是李約瑟博士的 hh中國的科學與 文明ii 第 3 卷數學。
3. 朱世傑卷積公式的解讀與傳播
錢寶琮[9] 1923 年在 hh學藝ii 上發表 “朱世傑垛積術廣義”, 首次將上文所稱的朱世傑卷 積公式推廣, 得到 (6) 式 (但不包括 (6) 式第二個等號右邊):
Xr|p|
1|p| · (n + 1 − r)|q|
1|q| =Xr|p+q|
1|p+q|
這裡 r|p|= r(r + 1) · · · (r + p − 1), 1|p|= p! 這是當時所用的記號。
湯天棟[10] 1926 年在 hh科學ii 上發表 “茭草形段羅草補注”, 用圖形和代數式表示 “茭草 形段” 部分的內容 (不涉及用組合符號表示的朱世傑卷積公式)。
錢先生 1929 年在 hh中國算學史ii 上編中再次引述, 1932 年再版[11]。
方淑姝[12] 1939 年在 hh數學雜誌ii 上發表 “朱世傑垛積術廣義”, 她使用的符號更接近現 代常用的組合符號, 提到了錢寶琮和湯天棟的成果。
現代科學史學科的奠基人、 著名科學史家喬治 · 薩頓 G. Sarton[13]較早將此公式介紹到 西方; 李約瑟的 hh中國的科學與文明ii (1953, 譯名 hh中國科學技術史ii) 第 3 卷數學[14] (王玲 參與) 中也引用了這一結果。
需要明確, 正是錢寶琮先生 1923 年首次將朱世傑的式 (6) 表述成現代組合形式發表在 hh學藝ii 上, 經喬治 · 薩頓、 李約瑟的介紹, 西方學者才對它有所瞭解[15]。至今在現代數學界流 傳開來, 錢先生發軔於先, 功不可沒。
但同時也須指出, 我國數學史界對此問題亦有不同的看法, 認為 “實際上這種類型的求和 問題已經超出了朱世傑著作中所有問題的範圍”, “這一不正確的論斷可能是受到了羅士琳的影 響。”[16]換言之, 朱世傑並無此式。 這涉及到羅士琳對 hh四元玉鑒ii 此題的解讀是否不著邊際或 無中生有, 是需要針對朱、 羅原著逐行、 逐句、 逐詞查證的。 本文的推證已見於前, 限於篇幅, 這 裡不能詳細說明; 如果感興趣, 這個問題是可以展開討論的。
4. Chu-Vandermonde 卷積公式的發展
Gauss 研究超比例級數時曾將類似的卷積關係推廣到複數的情形[17], 被稱為 Gauss- Vandermonde 公式。
Rothe(1793)-Hagen(1891)-Gould(1956) 的卷積型粧等式可寫成[18,20]
n
X
k=0
x x+kz
x+kz k
y
y+(n−k)z
y+(n−k)z n − k
= x + y x+y+nz
x+y+nz n
(9)
其中, Hagen 的粧等式為[19]
n
X
k=0
x + kz k
y − kz n − k
x
x + kz =x + y n
(10)
顯然 (10) 式中令 k = 0 便得出 (7) 式。 Mohanty-Handa (1969)[21]曾將 (9) 推廣成多 重和形式。 又, Hagen[19] 還將 (9) 發展成如下形式
n
X
k=0
x+kz k
y−kz n − k
p + qk
(x+kz)(y−kz) = p(x+y−nz)+nxq x(x + y)(y − nz)
x+y n
(11) 特別當 p = y, q = −z 時, 便有
n
X
k=0
x + kz k
y − kz n − k
1
x + kz = 1 x
x + y n
(12) 事實上可以證明, (9) 與 (11) 是等價的, 即兩者可互推。 當然它們均為朱−范公式的推廣。
注意到組合數學界同感興趣的事實——歷來公認為 Rothe-Hagen-Gould 卷積型粧等式 是朱−范公式的 “非平凡推廣” (non-trivial extension), 這裡卻有必要說明, 只須利用 Hsu- Gould 反演[22]公式, 即可證明 (11) 與 (9) 都可由朱−范公式推證。 這說明 (9) 與 (11) 實質 上仍與原始的朱−范公式等價 (詳見文獻[18][22])。
謹以此文紀念 hh四元玉鑒ii 成書 705 周年。
致謝: 徐利治 (Leetsch C. Hsu) 教授和羅伊德 (Keith Lloyd) 教授為本文提供了參考 文獻, 特致謝忱。
參考文獻
1. (元) 朱世傑, hh四元玉鑒ii (1303), 光緒二年 (1876) 丁取忠校刊本。 任繼愈主編, hh中國科學技術 典籍通匯ii 數學卷一, 河南教育出版社, 1994 年。
“茭草形段” 第四題為: 今有茭草八千五百六十八束, 欲令撒星更落一形垛之, 問底子幾何? 答曰: 一 十四束。 術曰: 立天元一為撒星更落一底子, 如積求之, 得一百二萬八千一百六十為益實, 二十四為 從方, 五十為從上廉, 三十五為從二廉, 一十為從三廉, 一為正隅。 四乘方開之, 合問。
hh四元玉鑒ii “果垛疊藏” 第六題為: 今有三角撒星更落一形果子積九百二十四個, 問底子幾何? 答 曰: 七個。 術曰: 立天元一為三角撒星更落一底子, 如積求之, 得六十六萬五千二百八十為益實, 一 百二十為從方, 二百七十四為從上廉, 二百二十五為從二廉, 八十五為從三廉, 一十五為從四廉, 一 為正隅。 五乘方開之, 合問。
2. Henry W. Gould, Combinatorial Identities, A standardized set of tables listing 500 binomial coefficient summations, Morgantown, W. Va., 1972.
3. (清) 羅士琳, hh四元玉鑒細草ii, 道光十五年 (1835) 揚州李棠寫本。 羅士琳的細草對朱世傑每道題 目增添了三行數表及詳解的 “術” 和 “草”。
4. Daniel I. A. Cohnen 著, 左孝淩等譯, hh組合理論的基本方法ii, 北大出版社, p.36, 1989。
5. 邵嘉裕, hh組合數學ii, 同濟大學出版社, p.14, 1991。
6. A. T. Vandermonde, M´em. Acad. Roy. Sciences Paris, p.489-498, 1772.
7. D. E. Knuth, The Art of Computer Programming, Vol. 1, Fundamental Algorithms, Third Edition (英文影印版), 清華大學出版社, p.407, 2002。
8. D. E. Knuth 克努特 (高德納) 著, 蘇運林譯, hh電腦程式設計藝術ii 第 1 卷基本演算法, 第 3 版, 國防工業出版社, p.55-56, 2002。
9. 錢寶琮, 朱世傑垛積術廣義, hh學藝ii 第 4 卷第 7 號, 1923; 另見 hh李嚴錢寶琮論文集ii 第 1 卷 hh古算考源ii, 遼寧教育出版社, p.45-54, 1998。
10. 湯天棟, 茭草形段羅草補注, hh科學ii 第 11 卷, p.1535-1558, 1926。
11. 錢寶琮, hh中國算學史ii 上編, p.134, 1932。
12. 方淑銖, 朱世傑垛積術廣義, hh數學雜誌ii 第 1 卷第 3 期, p.94-101, 1939.
13. G. Sarton, Introduction to the History of Science, Vol.3, p.701.
14. J. Needham and Wang Ling, Science and Civilization in China, Vol.3, Mathematics, (10) Series and Progressions, Cambridge at the University Press, p.138-139, 1959.
15. Ronald L. Graham, Donald E. Knuth and Oren Patashnik, Concrete Mathematics, Addison-Wesley, Reading MA, p.169-170, 1989.
16. 杜石然, 朱世傑研究, 見錢寶琮等者 hh宋元數學史論文集ii, 科學出版社, p.191。
17. Ronald L. Graham, Donald E. Knuth and Oren Patashnik, Concrete Mathematics, Addison-Wesley, Reading MA, p.212, 1989 (in a section on hypergeometric functions).
18. Fangjun Hsu and Leetsch C. Hsu, A unified treatment of a class of combinatorial sum, Discrete Mathematics, Vol.90, 1991, North Holland, p.191-197.
19. Johann G. Hagen, Synopsisder Hoeheren Mathematik, Vol. I, Arithmetische und al- gebraische analyse. Chapter III: “Theorie der combinationen”, Verlag Felix L. Dames, Berlin, p.55-68, 1891.
20. H. W. Gould, Some generalizations of Vandermonde’s convolution, Amer. Math. Mon- thly, Vol.63, p.84-91, 1956.
21. S. G. Mohanty and B. R. Handa, Extensions of Vandermonde type convolutions with serveral Summations and their applications, Canad. Math. Bull. Vol.12, p.45-62, 1969.
22. H. W. Gould and L. C. Hsu, Some new inverse series relations, Duke Math. J., Vol.40, p.885-891, 1973.
—本文作者現任教浙江大學科學技術與文化研究所—
• 緣 起: 為積極發掘並培養青少年數學人才, 透過專題研究, 培養新一代中學生的 數學素養, 引發青年人探索知識的興趣及提升學術水準, 故設立本數學獎。
• 參賽資格
(1) 現就讀高中職或以下之學生, 以個人身份參賽。
(2) 需有現任教高中職或以下之老師一人為指導老師。
• 參賽規定與時間
參賽的數學作品可以是純粹的數學研究, 也可以是數學應用在其他領域的綜合研 究。 唯必須具有原創性, 並且不可在其它科學展覽或數學競賽中曾參賽獲獎 (參 賽未獲獎者不在此限), 作品可以使用中文或英文書寫。 每年三月三十一日前於 網路報名; 五月三十一日前將參賽作品之完整 pdf 電子檔寄至台灣大學數學系。
• 評審方式
第一階段為初賽 (六月十六日至六月三十日), 進行書面審查, 選出若干隊伍; 第 二階段為決賽 (七月十七日), 進行口頭報告。
• 獎勵辦法
(1) 金牌獎 (不超過 1 名, 得從缺)
得獎學生可獲得證書、 獎座及研究補助金新台幣 60,000 元。 之後, 得獎者 若就讀國內外大學數學系或應用數學系, 並可獲得四年獎學金, 每年新台幣 120,000 元。 指導老師可獲得證書、 獎座及研究補助金新台幣 30,000 元。
(2) 銀牌獎 (不超過 2 名, 得從缺)
得獎學生可獲得證書、 獎座及研究補助金新台幣 40,000 元。 之後, 得獎者 若就讀國內外大學數學系或應用數學系, 並可獲得四年獎學金, 每年新台幣 60,000 元。 指導老師可獲得證書、 獎座及研究補助金新台幣 20,000 元。
(3) 佳作獎
得獎學生可獲得證書、 獎座及研究補助金新台幣 20,000 元。 指導老師可獲 得證書、 獎座及研究補助金新台幣 10,000 元。
(4) 團體獎
同一學校有兩件 (含) 以上作品獲獎則可獲頒團體獎座。 獲獎學校可優先與 台大數學系合作, 於學期中安排教授至該校進行數學教學研究之推廣活動。
• 聯絡方式
國立台灣大學數學系, 活動網址 http://www.math.ntu.edu.tw/∼shing tung/
聯絡電話:(02)3366-9902, 傳真:(02)3366-4471,mail: [email protected]
