和式【具體數學】

和式和遞歸式

和式很多情況下可以轉化爲遞歸式然後進行求解,比如:Sn=k=0nakS_n=\sum^n_{k=0}a_k,就等價於遞歸式S0=a0,Sn=Sn1+anS_0=a_0,S_n=S_{n-1}+a_n。迴歸上一章求解,雖然不一定簡單,但是將問題轉化爲已經解決的問題,是一定有效的方法。


回一下上一章講的那個“成套方法”。假設現在有遞歸式R0=αRn=Rn1+β+γnR_0=\alpha\\ R_n=R_{n-1}+\beta+\gamma n
我們將解寫成Rn=A(n)α+B(n)β+C(n)γR_n=A(n)\alpha+B(n)\beta+C(n)\gamma
α=1,β=0,γ=0\alpha=1,\beta=0,\gamma=0,得到A(n)=1A(n)=1
α=0,β=1,γ=0\alpha=0,\beta=1,\gamma=0,得到B(n)=nB(n)=n
α=0,β=0,γ=1\alpha=0,\beta=0,\gamma=1,得到C(n)=n2+n2C(n)=\frac{n^2+n}{2}
所以得到:Rn=α+nβ+n2+n2γR_n=\alpha+n\beta+\frac{n^2+n}{2}\gamma


迴歸到和式,如果我們要計算k=0n(a+bk)\sum_{k=0}^n(a+bk),那我們就可以帶入上面的式子,其中α=β=a,γ=b\alpha=\beta=a,\gamma=b。得到Rn=(n+1)a+n2+n2bR_n=(n+1)a+\frac{n^2+n}{2}b


這裏舉的例子是和式轉化爲遞歸式,反過來也是可以的,各有各的優點。比如上一章的河內塔問題。
我們這裏研究一下這個式子:anTn=bnTn1+cna_nT_n=b_nT_{n-1}+c_n
我們引入一個求和因子sns_n,原式變成snanTn=snbnTn1+sncns_na_nT_n=s_nb_nT_{n-1}+s_nc_n
我們之所以要引入求個因子,是想式子化簡爲樸素的遞歸形式。
我們使得sns_n滿足這個snbn=sn1an1s_nb_n=s_{n-1}a_{n-1}條件。
那我們記Sn=snanTnS_n=s_na_nT_n,就可以將SnS_n化簡爲Sn=Sn1+sncnS_n=S_{n-1}+s_nc_n
那我們就得到Sn=s0a0T0+k=1nskck=s1b1T0+k=1nskckS_n=s_0a_0T_0+\sum_{k=1}^ns_kc_k\\=s_1b_1T_0+\sum_{k=1}^ns_kc_k
所以Tn=1snan(s1b1T0+k=1nskck)T_n=\frac{1}{s_na_n}(s_1b_1T_0+\sum_{k=1}^ns_kc_k)
那麼現在這個式子看似解決了,但是ss該是什麼呢?
化簡約束式子得到sn=sn1an1bns_n=s_{n-1}\frac{a_{n-1}}{b_n}
得到sn=an1an2...a1bnbn1...b2ks_n=\frac{a_{n-1}a_{n-2}...a_1}{b_nb_{n-1}...b_2}*k
現在這個方法就完整了,注意這個方法並不是萬能的,是有限制的。


舉個例子,我們研究快排的平均次數。容易得知有下列遞歸式C0=C1=0Cn=n+1+2nk=0n1CkC_0=C_1=0\\C_n=n+1+\frac{2}{n}\sum_{k=0}^{n-1}C_k
這個遞歸式看着是可怕的,相信好多人是沒有嘗試解決過看着如此複雜的問題。
也許我們的工作的第一步目標是清晰的,一般會先消去最可怕的\sum符號。
Cn=n+1+2nk=0n1CknCn=n2+n+2k=0n1Ck(n1)Cn1=(n1)2+(n1)+2k=0n2CknCn(n1)Cn1=2n+2Cn1nCn=(n+1)Cn1+2nC_n=n+1+\frac{2}{n}\sum_{k=0}^{n-1}C_k\\nC_n=n^2+n+2\sum_{k=0}^{n-1}C_k\\(n-1)C_{n-1}=(n-1)^2+(n-1)+2\sum_{k=0}^{n-2}C_k\\nC_n-(n-1)C_{n-1}=2n+2C_{n-1}\\nC_n=(n+1)C_{n-1}+2n
所以,式子化簡爲下面的形式:C0=C1=0C2=3nCn=(n+1)Cn1+2nC_0=C_1=0\\C_2=3\\nC_n=(n+1)C_{n-1}+2n
帶入上面推出的結論,其中:an=nbn=n+1cn=2n2[n=1]+2[n=2]sn=an1an2...a1bnbn1...b1=2n2+na_n=n\\b_n=n+1\\c_n=2n-2[n=1]+2[n=2]\\s_n=\frac{a_{n-1}a_{n-2}...a_1}{b_nb_{n-1}...b_1}=\frac{2}{n^2+n}
所以CnC_n的解爲Cn=1snan(s1b1C0+k=1nskck)=2(n+1)k=1n1k+123(n+1)C_n=\frac{1}{s_na_n}(s_1b_1C_0+\sum_{k=1}^ns_kc_k)\\=2(n+1)\sum_{k=1}^n\frac{1}{k+1}-\frac{2}{3}(n+1)
我們記Hn=k=1n1kH_n=\sum_{k=1}^n\frac{1}{k},則原式等於:Cn=2(n+1)Hn83n23C_n=2(n+1)H_n-\frac{8}{3}n-\frac{2}{3}



