Martian148's blog Martian148's blog
首页
  • ICPC 算法笔记
  • ICPC 算法题解
  • 体系结构
  • 高等数学
  • 线性代数
  • 概率论与数理统计
  • 具体数学
  • Martian148的奇思妙想
  • 游记
  • 通识课笔记
关于
  • useful 网站
  • 友情链接
  • 分类
  • 归档

Martian148

一只热爱文科的理科生
首页
  • ICPC 算法笔记
  • ICPC 算法题解
  • 体系结构
  • 高等数学
  • 线性代数
  • 概率论与数理统计
  • 具体数学
  • Martian148的奇思妙想
  • 游记
  • 通识课笔记
关于
  • useful 网站
  • 友情链接
  • 分类
  • 归档
  • 高等数学

  • 线性代数

  • 概率论与数理统计

  • 具体数学

    • 第一章:递归问题
    • 第二章:和式
      • 记号
      • 和式和递归式的转化
      • 和式的变换法则
        • 扰动法
      • 多重和式
      • 一般性方法
        • 方法0:查找公式
        • 方法1:数学归纳法
        • 方法2:扰动法
        • 方法3:建立成套方法
        • 方法4:用积分替换和式
        • 方法5:展开和放缩
      • 有限微积分
      • 无限和式
    • 第三章:整值函数
  • 数学建模

  • 数学
  • 具体数学
martian148
2024-08-03
目录

第二章:和式

# 记号

求和的符号有两种形式

第一种是确定界限的形式,也叫封闭形式,例如:∑k=1nak\sum\limits_{k=1}^n a_kk=1∑n​ak​

第二种叫做一般形式,就是把一个或者多个条件写在 ∑\sum∑ 符号的下面,例如刚刚的例子可以写成 ∑1≤k≤nak\sum\limits_{1\le k \le n} a_k1≤k≤n∑​ak​

# 和式和递归式的转化

和式和递归式之间可以相互转化,

和式转化成递归式

Sn=∑anS_n=\sum a_nSn​=∑an​ 可以转化为 Sn=Sn−1+an,S1=a1S_n=S_{n-1}+a_n,S_1=a_1Sn​=Sn−1​+an​,S1​=a1​,后面这一项就是一个递归式子

例如:我们这里有一个递归式

R0=αRn=Rn−1+β+γn,n>0 \begin{aligned} R_0=\alpha\\ R_n=R_{n-1}+\beta+\gamma n, n>0 \end{aligned} R0​=αRn​=Rn−1​+β+γn,n>0​

这个可以写成通解的形式:Rn=A(n)α+B(n)β+C(n)γR_n=A(n)\alpha+B(n)\beta+C(n)\gammaRn​=A(n)α+B(n)β+C(n)γ

由成套方法得出 A(n)=1,B(n)=n,C(n)=(n2+n)/2A(n)=1,B(n)=n,C(n)=(n^2+n)/2A(n)=1,B(n)=n,C(n)=(n2+n)/2

我们需要计算 ∑k=0n(a+bk)\sum\limits_{k=0}^n(a+bk)k=0∑n​(a+bk),只需要把 a=α,a=β,b=γa=\alpha, a=\beta, b=\gammaa=α,a=β,b=γ 代入即可

递归式转化成和式

T0=0Tn=2Tn−1+1,n>0 \begin{aligned} T_0&=0\\ T_n&=2T_{n-1}+1, n>0 \end{aligned} T0​Tn​​=0=2Tn−1​+1,n>0​

两边除以 2n2^n2n 得到

T0/20=0Tn/2n=Tn−1/2n−1+1/2n,n>0 \begin{aligned} T_0/2^0&=0\\ T_n/2^n&=T_{n-1}/2^{n-1}+1/2^n, n>0 \end{aligned} T0​/20Tn​/2n​=0=Tn−1​/2n−1+1/2n,n>0​

设 Sn=Tn/2nS_n=T_n/2^nSn​=Tn​/2n

S0=0Sn=Sn−1+1/2n,n>0 \begin{aligned} S_0&=0\\ S_n&=S_{n-1}+1/2^n, n>0 \end{aligned} S0​Sn​​=0=Sn−1​+1/2n,n>0​

由此得到 Sn=∑k=1n12k=1−(12)nS_n=\sum\limits_{k=1}^n \frac{1}{2}^k = 1-(\frac{1}{2})^nSn​=k=1∑n​21​k=1−(21​)n,所以 Tn=2n−1T_n=2^n-1Tn​=2n−1

# 和式的变换法则

经典的三条性质,注意这里的 KKK 是一个集合

∑k∈Kcak=c∑k∈Kak(分配律) \sum_{k \in K} c a_k = c \sum_{k \in K} a_k \quad \text{(分配律)} k∈K∑​cak​=ck∈K∑​ak​(分配律)
∑k∈K(ak+bk)=∑k∈Kak+∑k∈Kbk(结合律) \sum_{k \in K} (a_k + b_k) = \sum_{k \in K} a_k + \sum_{k \in K} b_k \quad \text{(结合律)} k∈K∑​(ak​+bk​)=k∈K∑​ak​+k∈K∑​bk​(结合律)
∑k∈Kak=∑p(k)∈Kap(k)(交换律) \sum_{k \in K} a_k = \sum_{p(k) \in K} a_{p(k)} \quad \text{(交换律)} k∈K∑​ak​=p(k)∈K∑​ap(k)​(交换律)

注意 p(k)p(k)p(k) 是一个双射函数,也可以理解为一个排列

例如:K={1,2,3,4}K=\{1,2,3,4\}K={1,2,3,4} 那么,p(k)={4,2,3,1}p(k)=\{4,2,3,1\}p(k)={4,2,3,1} 就是一个排列

当然 p(k)p(k)p(k) 可以多出几个元素例如上面那个例子 p(k)={4,2,3,1,5,6}p(k)=\{4,2,3,1,5,6\}p(k)={4,2,3,1,5,6} 也是可以的,5,65,65,6 不在 KKK 中

接下来看一个例子:S=∑0≤k≤n(a+bk)S=\sum\limits_{0\le k \le n}(a+bk)S=0≤k≤n∑​(a+bk)

利用交换率,用 n−kn-kn−k 代替 kkk,得到

