返朴

1ᵏ+2ᵏ+ ⋯ + nᵏ = ?

加星标,才能不错过每日推送!方法见文末插图

这里我们介绍一种基于分部积分法的处理方法,不仅能获得连续自然数幂次的一般求和公式,还能附带证明一个有趣的猜想。说不定读完本文后,聪明的你也会灵感迸发,开辟一个与众不同的新解法。

撰文 | 朱慧坚(玉林师范学院数学与统计学院副教授)、丁玖(美国南密西西比大学数学系教授)
读者朋友看到本文标题,或许会不以为意地发笑一声:这还不简单,不就是
Image
等等吗?您甚至可能还要说下去:上面第一个公式连小学一年级的学生都会求;“数学王子”高斯(Carl Friedrich Gauss,1777-1855)在年幼时就巧妙地俘获了第二个等式;至少在我们读大学一年级期间,学到函数在区间上的定积分时,老师在演示用某种特殊的黎曼和的极限计算值的过程中,就告诉过听课的学生上述第三个公式。这个方法是由古希腊大数学家阿基米德(Archimedes,c.287 – c. 212 BC)首创的。
对,您说的大致不错,如果再去查一查微积分教科书,或许还可以继续看到公式
Image
等等,等等。但是,如果有人请教您等于什么,您能马上告诉他相应的公式,或指导他如何下手吗?更进一步地,在仔细察看了上面的六个代数恒等式后,您可能会突然发现,平方和表达式恰好是四次方和表达式的一个因子,立方和恰好是五次方和的一个因子。这时,您肯定在纳闷这仅仅是巧合,还是说对六次和七次幂,甚至更高次幂的求和公式也都是如此。好吧,我们再瞧一瞧接下来的三个公式:
Image
结果发现同样的模式依然存在。这似乎是个有意义的观察,勾起我们的好奇心:对一般的正奇数或正偶数,从到的连续自然数次幂之和的表达式是否分别具有上面所说的因子?
在全世界流传百年之久的传奇故事是这样说的:高斯年幼时,他的老师为了整治班内淘气学童,要求学生计算。聪明的小高斯将其反方向写成个同一个数,该数乘以再除以,便在片刻之间算出答案
Image
这个数学传奇激励着天下俊杰对各种形式的自然数群体相加、相乘寻觅规律,破解谜团。作为一个佳例,从古到今不知有多少人热衷于发明新的方法,回答本文标题所提的数学问题。为了建立一个既封闭又悦眼的求和公式,所涉及到的数学工具实在是五彩缤纷,思路和技巧各种各样,又“条条大道通罗马”,显示出不同数学分支之间的密切相关性。在这些手段中,有的是基于组合数学的生成函数概念,有的是运用线性差分方程的理论,也有的甚至借用了初等微分学里的洛必达法则。几何中最有名的勾股定理据说已有约五百种证明,但我们还不清楚现在有多少种办法能够导出的表达式。说不定读完本文后,有个聪明而好奇心极强的读者会滋生灵感,磨砺出一把新斧头,开辟一个与众不同的解法。这里,我们介绍一种基于分部积分法的处理方法,它历史悠久,不仅能获得连续自然数幂次的一般求和公式,而且还能证明上述关于和式因子结构的猜想为真。它来自于数值积分中一个著名求和公式对单项式函数的巧妙应用,曾被作者用于证明同事提出过的上述猜想,那时作者正为大学高年级学生和研究生开设一门数值分析课程,恰好在课堂上讲授到那个积分近似公式。
Image

欧拉-麦克劳林求和公式