和式的處理

我們知道,有些等價的和式,不同的寫法運算量是不一樣的。這裏就需要一些和式的運算法則。
kKcak=ckKak    ()\sum_{k\in K}ca_k=c\sum_{k\in K}a_k\ \ \ \ (分配率)kK(ak+bk)=kKak+kKbk() \sum_{k\in K}(a_k+b_k)=\sum_{k\in K}a_k+\sum_{k\in K}b_k(結合律)p(k)Kak=kKap(k)    () \sum_{p(k)\in K}a_k=\sum_{k\in K}a_{p(k)}\ \ \ \ (交換律)
擾動法,例如Sn=k=0naxkS_n=\sum_{k=0}^nax^k,那麼有Sn+1=Sn+axn+1=a+xSnS_{n+1}=S_n+ax^{n+1}=a+xS_n所以,Sn=aaxn+11xS_n=\frac{a-ax^{n+1}}{1-x},這就是擾動法。


現在,我們來考慮難一點的和式,Sn=k=0nk2kS_n=\sum_{k=0}^nk*2^k
現在,嘗試尋找Sn+1S_{n+1}SnS_n的關係,Sn+1=k=0n+1k2k=Sn+(n+1)2n+1S_{n+1}=\sum_{k=0}^{n+1}k*2^k=S_n+(n+1)*2^{n+1}Sn+1=k=0n+1k2k=k=0n(k+1)2k+1=k=0nk2k+1+k=0n2k+1=2Sn+2n+22S_{n+1}=\sum_{k=0}^{n+1}k*2^k=\sum_{k=0}^{n}(k+1)*2^{k+1}=\sum_{k=0}^nk*2^{k+1}+\sum_{k=0}^n2^{k+1}=2*S_n+2^{n+2}-2
兩個式子聯立,得:Sn+(n+1)2n+1=2Sn+2n+22Sn=(n+1)2n+12n+2+2Sn=(n1)2n+1+2S_n+(n+1)*2^{n+1}=2*S_n+2^{n+2}-2\\S_n=(n+1)*2^{n+1}-2^{n+2}+2\\S_n=(n-1)*2^{n+1}+2
那麼,我們現在將這個式子平凡化,求Sn=k=0nkxkS_n=\sum_{k=0}^nk*x^k
同上推導,在演示一遍,就是:Sn+1=k=0n+1kxk=Sn+(n+1)xn+1S_{n+1}=\sum_{k=0}^{n+1}k*x^k=S_n+(n+1)*x^{n+1}Sn+1=k=0n+1kxk=k=0n(k+1)xk+1=k=0nkxk+1+k=0nxk+1=xSn+xxn+21xS_{n+1}=\sum_{k=0}^{n+1}k*x^k=\sum_{k=0}^n(k+1)*x^{k+1}=\sum_{k=0}^nk*x^{k+1}+\sum_{k=0}^nx^{k+1}=x*S_n+\frac{x-x^{n+2}}{1-x}
聯立得:Sn+(n+1)xn+1=xSn+xxn+21xS_n+(n+1)*x^{n+1}=x*S_n+\frac{x-x^{n+2}}{1-x}(x1)Sn=(n+1)xn+1xxn+21x(x-1)*S_n=(n+1)*x^{n+1}-\frac{x-x^{n+2}}{1-x}(x1)Sn=(n+1)xn+1nxn+2x1x(x-1)*S_n=\frac{(n+1)*x^{n+1}-n*x^{n+2}-x}{1-x}Sn=x+nxn+2(n+1)xn+1(1x)2S_n=\frac{x+n*x^{n+2}-(n+1)*x^{n+1}}{(1-x)^2}



多重和式

多重和式可以交換求和次序,例如:
ijai,j[P(i,j)]=jiai,j[P(i,j)]\sum_i\sum_ja_{i,j}[P(i,j)]=\sum_j\sum_ia_{i,j}[P(i,j)]


下面,我們考慮求得SΔ=i=1nj=inaiajS_{\Delta 上}=\sum_{i=1}^n\sum_{j=i}^na_i*a_j
由於SΔ=i=1nj=inaiaj=i=1nj=inajai=SΔS_{\Delta 上}=\sum_{i=1}^n\sum_{j=i}^na_i*a_j=\sum_{i=1}^n\sum_{j=i}^na_j*a_i=S_{\Delta 下}
2SΔ=SΔ+SΔ=i=1nj=1naiaj+i=1nai2=(i=1nai)2+i=1nai22*S_{\Delta 上}=S_{\Delta 上}+S_{\Delta 下}=\sum_{i=1}^n\sum_{j=1}^na_i*a_j+\sum_{i=1}^na_i^2=(\sum_{i=1}^na_i)^2+\sum_{i=1}^na_i^2
完成。