S=∑0≤n−k≤n(a+b(n−k))=∑0≤k≤n(a+bn−bk) S=\sum\limits_{0\le n-k\le n}(a+b(n-k))=\sum\limits_{0\le k \le n}(a+bn-bk) S=0≤n−k≤n∑​(a+b(n−k))=0≤k≤n∑​(a+bn−bk)

利用结合律,把这两个方程相加,得到

2S=∑0≤k≤n(a+bk+a+bn−bk)=∑0≤k≤n(2a+bn)=(n+1)(2a+bn) 2S=\sum\limits_{0\le k \le n}(a+bk+a+bn-bk)=\sum\limits_{0\le k \le n}(2a+bn)=(n+1)(2a+bn) 2S=0≤k≤n∑​(a+bk+a+bn−bk)=0≤k≤n∑​(2a+bn)=(n+1)(2a+bn)

所以

S=(n+1)(2a+bn)2 S=\frac{(n+1)(2a+bn)}{2} S=2(n+1)(2a+bn)​

# 扰动法

和式中有一种非常厉害的方法叫扰动法,就是把一项从和式中去除出去,然后尝试把剩下的项变成 SnS_nSn​ 的形式,从而解出 SnS_nSn​

例如:Sn=∑0≤k≤naxkS_n=\sum\limits_{0\le k \le n}ax^kSn​=0≤k≤n∑​axk,使用扰动法的基本套路

Sn+axn+1=ax0+∑1≤k≤n+1axk=ax0+∑0≤k≤naxk+1=ax0+xSn S_n+ax^{n+1}=ax^0+\sum\limits_{1\le k\le n+1} ax^{k}=ax^0+\sum\limits_{0\le k\le n}ax^{k+1}=ax^0+xS_n Sn​+axn+1=ax0+1≤k≤n+1∑​axk=ax0+0≤k≤n∑​axk+1=ax0+xSn​

等式两边同时出现了 SnS_nSn​,解出 Sn=a−axn+11−xS_n=\frac{a-ax^{n+1}}{1-x}Sn​=1−xa−axn+1​

考虑另外一个例子:Sn=∑0≤k≤nk2kS_n=\sum\limits_{0\le k \le n}k2^kSn​=0≤k≤n∑​k2k,考虑扰动法

Sn+(n+1)2n+1=∑0≤k≤n(k+1)2k+1=2∑0≤k≤nk2k+∑0≤k≤n2k+1=2Sn+2n+2−2 \begin{aligned} S_n+(n+1)2^{n+1}&=\sum\limits_{0\le k \le n}(k+1)2^{k+1}\\ &=2\sum\limits_{0\le k\le n}k2^k+\sum\limits_{0\le k\le n}2^{k+1}\\ &=2S_n+2^{n+2}-2 \end{aligned} Sn​+(n+1)2n+1​=0≤k≤n∑​(k+1)2k+1=20≤k≤n∑​k2k+0≤k≤n∑​2k+1=2Sn​+2n+2−2​

解出 Sn=2n+1(n−1)+2S_n=2^{n+1}(n-1)+2Sn​=2n+1(n−1)+2

设未知数 xxx,能求出 Sn=∑0≤k≤nkxkS_n=\sum\limits_{0\le k \le n} kx^{k}Sn​=0≤k≤n∑​kxk 的通解:

Sn=x−xn+1(n+1)+nxn+2(x−1)2,x≠1 S_n=\frac{x-x^{n+1}(n+1)+nx^{n+2}}{(x-1)^2},x\ne 1 Sn​=(x−1)2x−xn+1(n+1)+nxn+2​,x​=1

我们也可以利用求导的方法,求出 Sn=∑0≤k≤nkxkS_n=\sum\limits_{0\le k \le n} kx^{k}Sn​=0≤k≤n∑​kxk 的通解:

已知:∑k=0nxk=1−xn+11−x\sum\limits_{k=0}^n x^k=\frac{1-x^{n+1}}{1-x}k=0∑n​xk=1−x1−xn+1​,两边同时求导

∑k=0nkxk−1=(1−x)(−(n+1)xn)+1−xn+1(1−x)2=1−(n+1)xn+nxn+1(1−x)2 \sum\limits_{k=0}^n kx^{k-1}=\frac{(1-x)(-(n+1)x^n)+1-x^{n+1}}{(1-x)^2}=\frac{1-(n+1)x^n+nx^{n+1}}{(1-x)^2} k=0∑n​kxk−1=(1−x)2(1−x)(−(n+1)xn)+1−xn+1​=(1−x)21−(n+1)xn+nxn+1​

两边同时乘以 xxx,得到

∑k=0nkxk=x−xn+1(n+1)+nxn+2(x−1)2 \sum\limits_{k=0}^n kx^k=\frac{x-x^{n+1}(n+1)+nx^{n+2}}{(x-1)^2} k=0∑n​kxk=(x−1)2x−xn+1(n+1)+nxn+2​

就能得到和上面扰动法一样的结果了

# 多重和式

多重和式对应于连续函数的积分,可以利用积分的一些思维来思考多重和式

对于自变量相互无关的情况很好理解

∑1≤j,k≤3ajbk=a1b1+a1b2+a1b3+a2b1+a2b2+a2b3+a3b1+a3b2+a3b3=a1(b1+b2+b3)+a2(b1+b2+b3)+a3(b1+b2+b3)=(a1+a2+a3)(b1+b2+b3)=(∑1≤j≤3aj)(∑1≤k≤3bk) \begin{aligned} \sum\limits_{1\le j,k\le 3}a_jb_k&=a_1b_1+a_1b_2+a_1b_3+a_2b_1+a_2b_2+a_2b_3+a_3b_1+a_3b_2+a_3b_3\\ &=a_1(b_1+b_2+b_3)+a_2(b_1+b_2+b_3)+a_3(b_1+b_2+b_3)\\ &=(a_1+a_2+a_3)(b_1+b_2+b_3)\\ &=\left(\sum\limits_{1\le j\le 3}a_j\right)\left(\sum\limits_{1\le k\le 3}b_k\right) \end{aligned} 1≤j,k≤3∑​aj​bk​​=a1​b1​+a1​b2​+a1​b3​+a2​b1​+a2​b2​+a2​b3​+a3​b1​+a3​b2​+a3​b3​=a1​(b1​+b2​+b3​)+a2​(b1​+b2​+b3​)+a3​(b1​+b2​+b3​)=(a1​+a2​+a3​)(b1​+b2​+b3​)=(1≤j≤3∑​aj​)(1≤k≤3∑​bk​)​