初等微积分里有关于不定积分和定积分的分部积分法,对于定积分,它指的是如下的等式:
Image
为叙述简单起见,假定两个函数和在区间上都具有连续的导函数。在(1)中,表达式image.png意指函数值和之差,即image.png。
学过初等微积分的读者都知道,不定积分的分部积分法就是函数导数的乘积法则
通过移项写成等价形式后,用不定积分的数学语言
而重新表述的。这和将陈述句“老布什总统是小布什总统的父亲”等价地写成“小布什总统是老布什总统的儿子”是同一个道理。至于定积分的分部积分法(1),它就是将上述不定积分的分部积分法同微积分的最重要定理——微积分学基本定理
; 一个原函数,即在区间结合后的自然产物。
在定积分的计算中,如果被积函数可以写成乘积的形式,并且上面(1)式右端的定积分值相较左端可以更容易地计算出来,那么此时的分部积分法算是用对了,这是定积分计算技巧中的通常思路。比如,要计算出定积分,如果令及,则,并且的一个原函数是。然后,(1)式给出
Image
上式右端的定积分是根据强大的微积分学基本定理赋值的:因为众所周知的一个原函数为,故该基本定理直接给出image.png反之,如果不恰当地选择及,则同一个公式(1)引来的仅仅是下面的等式
Image
这就与计算初衷南辕北辙,因为等式右边定积分中的被积函数看上去更加复杂,无法直接知道它的原函数是什么。这充分说明,在运用分部积分法之前,一定要聪明地选对公式(1)左端中的这两个因子函数和。为了本文要讲述的主题,我们先请读者注意,在(1)式中,只需是的一个原函数,因为根据微分学基本定理(即教科书里通常所称的拉格朗日中值定理),同一个连续函数的所有原函数之间只相差某个常数,因此该的选用有一个自由度,这给下面构造“欧拉-麦克劳林求和公式”提供了极大的方便。
现在,我们就“聪明地”选择分部积分法中的和,并反复运用(1)式,建立一种逼近定积分的数值积分方法——欧拉-麦克劳林求和公式。设想是一个光滑函数,意思是它直到某一个所需阶数的所有导函数在定义域上都是逐点连续的,其定义域为闭区间。我们的目的是要近似计算定积分。作为推导过程的主要步骤,我们先让
令。将定积分image.png(1)式左边的形式,其中我们如此选取一次多项式中的积分常数,以期满足额外条件image.png。这样,
Image
注意是区间的奇函数。运用一次分部积分法(1),并将上面的表达式代入,便有 
Image
同理,将上述等式右端的定积分写成image.png的形式,其中选取二次多项式
Image
这时,成了,此外还有。再用一次公式(1),便知
Image
Image
将上面的等式替换进前面一个等式里,我们就得到如下计算定积分image.png的一个公式:
Image
换言之,如果函数是二阶连续可微的,则其定积分可以被化为它的二阶导函数和二次函数乘积的定积分,尽管比起原先的定积分,这后一个定积分看上去更加复杂些,因为被积函数多了一个二次多项式的因子。然而,如果我们略微耐心一点,很快就会发现事实上这是一个聪明之举。
性急的读者或许从上式右端第二项的结构预测:它后面的积分可以继续用分部积分法展开,写成类似于它和第三项之和的形式,其中改为另一个数,函数在和一阶导数值分别升格为二阶导数值,其后的定积分中函数的二阶导函数和的乘积则进化成的三阶导函数和的乘积。您的思路当然是正确的,但我们最好还是从容不迫地再做一次分部积分法,或许能看到使人愉悦的新情况呢。好了,现在让我们求出三次多项式
Image
注意这时的奇函数满足,故有
Image
对上式再来一次分部积分法,其中所需的及image.png,
Image
从而的偶函数满足。综上所述,我们就获得了如下的等式:
Image
上述公式通过平移变换就可推广成计算定积分的一个实用公式。
为定义在区间上的一个光滑函数,则它与平移函数的复合函数上的光滑函数。将(2)式用于函数,结果就是
Image
现在令光滑函数的定义域为。对上面的等式求和,其中,再经过化简,我们便得到如下的一个等式
其中
遵循同样的思想,将上式中的函数的三阶导数值提高到任意奇数,我们便获得下列的一般性公式:
Image
其中。公式(3)就是数值积分领域里著名的欧拉-麦克劳林求和公式。等式的左端是对定积分用步长为的梯形法则得到的近似,等式右端的第二项则提供了某种数值“矫正”。注意到若被积函数的阶导函数恒为,则等式右边最末端的误差项就消失了。众所周知,多项式就是具有这个性质的好函数,该性质使得(3)是建立前个自然数同幂次求和公式的一个理想工具。
历史上,欧拉-麦克劳林求和公式于1732年由瑞士数学家欧拉(Leonhard Euler,1707-1783)发现,但未及时发表。1736年他与苏格兰数学家斯特林(James Sterling,1692-1770)通信时得知苏格兰数学家麦克劳林(Colin Maclaurin,1698-1746)对此获得更一般的结果后,放弃了自己的优先权。两人各自的公式分别发表于1738年和1742年。后人将此用途广泛的等式以他们的名字命名。公式推导中引进的这一族多项式被标准化后,与比他们更早半个世纪的瑞士数学家雅各布·伯努利(Jacob Bernoulli,1654-1705)发生了联系,而与后者名字也分不开的一列数字则对本文的主题贡献良多。
瑞士巴塞尔的伯努利家族声名显赫,学者云集,一门三代出了八位数学家。雅各布·伯努利是第一代中的“大哥”,他探索过的众多数学对象包括无穷级数和一类后来冠以他姓的微分方程,概率论中的大数定律是他的一大杰作,而银行家们最感兴趣的连续复利问题由他解决,导致绝对常数的出世(代表此数的通用字母e则来自于欧拉)。
雅各布·伯努利为自己设计的墓碑也蕴含了数学巧思,他希望在墓碑上刻一条数学意义丰富的对数螺线,根据雅各布本人的解释,这根自相似螺线“可以象征逆境中的坚韧和恒心,也可以象征人体在经历所有变化甚至死后,也能恢复到精确完美的状态。”可惜在他去世后,或许是雕刻墓碑的工匠误解了指示,错误刻下了不相干的阿基米德螺线。
Image
雅各布·伯努利的墓碑,螺线周围的拉丁文Eadem mutata resurgo意为“纵使变化,依然故我”|Wikipedia
雅各布的二弟约翰(Johann Bernoulli,1667-1748)微积分本领超越常人,常与大哥争强好胜,如比赛求解“最速降线问题”,甚至也曾向牛顿叫板,结果是后者提速了变分学的发展。约翰不仅培养和提携了少年欧拉,而且伯努利家族在他这一支涌现出的数学家最多。雅各布的大弟尼古拉斯是个画家,其子尼古拉斯一世(Nicolaus I Bernoulli,1687-1759)在伯父的教导下成长为数学家,对微分方程等学科多有建树。约翰的儿子丹尼尔(Daniel Bernoulli,1700-1782)对概率论有开创性研究,也是欧拉的同事和亲密战友。在伯努利家族中,上述几位是我们在大学教科书中见过面的杰出代表。第三代中的数学家约翰三世(Johann III Bernoulli,1744-1807)则是个多才多艺的罕见神童。
Image
伯努利多项式与伯努利数
上一节中出现的一列多项式。由它起步的所有后继多项式都通过下述两个性质被唯一地确定:对
Image
后面将要用到的伯努利多项式与上述多项式只相差一个常数因子。
对所有的非负整数,伯努利多项式由定义。在的前面乘上阶乘数,是为了使多项式中最高项的系数为,故所有的伯努利多项式都是首一多项式。
由的定义可知,伯努利多项式满足两个等式
Image
由上一节中给出的,它们各自所对应的伯努利多项式分别为:
Image
从上面可以看出伯努利多项式的一般模式:当为偶数时,偶函数,并且当当为奇数时,的奇函数,且。此外对所有都成立,而且当为大于的奇数时,
伯努利多项式中的常数项在本文中最有用,它们被称为伯努利数,记作
Image
根据上一段中的最后一句,当为奇数时。最前的两个伯努利数是面的几个偶数下标的伯努利数为
Image
显见它们是正负交错的。
现在我们用数学归纳法证明出伯努利多项式在标准的单项式基底下的一般表示形式。
引理1. 阶伯努利多项式具有表达式
Image
其中个东西放成一堆所有可能的组合个数。
证明.直接验证可知(4)式当时成立。现设当时(4)式为真,则当时,根据微积分学基本定理和等式,以及归纳假设,有
Image
上面推导中的倒数第三个等式来自于求和下标从到的平移变换。这就证明了(4)式对。证毕。
注1.因为,(4)式也可写为。
作为引理1的直接推论,我们有关于伯努利数的如下关系:
对所有的自然数,伯努利数满足递推关系
Image
证明.由于对大于的都满足,故有
Image
注2.因为(5)式两边都有同一项,它的等价形式是
Image
像一切有名的数那样,伯努利数也可以另法定义,其中普鲁士数论学家萨尔舒茨(Louis Saalschütz,1835-1913)于1893年列出的
Image
在表达上仅需加减乘除,在视觉上最直接明了。在大学教科书中,伯努利数则常常被定义为实解析函数
Image
更“高级”的定义涉及到著名的黎曼-zeta函数
Image
从上式可以推导出偶数下标伯努利数的一个显式表示
Image
由此公式以及zeta函数的性质,就可以严格证明所有偶数下标的伯努利数的符号是正负相间的,即对任意正整数,如同本节前面我们从到所看到的那样。
伯努利多项式还有一个在下一节里要用到的性质,我们也将它写成引理形式,并同引理1一样通过数学归纳法证之。
引理2. 伯努利多项式满足恒等式
Image
证明.直接计算可知上式对是对头的。设当时(6)式为真。令,则
Image
由归纳假设,我们得到微分恒等式,两边直接积分后有,而初始值条件则保证了积分常数,这就证明了(6)式对也成立。因此引理2得证。
Image
计算
我们现在回到欧拉-麦克劳林求和公式,希望将它用于单项式函数,其中为一自然数。先试验一下这个特殊情形,让自己热一热身,然后考虑一般情形,获取一个普适公式。由于的五阶及以上导函数恒等于0,这时在欧拉-麦克劳林求和公式(3)中取,得
Image
从而有
Image
它包含了的和作为一个因子。
现在考虑最一般的情形。假设注意到正整数无论为一奇数还是偶数,都有欧拉-麦克劳林求和公式(3)直接给出
Image
因而,
Image
这样就得出了前个自然数次幂的一个用偶数下标伯努利数表达的求和公式。
为了下一节的需要,我们需要对上式稍加变形,使得表达式更加美观实用,并能帮助解决一个小猜想。首先将下标为奇数的伯努利数也插进上面最后一行的求和记号中,并利用事实并且对成立。这样,
Image
如果我们将改为它的相反数,且将字母换写成花体image.png,而保持其他的伯努利数不变但另记成花体
Image
等式(8)传统上称为福尔哈伯公式(Faulhaber's formula),虽然更合理的名称应该是伯努利公式。我们还是固守原先的伯努利数,对(8)式中的求和下标做变量替换以获得另一个等价公式和某个定积分表达式。显见,当老下标从递增到时,新下标从递减到,因而(8)式转变成
Image
由前面引理1中伯努利多项式的表达式(4)知道,,加上引理2中的恒等式(6)推出image.png,故最终有
Image
这就给出了对任何正整数都成立的最简洁表达式(10)。读者可见,微积分之美就浓缩在表达自然数幂和这最后的定积分中!
(7)-(10)中的任何一个式子都回答了本文标题的提问,它给出了自然数从到的次幂求和的一个显式公式,这个表达式是的次多项式,其系数中的因子包含了伯努利数和组合数。
数学推导上我们休息片刻,轻松愉快地回顾一下这个公式出现之前的千年数学史。在两千年前的古希腊,毕达哥拉斯(Pythagoras,c. 572–c. 497 BC)、阿基米德等人就知道怎样求出前个自然数的和、平方和甚至立方和。在东方,一千五百年前的印度学者阿耶波多(Aryabhata,476–550)和一千年前的波斯数学家卡拉吉(Al-Karaji,c. 953–c. 1029)及伊斯兰黄金时代的学者伊本·海瑟姆(Ibn al-Haytham,c. 965–c. 1040)等也考虑过这些和式。中国古代数学家对此也有卓越的贡献,如南宋的杨辉(c.1238-c. 1298)和元代的朱世杰(1249-1314)等都深入研究过高阶等差级数,这与求解自然数幂和有异曲同工之妙。六十年前人民教育出版社出版的华罗庚教授《从杨辉三角谈起》这本小书,向中学生数学爱好者提供了对高阶等差级数用差分概念求和的思想,让少年时代的本文作者读得兴奋不已。
系统研究自然数幂次求和问题的首批人马迟至十六、十七世纪才降生,他们包括英国的哈里奥特(Thomas Harriot,1560-1621)、德国的福尔哈伯(JohannFaulhaber,1580–1635)以及法国的费马(Pierre de Fermat,1601–1665)和帕斯卡(Blaise Pascal,1623–1662)。哈里奥特是第一个用符号写出求和公式的人,然而他的公式止步于福尔哈伯不厌其烦地求得了直到的公式,发表在他于1631年出版的著作Academia Algebrae中,不过他并没有获得一般的求和公式。尽管如此,后人依然有时将(8)式称为福尔哈伯公式,以表彰他对此付出的辛勤劳动。
伯努利在历史上首次意识到数列 ,……的重要性,因为该序列为所有幂和提供了统一的公式。他在发现和的一般模式后百般欣喜,写下评述:“借助这张表,我花了不到半刻钟就找到了答案,将前1000个数字的10次方相加,得到的和是 91,409,924,241,424,243,424,241,924,242,500。”该成果于1713年被收进他身后出版的概率论著作Ars Conjectandi中,法国数学家棣莫弗(Abraham de Moivre,1667-1754)则提议将求和公式(8)中的那些常数,……命名为伯努利数。
Image
和式中的因子
现在,我们似乎可以解答本文开头的疑问了:为何当为奇数时,有因子为偶数时,有因子image.png?我们能对所有的奇数和偶数证明出同样的结论吗?
答案是肯定的,事实上,这就是用伯努利多项式来求出连续自然数同次幂和公式的一个有意义的副产品。为了下面的推导方便起见,我们将和数+记为。当,公式(9)给出
Image
而当,则有
Image
从而对偶数幂的情形,比值可以表达为
Image
由于,并且,故存在多项式
Image
若我们能够证明就是的一个因子。的确如此,这是由的奇函数这一事实所保证的。,因此关于前个自然数偶数幂和式的因子猜想是正确的,即
Image
对奇数幂的情形,我们有
Image
因为,我们只需要证明多项式也有两个根和,且重数至少为。证明如下:首先由伯努利数的定义,
Image
再从得
Image
最后,通过函数微分恒等式
Image
和数字等式,我们发现
Image
这就证明了是的一个因子,因此存在多项式
Image
最后,作为特例,取
Image
Image

相关阅读

1  发散级数怎样求和?

2 那些让学术名流刻骨铭心的数学题

3 所有自然数之和是-1/12?它在物理学中还有特别的应用?丨众妙之门

4 看数学大神欧拉如何处理微积分?“非常简单!”

近期推荐

1  “美国已不是过去的美国”,特朗普挥向高等教育的大棒会带来什么?

2  巴黎高师数学系4年仅招16名国际生,7人来自这个班

3  关于量子计算,我们仍不知道它到底能做什么

4  对话安德鲁·怀尔斯:寻找问题是数学家最关键的能力

5  走进“熵”的心脏:系综——吉布斯理论追溯

特 别 提 示

1. 进入『返朴』微信公众号底部菜单“精品专栏“,可查阅不同主题系列科普文章。

2. 『返朴』提供按月检索文章功能。关注公众号,回复四位数组成的年份+月份,如“1903”,可获取2019年3月的文章索引,以此类推。

版权说明:欢迎个人转发,任何形式的媒体或机构未经授权,不得转载和摘编。转载授权请在「返朴」微信公众号内联系后台。

找不到《返朴》了?快加星标!!
Image
Image

长按下方图片关注「返朴」,查看更多历史文章

Image
微信实行乱序推送,常点“在看”,可防失联
Image