下面繼續看這個和式:Sn=1i<jn(aiaj)(bibj)S_n=\sum_{1\le i<j\le n}(a_i-a_j)(b_i-b_j)
2Sn=1i,jn(aiaj)(bibj)1i=jn(aiaj)(bibj)=1i,jnaibi1i,jnaibj1i,jnajbi+1i,jnajbj1i=jn(aiaj)(bibj)=2ni=1naibi2i=1naii=1nbi1i=jn(aiaj)(bibj)2*S_n=\sum_{1\le i,j\le n}(a_i-a_j)(b_i-b_j)-\sum_{1\le i=j\le n}(a_i-a_j)(b_i-b_j)\\=\sum_{1\le i,j\le n}a_ib_i-\sum_{1\le i,j\le n}a_ib_j-\sum_{1\le i,j\le n}a_jb_i+\sum_{1\le i,j\le n}a_jb_j-\sum_{1\le i=j\le n}(a_i-a_j)(b_i-b_j)\\=2n*\sum_{i=1}^na_ib_i-2*\sum_{i=1}^na_i\sum_{i=1}^nb_i-\sum_{1\le i=j\le n}(a_i-a_j)(b_i-b_j)Sn=ni=1naibii=1naii=1nbiS_n=n*\sum_{i=1}^na_ib_i-\sum_{i=1}^na_i\sum_{i=1}^nb_i
所以,移項得:i=1naii=1nbi=ni=1naibi1i<jn(aiaj)(bibj)\sum_{i=1}^na_i\sum_{i=1}^nb_i=n*\sum_{i=1}^na_ib_i-\sum_{1\le i<j\le n}(a_i-a_j)(b_i-b_j)


那麼研究這個和式Sn=1i<jn1jiS_n=\sum_{1\le i<j\le n}\frac{1}{j-i}
和式展開:Sn=j=1ni=1j11ji=j=1nji=1j11i=j=1ni=1j11i=j=1nHj1=j=1n1HjS_n=\sum_{j=1}^n\sum_{i=1}^{j-1}\frac{1}{j-i}\\=\sum_{j=1}^n\sum_{j-i=1}^{j-1}\frac{1}{i}\\=\sum_{j=1}^n\sum_{i=1}^{j-1}\frac{1}{i}\\=\sum_{j=1}^nH_{j-1}\\=\sum_{j=1}^{n-1}H_j
但是這樣還是2\sum^2的複雜程度。我們這個思路還是不是太優。這是很正常的,我們需要嘗試另一種方式求解(?怎麼嘗試另一種思路呢?上面是消除jj,下面我們消除ii)。
1i<jn1ji=1j<j+in1j=j=1ni=1nj1j=j=1nnjj=j=1nnjn=nHnn\sum_{1\le i<j\le n}\frac{1}{j-i}\\=\sum_{1\le j<j+i\le n}\frac{1}{j}\\=\sum_{j=1}^n\sum_{i=1}^{n-j}\frac{1}{j}\\=\sum_{j=1}^n\frac{n-j}{j}\\=\sum_{j=1}^n\frac{n}{j}-n\\=n*H_n-n
我們將上面的錯誤結論和現在的正確結論結合起來,就得到一個很有趣的式子。
i=1n1Hi=nHnn\sum_{i=1}^{n-1}H_i=n*H_n-n
這個式子看起來是如此的美妙。



一般性方法

推導了那麼多,我們現在來研究研究,遇到一個式子應該怎麼做,來求出他的封閉形式的解。
比如我們現在遇到了這樣的和式Sn=k=1nk2S_n=\sum_{k=1}^nk^2
我們遇到這樣的和式,一般都是要列一下前幾項。
n012345678n201491625364964Sn01514305591140204\begin{array}{c|lcr} n & 0 & 1 & 2 & 3 & 4 & 5 & 6 & 7 & 8 \\ \hline n^2 & 0 & 1 & 4 & 9 & 16 & 25 & 36 & 49 & 64\\ \hline S_n & 0 & 1 & 5 & 14 &30 & 55 & 91 & 140 & 204 \\ \end{array}


方法一、OEIS

首先,我們可以查表,查出這個函數是什麼。我們可以登錄OEIS進行查找。
然後我們就可以搜索到Sn=n(n+1)(2n+1)6S_n=\frac{n(n+1)(2n+1)}{6}
這樣的結果。


方法二、猜測答案,然後歸納法驗證

有些數列數列是特別容易就猜到答案了,我們需要在這基礎上進行驗證。
現在假設我們根據前幾項猜出來了通項公式是Sn=n(n+1)(2n+1)6S_n=\frac{n(n+1)(2n+1)}{6}
那麼我們怎麼證明他的正確性呢?
Sn=Sn1+n2=(n1)n(2n1)+6n26=2n3+3n2+n6=n(n+1)(2n+1)6S_n=S_{n-1}+n^2\\=\frac{(n-1)n(2n-1)+6n^2}{6}\\=\frac{2n^3+3n^2+n}{6}\\=\frac{n(n+1)(2n+1)}{6}
又由於S0S_0也成立,所以得證。


方法三、擾動法