一个多重和式可以转化成两个和式相乘的形式,前提是 J=KJ=KJ=K,其中 j∈J,k∈Kj\in J,k\in Kj∈J,k∈K,这一点在级数相乘的时候极为有效

在具体求多重和式的时候,往往是一层一层往外求,理论上来说,先固定住哪个自变量都是无所谓的,所以我们可以交换求和符号的位置

设 P(j,k)P(j,k)P(j,k) 返回的是 j,kj,kj,k 是否满足某种性质,那么有:

∑j∈J∑k∈Kaj,k[P(j,k)]=∑k∈K∑j∈Jaj,k[P(j,k)] \sum\limits_{j\in J}\sum\limits_{k\in K}a_{j,k}[P(j,k)]=\sum\limits_{k\in K}\sum\limits_{j\in J}a_{j,k}[P(j,k)] j∈J∑​k∈K∑​aj,k​[P(j,k)]=k∈K∑​j∈J∑​aj,k​[P(j,k)]

讨论自变量之间有某些限制

一种比较常见的方式:1≤j≤k≤n1\le j \le k \le n1≤j≤k≤n,我们可以通过画图找出需要求和的元素

image.png|512

可以看到对角线以及右上角的那一块都是需要求和的

另外一种常见的方式 1≤j,k≤n,j+k=n1\le j,k\le n,j+k=n1≤j,k≤n,j+k=n,画出图来是按对角线求和

image.png

这种求和方法有一种自己的名字,叫卷积(Convolution)

例一:求 ∑1≤j≤k≤najak\sum\limits_{1\le j\le k\le n}a_ja_k1≤j≤k≤n∑​aj​ak​

image.png

可以看到要求部分是右上三角形

由于图标的对称性,可以得到 S左下=S右上S_{\text{左下}}=S_{\text{右上}}S左下​=S右上​

由容斥原理:S左下+S右上=S全部−S对角线S_{\text{左下}}+S_{\text{右上}}=S_{\text{全部}}-S_{\text{对角线}}S左下​+S右上​=S全部​−S对角线​

于是就能得我们想要的答案了

S右上=∑1≤j≤k≤najak=12((∑k=1nak)2+∑k=1nak2) S_{\text{右上}}=\sum\limits_{1\le j\le k\le n}a_ja_k=\frac{1}{2}\left(\left(\sum\limits_{k=1}^na_k\right)^2+\sum\limits_{k=1}^n a_k^2\right) S右上​=1≤j≤k≤n∑​aj​ak​=21​⎝⎜⎛​(k=1∑n​ak​)2+k=1∑n​ak2​⎠⎟⎞​

例二:求 S=∑1≤j<k≤n(ak−aj)(bk−bj)S=\sum\limits_{1\le j < k \le n}(a_k-a_j)(b_k-b_j)S=1≤j<k≤n∑​(ak​−aj​)(bk​−bj​)

交换 j,kj,kj,k 仍然由对称性:

S=∑1≤k<j≤n(aj−ak)(bj−bk)=∑1≤k<j≤n(ak−aj)(bk−bj) S=\sum\limits_{1\le k< j\le n}(a_j-a_k)(b_j-b_k)=\sum\limits_{1\le k< j\le n}(a_k-a_j)(b_k-b_j) S=1≤k<j≤n∑​(aj​−ak​)(bj​−bk​)=1≤k<j≤n∑​(ak​−aj​)(bk​−bj​)

所以得到答案

2S=∑1≤j,k≤n(aj−ak)(bj−bk)−∑1≤j=k≤n(aj−ak)(bj−bk) 2S=\sum\limits_{1\le j, k\le n}(a_j-a_k)(b_j-b_k)-\sum\limits_{1\le j=k\le n}(a_j-a_k)(b_j-b_k) 2S=1≤j,k≤n∑​(aj​−ak​)(bj​−bk​)−1≤j=k≤n∑​(aj​−ak​)(bj​−bk​)

显然,第二个和式等于 000,把第一个和式展开

2S=∑1≤j,k≤najbj−∑1≤j,k≤najbk−∑1≤j,k≤nakbj+∑1≤j,k≤nakbk=2n∑1≤k≤nakbk−2(∑k=1nak)(∑k=1nbk) \begin{aligned} 2S&=\sum\limits_{1\le j, k\le n}a_jb_j-\sum\limits_{1\le j, k\le n}a_jb_k-\sum\limits_{1\le j, k\le n}a_kb_j+\sum\limits_{1\le j, k\le n}a_kb_k\\ &=2n\sum\limits_{1\le k\le n}a_kb_k-2\left(\sum\limits_{k=1}^n a_k\right)\left(\sum\limits_{k=1}^n b_k\right) \end{aligned} 2S​=1≤j,k≤n∑​aj​bj​−1≤j,k≤n∑​aj​bk​−1≤j,k≤n∑​ak​bj​+1≤j,k≤n∑​ak​bk​=2n1≤k≤n∑​ak​bk​−2(k=1∑n​ak​)(k=1∑n​bk​)​

所以就得出了 S=n∑1≤k≤nakbk−(∑k=1nak)(∑k=1nbk)S=n\sum\limits_{1\le k\le n}a_kb_k-\left(\sum\limits_{k=1}^n a_k\right)\left(\sum\limits_{k=1}^n b_k\right)S=n1≤k≤n∑​ak​bk​−(k=1∑n​ak​)(k=1∑n​bk​)

  • 若 {a},{b}\{a\},\{b\}{a},{b} 都是升序的,那么显然 S≥0S\ge 0S≥0,所以 n∑1≤k≤nakbk−(∑k=1nak)(∑k=1nbk)≥nn\sum\limits_{1\le k\le n}a_kb_k-\left(\sum\limits_{k=1}^n a_k\right)\left(\sum\limits_{k=1}^n b_k\right) \ge nn1≤k≤n∑​ak​bk​−(k=1∑n​ak​)(k=1∑n​bk​)≥n ,可以推出不等式:n∑1≤k≤nakbk≥(∑k=1nak)(∑k=1nbk)n\sum\limits_{1\le k\le n}a_kb_k\ge \left(\sum\limits_{k=1}^n a_k\right)\left(\sum\limits_{k=1}^n b_k\right)n1≤k≤n∑​ak​bk​≥(k=1∑n​ak​)(k=1∑n​bk​)

  • 若 {a}\{a\}{a} 是升序的,{b}\{b\}{b} 是降序的,那么显然 S≤0S\le 0S≤0,所以 n∑1≤k≤nakbk−(∑k=1nak)(∑k=1nbk)≤nn\sum\limits_{1\le k\le n}a_kb_k-\left(\sum\limits_{k=1}^n a_k\right)\left(\sum\limits_{k=1}^n b_k\right) \le nn1≤k≤n∑​ak​bk​−(k=1∑n​ak​)(k=1∑n​bk​)≤n ,可以推出不等式:n∑1≤k≤nakbk≤(∑k=1nak)(∑k=1nbk)n\sum\limits_{1\le k\le n}a_kb_k\le \left(\sum\limits_{k=1}^n a_k\right)\left(\sum\limits_{k=1}^n b_k\right)n1≤k≤n∑​ak​bk​≤(k=1∑n​ak​)(k=1∑n​bk​)

