超幾何数列の基礎5：WZ method

はじめに

　この記事では 前回の記事 に引き続き超幾何数列の基本事項についてまとめていきます。
　前回の記事ではproperな超幾何数列$F(n,k)$に対して
$$\sum _{i=0}^I{p}_i(n)F(n+i,k)=G(n,k+1)-G(n,k)$$
というタイプのcreative telescopingができることを示しました。そして今回の記事ではこの最たる例として
$$F(n+1,k)-F(n,k)=G(n,k+1)-G(n,k)$$
というタイプのcreative telescopingができる場合に関する手法Wilf-Zeilberger method、略してWZ methodについて解説していきます。

導入

　いま(properな)超幾何数列$A(n),F(n,k)$についての等式
$$\sum _{k}F(n,k)=A(n)$$
があったとしましょう。このときこの右辺は一階の漸化式
$$f(n+1)-\frac{A(n+1)}{A(n)}f(n)=0$$
を満たすので、Zeilberger's methodによって得られる
$$\sum _{i=0}^I{p}_i(n)F(n+i,k)=G(n,k+1)-G(n,k)$$
という$F(n,k)$に関する漸化式も
$$F(n+1,k)-\frac{A(n+1)}{A(n)}F(n,k)=G(n,k+1)-G(n,k)$$
のような形になると期待できます。
　そして、驚くべきことに、実際それが多くの場合において成り立つということが経験的な事実として知られています。特にこのとき
$$\bar{F}(n,k)=\frac{F(n,k)}{A(n)},\bar{G}(n,k)=\frac{G(n,k)}{A(n+1)}$$
とおくと
$$\bar{F}(n+1,k)-\bar{F}(n,k)=\bar{G}(n,k+1)-\bar{G}(n,k)$$
という等式が成り立つことに注意しましょう。

WZ method

　このように
$$F(n+1,k)-F(n,k)=G(n,k+1)-G(n,k)$$
を満たすような超幾何数列の組$F,G$を構成することで
$$\sum _{k}F(n,k)=1$$
のような等式を示す手法のことをWZ methodと言います。

　ちなみに$F(n,k)$に対してパートナーとなる$G(n,k)$のことをWZ-mateと言うことがあります。

WZ-pairの構成法

　いま与えられた超幾何数列$F(n,k)$に対してWZ methodが適用できるのかどうか、そしてそのときWZ-mate$G(n,k)$はどうやって求められるのかについては
$$A(k)=F(n+1,k)-F(n,k)$$
に対して Gosper's algorithm を実行するだけで解決します。
　また上での議論を踏まえると与えられた超幾何数列$F(n,k)$に対してあるWZ-pair$\bar{F},\bar{G}$を構成する方法として次のようなものが考えられます。

　またWZ-pairの構成に関するより一般的な方法としてMarkov-WZ methodというものも知られていますがそれについては 次回の記事 にて詳しく解説しようと思います。

計算例

応用法1：級数の求値

　上でも言及したようにWZ methodの基本的な使い方として
$$f(n)=\sum _{k}F(n,k)=\mathrm{Const}.$$
というタイプの等式を示すことが挙げられます。
　このような等式は
$$F(n+1,k)-F(n,k)=G(n,k+1)-G(n,k)$$
の両辺を$k$について足し合わせたとき、右辺の和が
$$\sum _{k=-\mathrm{\infty }}^{\mathrm{\infty }}(G(n,k+1)-G(n,k))=G(n,\mathrm{\infty })-G(n,-\mathrm{\infty })=0$$
となることを仮定すれば
$$f(n+1)-f(n)=0$$
が成り立つことから導けます。

Carlsonの定理

