この記事では 前回の記事 に引き続き超幾何数列の基本事項についてまとめていきます。
さて
前回の記事
では手計算向けのGosper's algorithmこと掃き出し法というものについて解説しました。
しかし掃き出し法はいつでも直感や手に優しいということはなく、例えば
$$H(n)=\frac1{(n+1)(n+m+1)}$$
のような数列の扱いには手を焼かすこととなります。
実際この手法においてはわざわざ分母を階乗に書き直して
$$H(n)=(n+2)(n+3)\cdots(n+m)\times\frac{n!}{(n+m+1)!}$$
と変形しなければならず、またこの数列に対応する$T(n)$は
$$T(n)=\frac{f_m(n)}{(n+1)(n+2)\cdots(n+m)}$$
という形で求めることになりますが、例えば$m=1,2,3,4,5$において
\begin{align}
f_1(n)&=-1\\
f_2(n)&=-\frac{2n+3}2\\
f_3(n)&=-\frac{3n^2+12n+11}3\\
f_4(n)&=-\frac{4n^3+30n^2+70n+50}4\\
f_5(n)&=-\frac{5n^4+60n^3+255n^2+450n+274}5
\end{align}
と求まるように非常に計算がダルそうであることが見て取れます。
さらにこのように掃き出し法にて求まる$T(n)$の表示は煩雑な見た目になるのに対し、実のところ$T(n)$は
$$T(n)=-\frac1m\l(\frac1{n+1}+\frac1{n+2}+\cdots+\frac1{n+m}\r)$$
という簡潔な表示を持っており、初めからこのような表示で$T(n)$を求めることができるならそれに越したことはありません。
そしてそんな近道を実現できるのがAbramov-Petkovšek reductionという手法になります。
Abramov-Petkovšek reduction、以下これを簡約化の手法と呼びます、とは大雑把に言うと$H(n)$を適当に
$$H(n)=(\mbox{有理関数})\times\frac{A(n)}{B(n)}$$
という形に分解したとき、部分分数分解と平行移動を繰り返し用いることで
$$H(n)=(\mbox{多項式})\times\frac{A(n)}{B(n)}+T'(n+1)-T'(n)$$
という掃き出し法で処理しやすい形に帰着させる手法となっています。
と言ってもやることは簡単でまず
$$H(n)=(\mbox{有理関数})\times\frac{A(n)}{B(n)}$$
と分解します。
この場合においては$A(n),B(n)$が
を満たすように取ることが最適となります。
そしてこの有理関数部分を部分分数分解したときに
$$\l(\frac{p(n)}{f(n+m)}+\frac{q(n)}{f(n)}\r)\times\frac{A(n)}{B(n)}$$
という因子があったとすると、この第一項を
\begin{align}
\frac{p(n)}{f(n+m)}\frac{A(n)}{B(n)}
&=\frac{p(n-m)}{f(n)}\frac{A(n-m)}{B(n-m)}+T'(n+1)-T'(n)\\
\bigg(T'(n)&=\sum^m_{k=1}\frac{p(n-k)}{f(n+m-k)}\frac{A(n-k)}{B(n-k)}\bigg)
\end{align}
と平行移動させる、あるいは第二項を
\begin{align}
\frac{q(n)}{f(n)}\frac{A(n)}{B(n)}
&=\frac{q(n+m)}{f(n+m)}\frac{A(n+m)}{B(n+m)}+T'(n+1)-T'(n)\\
\bigg(T'(n)&=-\sum^{m-1}_{k=0}\frac{q(n+k)}{f(n+k)}\frac{A(n+k)}{B(n+k)}\bigg)
\end{align}
と平行移動させることで分母を揃える、というのが主な方針となります。
ただこのとき
$$\frac{p(n-m)}{f(n)}\frac{A(n-m)}{B(n-m)}+\frac{q(n)}{f(n)}\frac{A(n)}{B(n)}$$
とか
$$\frac{p(n)}{f(n+m)}\frac{A(n)}{B(n)}+\frac{q(n+m)}{f(n+m)}\frac{A(n+m)}{B(n+m)}$$
とかをそのまま足し合わせても却って煩雑になるだけなので、これを
$$\frac{(\mbox{多項式})}{f(n)}\frac{A(n)}{B(n)}\quad\mbox{とか}\quad
\frac{(\mbox{多項式})}{f(n+m)}\frac{A(n)}{B(n)}$$
という形に帰着させられるように一工夫する必要があり、その一工夫こそがこの手法の肝と言えます。
具体的には平行移動する前に部分分数分解を上手く用いるで
$$\frac{p(n)}{f(n+m)}\frac{A(n)}{B(n)}
=\frac{P(n)}{f(n+m)}\frac{A(n+m)}{B(n+m)}+R(n)\frac{A(n)}{B(n)}$$
とか
$$\frac{q(n)}{f(n)}\frac{A(n)}{B(n)}
=\frac{Q(n)}{f(n)}\frac{A(n-m)}{B(n-m)}+S(n)\frac{A(n)}{B(n)}$$
という形に変形しておくことができ、この第一項を平行移動させることで
\begin{align}
\l(\frac{p(n)}{f(n+m)}+\frac{q(n)}{f(n)}\r)\frac{A(n)}{B(n)}
&=\l(\frac{P(n-m)+q(n)}{f(n)}+R(n)\r)\frac{A(n)}{B(n)}+T'(n+1)-T'(n)\\
&=\l(\frac{p(n)+Q(n+m)}{f(n+m)}+S(n)\r)\frac{A(n)}{B(n)}+T''(n+1)-T''(n)
\end{align}
のように計算することができます。
なおその変形の計算にはいくつかの方法があるので詳しくは以下の折り畳みをご覧ください。ただこれを見てもよくわかんないと思うので、計算例を見て学んだ方が早いかもしれません。
$A(n)$や$B(n)$から平行移動させたいだけの因子を取り出して部分分数分解することで
\begin{align}
\frac{p(n)}{f(n+m)}\frac{A(n)}{B(n)}
&=\frac{p(n)}{f(n+m)a(n)a(n+1)\cdots a(n+m-1)}\frac{A(n+m)}{B(n)}\\
&=\l(\frac{P(n)}{f(n+m)}+\frac{R(n)}{a(n)a(n+1)\cdots a(n+m-1)}\r)\frac{A(n+m)}{B(n)}\\
&=\frac{P(n)b(n)b(n+1)\cdots b(n+m-1)}{f(n+m)}\frac{A(n+m)}{B(n+m)}+R(n)\frac{A(n)}{B(n)}
\end{align}
とか
\begin{align}
\frac{q(n)}{f(n)}\frac{A(n)}{B(n)}
&=\frac{q(n)}{f(n)b(n-1)b(n-2)\cdots b(n-m)}\frac{A(n)}{B(n-m)}\\
&=\l(\frac{Q(n)}{f(n)}+\frac{S(n)}{b(n-1)b(n-2)\cdots b(n-m)}\r)\frac{A(n)}{B(n-m)}\\
&=\frac{Q(n)a(n-1)a(n-2)\cdots a(n-m)}{f(n)}\frac{A(n-m)}{B(n-m)}+S(n)\frac{A(n)}{B(n)}
\end{align}
と変形する。
・メリット:方法2,3と違って一遍に平行移動できるので考えるのが楽。
・デメリット:分母がデカくなるので部分分数分解の計算が面倒。しばしば最後の式から分子の次数下げをして$R(n)$や$S(n)$に組み込む必要があって二度手間感がある。
方法1を用いて一つずつ平行移動していき
\begin{align}
\frac{P_0(n)}{f(n+m)}\frac{A(n)}{B(n)}
&=\frac{P_1(n)}{f(n+m)}\frac{A(n+1)}{B(n+1)}+R_1(n)\frac{A(n)}{B(n)}\\
\frac{P_1(n-1)}{f(n+m-1)}\frac{A(n)}{B(n)}
&=\frac{P_2(n)}{f(n+m-1)}\frac{A(n+1)}{B(n+1)}+R_2(n)\frac{A(n)}{B(n)}\\
&\ \ \vdots\\
\frac{P_{m-1}(n-1)}{f(n+1)}\frac{A(n)}{B(n)}
&=\frac{P_m(n)}{f(n+1)}\frac{A(n+1)}{B(n+1)}+R_m(n)\frac{A(n)}{B(n)}
\end{align}
と変形することで
\begin{align}
\frac{P_0(n)}{f(n+m)}\frac{A(n)}{B(n)}
=\frac{P_m(n-1)}{f(n)}\frac{A(n)}{B(n)}+R(n)\frac{A(n)}{B(n)}+T'(n+1)-T'(n)\\
\l(R(n)=\sum^m_{k=1}R_k(n),\quad
T'(n)=\sum^m_{k=1}\frac{P_k(n-1)}{f(n+m-k)}\frac{A(n)}{B(n)}\r)
\end{align}
とか、同様にして
\begin{align}
\frac{Q_0(n)}{f(n)}\frac{A(n)}{B(n)}
=\frac{Q_m(n+1)}{f(n+m)}\frac{A(n)}{B(n)}+S(n)\frac{A(n)}{B(n)}+T'(n+1)-T'(n)\\
\l(S(n)=\sum^m_{k=1}S_k(n),\quad
T'(n)=-\sum^m_{k=1}\frac{Q_k(n)}{f(n+k-1)}\frac{A(n-1)}{B(n-1)}\r)
\end{align}
と平行移動させる。
・メリット:一回一回の計算は方法1より楽。$R(n)$や$S(n)$の次数が大きくなり過ぎないため後が楽。
・デメリット:$m$回も部分分数分解する必要がありダルいし計算間違いを起こしやすい
概ね方法2と同様だが、最初から$A(n)$と$B(n)$の両方から因子を取り出して
\begin{align}
\frac{P_0(n)}{f(n+m)}\frac{A(n)}{B(n)}
&=\frac{P_0(n)b(n)}{f(n+m)a(n)}\frac{A(n+1)}{B(n+1)}\\
&=\l(\frac{P_1(n)}{f(n+m)}+\frac{R_1(n)}{a(n)}\r)\frac{A(n+1)}{B(n+1)}\\
&=\frac{P_1(n)}{f(n+m-1)}\frac{A(n+1)}{B(n+1)}+R_1(n)\frac{A(n)}{B(n+1)}
\end{align}
とか
\begin{align}
\frac{Q_0(n)}{f(n)}\frac{A(n)}{B(n)}
&=\frac{Q_0(n)a(n-1)}{f(n)b(n-1)}\frac{A(n-1)}{B(n-1)}\\
&=\l(\frac{Q_1(n)}{f(n)}+\frac{S_1(n)}{b(n-1)}\r)\frac{A(n-1)}{B(n-1)}\\
&=\frac{Q_1(n)}{f(n)}\frac{A(n-1)}{B(n-1)}+S_1(n)\frac{A(n-1)}{B(n)}
\end{align}
のように部分分数分解していくことで
$$\frac{P_0(n)}{f(n+m)}\frac{A(n)}{B(n)}
=\frac{P_m(n-1)}{f(n)}\frac{A(n)}{B(n)}+R(n)\frac{A(n)}{B(n+1)}+T'(n+1)-T'(n)$$
とか
$$\frac{Q_0(n)}{f(n)}\frac{A(n)}{B(n)}
=\frac{Q_m(n+1)}{f(n+m)}\frac{A(n)}{B(n)}+S(n)\frac{A(n-1)}{B(n)}+T'(n+1)-T'(n)$$
という形に帰着させる。
・メリット:方法1,2における分子の次数下げという二度手間が必要ない。
・デメリット:方法2と同じくダルい。方法1のように一遍に平行移動させる手法に拡張できない(多分)。$A$と$B$の引数がズレるのが厄介。
この手法の便利なところは$H(n)$がGosper総和可能であれば、このように分母が$f(n)$に関する項たちを揃えたときにそれらの項が打ち消し合う、つまり
\begin{align}
&\l(\frac{p_0(n)}{f(n)}+\frac{p_1(n)}{f(n+1)}+\cdots+\frac{p_m(n)}{f(n+m)}\r)\times\frac{A(n)}{B(n)}\\
={}&(\mbox{多項式})\times\frac{A(n)}{B(n)}+T'(n+1)-T'(n)
\end{align}
という形に帰着され、したがって$H(n)$は
$$H(n)=(\mbox{多項式})\times\frac{A(n)}{B(n)}+T'(n+1)-T'(n)$$
という形に簡約化されることにあります。
こうなってしまえば後はこの第一項について掃き出し法を適用することで$T(n)$を求めることができるわけです。
そして逆に言えば分母が$f(n)$に関する項を揃えても$f(n)$が消えず
\begin{align}
&\l(\frac{p_0(n)}{f(n)}+\frac{p_1(n)}{f(n+1)}+\cdots+\frac{p_m(n)}{f(n+m)}\r)\times\frac{A(n)}{B(n)}\\
={}&\frac{(\mbox{多項式})}{f(n)}\times\frac{A(n)}{B(n)}+T'(n+1)-T'(n)
\end{align}
という形にしかならないとき、それだけで$H(n)$はGosper総和不可能であることが言えてしまいます。
このことについては次のような主張によって保証されることとなります。
超幾何数列
$$H(n)=\frac{p(n)}{q(n)}\frac{A(n)}{B(n)}$$
が以下の条件
を満しているとき$H(n)$はGosper総和不可能である。
この条件(iii)は分母$q(n)$のある因子$\a(n)$が$A(n)$を"割り切らない"ことを意味しており、条件(iv)は$q(n)$のある因子$\b(n)$によって$B(n)$が"延長されない"ことを意味しています。
例えば
$$H(n)=\frac1{n-\frac23}\frac{(3n)!}{n!}$$
という表示は
$$H(n)=9n(3n-1)\cdot\frac{(3n-2)!}{n!}$$
と約分できるため条件(iii)を満たさず、また例えば
$$H(n)=\frac1{(n+4)n!}$$
という表示は分母の$n!$を
$$H(n)=\frac{(n+1)(n+2)(n+3)}{(n+4)!}$$
と延長できるため条件(iv)を満たしません。
例えば定理の条件を満たす数列
$$H(n)=\frac1{n(2n+1)}\frac{(\frac12)_n}{n!}$$
を考えたとき、掃き出し法においてこれは
$$H(n)=1\times\frac{(\frac12)_n(n-1)!}{(\frac12)_{n+1}n!}\frac{(\frac12)_n}{n!}$$
と分解しなければならず、したがって求める$T(n)$は
\begin{align}
T(n)
&=f(n)\frac{(\frac12)_n(n-1)!}{(\frac12)_{(n-1)+1}(n-1)!}\frac{(\frac12)_n}{(n-1)!}\\
&=f(n)\frac{(\frac12)_n}{(n-1)!}
\end{align}
という形になりますが、どのような多項式$f(n)$を持ってきても
\begin{align}
T(n+1)-T(n)
&=((n+\tfrac12)f(n+1)-nf(n))\frac{(\frac12)_n}{n!}\\
&\neq\frac1{n(2n+1)}\frac{(\frac12)_n}{n!}
\end{align}
となることから$T(n)$は求まらない、つまり$H(n)$はGosper総和不可能であることが示されます。
これを一般化したのが以下の証明となります。
条件(ii)-(iv)より掃き出し法において$H(n)$は
$$Q(n)=\prod^n_{k=0}q(k)$$
なる超幾何数列$Q(n)$を用いて
$$A(n)=p(n)\frac{A(n)Q(n-1)}{B(n)Q(n)}$$
と分解しなければならないことがわかる。
したがって求める$T(n)$は
\begin{align}
T(n)
&=f(n)\frac{A(n)Q(n-1)}{B(n-1)Q(n-1)}\\
&=f(n)\frac{A(n)}{B(n-1)}
\end{align}
という形になるが、もし$T(n)$が求まれば
\begin{align}
H(n)&=T(n+1)-T(n)\\
&=(\mbox{多項式})\times\frac{A(n)}{B(n)}
\end{align}
と表せることになり条件(i)に矛盾。
よって$H(n)$はGosper総和不可能でなければならない。
さて訳のわからん文字式をこねくり回すのはこの辺にして、あとは計算あるのみです。
$$\sum^n_{k=0}\frac1{k^2+\sqrt5k-1}$$
を求めよ。
$$\frac1{n^2+\sqrt5n-1}
=\frac13\l(\frac1{n+\frac{\sqrt5-3}2}-\frac1{n+\frac{\sqrt5+3}2}\r)$$
と部分分数分解できるので
\begin{align}
\frac1{n+\frac{\sqrt5-3}2}
&=\frac1{n+\frac{\sqrt5+3}2}+T(n+1)-T(n)\\
\bigg(T(n)&=-\frac1{n+\frac{\sqrt5-3}2}-\frac1{n+\frac{\sqrt5-1}2}-\frac1{n+\frac{\sqrt5+1}2}\bigg)
\end{align}
と平行移動させることで
\begin{align}
\sum^n_{k=0}\frac1{k^2+\sqrt5k-1}
&=\frac{T(n+1)-T(0)}3\\
&=-\frac13\l(\frac1{n+\frac{\sqrt5-1}2}+\frac1{n+\frac{\sqrt5+1}2}+\frac1{n+\frac{\sqrt5+3}2}+\frac{3-\sqrt5}2\r)
\end{align}
を得る。
$$\sum^n_{k=0}\frac{3k-1}{(k+\frac13)(k+\frac43)}\frac{(\frac12)_k}{(1)_k}$$
を求めよ。
$$\frac{3n-1}{(n+\frac13)(n+\frac43)}
=\frac5{n+\frac43}-\frac2{n+\frac13}$$
と部分分数分解でき、また
\begin{align}
\frac5{n+\frac43}\frac{(\frac12)_n}{(1)_n}
&=\frac5{(n+\frac43)(n+\frac12)}\frac{(\frac12)_{n+1}}{(1)_n}\\
&=6\l(\frac1{n+\frac12}-\frac1{n+\frac43}\r)\frac{(\frac12)_{n+1}}{(1)_n}\\
&=6\frac{(\frac12)_n}{(1)_n}-\frac{6(n+1)}{n+\frac43}\frac{(\frac12)_{n+1}}{(1)_{n+1}}
\end{align}
と変形できるので
\begin{align}
\frac{3n-1}{(n+\frac13)(n+\frac43)}\frac{(\frac12)_n}{(1)_n}
&=\l(6-\frac{6n}{n+\frac13}-\frac2{n+\frac13}\r)\frac{(\frac12)_n}{(1)_n}+T(n+1)-T(n)\\
&=T(n+1)-T(n)\\
\bigg(T(n)
&=-\frac{6n}{n+\frac13}\frac{(\frac12)_n}{(1)_n}\bigg)
\end{align}
つまり
$$\sum^n_{k=0}\frac{3k-1}{(k+\frac13)(k+\frac43)}\frac{(\frac12)_k}{(1)_k}
=-\frac6{n+\frac43}\frac{(\frac12)_{n+1}}{(1)_n}$$
を得る。
$$\sum^{n-1}_{k=2}\frac{(k-2)!}{3^k(k^2-11k+27)(k^2-9k+17)}$$
をいい感じに変形せよ。
$$\frac1{(n^2-11n+27)(n^2-9n+17)}
=\frac16\l(\frac{n-6}{n^2-11n+27}-\frac{n-4}{n^2-9n+17}\r)$$
と部分分数分解でき、また
\begin{align}
\frac{n-4}{n^2-9n+17}\frac{(n-2)!}{3^n}
&=\frac{n-4}{(n^2-9n+17)(n-1)}\frac{(n-1)!}{3^n}\\
&=\frac13\l(\frac{n-5}{n^2-9n+17}-\frac1{n-1}\r)\frac{(n-1)!}{3^n}\\
&=\frac{n-6}{n^2-11n+27}\frac{(n-2)!}{3^n}-\frac{(n-2)!}{3^{n+1}}+T(n+1)-T(n)\\
\bigg(T(n)&=\frac{n-6}{n^2-11n+27}\frac{(n-2)!}{3^n}\bigg)
\end{align}
と変形できるので
\begin{align}
\sum^{n-1}_{k=2}\frac{(k-2)!}{3^k(k^2-11k+27)(k^2-9k+17)}
&=\frac{T(2)-T(n)}6+\frac16\sum^{n-2}_{k=2}\frac{(k-2)!}{3^{k+1}}\\
&=-\frac2{243}-\frac{n-6}{6(n^2-11n+27)}\frac{(n-2)!}{3^n}+\frac1{18}\sum^{n-2}_{k=2}\frac{(k-2)!}{3^k}
\end{align}
を得る。
ちなみに簡約化の手法はその特性から次のような無限和を求めるのにも応用することができます。
$$\sum^\infty_{n=0}\frac1{(n^4+n^2+1)n!}$$
を求めよ。
$$\sum^\infty_{n=0}\frac1{n+2}\frac{(\frac12)_n}{(1)_n}x^n$$
を求めよ。
$$\frac1{n+2}\frac{(\frac12)_n}{(1)_n}x^n=(n+1)\frac{(\frac12)_n}{(1)_{n+2}}x^n$$
に注意すると
$$\frac1{n+2}\frac{(\frac12)_n}{(1)_n}x^n
=f(n)\frac{(\frac12)_{n+2}}{(1)_{n+2}}x^{n+2}+g(n)\frac{(\frac12)_n}{(1)_n}x^n$$
のような形に変形できることが期待できる。
実際
\begin{align}
\frac1{n+2}\frac{(\frac12)_n}{(1)_n}x^n
&=\frac23\l(\frac1{n+\frac12}-\frac1{n+2}\r)\frac{(\frac12)_{n+1}}{(1)_n}x^n\\
&=\frac23\l(\frac{(\frac12)_n}{(1)_n}x^n-\frac{n+1}{n+2}\frac{(\frac12)_{n+1}}{(1)_{n+1}}x^n\r)\\
\frac n{n+1}\frac{(\frac12)_n}{(1)_n}x^{n-1}
&=\l(\frac2{n+1}-\frac1{n+\frac12}\r)\frac{(\frac12)_{n+1}}{(1)_n}x^{n-1}\\
&=2\frac{(\frac12)_{n+1}}{(1)_{n+1}}x^{n-1}-\frac{(\frac12)_n}{(1)_n}x^{n-1}
\end{align}
と変形できるので
\begin{align}
\frac1{n+2}\frac{(\frac12)_n}{(1)_n}x^n
&=\frac23\l(x^2-(2-x)\r)\frac{(\frac12)_n}{(1)_n}x^{n-2}+T(n+1)-T(n)\\
&=\frac23(x+2)(x-1)\frac{(\frac12)_n}{(1)_n}x^{n-2}+T(n+1)-T(n)\\
\Bigg(T(n)&=-\frac23\bigg(\frac n{n+1}\frac{(\frac12)_n}{(1)_n}x^{n-1}+2\frac{(\frac12)_n}{(1)_n}x^{n-2}\bigg)\Bigg)
\end{align}
と表せる。
よって
\begin{align}
\sum^\infty_{n=0}\frac1{n+2}\frac{(\frac12)_n}{(1)_n}x^n
&=-T(0)+\frac23(x+2)(x-1)\sum^\infty_{n=0}\frac{(\frac12)_n}{(1)_n}x^{n-2}\\
&=\frac43\frac1{x^2}-\frac23(2+x)(1-x)\frac1{x^2\sqrt{1-x}}\\
&=\frac23\frac{2-(2+x)\sqrt{1-x}}{x^2}
\end{align}
を得る。
また 第八回の記事 にて紹介したように、簡約化の手法は次のような数列の漸化式を求めるのにも応用できます。ただ以下の計算例を見てもわかるように手計算で実行するのはあまり現実的ではありません。
$$f(n)=\sum^n_{k=0}\binom nk^3$$
の満たす漸化式を求めよ。
$$F(n,k)=\binom nk^3=\l((-1)^k\frac{(-n)_k}{(1)_k}\r)^3$$
とおいたとき
\begin{align}
F(n+1,k)
&=\l(\frac{n+1}{n+1-k}\r)^3\times F(n,k)\\
&=\frac{(n+1)^3}{k^3(n+1-k)^3}\times(-1)^k\frac{(-n)_k^3}{(1)_{k-1}^3}\\
&=\bigg(\frac1{k^3}+\frac3{(n+1)k^2}+\frac6{(n+1)^2k}\\
&\qquad+\frac1{(n+1-k)^3}+\frac3{(n+1)(n+1-k)^2}+\frac6{(n+1)^2(n+1-k)}\bigg)
\times(-1)^k\frac{(-n)_k^3}{(1)_{k-1}^3}\\
&=\bigg(1+\frac{3k}{n+1}+\frac{6k^2}{(n+1)^2}\\
&\qquad{}+1+\frac{3(n-k)}{n+1}+\frac{6(n-k)^2}{(n+1)^2}\bigg)
\times(-1)^k\frac{(-n)_k^3}{(1)_k^3}+T(n,k+1)-T(n,k)\\
&=\l(\frac{5n+2}{n+1}+\frac{6(k^2+(n-k)^2)}{(n+1)^2}\r)F(n,k)+T(n,k+1)-T(n,k)
\end{align}
および
\begin{align}
F(n+2,k)
&=\l(\frac{5n+7}{n+2}+\frac{6(k^2+(n+1-k)^2)}{(n+2)^2}\r)F(n+1,k)+T(n+1,k+1)-T(n+1,k)
\times(-1)^k\frac{(-n)_k^3}{(1)_{k-1}^3}\\
(k^2+(n+1-k)^2)F(n+1,k)
&=\bigg(\frac{(n+1)^2}{k^3}+\frac{n+1}{k^2}+\frac2k\\
&\qquad{}+\frac{(n+1)^2}{(n+1-k)^3}+\frac{n+1}{(n+1-k)^2}+\frac2{n+1-k}\bigg)
\times(-1)^k\frac{(-n)_k^3}{(1)_{k-1}^3}\\
&=((n+1)^2+(n+1)k+2k^2\\
&\qquad{}+(n+1)^2+(n+1)(n-k)+2(n-k)^2)\times(-1)^k\frac{(-n)_k^3}{(1)_k^3}+T'(n,k+1)-T(n,k)\\\
&=((3n+2)(n+1)+2(k^2+(n-k)^2))F(n,k)+T'(n,k+1)-T(n,k)
\end{align}
と変形できることに注意する。
いま$f(n)$が三項間漸化式を満たすと期待し、未知の多項式$p_1(n),p_2(n),p_3(n)$についての方程式
\begin{align}
&p_0(n)F(n,k)+p_1(n)F(n+1,k)+p_2(n)F(n+2,k)\\
&=\Bigg(p_0(n)+p_1(n)\l(\frac{5n+2}{n+1}+\frac{6(k^2+(n-k)^2)}{(n+1)^2}\r)\\
&\qquad{}+p_2(n)\frac{5n+7}{n+2}\l(\frac{5n+2}{n+1}+\frac{6(k^2+(n-k)^2)}{(n+1)^2}\r)\\
&\qquad{}+p_2(n)\frac6{(n+2)^2}\l((3n+2)(n+1)+2(k^2+(n-k)^2)\r)\Bigg)F(n,k)\\
&\qquad{}+T'''(n,k+1)-T'''(n,k)\\
&=T'''(n,k+1)-T'''(n,k)\\
\end{align}
を考える。$p_2(n)=(n+2)^2p_3(n)$とおくとこれは二本の方程式
\begin{align}
0&=p_1(n)+(n+2)(5n+7)p_3(n)+2(n+1)^2p_3(n)\\
0&=p_0(n)+\frac{5n+2}{n+1}(p_1(n)+(n+2)(5n+7)p_3(n))+6(3n+2)(n+1)p_3(n)\\
&=p_0(n)-2(5n+2)(n+1)p_3(n)+6(3n+2)(n+1)p_3(n)
\end{align}
に整理できるので$p_3(n)=1$とすることで
\begin{align}
p_1(n)
&=-(n+2)(5n+7)-2(n+1)^2\\
&=-(7n^2+21n+16)\\
p_0(n)&=2(5n+2)(n+1)-6(3n+2)(n+1)\\
&=-8(n+1)^2
\end{align}
と求まり、以上より
$$(n+2)^2f(n+2)-(7n^2+21n+16)f(n+1)-8(n+1)^2f(n)=0$$
を得る。