这两个不等式叫做 切比雪夫单调不等式

例三

Sn=∑1≤j<k≤n1k−j S_n=\sum\limits_{1\le j<k\le n}\frac{1}{k-j} Sn​=1≤j<k≤n∑​k−j1​

我们可以用 k−jk-jk−j 替换 jjj,然后固定 kkk

Sn=∑1≤k≤n∑1≤j≤k1k−j=∑1≤k≤n∑0≤j≤k−11j=∑1≤k≤nHk−1=∑0≤k≤n−1Hk \begin{aligned} S_n&=\sum\limits_{1\le k \le n}\sum\limits_{1\le j\le k}\frac{1}{k-j}\\ &=\sum\limits_{1\le k \le n}\sum\limits_{0\le j\le k-1}\frac{1}{j}\\ &=\sum\limits_{1\le k \le n}H_{k-1}\\ &=\sum\limits_{0\le k \le n-1}H_{k}\\ \end{aligned} Sn​​=1≤k≤n∑​1≤j≤k∑​k−j1​=1≤k≤n∑​0≤j≤k−1∑​j1​=1≤k≤n∑​Hk−1​=0≤k≤n−1∑​Hk​​

其中 HkH_kHk​ 是调和级数,我们对调和级数有公式,但是对于调和级数求和没有公式,说明固定 kkk 走不通考虑固定 jjj

Sn=∑1≤j<n∑j<k≤n1k−j=∑1≤j≤n∑j<k+j≤n1k+j−j=∑1≤j≤n∑0<k≤n−j1k=∑1≤j≤nHn−j=∑0≤j≤n−1Hj \begin{aligned} S_n&=\sum\limits_{1\le j < n}\ \sum\limits_{j < k\le n}\frac{1}{k-j}\\ &=\sum\limits_{1\le j \le n}\ \sum\limits_{j<k+j\le n}\frac{1}{k+j-j}\\ &=\sum\limits_{1\le j \le n}\ \sum\limits_{0<k\le n-j}\frac{1}{k}\\ &=\sum\limits_{1\le j \le n}H_{n-j}\\ &=\sum\limits_{0\le j \le n-1}H_{j}\\ \end{aligned} Sn​​=1≤j<n∑​ j<k≤n∑​k−j1​=1≤j≤n∑​ j<k+j≤n∑​k+j−j1​=1≤j≤n∑​ 0<k≤n−j∑​k1​=1≤j≤n∑​Hn−j​=0≤j≤n−1∑​Hj​​

还是回到了调和级数求和的问题,考虑按对角线求和

Sn=∑1≤j<n∑j<k≤n1k−j=∑1≤j≤n∑j<k+j≤n1k+j−j=∑1≤j≤n∑0<k≤n−j1k=∑0<k<n∑1≤j≤n−k1k=∑0<k<nn−kk=∑0<k<n(nk−1)=n∑1≤k≤n1k−nn−(n−1)=nHn−n \begin{aligned} S_n&=\sum\limits_{1\le j < n}\ \sum\limits_{j < k\le n}\frac{1}{k-j}\\ &=\sum\limits_{1\le j \le n}\ \sum\limits_{j<k+j\le n}\frac{1}{k+j-j}\\ &=\sum\limits_{1\le j \le n}\ \sum\limits_{0<k\le n-j}\frac{1}{k}\\ &=\sum\limits_{0<k<n}\ \sum\limits_{1\le j\le n-k}\frac{1}{k}\\ &=\sum\limits_{0<k<n}\frac{n-k}{k}=\sum\limits_{0<k<n}(\frac{n}{k}-1)\\ &=n\sum\limits_{1\le k\le n}\frac{1}{k}-\frac{n}{n}-(n-1)\\ &=nH_n-n \end{aligned} Sn​​=1≤j<n∑​ j<k≤n∑​k−j1​=1≤j≤n∑​ j<k+j≤n∑​k+j−j1​=1≤j≤n∑​ 0<k≤n−j∑​k1​=0<k<n∑​ 1≤j≤n−k∑​k1​=0<k<n∑​kn−k​=0<k<n∑​(kn​−1)=n1≤k≤n∑​k1​−nn​−(n−1)=nHn​−n​

我们还可以从几何的层面来理解这个式子:

image.png|421

刚开始的几次我们按照行和列求和,都会得到一个调和级数求和的式子,最后一种方式是按照对角线,这里得到 31+22+13\frac{3}{1}+\frac{2}{2}+\frac{1}{3}13​+22​+31​

# 一般性方法

这一节基本上就是总结性质的,把前面介绍的方法综合运用一下

要求一个立方和 Sn=∑k=0nk2,n≥0S_n=\sum\limits_{k=0}^n k^2, \ n \ge 0Sn​=k=0∑n​k2, n≥0

# 方法0:查找公式

OEIS (opens new window) 能帮我们很好的查找公式

得到的答案是 Sn=n(n+1)(2n+1)6S_n=\frac{n(n+1)(2n+1)}{6}Sn​=6n(n+1)(2n+1)​

# 方法1:数学归纳法

已知:Sn=n(n+12)(n+1)3S_n=\frac{n(n+\frac{1}{2})(n+1)}{3}Sn​=3n(n+21​)(n+1)​