　ところでproperな超幾何数列とは
$$F(n,k)=P(n,k)\frac{\prod _{j=1}^s\mathrm{\Gamma }(a_{j}n+b_{j}k+c_j)}{\prod _{j=1}^t\mathrm{\Gamma }(a_j^{\prime }n+b_j^{\prime }k+c_j^{\prime })}{x}^n{y}^k$$
と表せるようなもののことを言うのでした。ここで適当な条件下においてこの右辺を$n\in \mathbb{C}$についての複素関数に拡張すると上の結果
$$\sum _{k}F(n,k)=\mathrm{Const}.(n=0,1,2,\dots )$$
から
$$\sum _{k}F(n,k)=\mathrm{Const}.(n\in \mathbb{C})$$
を導くことができます。
　これはCarlsonの定理から直ちに従います。

　いま$n$を複素数とできることによって特殊値が確認しやすくなるという嬉しさが発生します。
　例えば$F(n,k)$が
$$\frac{1}{\mathrm{\Gamma }(n-k+c)}(c\notin \mathbb{Z})$$
という因子を持つとき、$n\in \mathbb{Z}$の範囲であればこれはうんともすんとも言いませんが、$n\in \mathbb{C}$の範囲であれば$n=1-c$とすることでこれは$k\ge 1$において$0$となるので
$$\sum _{k=0}^{\mathrm{\infty }}F(1-c,k)=F(1-c,0)$$
と求まり、したがって任意の$n\in \mathbb{C}$に対して
$$\sum _{k=0}^{\mathrm{\infty }}F(n,k)=F(1-c,0)$$
が成り立つ、といった議論をすることができます。

応用法2：級数の変形

　またWZ methodのもう一つの使い方として級数同士の等式
$$\sum _{n}G(n,k)=\sum _n{G}^{\prime }(n,k)$$
を量産するということが挙げられます。

離散版Greenの定理

　これにはWZ-pairのなす差分形式というものがいわゆる"渦なし"であることが重要となります。

　このことは渦なし場
$$\mathit{r}=(v(x,y),u(x,y))(\frac{\partial u}{\partial x}=\frac{\partial v}{\partial y})$$
においてGreenの定理より
$${\oint }_C\mathit{r}\cdot d\mathit{x}={\oint }_C(vdx+udy)=0$$
が成り立つことの類似となっています。
　ちなみにより一般的な事実として 離散版Stokesの定理 というものも知られています。

経路の取り換え

　いま適当な閉経路における和分を考えることで例えば次のような変換公式が得られます。

　またこの階段の各辺の長さを変えることでこの公式は次のように一般化できます。

　例えば
$$\begin{array}{rl}\sum _{n=0}^{\mathrm{\infty }}G(n,0)& =\sum _{n=0}^{\mathrm{\infty }}(F(n,n)+G(n,n+1))\\ & =\sum _{n=0}^{\mathrm{\infty }}(F(n,2n)+F(n,2n+1)+G(n,2n+2))\\ & =\sum _{n=0}^{\mathrm{\infty }}(F(n,3n)+F(n,3n+1)+F(n,3n+2)+G(n,3n+3))\\ & =\sum _{n=0}^{\mathrm{\infty }}(F(2n,2n)+F(2n,2n+1)+G(2n,2n+2)+G(2n+1,2n+2))\end{array}$$
のような式が成り立ちます。
　他にも様々な経路の取り方を考えることで色々な変換公式が得られます。例えば次のような有限和の等式なども興味深いですね。

具体例

　さて御託を並べるのはこのくらいにして、ここからは実際の応用例について見ていくこととしましょう。

ラマヌジャンの円周率公式

　まずはラマヌジャンの円周率公式、の中でも簡単なもの
$$\frac{2}{\pi }=\sum _{k=0}^{\mathrm{\infty }}(-1{)}^k(4k+1)\frac{(\frac{1}{2}{)}_k^3}{(1{)}_k^3}$$
を証明してみましょう。
　これはEkhad-Zeilbergerによって次のように一般化されています。

　ちなみにその他の円周率公式にも同様の一般化が得られており、それらの公式については この記事 にてそれとなくまとめています。

$\zeta (2),\zeta (3)$の加速級数

