自然數(shù)的等冪和——伯努利數(shù)

高斯幼年時的老師為了刁難學(xué)生,在黑板上寫下了這個式子
這項(xiàng)麻煩的工作讓全班同學(xué)都在為它忙著,但此時高斯脫口而出了答案,等于5050,讓全班同學(xué)都驚呆了
這是一個廣為流傳的故事,經(jīng)過多次修改后產(chǎn)生了許多版本,但事實(shí)上沒人知道高斯當(dāng)時到底是用的什么方法,
不過我們今天不僅僅只解決這一個式子
其實(shí)這個故事之前有一個小插曲,17世紀(jì)時雅各布伯努利(Jacob Bernoulli)他的《猜度術(shù)》中說他能發(fā)現(xiàn)1到1000的十次方之和,這可以在七分半內(nèi)解決,最后的結(jié)果是91409924241424243424241924242500
這時你也許會想:三十二位數(shù),什么鬼?七分半!手算?
當(dāng)然不可能是手算的,為了弄清楚他的方法,我們進(jìn)入今天的主題

自然數(shù)等冪和
首先我們考慮
找到它的通項(xiàng),很簡單,將重新排列
兩式相加,其中每一項(xiàng)都等于n+1,且一共n項(xiàng),即
接下來提升難度
要找到它的通項(xiàng),乍一看有些不知所措,不過我們注意到一次冪的和通項(xiàng)為二次的多項(xiàng)式,不妨從三次方來考慮:
將后面的括號展開
? ? ? ? ? ? ? ? ? ? ? ? ? ? ?
對x從1加到n,上式變?yōu)?/p>
根據(jù)上面的結(jié)果就能得到
接下來難度繼續(xù)增加,
沿用二次冪的思路,將
中的括號用牛頓二項(xiàng)式定理展開
其中為二項(xiàng)式系數(shù),
對x從1加到n,左側(cè)的差分就變?yōu)?img type="latex" class="latex" src="http://api.bilibili.com/x/web-frontend/mathjax/tex?formula=n%5E%7Bk%2B1%7D" alt="n%5E%7Bk%2B1%7D">,而右邊則可交換求和次序,上式變?yōu)?/p>
于是,我們可以得到自然數(shù)等冪和的遞推關(guān)系式
伯努利數(shù)
由上可知前n個自然數(shù)的k次冪和的通項(xiàng)為一k+1次多項(xiàng)式
不過這個遞推表達(dá)式有些復(fù)雜了,于是我們來試著簡化它一下
不妨來看一看這些多項(xiàng)式中有沒有什么規(guī)律,這里我已經(jīng)計算了到
以降冪次排列,當(dāng)中第二列的系數(shù)和后面偶數(shù)列的系數(shù)0格外顯眼,于是我將它標(biāo)上了紅色
首先由前面的結(jié)論我們知道最高次項(xiàng)的系數(shù)為
,又大概知道了所有偶數(shù)列的規(guī)律,接下來看第三列,當(dāng)中
引起了我們注意,于是我們發(fā)現(xiàn)似乎
第三列系數(shù)的規(guī)律是
,但其他列的規(guī)律就有點(diǎn)不好找了,那么我們來看看Bernoulli的想法吧
Bernoulli猜測自然數(shù)等冪和按降次冪排列的第n列取決于該列第一個數(shù),后來人們將這些數(shù)命名為伯努利數(shù)(Bernoulli?numbers),這里將第n個伯努利數(shù)記為??,且?
嗯,這跟我們發(fā)現(xiàn)的規(guī)律有一丟丟相似,下面我們就從伯努利的思路出發(fā)吧
首先我們應(yīng)該找出一個能夠計算伯努利數(shù)比較方便的式子,根據(jù)
對n取導(dǎo)數(shù),可以得到
這顯然是一個十分弱智的行為,但如果將n=0代入,其實(shí)就是第u個伯努利數(shù)
再根據(jù)就能得到伯努利數(shù)的一種計算方法了
有了伯努利數(shù)的計算方法,就可以計算前幾個伯努利數(shù)了,分別是
對上式求“高階導(dǎo)”并代入n=0亦能得到自然數(shù)等冪和的通項(xiàng)中其他系數(shù)間的關(guān)系

接下來我們回歸主題——等冪和
引入??的生成函數(shù):
由于的遞推關(guān)系比較復(fù)雜,我們直接將其定義代入
注意到其中綠色部分為零點(diǎn)處的Taylor級數(shù)展開,再根據(jù)等比數(shù)列求和公式,有
這樣就能用來計算
了:
當(dāng)然,對這個分式求導(dǎo)是個麻煩的工作,我們來尋找一種更簡便的方法
在生成函數(shù)中對n取“導(dǎo)數(shù)”,并代入n=0
于是我們又得到了伯努利數(shù)的生成函數(shù)定義,將它代入到的生成函數(shù)中
再將在零點(diǎn)處展開為Taylor級數(shù)
由柯西乘積公式,可得
對比的定義中的系數(shù),有
最終,我們將二維的自然數(shù)等冪和用一維的伯努利數(shù)表示了出來,實(shí)現(xiàn)了“降維打擊”

第二類伯努利數(shù)
前面所介紹的是第一類伯努利數(shù)(也叫“原始的伯努利數(shù)”)
在上面的公式中,令n=1,因,有
取k為任一大于1的奇數(shù),又有
將兩式作差,偶數(shù)項(xiàng)就全被減掉了,
我們已經(jīng)計算了,則
代入到上式中,
因?yàn)樗?strong>二項(xiàng)式系數(shù)均不為零,又根據(jù)k可取任意大于1的奇數(shù),我們得到
根據(jù)伯努利數(shù)的生成函數(shù)定義
用替換
,
當(dāng)中所有偶數(shù)項(xiàng)均與伯努利數(shù)相同,而又由于伯努利數(shù)中除1外所有奇數(shù)項(xiàng)都等于0,所以其實(shí)與
的Maclaurin級數(shù)展開中只有一次項(xiàng)系數(shù)不同
就是第二類伯努利數(shù),只需在第一類伯努利數(shù)中將
換成
就可以得到,雖然實(shí)際上它跟第一類伯努利數(shù)就只有一個數(shù)的區(qū)別,但這一個數(shù)的區(qū)別卻使第二類伯努利數(shù)在應(yīng)用中方便許多,因此許多地方用的都是第二類伯努利數(shù)

本期只介紹了Bernoulli numbers的起源——等冪和問題,其實(shí)這只是它的應(yīng)用中的冰山一角,它與Riemann?Zeta函數(shù)偶數(shù)、負(fù)奇數(shù)處的值密切相關(guān),正切函數(shù)的麥克勞林展開中也有它,此外還有許多地方都能看到它的影子
Bernoulli numbers因此視為數(shù)學(xué)中最重要的數(shù)列之一