显然 S0=0S_0=0S0​=0 成立,假设 n>0n>0n>0 时,Sn−1=(n−1)(n−12)n3S_{n-1}=\frac{(n-1)(n-\frac{1}{2})n}{3}Sn−1​=3(n−1)(n−21​)n​ 成立,有:Sn=Sn−1+SS_n=S_{n-1}+SSn​=Sn−1​+S

3Sn=(n−1)(n−12)(n)+3n2=n3+32n2+12n=n(n+1)(n+12) \begin{aligned} 3S_n&=(n-1)(n-\frac{1}{2})(n)+3n^2 \\ &=n^3+\frac{3}{2}n^2+\frac{1}{2}n \\ &=n(n+1)(n+\frac{1}{2}) \end{aligned} 3Sn​​=(n−1)(n−21​)(n)+3n2=n3+23​n2+21​n=n(n+1)(n+21​)​

所以 Sn=n(n+12)(n+1)3S_n=\frac{n(n+\frac{1}{2})(n+1)}{3}Sn​=3n(n+21​)(n+1)​ 对所有 n≥0n\ge 0n≥0 都成立

# 方法2:扰动法

Sn+(n+1)2=∑k=0n(k+1)2=∑k=0n(k2+2k+1)=∑k=0nk2+2∑k=0nk+∑k=0n1=Sn+2∑k=0nk+n+1 \begin{aligned} S_n+(n+1)^2&=\sum\limits_{k=0}^{n}(k+1)^2=\sum\limits_{k=0}^n(k^2+2k+1)\\ &=\sum\limits_{k=0}^{n}k^2+2\sum\limits_{k=0}^nk+\sum\limits_{k=0}^n1\\ &=S_n+2\sum\limits_{k=0}^nk+n+1 \end{aligned} Sn​+(n+1)2​=k=0∑n​(k+1)2=k=0∑n​(k2+2k+1)=k=0∑n​k2+2k=0∑n​k+k=0∑n​1=Sn​+2k=0∑n​k+n+1​

惊人的发现,我们需要求的 SnS_nSn​ 被消掉了,留下一个 ∑k=0nk=(n+1)2−(n+1)2=n(n+1)2\sum\limits_{k=0}^nk=\frac{(n+1)^2-(n+1)}{2}=\frac{n(n+1)}{2}k=0∑n​k=2(n+1)2−(n+1)​=2n(n+1)​

我们猜想能不能用立方和来把平方和求出来,设 Tn=∑k=0nk3T_n=\sum\limits_{k=0}^n k^3Tn​=k=0∑n​k3

Tn+(n+1)3=∑k=0n(k+1)3=∑k=0n(k3+3k2+3k+1)=∑k=0nk3+3∑k=0nk2+3∑k=0nk+∑k=0n1=Tn+3Sn+3n(n+1)2+n+1 \begin{aligned} T_n+(n+1)^3&=\sum\limits_{k=0}^n(k+1)^3=\sum\limits_{k=0}^n(k^3+3k^2+3k+1)\\ &=\sum\limits_{k=0}^n k^3+3\sum\limits_{k=0}^n k^2+3\sum\limits_{k=0}^n k+\sum\limits_{k=0}^n 1\\ &=T_n+3S_n+3\frac{n(n+1)}{2}+n+1 \end{aligned} Tn​+(n+1)3​=k=0∑n​(k+1)3=k=0∑n​(k3+3k2+3k+1)=k=0∑n​k3+3k=0∑n​k2+3k=0∑n​k+k=0∑n​1=Tn​+3Sn​+32n(n+1)​+n+1​

TnT_nTn​ 被消掉了,留下我们需要求的 SnS_nSn​,这样 3Sn=(n+1)3−3n(n+1)2−(n+1)3S_n=(n+1)^3-3\frac{n(n+1)}{2}-(n+1)3Sn​=(n+1)3−32n(n+1)​−(n+1)

所以 Sn=n(n+1)(n+12)3S_n=\frac{n(n+1)(n+\frac{1}{2})}{3}Sn​=3n(n+1)(n+21​)​

# 方法3:建立成套方法

还是老规矩,建立成套方法,设:

R0=α,n=0Rn=Rn−1+β+γn+δn2,n>0 \begin{aligned} R_0&=\alpha, n=0\\ R_n&=R_{n-1}+\beta+\gamma n+\delta n^2, n>0 \end{aligned} R0​Rn​​=α,n=0=Rn−1​+β+γn+δn2,n>0​

解的一般形式就是:Rn=A(n)α+B(n)β+C(n)γ+D(n)δR_n=A(n)\alpha+B(n)\beta+C(n)\gamma+D(n)\deltaRn​=A(n)α+B(n)β+C(n)γ+D(n)δ

我们已知,δ=0\delta=0δ=0 时,A(n)=1,B(n)=n,C(n)=(n2+n)/2A(n)=1,B(n)=n,C(n)=(n^2+n)/2A(n)=1,B(n)=n,C(n)=(n2+n)/2

现在要求 D(n)D(n)D(n) 和他们的关系,设 Rn=n3R_n=n^3Rn​=n3,有 n3=(n−1)3+β+γn+δn2n^3=(n-1)^3+\beta+\gamma n+\delta n^2n3=(n−1)3+β+γn+δn2

可得:α=0\alpha=0α=0,β=1\beta=1β=1,γ=−3\gamma=-3γ=−3,δ=3\delta=3δ=3,代入:3D(n)−3C(n)+B(n)=n33D(n)-3C(n)+B(n)=n^33D(n)−3C(n)+B(n)=n3

就把 D(n)D(n)D(n) 接出来了,D(n)=n3−3n(n+1)2+n=n(n+12)(n+1)3D(n)=n^3-3\frac{n(n+1)}{2}+n=\frac{n(n+\frac{1}{2})(n+1)}{3}D(n)=n3−32n(n+1)​+n=3n(n+21​)(n+1)​

我们需要求的和式 SnS_nSn​ 转化成递归式子,只需要令 α=β=γ=0\alpha=\beta=\gamma=0α=β=γ=0 以及 δ=1\delta=1δ=1

则 Sn=D(n)S_n=D(n)Sn​=D(n)

# 方法4:用积分替换和式

我们已知 ∫0nx2dx=n33\int_{0}^{n}x^2dx=\frac{n^3}{3}∫0n​x2dx=3n3​,显然,和式和积分之间的差距并不差,只是差一个误差