　次にApéryによる$\zeta (3)$の無理性の証明に使われたことでも有名な公式

をWZ methodによって証明してみましょう。

　また類似のWZ-pairを用いることで$\zeta (2)$に関する次のような公式が得られます。

　さらに上の問題2として求めたWZ-pair
$$\begin{array}{rl}H(n,k)& ={\left(\frac{n!k!}{(n+k)!}\right)}^2\\ F(n,k)& =\frac{2(2n-1)}{n^3(\genfrac{}{}{0ex}{}{2n}{n})}H(n,k)\\ G(n,k)& =\frac{2k+3n}{n^3(\genfrac{}{}{0ex}{}{2n}{n})}H(n,k)\end{array}$$
を用いると次のような公式が導かれます。

おわりに

　以上がWZ methodの概要となります。
　見ての通りWZ methodは自由変数のある与えられた級数
$$\sum _{k=0}^{\mathrm{\infty }}F(n,k),\sum _{n=0}^{\mathrm{\infty }}G(n,k)$$
の値を求めたり加速級数に変形したりするのには使えますが、自由変数のない級数
$$\sum _{k=0}^{\mathrm{\infty }}A(k)$$
に対して適当なWZ-pairを探し出し、その性質を調べようとするのには向いていないと思います。
　ただそのような点に目を瞑れば、WZ methodは主に超幾何関数にまつわる多くの等式を(機械的に)示す上で非常に重要な手法であり、また無作為に取ってきた超幾何数列$H$に対してWZ-pair$F,G$を構成してみるだけでも色々と非自明な等式が現れて面白いので、皆さんも色々なWZ methodの例に触れて遊んでみてはいかがでしょうか。
　私もWZ methodについてそこまで詳しいわけではないので、また面白い話を見かけたら続きの記事でも書こうと思います。とりあえず個人的に掲げていた「WZ methodについて解説する」という目標は達成できたので今回のシリーズはこんなところで。
　では。

おまけ：ポテンシャル

　ところでWZ form
$$\omega =F(n,k)\delta k+G(n,k)\delta n$$
とはいわゆる"渦なし"の差分形式であったのでこれに関するポテンシャルという概念を考えることができます。

　特に経路$C:(0,0)\to (n,k)$における不定和分が定める関数
$$c(n,k)=\sum _C\omega$$
は$\omega$のポテンシャルの一つとなります。

WZ-pairの変換

　また逆に任意に数列$c(n,k)$を持ってきたとき、その"勾配"
$$\begin{array}{rl}F(n,k)& =c(n,k+1)-c(n,k)\\ G(n,k)& =c(n+1,k)-c(n,k)\end{array}$$
を取るとこれらは
$$F(n+1,k)-F(n,k)=G(n,k+1)-G(n,k)$$
を満たします。
　特に$c(n,k)$をあるWZ form$\omega$のポテンシャルとし
$$c^{\prime }(n,k)=c(pn+qk,rn+sk),(p,q,r,s\in \mathbb{Z})$$
という関数の"勾配"を考えることで新たなWZ-pairを構成することができます。
　例えば
$$c^{\prime }(n,k)=c(n,n+k)$$
とおくと
$$\begin{array}{rl}{c}^{\prime }(n,k+1)-c^{\prime }(n,k)& =c(n,n+k+1)-c(n,n+k)\\ & =F(n,n+k)\\ c^{\prime }(n+1,k)-c^{\prime }(n,k)& =c(n+1,n+k+1)-(c(n+1,n+k)-c(n+1,n+k))-c(n,n+k)\\ & =F(n+1,n+k)+G(n,n+k)\end{array}$$
より
$$\begin{array}{rl}{F}_1(n,k)& =F(n,n+k)\\ G_1(n,k)& =F(n+1,n+k)+G(n,n+k)\end{array}$$
というWZ-pairが得られます。
　だから何だという感じではありますが、たまに役に立つので豆知識程度に覚えておくといいと思います。