我們加上最有一項。
Sn+(n+1)2=k=0n(k+1)2=k=0nk2+2k=0nk+k=0n1=Sn+2k=0nk+n+1S_n+(n+1)^2=\sum_{k=0}^n(k+1)^2\\=\sum_{k=0}^nk^2+2*\sum_{k=0}^nk+\sum_{k=0}^n1\\=S_n+2*\sum_{k=0}^nk+n+1
然鵝,左右兩邊的SnS_n消去了,得到了一個毫無卵用的恆等式。
但是,我們發現,這樣能消去一維,那麼我們將初始的式子設爲三維的不就可以了嗎?
Tn=k=0nk3T_n=\sum_{k=0}^nk^3
Tn+(n+1)3=k=0n(n+1)3=k=0nk3+3k=0nk2+3k=0nk+k=0n1=Tn+3Sn+3n(1+n)2+n+1T_n+(n+1)^3=\sum_{k=0}^n(n+1)^3\\=\sum_{k=0}^nk^3+3\sum_{k=0}^nk^2+3\sum_{k=0}^nk+\sum_{k=0}^n1\\=T_n+3S_n+\frac{3n(1+n)}{2}+n+1
3Sn=(n+1)33n2+3n2n13S_n=(n+1)^3-\frac{3n^2+3n}{2}-n-1Sn=2n3+3n2+n6S_n=\frac{2n^3+3n^2+n}{6}Sn=n(n+1)(2n+1)6S_n=\frac{n(n+1)(2n+1)}{6}


方法四、建立成套方法