我们设这个误差为 En=Sn−n33E_n=S_n-\frac{n^3}{3}En​=Sn​−3n3​,通过 SnS_nSn​ 的递归式,我们能得出 Sn=Sn−1+n2S_n=S_{n-1}+n^2Sn​=Sn−1​+n2

En=Sn−n33=Sn−1+n2−n33=En−1+(n−1)33+n2−n33=En−1+n−13 \begin{aligned} E_n&=S_n-\frac{n^3}{3}=S_{n-1}+n^2-\frac{n^3}{3}\\ &=E_{n-1}+\frac{(n-1)^3}{3}+n^2-\frac{n^3}{3}\\ &=E_{n-1}+n-\frac{1}{3} \end{aligned} En​​=Sn​−3n3​=Sn−1​+n2−3n3​=En−1​+3(n−1)3​+n2−3n3​=En−1​+n−31​​

于是,我们就得到了 EnE_nEn​ 的递推式,En=En−1+n−13E_n=E_{n-1}+n-\frac{1}{3}En​=En−1​+n−31​,E0=0E_0=0E0​=0

很容易算出 EnE_nEn​ 的通项公式为 En=n(n+1)2−n3E_n=\frac{n(n+1)}{2}-\frac{n}{3}En​=2n(n+1)​−3n​,所以 Sn=En+n33=n(n+1)(n+12)3S_n=E_n+\frac{n^3}{3}=\frac{n(n+1)(n+\frac{1}{2})}{3}Sn​=En​+3n3​=3n(n+1)(n+21​)​

# 方法5:展开和放缩

这种方法极具技巧性,把一个一重和式转化成二重和式

Sn=∑k=0nk2=∑1≤j≤k≤nk=∑1≤j≤n∑j≤k≤nk=∑1≤j≤n(j+n2)(n−j+1)=12∑1≤j≤n(n(n+1)+j−j2)=12n2(n+1)+14n(n+1)−12Sn \begin{aligned} S_n&=\sum\limits_{k=0}^n k^2=\sum\limits_{1\le j \le k \le n} k \\ &=\sum\limits_{1\le j \le n}\ \sum\limits_{j\le k\le n} k\\ &=\sum\limits_{1\le j \le n} (\frac{j+n}{2})(n-j+1)\\ &=\frac{1}{2}\sum\limits_{1\le j \le n} \left(n(n+1)+j-j^2\right)\\ &=\frac{1}{2}n^2(n+1)+\frac{1}{4}n(n+1)-\frac{1}{2} S_n\\ \end{aligned} Sn​​=k=0∑n​k2=1≤j≤k≤n∑​k=1≤j≤n∑​ j≤k≤n∑​k=1≤j≤n∑​(2j+n​)(n−j+1)=21​1≤j≤n∑​(n(n+1)+j−j2)=21​n2(n+1)+41​n(n+1)−21​Sn​​

这样等式两边都出现了 SnS_nSn​,解出 Sn=n(n+1)(2n+1)6S_n=\frac{n(n+1)(2n+1)}{6}Sn​=6n(n+1)(2n+1)​

# 有限微积分

我们类似于积分和微分的思维,在离散域上定义有限微积分,定义差分算子 Δf(x)\Delta f(x)Δf(x) 是 (f(x+h)−f(x))/h(f(x+h)-f(x))/h(f(x+h)−f(x))/h 当 h=1h=1h=1 时的值

还需要定义两种特殊的次幂:

  • 下降阶乘幂 xm‾=x(x−1)⋯(x−m+1),m≥0x^{\underline{m}}=x(x-1)\cdots(x-m+1),\ m \ge 0xm​=x(x−1)⋯(x−m+1), m≥0
  • 上升阶乘幂 xm‾=x(x+1)⋯(x+m−1),m≥0x_{\overline{m}}=x(x+1)\cdots(x+m-1),\ m \ge 0xm​=x(x+1)⋯(x+m−1), m≥0

我们为什么要定义这样奇特的形式呢,因为这样求差分的时候会和微分有很多共同点

Δxm‾=(x+1)m‾−xm‾=(x+1)x(x−1)⋯(x−m+2)−x(x−1)(x−2)⋯(x−m+1)=(x+1−(x−m+1))x(x−1)(x−2)⋯(x−m+2)=mx(x−1)(x−2)⋯(x−m+2)=mxm−1‾ \begin{aligned} \Delta x^{\underline{m}}&=(x+1)^{\underline{m}}-x^{\underline{m}}\\ &=(x+1)x(x-1)\cdots(x-m+2)-x(x-1)(x-2)\cdots(x-m+1)\\ &=(x+1-(x-m+1))x(x-1)(x-2)\cdots(x-m+2)\\ &=mx(x-1)(x-2)\cdots(x-m+2)\\ &=m x^{\underline{m-1}} \end{aligned} Δxm​​=(x+1)m​−xm​=(x+1)x(x−1)⋯(x−m+2)−x(x−1)(x−2)⋯(x−m+1)=(x+1−(x−m+1))x(x−1)(x−2)⋯(x−m+2)=mx(x−1)(x−2)⋯(x−m+2)=mxm−1​​

有限微积分也存在类似于微积分基本定理的东西

g(x)=Δf(x)⟺∑g(x)δx=f(x)+C g(x)=\Delta f(x)\Longleftrightarrow \sum g(x)\delta x=f(x)+C g(x)=Δf(x)⟺∑g(x)δx=f(x)+C

这里的 ∑g(x)δx\sum g(x)\delta x∑g(x)δx 是 g(x)g(x)g(x) 的不定和式,这里的 CCC 可以是周期为 111 的任意函数

大名鼎鼎的牛顿-莱布尼茨公式在有限微积分中也有体现 ∑abg(x)δx=f(x)∣ab=f(b)−f(a)\sum_a^b g(x) \delta x = f(x) \bigg|_a^b = f(b) - f(a)∑ab​g(x)δx=f(x)∣∣∣∣∣​ab​=f(b)−f(a)

思考 ∑abg(x)δx=f(x)∣ab=f(b)−f(a)\sum_a^b g(x) \delta x = f(x) \bigg|_a^b = f(b) - f(a)∑ab​g(x)δx=f(x)∣∣∣∣∣​ab​=f(b)−f(a) 当,上下限相同时:

