想必我们都听说过Gauss小学的时候便快速算出了的答案5050,这让我们不禁好奇这个问题应该如何计算或者说如果将100改成变量n应该怎么办呢?
事实上,我们可以这么做:
仔细观察,这第二个公式似乎只是第一个公式的逆序但是让我们观察一下 ,第一个公式的第i项应该是i,而第二个公式的第i项应该是n+1-i,如果上下第i项相加应该都是n+1,而i从1到n,因此,我们知道,带入100正是5050,这让我们不禁感到神奇,原本需要加100次,现在我们只需要一个简单的公式就能算出来!
这让我不禁好奇这些东西是不是也可以像上面一样,直接化简为一个公式,答案当然是Yes!
新的办法
如果我们采用求这样子采用逆序求和我们会发现错的一塌糊涂,所以我们必须采用新的办法。但是在这之前我们需要先规定一些新的记法:
.容易知道
,
。让我们看一种更有启发式的办法来证明它们,首先我们发现
,如果将上面左侧的全部加起来右侧也加起来得到
,看上去还不错,让我们再试试
,这次我们用
重复上面的流程,容易发现
,咦?似乎有点门道,让我们移项化简,容易得到
。
这次直接让我们跳到x次幂使用这种方法,。(注意这里的
是组合数)还是熟悉的配方还是原来的味道,左右侧相加,我们得到了
.我们发现
依赖于
那么我们只要化简上式就可以得到
的递推式,不过用这种办法算的话似乎有点慢并且出错的概率极大(因为
依赖于
,所以只要有一个式子出错后面的计算便会前功尽弃),我曾经用这种办法算到了
计算量就已经非常的复杂,而且每次算出一个式子我就会代入三四个数字验证对错。不过现在我发现了一种投机取巧的方式,只要解出一组线性方程我们就可以得到
。
不妨让我们看几个例子吧!
我们发现当下标为i时,的最高次幂似乎是i+1次幂。不过毕竟只是针对前四项成立,如何证明他对第x项也成立呢?
我们可以用强归纳法:假设对于小于n的数字i而言,的最高次幂是i+1。
而我们已经证明,左侧多项式最高次幂显然是x+1,我们可以再假设
最高次幂不是x+1:
,
最高次幂大于x+1,如果是这样就会得到x+1次幂函数等于大于x+1次幂的函数,这显然是荒谬的!
,
最高次幂小于x+1,如果是这样就会得到x+1次幂函数等于小于x+1次幂的函数,因为我们的假设是对于小于n的数字i而言,
的最高次幂是i+1,那么左侧最高次幂不是
的最高次幂就是
的最高次幂,但是无论哪个都会小于x+1.
所以,最高次幂必然是x+1。
因此,我们知道,带入n=0,我们知道
=0,因此
,因此
,如果我们解出各个系数我们也就解出了
的公式,听上去似乎很不错,但是该怎么做呢?
当我们解出这个方程的时候各个系数也就一目了然了,对于这个方程我们可以使用matlab进行求解。具体代码实现大家可以自行尝试!