之前我們研究了遞歸式{R0=α,Rn=Rn1+β+γn+δn2.\left \{ \begin{array}{c} R_0=\alpha,\\ R_n=R_{n-1}+\beta+\gamma n+\delta n^2. \end{array} \right .
我們可以說這樣的遞歸式的一般解是這樣的:Rn=A(n)α+B(n)β+C(n)γ+D(n)δR_n=A(n)\alpha+B(n)\beta+C(n)\gamma+D(n)\delta
現在,我們可以確信的是,D(n)=Sn,C(n)=k=1nk,B(n)=k=1n1D(n)=S_n,C(n)=\sum_{k=1}^nk,B(n)=\sum_{k=1}^n1,這是不容置否的。
我們現在令Rn=n3R_n=n^3,那麼α=0,β=1,γ=3,δ=3\alpha=0,\beta=1,\gamma=-3,\delta=3
所以:B(n)3C(n)+3D(n)=n3B(n)-3C(n)+3D(n)=n^3
Sn=D(n)=n3+3C(n)B(n)3=n3+3k=0nkn3=n3+3n(1+n)2n3=2n3+3n2+n6=n(n+1)(2n+1)6S_n=D(n)=\frac{n^3+3C(n)-B(n)}{3}\\=\frac{n^3+3\sum_{k=0}^nk-n}{3}\\=\frac{n^3+\frac{3n(1+n)}{2}-n}{3}\\=\frac{2n^3+3n^2+n}{6}\\=\frac{n(n+1)(2n+1)}{6}


方法五、用積分代替和式

我們用f(x)=x2f(x)=x^2這一曲線的面積近似代替k=0nk2\sum_{k=0}^nk^2,求得:
k=0nk20nx2 dx=n33\sum_{k=0}^nk^2\approx\int_0^nx^2\ dx=\frac{n^3}{3}
但是這是約等,我們記Sn=k=0nk2=n33+EnS_n=\sum_{k=0}^nk^2=\frac{n^3}{3}+E_n
En=Snn33=Sn1+n2n33=(n1)33+n2n33+En1E_n=S_n-\frac{n^3}{3}=S_{n-1}+n^2-\frac{n^3}{3}=\frac{(n-1)^3}{3}+n^2-\frac{n^3}{3}+E_{n-1}En=En1+n13E_n=E_{n-1}+n-\frac{1}{3}
然後EnE_n就可以借用上面的方法進行求解了。


方法六、展開和收縮

Sn=k=0nk2=i=0nj=inj=i=0ni+n2(ni+1)=i=0nn(n+1)i2+i2=n3+n22+n(n+1)4Sn2S_n=\sum_{k=0}^nk^2=\sum_{i=0}^n\sum_{j=i}^nj\\=\sum_{i=0}^n\frac{i+n}{2}*(n-i+1)\\=\frac{\sum_{i=0}^nn(n+1)-i^2+i}{2}\\=\frac{n^3+n^2}{2}+\frac{n(n+1)}{4}-\frac{S_n}{2}
3Sn2=2n3+3n2+n4\frac{3*S_n}{2}=\frac{2n^3+3n^2+n}{4}
Sn=n(n+1)(2n+1)6S_n=\frac{n(n+1)(2n+1)}{6}


小結

方法實在是特別多,後面也會遇到更精妙的算法。後面再補充。



有限微積分和無限微積分

什麼是有限微積分?什麼是無限微積分?
無限微積分是由微分算子DD所定義的性質。Df(x)=limh0f(x+h)f(x)hDf(x)=\lim_{h\to0}\frac{f(x+h)-f(x)}{h}
有限微積分是由差分算子Δ\Delta所定義的性質。Δf(x)=f(x+1)f(x)\Delta f(x)=f(x+1)-f(x)
我們知道,設f(x)=xmf(x)=x^m,則Df(x)=mxm1Df(x)=mx^{m-1}
但是,Δ\Delta算子則沒有這樣的性質。
那麼,我們就規定與之相對應的運算。
我們記xm=x(x1)(x2)...(xm+1)x^{\underline m}=x(x-1)(x-2)...(x-m+1)xm=x(x+1)(x+2)...(x+m1)x^{\overline m}=x(x+1)(x+2)...(x+m-1)
這個運算時和階乘有一點點聯繫的,n!=nn=1nn!=n^{\underline n}=1^{\overline n}
那麼Δ(xm)=(x+1)mxm=mxm1\Delta(x^{\underline m})=(x+1)^{\underline m}-x^{\underline m}=m*x^{\underline {m-1}}


下面介紹算子Δ\DeltaDD的逆運算。
DD運算的逆運算就是\intg(x)=D f(x)g(x) dx=f(x)+Cg(x)=D\ f(x) \Leftarrow\Rightarrow\int g(x)\ dx=f(x)+C
Δ\Delta運算的逆運算就是\sumg(x)=Δf(x)g(x)δx=f(x)+Cg(x)=\Delta f(x)\Leftarrow\Rightarrow\sum g(x)\delta x=f(x)+C
注意,Δ\Delta運算的CC並不一定是常數,也有可能是函數,只要滿足C(x)=C(x+1)C(x)=C(x+1)就好了。
那麼知道這一點,如果我們要求k=ab1g(x)\sum_{k=a}^{b-1}g(x),那麼我們僅僅需要找出g(x)=f(x+1)f(x)g(x)=f(x+1)-f(x),那麼答案就等於f(b)f(a)f(b)-f(a)
我們知道ab+bc=ac\int_a^b+\int_b^c=\int_a^c,同樣的ab+bc=ac\sum_a^b+\sum_b^c=\sum_a^c
基本積分規則:
0nxm=xm+1m+10n=nm+1m+1\int_0^nx^m=\frac{x^{m+1}}{m+1}|_0^n=\frac{n^{m+1}}{m+1}
x=0n1xm=xm+1m+10n=nm+1m+1\sum_{x=0}^{n-1}x^{\underline m}=\frac{x^{\underline{m+1}}}{m+1}|_0^n=\frac{n^{\underline{m+1}}}{m+1}
那麼,根據這個式子,就可以輕鬆地推出1+2+...+n1+2+...+n的結果:
x=0nx=x220n+1=(n+1)22=n(n+1)2\sum_{x=0}^nx=\frac{x^{\underline2}}{2}|_0^{n+1}=\frac{(n+1)^{\underline{2}}}{2}=\frac{n(n+1)}2
現在,我們嘗試推導12+22+...+n21^2+2^2+...+n^2的通項公式:
x=0nx2=x=0nx2+x=0nx=(n+1)33+(n+1)22=2n3+3n2+n6=n(n+1)(2n+1)6\sum_{x=0}^nx^2=\frac{}{}\sum_{x=0}^nx^{\underline2}+\sum_{x=0}^nx=\frac{(n+1)^{\underline3}}{3}+\frac{(n+1)^{\underline{2}}}{2}= \frac{2n^3+3n^2+n}{6}=\frac{n(n+1)(2n+1)}{6}
甚至,這個式子可以計算立方,只要我們求出來x3=x3+3x2+xx^3=x^{\underline3}+3*x^{\underline2}+x
所以:
x=0nx3=(n+1)44+(n+1)3+(n+1)22=n2(1+n)24\sum_{x=0}^nx^3=\frac{(n+1)^{\underline4}}4+(n+1)^{\underline3}+\frac{(n+1)^{\underline2}}{2}=\frac{n^2(1+n)^2}{4}
而且,下降冪滿足分配率,比如:
(x+y)2=x2+2xy+y2(x+y)^{\underline 2}=x^{\underline 2}+2xy+y^{\underline 2}
還有其他的一些基本的性質:
xm+n=xm(xm)nx^{\underline {m+n}}=x^{\underline m}(x-m)^{\underline n}


到這,我們開始回想下剛纔的式子
x=ab1xm δx=xm+1m+1ab\sum_{x=a}^{b-1}x^{\underline m}\ \delta x=\frac{x^{\underline{m+1}}}{m+1}|_a^b
這個式子是一直成立的,但是,當m=1m=-1時,這個式子會變成什麼。
首先,我們先探討下x1x^{\underline{-1}}是什麼?
觀察x+x^{\underline{+}}的規律,得出xm+1=xm(xm)x^{\underline{m+1}}=x^{\underline m}*(x-m)
所以,爲了保持好的性質,xm=1(x+1)(x+2)...(x+m)x^{\underline{-m}}=\frac{1}{(x+1)(x+2)...(x+m)}
經過檢驗,xx^{\underline -}x+x^{\underline +}沒有什麼不同。
所以x1=1x+1x^{\underline{-1}}=\frac{1}{x+1}
我們設:f(x)=x=ab1xm δxf(x)=\sum_{x=a}^{b-1}x^{\underline m}\ \delta x
所以:Δf(x)=f(x+1)f(x)=1x+1\Delta f(x)=f(x+1)-f(x)=\frac{1}{x+1}
所以:f(x)=Hxf(x)=H_x
那麼,我們的公式就得到了拓展:
x=ab1xm δx={xm+1m+1ab,if m̸=1HbHa,if m=1\sum_{x=a}^{b-1}x^{\underline m}\ \delta x= \begin{cases} \frac{x^{\underline{m+1}}}{m+1}|_a^b, &if\ m\not=-1\\[5ex] H_b-H_a, &if\ m=-1 \end{cases}
這個其實和\int運算是很像的。
abxm dx={xm+1m+1ab,if m̸=1ln(b)ln(a),if m=1\int_a^bx^m\ d x= \begin{cases} \frac{x^{m+1}}{m+1}|_a^b, &if\ m\not=-1\\[5ex] ln(b)-ln(a), &if\ m=-1 \end{cases}
其實HxH_x就是對ln(x)ln(x)的離散模擬。


我們找到了D ex=exD\ e^x=e^x,那麼是否存在Δ f(x)=f(x)\Delta\ f(x)=f(x)呢?
這其實是很好找的:Δ f(x)=f(x+1)f(x)=f(x)f(x+1)=2f(x)\Delta\ f(x)=f(x+1)-f(x)=f(x)\\[2ex]f(x+1)=2f(x)
只要滿足這個式子就好了。


Δ(cx)=cx+1cx=(c1)cx\Delta(c^x)=c^{x+1}-c^x=(c-1)c^x
ax<bcx=x=ab1cx δx=cxc1ab=cbcac1\sum_{a\le x<b}c^x=\sum_{x=a}^{b-1}c^x\ \delta x=\frac{c^x}{c-1}|_a^b=\frac{c^b-c^a}{c-1}



習題


記號k=40qk\sum_{k=4}^0q_k的含義是什麼?

這個其實是沒有嚴格的定義的,你可以認爲是0。也可以認爲是q0+q1+q2+q3+q4。也可以認爲是-q4-q3-q2-q1-q0。


化簡表達式x([x>0][x<0])x*([x>0]-[x<0])
x([x>0][x<0])={x,if x0x,if x<0=xx*([x>0]-[x<0])=\begin{cases} x, &if\ x\ge0\\[2ex] -x, &if\ x<0 \end{cases}=|x|


將和式k=05ak\sum_{k=0}^5a_kk2=05ak2\sum_{k^2=0}^5a_{k^2}完全展開以證實你對記號\sum的理解。
k=05ak=a0+a1+a2+a3+a4+a5\sum_{k=0}^5a_k=a_0+a_1+a_2+a3+a_4+a_5
k2=05ak2=a(2)2+a(1)2+a0+a12+a22=a4+a1+a0+a1+a4\sum_{k^2=0}^5a_{k^2}=a_{(-2)^2}+a_{(-1)^2}+a_{0}+a_{1^2}+a_{2^2}=a_4+a_1+a_0+a_1+a_4


將三重和式1i<j<k4aijk\sum_{1\le i<j<k\le4}a_{ijk}表示成三個重疊的和式
a、先對kk,再對jj,再對ii求和。
1i<j<k4aijk=i=14j=i+14k=j+14aijk=a123+a124+a134+a234\sum_{1\le i<j<k\le4}a_{ijk}=\sum_{i=1}^4\sum_{j=i+1}^4\sum_{k=j+1}^4a_{ijk}=a_{123}+a_{124}+a_{134}+a_{234}
b、先對ii,再對jj,再對kk求和。
1i<j<k4aijk=k=14j=1k1i=1j1aijk=a123+a124+a134+a234\sum_{1\le i<j<k\le4}a_{ijk}=\sum_{k=1}^4\sum_{j=1}^{k-1}\sum_{i=1}^{j-1}a_{ijk}=a_{123}+a_{124}+a_{134}+a_{234}


下面的推導何處有錯?(j=1naj)(k=1n1ak)=j=1nk=1najak=k=1nk=1nakak=n2(\sum_{j=1}^na_j)(\sum_{k=1}^n\frac{1}{a_k})=\sum_{j=1}^n\sum_{k=1}^n\frac{a_j}{a_k}=\sum_{k=1}^n\sum_{k=1}^n\frac{a_k}{a_k}=n^2

第二步到第三步錯了。


作爲jjnn的函數,k[1jkn]\sum_k[1\le j\le k\le n]的值是什麼?
k[1jkn]={0,if (j<0)  (n<1)  (n<j)nj+1,else\sum_k[1\le j\le k\le n]=\begin{cases} 0, &if\ (j<0)\ ||\ (n<1)\ ||\ (n<j)\\[2ex] n-j+1,&else \end{cases}


f(x)=f(x)f(x1)\nabla f(x)=f(x)-f(x-1)(xm)\nabla(x^{\overline{m}})是什麼?
(xm)=xm(x1)m=mxm1\nabla(x^{\overline m})=x^{\overline m}-(x-1)^{\overline m}=m*x^{\overline{m-1}}


mm是給定的整數時,0m0^{\underline m}的值是什麼?
0m={0,if m01m!,if m<00^{\underline m}=\begin{cases} 0, &if\ m\ge0\\[2ex] \frac{1}{|m|!},&if\ m<0 \end{cases}


證明k=0n1(ak+1ak)bk=anbna0b0k=0n1ak+1(bk+1bk)\sum_{k=0}^{n-1}(a_{k+1}-a_k)b_k=a_nb_n-a_0b_0-\sum_{k=0}^{n-1}a_{k+1}(b_{k+1}-b_k)
這裏需要先證明一個引理:
Δ(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)=v(x+1)Δu(x)+u(x)Δv(x)\Delta(u(x)v(x))=u(x+1)v(x+1)-u(x)v(x)\\[2ex]=u(x+1)v(x+1)-u(x)v(x+1)+u(x)v(x+1)-u(x)v(x)\\[2ex]=v(x+1)\Delta u(x)+u(x)\Delta v(x)
這裏引入一個算子EE,表名Ef(x)=f(x+1)Ef(x)=f(x+1)
Δ(uv)=uΔv+EvΔuuΔv=Δ(uv)EvΔuuΔv=uvEvΔu\Delta(uv)=u\Delta v+Ev\Delta u\\[2ex]u\Delta v=\Delta(uv)-Ev\Delta u\\[2ex]\sum u\Delta v=uv-\sum Ev\Delta u
同樣,定積分也是這樣:
k=0n1uΔv=uv0nk=0n1EvΔu\sum_{k=0}^{n-1}u\Delta v=uv|_0^n-\sum_{k=0}^{n-1}Ev\Delta u
那麼,答案就顯而易見了:
k=0n1(ak+1ak)bk=k=0n1bkΔak=anbna0b0k=0n1ak+1(bk+1bk)\sum_{k=0}^{n-1}(a_{k+1}-a_k)b_k=\sum_{k=0}^{n-1}b_k\Delta a_k=a_nb_n-a_0b_0-\sum_{k=0}^{n-1}a_{k+1}(b_{k+1}-b_k)


證明:只要c是一個整數,函數p(k)=k+(1)kcp(k)=k+(-1)^kc就是對所有整數的全排列。
p(k)=np(k)=n,則:
n=k+(1)kcn+c=k+((1)k+1)cn=k+(-1)^kc\\[2ex]n+c=k+((-1)^k+1)c
因爲((1)k+1)((-1)^k+1)是偶數,所以:
(1)n+c=(1)k  &&  k=n(1)kc(-1)^{n+c}=(-1)^k\ \ \&\&\ \ k=n-(-1)^kc
所以,由kk能夠唯一的推出nn,而且由於奇偶性,不可能重複,所以就是全排列。


利用成套方法求k=0n(1)kk2\sum_{k=0}^n(-1)^kk^2的封閉形式。
設:
{R0=α,Rn=Rn1+(1)n(β+nγ+n2δ)\begin{cases} R_0=\alpha,\\[5ex] R_n=R_{n-1}+(-1)^n(\beta+n\gamma+n^2\delta) \end{cases}
則:R(n)=A(n)α+B(n)β+C(n)γ+D(n)δR(n)=A(n)\alpha+B(n)\beta+C(n)\gamma+D(n)\delta
易知,這時候我們要求的答案就蘊含在D(n)D(n)裏面。
考慮Rn=(1)nnR_n=(-1)^nn時,研究A(n)B(n)C(n)A(n)、B(n)、C(n)的關係:2C(n)B(n)=(1)nn2C(n)-B(n)=(-1)^nn
考慮Rn=(1)nn2R_n=(-1)^nn^2時,研究A(n)B(n)C(n)D(n)A(n)、B(n)、C(n)、D(n)的關係:2D(n)2C(n)+B(n)=(1)nn22D(n)-2C(n)+B(n)=(-1)^nn^2
所以,k=0n(1)kk2=D(n)=(1)n(n2+n)2\sum_{k=0}^n(-1)^kk^2=D(n)=\frac{(-1)^n(n^2+n)}{2}

//mathematics代碼
Simplify[Sum[(-1)^k*k^2, {k, 0, n}]]

k=1nk2k\sum_{k=1}^nk*2^k重新改寫成多重和式1jkn2k\sum_{1\le j\le k\le n}2^k的形式來對他進行計算。
k=1nk2k=j=1nk=jn2k=2n+1(n1)+2\sum_{k=1}^nk*2^k=\sum_{j=1}^n\sum_{k=j}^n2^k=2^{n+1}(n-1)+2

//Mathematics代碼
FullSimplify[Sum[k*2^k, {k, 1, n}]]

用正文中的方法5來計算Tn=k=1nk3T_n=\sum_{k=1}^nk^3:首先記Tn+Sn=21jknjkT_n+S_n=2\sum_{1\le j\le k\le n}jk,然後應用。
Tn+Sn=k=1nk2(k+1)=2j=1nk=jnjk=j=1nj2+(j=1nj)2=Sn+(n(n+1))24Tn=n2(n+1)24T_n+S_n=\sum_{k=1}^nk^2(k+1)=2*\sum_{j=1}^n\sum_{k=j}^njk=\sum_{j=1}^nj^2+(\sum_{j=1}^nj)^2=S_n+\frac{(n(n+1))^2}{4}\\[2ex]T_n=\frac{n^2(n+1)^2}{4}

//mathematics代碼
FullSimplify[Sum[k^3, {k, 1, n}]]

證明xm/(xn)m=xn/(xm)nx^{\underline m}/(x-n)^{\underline m}=x^{\underline n}/(x-m)^{\underline n},除非其中有一個分母爲0。
xm(xm)n=xm+n=xn(xn)mx^{\underline m}(x-m)^{\underline n}=x^{\underline{m+n}}=x^{\underline n}(x-n)^{\underline{m}}


證明,對於所有的整數mm,下面的公式可以用來在上升階乘冪與下降階乘冪之間進行轉換:
xm=(1)m(x)m=(x+m1)m=1(x1)mx^{\overline m}=(-1)^m(-x)^{\underline m}=(x+m-1)^{\underline m}=\frac{1}{(x-1)^{\underline{-m}}}
xm=(1)m(x)m=(xm+1)m=1(x+1)mx^{\underline m}=(-1)^m(-x)^{\overline m}=(x-m+1)^{\overline m}=\frac{1}{(x+1)^{\overline{-m}}}


z\Re_zz\Im_z是複數zz的實部和虛部,其絕對值z|z|(z)2+(z)2\sqrt{(\Re_z)^2+(\Im_z)^2},當實值的和式kKak\sum_{k\in K}\Re_{a_k}kKak\sum_{k\in K}\Im_{a_k}兩者都絕對收斂時,就說複數項aka_k組成的和式kKak\sum_{k\in K}a_k絕對收斂。證明:kKak\sum_{k\in K}a_k絕對收斂,當且僅當存在一個有界常數BB,使得對所有的有限子集FKF\subseteq K,都有kFakB\sum_{k\in F}|a_k|\le B
沒看懂題目。


利用求和因子來求解遞歸式T0=5;2Tn=nTn1+3n!,n>0T_0=5;2T_n=nT_{n-1}+3*n!,n>0
2Tn=nTn1+3n!2nTnn!=2n1Tn1(n1)!+322n2T_n=nT_{n-1}+3*n!\\[2ex]\frac{2^nT_n}{n!}=\frac{2^{n-1}T_{n-1}}{(n-1)!}+\frac{3}{2}2^n
Gn=2nTnn!G_n=\frac{2^nT_n}{n!},那麼GnG_n的遞歸式爲:
Gn={5,if n=0,Gn1+322n,if n>0.G_n= \begin{cases} 5, &if\ n=0,\\[5ex] G_{n-1}+\frac{3}{2}2^n, &if\ n>0. \end{cases}
所以:Gn=5+32i=1n2i=32n+2G_n=5+\frac32\sum_{i=1}^n2^i=3*2^n+2
所以:Tn=3n!+2n!2nT_n=3*n!+\frac{2*n!}{2^n}


試用擾動法計算k=0nkHk\sum_{k=0}^nkH_k,不過改爲推導出k=0nHk\sum_{k=0}^nH_k的值。
這一道題根據mathematics秒出答案。
k=0nkHk=n(n+1)(2Hn+11)4k=0nHk=(n+1)(Hn+11)\sum_{k=0}^nkH_k=\frac{n(n+1)(2H_{n+1}-1)}{4}\\[2ex]\sum_{k=0}^nH_k=(n+1)(H_{n+1}-1)

Simplify[Sum[k*HarmonicNumber[k], {k, 0, n}]]
Simplify[Sum[HarmonicNumber[k], {k, 0, n}]]

不過,還是要親自計算一下,記Sn=k=0nkHkS_n=\sum_{k=0}^nkH_k
Sn+(n+1)Hn+1=Sn+k=0nHk+n+1k=0nHk=(n+1)(Hn+11)S_n+(n+1)H_{n+1}=S_n+\sum_{k=0}^nH_k+n+1\\[2ex]\sum_{k=0}^nH_k=(n+1)(H_{n+1}-1)


假設n0n\ge 0,用擾動法計算和式Sn=k=0n(1)nkS_n=\sum_{k=0}^n(-1)^{n-k}Tn=k=0n(1)nkkT_n=\sum_{k=0}^n(-1)^{n-k}k以及Un=k=0n(1)nkk2U_n=\sum_{k=0}^n(-1)^{n-k}k^2
Sn=k=0n(1)nk=(1+(1)n)(1)n2Tn=k=0n(1)nkk=(1)n2(1)nn+(1)n14Un=k=0n(1)nkk2=(1)2nn(n+1)2S_n=\sum_{k=0}^n(-1)^{n-k}=\frac{(1+(-1)^n)(-1)^n}{2}\\[2ex]T_n=\sum_{k=0}^n(-1)^{n-k}k=(-1)^n\frac{2(-1)^nn+(-1)^n-1}{4}\\[2ex]U_n=\sum_{k=0}^n(-1)^{n-k}k^2=\frac{(-1)^{2n}n(n+1)}{2}

mathematics代碼:
Sum[(-1)^(n - k), {k, 0, n}]
Sum[(-1)^(n - k) k, {k, 0, n}]
Sum[(-1)^(n - k) k^2, {k, 0, n}]

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章