∑aag(x)δx=f(x)∣aa=f(a)−f(a)=0 \sum_a^a g(x) \delta x = f(x) \bigg|_a^a = f(a) - f(a) = 0 a∑a​g(x)δx=f(x)∣∣∣∣∣​aa​=f(a)−f(a)=0

当上下限差 111 时:

∑aa+1g(x)δx=f(x)∣aa+1=f(a+1)−f(a)=Δf(x)=g(x) \sum_a^{a+1} g(x) \delta x = f(x) \bigg|_a^{a+1} = f(a+1) - f(a)=\Delta f(x)=g(x) a∑a+1​g(x)δx=f(x)∣∣∣∣∣​aa+1​=f(a+1)−f(a)=Δf(x)=g(x)

着预示着和传统的求和符号有点不同,可以理解为 ∑abg(x)δx=∑k=ab−1g(k)=∑a≤k<bg(k),b≥a\sum_{a}^b g(x)\delta x=\sum\limits_{k=a}^{b-1}g(k)=\sum\limits_{a\le k<b}g(k), b\ge a∑ab​g(x)δx=k=a∑b−1​g(k)=a≤k<b∑​g(k),b≥a

也就是说,如果存在能找到一个 f(x)f(x)f(x) 使得 f(x+1)−f(x)=g(x)f(x+1)-f(x)=g(x)f(x+1)−f(x)=g(x),那么理论上来说就可以求 ∑g(x)\sum g(x)∑g(x) 这类的和了

我们尝试去求类似于原函数的东西 f(x)f(x)f(x)

比如:g(x)=xm‾g(x)=x^{\underline{m}}g(x)=xm​,f(x)=xm+1‾m+1f(x)=\frac{x^{\underline{m+1}}}{m+1}f(x)=m+1xm+1​​

我们用前面的例子,求 ∑0≤k<nk2\sum\limits_{0\le k < n} k^20≤k<n∑​k2 已知 k2=k2‾+k1‾k^2=k^{\underline{2}}+k^{\underline{1}}k2=k2​+k1​,所以

∑0≤k<nk2=n3‾3+n2‾2=n(n−1)(n−2+32)=13n(n−12)(n−1) \sum\limits_{0\le k < n} k^2=\frac{n^{\underline{3}}}{3}+\frac{n^{\underline{2}}}{2}=n(n-1)\left(n-2+\frac{3}{2}\right)=\frac{1}{3}n\left(n-\frac{1}{2}\right)(n-1) 0≤k<n∑​k2=3n3​​+2n2​​=n(n−1)(n−2+23​)=31​n(n−21​)(n−1)

利用 n+1n+1n+1 替换 nnn 就可以得到 SnS_nSn​ 了

尝试继续推广无限微积分,考虑下降次幂为负数的情况,xm‾,m<0x^{\underline{m}},\ m<0xm​, m<0

观察:

  • x3‾=x(x−1)(x−2)x^{\underline{3}}=x(x-1)(x-2)x3​=x(x−1)(x−2),x2‾=x(x−1)x^{\underline{2}}=x(x-1)x2​=x(x−1) 两个相除得出 x−2x-2x−2
  • x2‾=x(x−1)x^{\underline{2}}=x(x-1)x2​=x(x−1),x1‾=xx^{\underline{1}}=xx1​=x 两个相除得出 x−1x-1x−1
  • x1‾=xx^{\underline{1}}=xx1​=x,x0‾=1x^{\underline{0}}=1x0​=1 两个相除得出 xxx

于是找规律得出 x−1‾=1x+1,x−2‾=1(x+1)(x+2)x^{\underline{-1}}=\frac{1}{x+1},x^{\underline{-2}}=\frac{1}{(x+1)(x+2)}x−1​=x+11​,x−2​=(x+1)(x+2)1​,x−m=1(x+1)(x+2)⋯(x+m)x^{-m}=\frac{1}{(x+1)(x+2)\cdots(x+m)}x−m=(x+1)(x+2)⋯(x+m)1​

通常的幂法则:xm+n=xmxnx^{m+n}=x^mx^nxm+n=xmxn,推广到有限微积分就变成了 xm+n‾=xm‾(x−m)n‾x^{\underline{m+n}}=x^{\underline{m}}(x-m)^{\underline{n}}xm+n​=xm​(x−m)n​

现在确认下降幂的差分性质在 m<0m<0m<0 是否成立也就是 Δxm‾=mxm−1‾\Delta x^{\underline{m}}=mx^{\underline{m-1}}Δxm​=mxm−1​

如果 m=−2m=-2m=−2,有

Δx−2‾=(x+1)−2‾−x−2‾=1(x+2)(x+3)−1(x+1)(x+2)=(x+1)−(x+3)(x+1)(x+2)(x+3)=−2(x+1)(x+2)(x+3)=−2x−3‾ \begin{aligned} \Delta x^{\underline{-2}}&=(x+1)^{\underline{-2}}-x^{\underline{-2}}\\ &=\frac{1}{(x+2)(x+3)}-\frac{1}{(x+1)(x+2)}\\ &=\frac{(x+1)-(x+3)}{(x+1)(x+2)(x+3)}\\ &=\frac{-2}{(x+1)(x+2)(x+3)}\\ &=-2x^{\underline{-3}} \end{aligned} Δx−2​​=(x+1)−2​−x−2​=(x+2)(x+3)1​−(x+1)(x+2)1​=(x+1)(x+2)(x+3)(x+1)−(x+3)​=(x+1)(x+2)(x+3)−2​=−2x−3​​

说明求和性质对 m<0m<0m<0 也成立,即 ∑abxm‾δx=xm+1‾m+1∣ab,m≠1\sum_a^b x^{\underline{m}}\delta x=\frac{x^{\underline{m+1}}}{m+1}\bigg|_a^b,\ m \ne 1∑ab​xm​δx=m+1xm+1​​∣∣∣∣∣​ab​, m​=1

现在考虑 m=1m=1m=1 的时候,对于微积分,我们又 ∫abx−1dx=ln⁡x∣ab\int_a^b x^{-1}dx=\ln x\bigg|_{a}^b∫ab​x−1dx=lnx∣∣∣∣∣​ab​

我们需要在有限微积分中找一个类似于 ln⁡x\ln xlnx 的函数 ,这个函数的差分为 1x+1\frac{1}{x+1}x+11​

很显然能得到这个函数就是 HxH_xHx​,于是就得到了求和的完整形式

∑abx−mδx={xm+1m+1∣ab,m≠−1,Hx∣ab,m=−1. \sum_{a}^{b} x^{-m} \delta x = \begin{cases} \left. \frac{x^{m+1}}{m+1} \right|_{a}^{b}, & m \neq -1, \\ \left. H_{x} \right|_{a}^{b}, & m = -1. \end{cases} a∑b​x−mδx=⎩⎪⎨⎪⎧​m+1xm+1​∣∣∣∣​ab​,Hx​∣ab​,​m​=−1,m=−1.​

类似的,我们需要找一个 exe^xex 类似物,根据定义 d(ex)=exd(e^x)=e^xd(ex)=ex,所以有 f(x+1)−f(x)=f(x)⟺f(x+1)=2f(x)f(x+1)-f(x)=f(x)\Longleftrightarrow f(x+1)=2f(x)f(x+1)−f(x)=f(x)⟺f(x+1)=2f(x),即 f(x)=2xf(x)=2^xf(x)=2x

cxc^xcx 的差分也相当简单,即对任意的 ccc 有

Δ(cx)=cx+1−cx=(c−1)cx. \Delta (c^x) = c^{x+1} - c^x = (c - 1) c^x. Δ(cx)=cx+1−cx=(c−1)cx.

那么,c≠1c\ne 1c​=1 时,cxc^xcx 的原函数就是 cx/(c−1)c^x/(c-1)cx/(c−1)

两个相乘的函数求微分:d(uv)=udv+vdud(uv)=udv+vdud(uv)=udv+vdu,两边同时积分得到分部积分法的公式:∫udv=uv−∫vdu\int udv=uv-\int vdu∫udv=uv−∫vdu

有限微积分也类似

Δ(u(x)v(x))=u(x+1)v(x+1)−u(x)v(x)=u(x+1)v(x+1)−u(x)v(x+1)+u(x)v(x+1)−u(x)v(x)=u(x)Δv(x)+v(x+1)Δu(x). \begin{aligned} \Delta \left( u(x)v(x) \right) &= u(x+1)v(x+1) - u(x)v(x) \\ &= u(x+1)v(x+1) - u(x)v(x+1) + u(x)v(x+1) - u(x)v(x) \\ &= u(x)\Delta v(x) + v(x+1)\Delta u(x). \end{aligned} Δ(u(x)v(x))​=u(x+1)v(x+1)−u(x)v(x)=u(x+1)v(x+1)−u(x)v(x+1)+u(x)v(x+1)−u(x)v(x)=u(x)Δv(x)+v(x+1)Δu(x).​

这里的 v(x+1)v(x+1)v(x+1) 看着很烦,所以定义一个位移算子 Ef(x)=f(x+1)\text{E} f(x)=f(x+1)Ef(x)=f(x+1)

所以乘积差分法则可以写为:

Δ(uv)=uΔv+EvΔu\Delta(uv)=u\Delta v+\text{E}v\Delta uΔ(uv)=uΔv+EvΔu

两边求和得:

∑uΔv=uv−∑EvΔu\sum u\Delta v=uv-\sum \text{E}v\Delta u∑uΔv=uv−∑EvΔu

有限微积分和无限微积分对应着离散和连续,有很多共通点,在下列表中给出

image.png|789

来看几个例子:

例一:

∑k=0nk2k=∑0n+1x2xδx=x2x−2x+1∣0n+1=((n+1)2n+1−2n+2)−(0×20−21)=(n−1)2n+1+2 \begin{aligned} \sum_{k=0}^{n} k 2^k &= \sum_{0}^{n+1} x 2^x \delta x \\ &= x 2^x - 2^{x+1} \bigg|_{0}^{n+1}\\ &= \left( (n+1) 2^{n+1} - 2^{n+2} \right) - (0 \times 2^0 - 2^1)\\ &= (n-1) 2^{n+1} + 2 \end{aligned} k=0∑n​k2k​=0∑n+1​x2xδx=x2x−2x+1∣∣∣∣∣​0n+1​=((n+1)2n+1−2n+2)−(0×20−21)=(n−1)2n+1+2​

例二: 求 ∑0≤k<nkHk\sum\limits_{0\le k < n} kH_k0≤k<n∑​kHk​

先求原函数,根据分部积分法则:

∑xHxδx=x2‾2Hx−∑(x+1)2‾xx−1‾δx=x2‾2Hx−12∑x1‾δx=x2‾2Hx−x2‾4+C \begin{aligned} \sum xH_x\delta x&=\frac{x^{\underline{2}}}{2}H_x-\sum\frac{(x+1)^{\underline{2}}}{x}x^{\underline{-1}}\delta x\\ &=\frac{x^{\underline{2}}}{2}H_x-\frac{1}{2}\sum x^{\underline{1}}\delta x\\ &=\frac{x^{\underline{2}}}{2}H_x-\frac{x^{\underline{2}}}{4}+C \end{aligned} ∑xHx​δx​=2x2​​Hx​−∑x(x+1)2​​x−1​δx=2x2​​Hx​−21​∑x1​δx=2x2​​Hx​−4x2​​+C​

然后把上下限代入:

∑0≤k<nkHk=∑0nxHxδx=n2‾2(Hn−12) \sum\limits_{0\le k <n} kH_k=\sum_{0}^{n}xH_x\delta x=\frac{n^{\underline{2}}}{2}\left(H_n-\frac{1}{2}\right) 0≤k<n∑​kHk​=0∑n​xHx​δx=2n2​​(Hn​−21​)

# 无限和式

引用大乌拉的一句话:

我们对无穷的东西几乎一无所知

这一章书上举了几个具体的例子,都在表达无穷和式似乎没有我们想象的这么简单

上次更新: 2025/04/08, 18:03:31
第一章:递归问题
第三章:整值函数

← 第一章:递归问题 第三章:整值函数→

最近更新
01
Java基础语法
05-26
02
开发环境配置
05-26
03
pink 老师 JavaScript 学习笔记
05-26
更多文章>
Theme by Vdoing | Copyright © 2024-2025 Martian148 | MIT License
  • 跟随系统
  • 浅色模式
  • 深色模式
  • 阅读模式