文献あり

Apéryの定理とBeukersの手法

433

はじめに

　この記事ではApéryの数列
$\begin{aligned} a_{n} & = \sum_{k = 0}^{n} {(\binom{n}{k})}^{2} {(\binom{n + k}{k})}^{2} c_{n, k} \\ b_{n} & = \sum_{k = 0}^{n} {(\binom{n}{k})}^{2} {(\binom{n + k}{k})}^{2} \\ (c_{n, k} & = \sum_{m = 1}^{n} \frac{1}{m^{3}} + \sum_{m = 1}^{k} \frac{(- 1)^{m - 1}}{2 m^{3} (\binom{n}{m}) (\binom{n + m}{m})}) \end{aligned}$
の積分表示
$\begin{aligned} a_{n} ζ (3) - b_{n} & = \frac{1}{2} \int_{0}^{1} \int_{0}^{1} \frac{- \log x y}{1 - x y} P_{n} (x) P_{n} (y) d x d y \\ (P_{n} (x) & = \frac{1}{n!} \frac{d^{n}}{d x^{n}} (x (1 - x))^{n}) \end{aligned}$
について解説していきます。

あらすじ

　1978年に証明されたApéryの定理には、その翌年Beukersによって別証明が与えられました。

Apéryの定理

　 $ζ (3)$ は無理数である。

　Beukersの手法でも証明の決め手となるものは Apéryの手法と同じであり、次の補題を満たすような整数列 $p_{n}, q_{n}$ を構成することにあります。

　実数 $β$ と整数列 $p_{n}, q_{n}$ に対し
$β \neq \frac{p_{n}}{q_{n}} (\forall n \in Z_{\geq 0}), lim_{n \to \infty} (q_{n} β - p_{n}) = 0$
が成り立つとき、 $β$ は無理数である。

　ある整数 $p, q$ を用いて $β = p / q$ と表せるとすると、十分大きいある $n$ に対し
$0 < | q_{n} p - p_{n} q | < 1$
が成り立つことになるが、これは $q_{n} p - p_{n} q$ が整数であることに矛盾。よって $β$ は無理数でなければならない。

　 Apéryの手法では
$\begin{aligned} a_{n} & = \sum_{k = 0}^{n} {(\binom{n}{k})}^{2} {(\binom{n + k}{k})}^{2} c_{n, k} \\ b_{n} & = \sum_{k = 0}^{n} {(\binom{n}{k})}^{2} {(\binom{n + k}{k})}^{2} \\ (c_{n, k} & = \sum_{m = 1}^{n} \frac{1}{m^{3}} + \sum_{m = 1}^{k} \frac{(- 1)^{m - 1}}{2 m^{3} (\binom{n}{m}) (\binom{n + m}{m})}) \end{aligned}$
という有理数列を用いたのに対し、Beukersの手法ではある多項式列 $P_{n} (x)$ に対し
$\frac{1}{2} \int_{0}^{1} \int_{0}^{1} \frac{- \log x y}{1 - x y} P_{n} (x) P_{n} (y) d x d y = b_{n}^{'} ζ (3) - a_{n}^{'}$
を満たすような有理数列 $a_{n}^{'}, b_{n}^{'}$ を用いることで所望の整数列 $p_{n}, q_{n}$ を構成しました。
　そして興味深いことに、実はこの $a_{n}, b_{n}$ および $a_{n}^{'}, b_{n}^{'}$ はそれぞれ同じ数列を定めているというのです。この記事ではこの事実について簡単に解説していきます。

Apéryの手法との比較

　ちなみに上の有理数列 $a_{n}, b_{n}$ からApéryの定理を導くにあたって大まかに

$d_{n} = lcm (1, 2, 3, \dots, n)$ とおくと、 $b_{n}$ および $2 d_{n}^{3} a_{n}$ は整数である。
ある定数 $C < e^{- 3}$ が存在し、十分大きい任意の $n$ に対し
$0 < | b_{n} ζ (3) - a_{n} | < C^{n}$
つまり
$0 < | 2 d_{n}^{3} b_{n} ζ (3) - 2 d_{n}^{3} a_{n} | < 2 d_{n}^{3} C^{n} \approx (e^{3} C)^{n} \to 0$
が成り立つ。

という二つの主張を示していくこととなります。
　これらの主張について Apéryの手法では

$a_{n}, b_{n}$ を定める $\sum$ の中身に対し、それぞれ $2 d_{n}^{3}$ を掛けると整数となることを確かめる。
$a_{n}, b_{n}$ が満たす漸化式のCasoratianや特性方程式を用いることで
$| b_{n} ζ (3) - a_{n} | < \frac{6 ζ (3)}{b_{n}} \approx \frac{1}{(\sqrt{2} + 1)^{4 n}}$
と評価する。

といったアプローチを取っていたのに対し、Beukersの手法では

$P_{n} (x) P_{n} (y)$ を単項式 $A x^{r} y^{s}$ の和に分解して考えたとき、各 $r, s$ に対しある整数 $l, m$ が存在して
$\int_{0}^{1} \int_{0}^{1} \frac{- \log x y}{1 - x y} x^{r} y^{s} d x d y = l ζ (3) + \frac{m}{d_{n}^{3}}$
が成り立つことを確かめる。
積分を適当に変形することで
$\begin{aligned} 2 (b_{n} ζ (3) - a_{n}) & = \int_{0}^{1} \int_{0}^{1} \int_{0}^{1} {(\frac{x (1 - x) y (1 - y) z (1 - z)}{1 - (1 - x y) z})}^{n} \frac{d x d y d z}{1 - (1 - x y) z} \\ \leq \frac{2 ζ (3)}{(\sqrt{2} + 1)^{4 n}} \end{aligned}$
と評価する。

というアプローチを取ることになります。
　このような類似点にも注目してみると面白いかもしれません。

Beukersによる無理性証明

　まずはBeukersによるApéryの定理の証明を完成させておきましょう。ついでに $ζ (2)$ の無理性も同時に証明していきます。

補題

Beukers-Hadjicostasの公式

　非負整数 $n, r, s$ に対し
$\int_{0}^{1} \int_{0}^{1} \frac{(- \log x y)^{n}}{1 - x y} x^{r} y^{s} d x d y = {\begin{cases} \frac{n!}{r - s} \sum_{k = s + 1}^{r} \frac{1}{k^{n + 1}} & (r > s) \\ (n + 1)! (ζ (n + 2) - \sum_{k = 1}^{r} \frac{1}{k^{n + 2}}) & (r = s) \end{cases}$
が成り立つ。

　この記事の例1として示した。

補題3

$d_{r} = lcm (1, 2, 3, \dots, r)$
とおくと、非負整数 $n, r, s$ に対し
$\int_{0}^{1} \int_{0}^{1} \frac{(- \log x y)^{n}}{1 - x y} x^{r} y^{s} d x d y \in {\begin{cases} \frac{1}{d_{r}^{n + 1}} Z & (r > s) \\ \frac{1}{d_{r}^{n + 2}} Z + ζ (n + 2) Z & (r = s) \end{cases}$
が成り立つ。

　整数係数の多項式列 $P_{n} (x)$ を
$\begin{aligned} P_{n} (x) & = \frac{1}{n!} \frac{d^{n}}{d x^{n}} (x (1 - x))^{n} \\ = \sum_{k = 0}^{n} (- 1)^{k} (\binom{n}{k}) (\binom{n + k}{k}) x^{k} \end{aligned}$
によって定めると、 $C^{n}$ 級関数 $f$ に対し
$\int_{0}^{1} f (x) P_{n} (x) d x = \frac{(- 1)^{n}}{n!} \int_{0}^{1} (x (1 - x))^{n} \frac{d^{n}}{d x^{n}} f (x) d x$
が成り立つ。

　部分積分を繰り返し行うことで $k = 0, 1, 2, \dots, n$ に対し
$\int_{0}^{1} f (x) P_{n} (x) d x = \frac{(- 1)^{k}}{n!} \int_{0}^{1} (\frac{d^{n - k}}{d x^{n - k}} (x (1 - x))^{n}) \frac{d^{k}}{d x^{k}} f (x) d x$
が成り立つことがわかる。

証明

$I_{n} = \int_{0}^{1} \int_{0}^{1} \frac{(1 - y)^{n}}{1 - x y} P_{n} (x) d x d y$
とおくと
$I_{n} = (- 1)^{n} \int_{0}^{1} \int_{0}^{1} \frac{(x (1 - x) y (1 - y))^{n}}{(1 - x y)^{n + 1}} d x d y$
特に
$0 < | I_{n} | < ζ (2) \cdot (\frac{\sqrt{5} - 1}{2})^{5 n}$
が成り立つ。

証明(長いので折り畳み)

　 $f (x) = 1 / (1 - x y)$ について補題4を適用することで
$I_{n} = (- 1)^{n} \int_{0}^{1} \int_{0}^{1} \frac{(x (1 - x) y (1 - y))^{n}}{(1 - x y)^{n + 1}} d x d y$
が成り立つ。
　いま
$F (x, y) = \frac{x (1 - x) y (1 - y)}{1 - x y}$
とおくと
$\begin{aligned} F_{x} (x, y) & = \frac{1 - 2 x + x^{2} y}{(1 - x y)^{2}} y (1 - y) \\ F_{y} (x, y) & = \frac{1 - 2 y + x y^{2}}{(1 - x y)^{2}} x (1 - x) \end{aligned}$
より方程式 $F_{x} (x, y) = F_{y} (x, y) = 0$ を解くと
$x = y, 1 - x - x^{2} = 0$
つまり $F (x, y)$ は $x, y \in [0, 1]$ において
$F (\frac{\sqrt{5} - 1}{2}, \frac{\sqrt{5} - 1}{2}) = (\frac{\sqrt{5} - 1}{2})^{5}$
を最大値に持つことがわかる。
　よって
$\begin{aligned} 0 < | I_{n} | & < (\frac{\sqrt{5} - 1}{2})^{5 n} \int_{0}^{1} \int_{0}^{1} \frac{d x d y}{1 - x y} \\ = (\frac{\sqrt{5} - 1}{2})^{5 n} ζ (2) \end{aligned}$
を得る。

$J_{n} = \int_{0}^{1} \int_{0}^{1} \frac{- \log x y}{1 - x y} P_{n} (x) P_{n} (y) d x d y$
とおくと
$J_{n} = \int_{0}^{1} \int_{0}^{1} \int_{0}^{1} \frac{(x (1 - x) y (1 - y) z (1 - z))^{n}}{(1 - (1 - x y) z)^{n + 1}} d x d y d z$
特に
$0 < J_{n} < 2 ζ (3) \cdot (\sqrt{2} - 1)^{4 n}$
が成り立つ。

証明(長いので折り畳み)

$\frac{- \log x y}{1 - x y} = \int_{0}^{1} \frac{d z}{1 - (1 - x y) z}$
と変形し、 $f (x) = 1 / (1 - (1 - x y) z)$ について補題4を適用することで
$J_{n} = \int_{0}^{1} \int_{0}^{1} \int_{0}^{1} \frac{(x (1 - x) y z)^{n}}{(1 - (1 - x y) z)^{n + 1}} P_{n} (y) d x d y d z$
が成り立つ。
　ここで $z$ について
$z = \frac{1 - w}{1 - (1 - x y) w} \leftrightarrow w = \frac{1 - z}{1 - (1 - x y) z}$
と変数変換すると
$1 - w = \frac{x y z}{1 - (1 - x y) z}, \frac{d z}{1 - (1 - x y) z} = - \frac{d w}{1 - (1 - x y) w}$
より
$J_{n} = \int_{0}^{1} \int_{0}^{1} \int_{0}^{1} \frac{((1 - x) (1 - w))^{n}}{1 - (1 - x y) w} P_{n} (y) d x d y d w$
が成り立ち、 $f (y) = 1 / (1 - (1 - x y) w)$ に対し補題4を適用することで
$J_{n} = \int_{0}^{1} \int_{0}^{1} \int_{0}^{1} \frac{(x (1 - x) y (1 - y) w (1 - w))^{n}}{(1 - (1 - x y) w)^{n + 1}} d x d y d w$
を得る。
　いま
$F (x, y, z) = \frac{x (1 - x) y (1 - y) z (1 - z)}{1 - (1 - x y) z}$
とおくと
$\begin{aligned} F_{x} (x, y, z) & = \frac{(1 - 2 x) (1 - z) - x^{2} y z}{(1 - (1 - x y) z)^{2}} y (1 - y) z (1 - z) \\ F_{y} (x, y, z) & = \frac{(1 - 2 y) (1 - z) - x y^{2} z}{(1 - (1 - x y) z)^{2}} x (1 - x) z (1 - z) \\ F_{z} (x, y, z) & = \frac{1 - 2 z + (1 - x y) z^{2}}{(1 - (1 - x y) z)^{2}} x (1 - x) y (1 - y) \end{aligned}$
より方程式 $F_{x} = F_{y} = F_{z} = 0$ を解くと
$x = y, z = \frac{1}{x + 1}, 1 - 2 x - x^{2} = 0$
つまり $F (x, y, z)$ は $x, y, z \in [0, 1]$ において
$F (\sqrt{2} - 1, \sqrt{2} - 1, \frac{1}{\sqrt{2}}) = (\sqrt{2} - 1)^{4}$
を最大値に持つことがわかる。
　よって
$\begin{aligned} 0 < J_{n} & < (\sqrt{2} - 1)^{4 n} \int_{0}^{1} \int_{0}^{1} \int_{0}^{1} \frac{d x d y d z}{1 - (1 - x y) z} \\ = (\sqrt{2} - 1)^{4 n} \cdot 2 ζ (3) \end{aligned}$
を得る。

　 $ζ (2), ζ (3)$ は無理数である。

　補題3およびその系よりある整数 $p_{n}, q_{n}, p_{n}^{'}, q_{n}^{'}$ が存在して
$I_{n} = \frac{q_{n} ζ (2) - p_{n}}{d_{n}^{2}}, J_{n} = \frac{q_{n}^{'} ζ (3) - p_{n}^{'}}{d_{n}^{3}}$
と表せる。
　また素数定理から十分大きい任意の $n$ に対して
$\begin{aligned} d_{n} & = \prod_{p ∣ n} p^{⌊ \log_{p} n ⌋} \\ \leq \prod_{p ∣ n} p^{\log_{p} n} = n^{π (n)} \\ \leq n^{(1 + ε) \frac{n}{\log n}} = e^{(1 + ε) n} < 3^{n} \end{aligned}$
と評価できることに注意すると
$\begin{aligned} 0 < | q_{n} ζ (2) - p_{n} | & < ζ (2) \cdot d_{n}^{2} {(\frac{\sqrt{5} - 1}{2})}^{5 n} < \frac{9^{n}}{11^{n}} \to 0 \\ 0 < q_{n}^{'} ζ (3) - p_{n}^{'} & < 2 ζ (3) \cdot d_{n}^{3} (\sqrt{2} - 1)^{4 n} < \frac{27^{n}}{33^{n}} \to 0 \end{aligned}$
が成り立つので補題2より $ζ (2), ζ (3)$ は無理数であることが示された。

Apéryの数列との関係

　上では
$\begin{aligned} I_{n} & = (- 1)^{n} \int_{0}^{1} \int_{0}^{1} \frac{(x (1 - x) y (1 - y))^{n}}{(1 - x y)^{n + 1}} d x d y \\ = \frac{q_{n} ζ (2) - p_{n}}{d_{n}^{2}} \\ J_{n} & = \int_{0}^{1} \int_{0}^{1} \int_{0}^{1} \frac{(x (1 - x) y (1 - y) z (1 - z))^{n}}{(1 - (1 - x y) z)^{n + 1}} d x d y d z \\ = \frac{q_{n}^{'} ζ (3) - p_{n}^{'}}{d_{n}^{3}} \end{aligned}$
によって定まる整数列 $p_{n}, q_{n}, p_{n}^{'}, q_{n}^{'}$ が $ζ (2), ζ (3)$ の無理性を導くことを示しました。
　ところでこれらの被積分関数のように、離散変数 $n$ と連続変数 $x = (x_{1}, x_{2}, \dots)$ を持つ関数 $f (n, x)$ であって、その階比や偏導関数との比
$\frac{f (n + 1, x)}{f (n, x)}, \frac{\partial_{x_{1}} f (n, x)}{f (n, x)}, \frac{\partial_{x_{2}} f (n, x)}{f (n, x)}, \dots$
が $n, x_{1}, x_{2}, \dots$ についての有理関数となるようなものについて、その適当な領域 $D$ 上での積分
$I_{n} = \int_{D} f (n, x) d x$
は何らかの多項式係数の線形漸化式
$\sum_{k = 0}^{m} p_{k} (n) I_{n + k} = 0$
を満たすことが知られています。

Almkvist-Zeilberger Algorithm

　より正確にはある $n, x_{1}, x_{2}, \dots$ についての有理関数
$r_{i} (n, x) (i = 1, 2, \dots)$
が存在して
$\begin{aligned} \sum_{k = 0}^{m} p_{j} (n) f (n, x) & = \sum_{i} \partial_{x_{i}} g_{i} (n, x) \\ (g_{i} (n, x) & = r_{i} (n, x) f_{i} (n, x)) \end{aligned}$
という等式が成り立つことが知られており、逐次積分によって $g_{i} (n, x)$ が消えるような領域 $D$ においてこれを積分することで
$\sum_{k = 0}^{m} p_{k} (n) I_{n + k} = 0$
という等式を得ることができるわけです。
　またこのような多項式 $p_{k} (n)$ と有理関数 $r_{i} (n, x)$ を機械的に求めるアルゴリズムとしてAlmkvist-Zeilberger Algorithmというものが知られています。連続変数が一個の場合についてはこの記事にて解説しており、また一般の場合については参考文献AZAなどを参照してください(気が向いたら解説記事を書こうと思います)。

　この事実を上の積分
$\begin{aligned} I_{n} & = (- 1)^{n} \int_{0}^{1} \int_{0}^{1} \frac{(x (1 - x) y (1 - y))^{n}}{(1 - x y)^{n + 1}} d x d y \\ J_{n} & = \int_{0}^{1} \int_{0}^{1} \int_{0}^{1} \frac{(x (1 - x) y (1 - y) z (1 - z))^{n}}{(1 - (1 - x y) z)^{n + 1}} d x d y d z \end{aligned}$
に対して適用すると以下の主張が導かれます。

　 $I_{n}, J_{n}$ は漸化式
$\begin{aligned} (n + 1)^{2} I_{n + 1} - (11 n^{2} + 11 n + 3) I_{n} - n^{2} I_{n - 1} & = 0 \\ (n + 1)^{3} J_{n + 1} - (2 n + 1) (17 n^{2} + 17 n + 5) J_{n} + n^{3} J_{n - 1} & = 0 \end{aligned}$
を満たす。

　このことは上のアルゴリズムなどを用いることで確かめられますが、手計算で示すには計算が煩雑なので詳しくは省略します。
　また
$I_{n} = b_{n}^{'} ζ (2) - a_{n}^{'}, J_{n} = 2 (b_{n} ζ (3) - a_{n})$
なる有理数列 $a_{n}^{'}, b_{n}^{'}, a_{n}, b_{n}$ を考えると、これらは $I_{n}, J_{n}$ と同じ漸化式を満たすことがわかります(実際、そうでなければ $I_{n}, J_{n}$ の満たす漸化式から例えば
$ζ (2) = \frac{(n + 1)^{2} a_{n + 1}^{'} - (11 n^{2} + 11 n + 3) a_{n}^{'} - n^{2} a_{n - 1}^{'}}{(n + 1)^{2} b_{n + 1}^{'} - (11 n^{2} + 11 n + 3) b_{n}^{'} - n^{2} b_{n - 1}^{'}} \in Q$
が成り立つことになり $ζ (2), ζ (3)$ が無理数であることに矛盾する)。したがって
$\begin{aligned} I_{0} & = ζ (2), & I_{1} & = 3 ζ (2) - 5 \\ \frac{1}{2} J_{0} & = ζ (3), & \frac{1}{2} J_{1} & = 5 ζ (3) - 6 \end{aligned}$
に注意すると $a_{n}^{'}, b_{n}^{'}, a_{n}, b_{n}$ は Apéryの手法において用いられた数列と同じ漸化式と初期値を満たすこととなり、以下の主張が得られます( $c_{n, k}^{'}$ の定義がこの記事で紹介したものとは $2$ 倍だけ異なることに注意する)。

$\begin{aligned} a_{n}^{'} & = \sum_{k = 0}^{n} {(\binom{n}{k})}^{2} (\binom{n + k}{k}) c_{n, k}^{'} \\ b_{n}^{'} & = \sum_{k = 0}^{n} {(\binom{n}{k})}^{2} (\binom{n + k}{k}) \\ (c_{n, k}^{'} & = 2 \sum_{m = 1}^{n} \frac{(- 1)^{m - 1}}{m^{2}} + \sum_{m = 1}^{k} \frac{(- 1)^{m + n - 1}}{m^{2} (\binom{n}{m}) (\binom{n + m}{m})}) \end{aligned}$
とおくと
$b_{n}^{'} ζ (2) - a_{n}^{'} = \int_{0}^{1} \int_{0}^{1} \frac{(1 - y)^{n}}{1 - x y} P_{n} (x) d x d y$
が成り立つ。

$\begin{aligned} a_{n} & = \sum_{k = 0}^{n} {(\binom{n}{k})}^{2} {(\binom{n + k}{k})}^{2} c_{n, k} \\ b_{n} & = \sum_{k = 0}^{n} {(\binom{n}{k})}^{2} {(\binom{n + k}{k})}^{2} \\ (c_{n, k} & = \sum_{m = 1}^{n} \frac{1}{m^{3}} + \sum_{m = 1}^{k} \frac{(- 1)^{m - 1}}{2 m^{3} (\binom{n}{m}) (\binom{n + m}{m})}) \end{aligned}$
とおくと
$b_{n} ζ (3) - a_{n} = \int_{0}^{1} \int_{0}^{1} \frac{- \log x y}{1 - x y} P_{n} (x) P_{n} (y) d x d y$
が成り立つ。

余談

　上では $I_{n}, J_{n}$ が満たす漸化式からこれらがApéryの数列と結びついていることを示しましたが、実は $ζ (2), ζ (3)$ の係数がApéryの数列となることは普通に積分を計算することでも確かめられます。
　実際
$P_{n} (x) = \sum_{k = 0}^{n} (- 1)^{k} (\binom{n}{k}) (\binom{n + k}{k}) x^{k}$
より
$\begin{aligned} (1 - y)^{n} P_{n} (x) & = \sum_{r, s = 0}^{n} (- 1)^{r + s} (\binom{n}{r}) (\binom{n + r}{r}) (\binom{n}{s}) x^{r} y^{s} \\ P_{n} (x) P_{n} (y) & = \sum_{r, s = 0}^{n} (- 1)^{r + s} (\binom{n}{r}) (\binom{n + r}{r}) (\binom{n}{s}) (\binom{n + s}{s}) x^{r} y^{s} \end{aligned}$
と展開できることと
$\int_{0}^{1} \int_{0}^{1} \frac{(- \log x y)^{n}}{1 - x y} x^{r} y^{s} d x d y = {\begin{aligned} 0 \cdot ζ (n + 2) + (有理数) & (r \neq s) \\ (n + 1)! ζ (n + 2) + (有理数) & (r = s) \end{aligned}$
であったことに注意すると
$\begin{aligned} I_{n} & = \sum_{k = 0}^{n} {(\binom{n}{k})}^{2} (\binom{n + k}{k}) \cdot ζ (2) + (有理数) \\ J_{n} & = \sum_{k = 0}^{n} {(\binom{n}{k})}^{2} {(\binom{n + k}{k})}^{2} \cdot ζ (3) + (有理数) \end{aligned}$
となることがわかります。

おわりに

　以上がBeukersによる $ζ (2), ζ (3)$ の無理性証明とそのApéryの手法との関係でした。
　とりあえず個人的に把握できていることについてはまとめてみましたが、やはりまだApéryやBeukersの手法にはなぜ上手くいくのかが解明されていないところが多く、まだまだ謎は深まるばかりですね。
　ちなみに他にも
$I_{n} = \int_{[0, 1]^{k}} f (x) (x_{1} (1 - x_{1}) x_{2} (1 - x_{2}) \dots x_{k} (1 - x_{k}) g (x))^{n} d x$
という形の積分を考えることで
$C = \int_{[0, 1]^{k}} f (x) d x$
という形の定数の無理性を示せないか、ということが期待されていますが未だ目ぼしい結果は得られていないように思います。
　ゆくゆくは同様の手法による $ζ (5)$ などの無理性証明も発見されていくのでしょうか。それともこれらの手法でできることには限りがあるのでしょうか。皆さんも興味があればApéryやBeukersの手法について考察してみてはいかがでしょうか。
　とりあえず今回はこんなところで。では。

おまけ： $e$ の無理性

　ちなみにBeukersの手法は $e$ や $\log$ の無理性を示すのにも応用することができます。以下ではそのことについて紹介していきます。
　まず
$I_{n} = \int_{0}^{1} P_{n} (x) e^{1 - x} d x = \frac{1}{n!} \int_{0}^{1} (x (1 - x))^{n} e^{1 - x} d x$
という積分を考えることによって $e$ は無理数であることを示していきましょう。

　非負整数 $k$ に対し
$\int_{0}^{1} x^{k} e^{1 - x} d x = k! e - \sum_{m = 0}^{k} \frac{k!}{m!}$
が成り立つ。特に
$I_{n} \in e Z + Z$
が成り立つ。

$\int_{0}^{1} x^{k} e^{1 - x} d x = {[- \sum_{m = 0}^{k} \frac{k!}{m!} x^{m} e^{1 - x}]}_{0}^{1} = k! e - \sum_{m = 0}^{k} \frac{k!}{m!}$
とわかる。

$0 < I_{n} < \frac{e - 1}{4^{n} n!}$
が成り立つ。

　 $f (x) = x (1 - x)$ は $0 \leq x \leq 1$ において $f (1 / 2) = 1 / 4$ を最大値に持つことから
$0 < I_{n} < \frac{1}{4^{n} n!} \int_{0}^{1} e^{1 - x} d x = \frac{e - 1}{4^{n} n!}$
を得る。

　 $e$ は無理数である。

　補題11からある整数列 $a_{n}, b_{n}$ が存在して
$I_{n} = b_{n} e - a_{n}$
と表せ、また補題12から
$0 < b_{n} e - a_{n} < \frac{e - 1}{4^{n} n!} \to 0$
が成り立つので補題2より $e$ は無理数であることが示された。

$I_{n}$ の満たす漸化式

　ちなみにこの $I_{n}$ は以下の漸化式を満たします。

　 $I_{n}$ は漸化式
$I_{n + 1} + 2 (2 n + 1) I_{n} - I_{n - 1} = 0$
を満たす。またその初項は
$I_{0} = e - 1, I_{1} = 3 - e$
である。

　部分積分によって
$I_{n + 1} = \frac{1}{(n + 1)!} \int_{0}^{1} (\frac{d^{2}}{d x^{2}} (x (1 - x))^{n + 1}) e^{1 - x} d x$
と変形できることと
$\frac{d^{2}}{d x^{2}} \frac{(x (1 - x))^{n + 1}}{(n + 1)!} = - 2 (2 n + 1) \frac{(x (1 - x))^{n}}{n!} + \frac{(x (1 - x))^{n - 1}}{(n - 1)!}$
が成り立つことからわかる。

　また積分を直接計算することでその一般項は以下のように求まることがわかります。

$\begin{aligned} a_{n} & = \sum_{k = 0}^{n} (- 1)^{k} \frac{(n + k)!}{k! (n - k)!} \sum_{m = 0}^{k} \frac{1}{m!} \\ b_{n} & = \sum_{k = 0}^{n} (- 1)^{k} \frac{(n + k)!}{k! (n - k)!} \end{aligned}$
とおくと
$b_{n} e - a_{n} = \int_{0}^{1} P_{n} (x) e^{1 - x} d x$
が成り立つ。

　なお $P_{n} (1 - x) = (- 1)^{n} P_{n} (x)$ から
$I_{n} = (- 1)^{n} \int_{0}^{1} P_{n} (x) e^{x} d x$
と変形すると
$\begin{aligned} a_{n} & = (- 1)^{n} \sum_{k = 0}^{n} \frac{(n + k)!}{k! (n - k)!} \\ b_{n} & = (- 1)^{n} \sum_{k = 0}^{n} \frac{(n + k)!}{k! (n - k)!} \sum_{m = 0}^{k} \frac{(- 1)^{m}}{m!} \end{aligned}$
とも表せます。

おまけ： $\log$ の無理性

　次に
$I_{n} = \int_{0}^{1} \frac{P_{n} (x)}{1 + c x} d x = c^{n} \int_{0}^{1} \frac{(x (1 - x))^{n}}{(1 + c x)^{n + 1}} d x$
という積分を考えることによって、特定の有理数 $c$ に対し $\log (1 + c)$ が無理数となることを示していきましょう。
　余談ですがこの積分は
$I_{n} = \int_{0}^{1} \int_{0}^{1} \frac{P_{n} (x y)}{(1 + c x y) (- \log x y)} d x d y$
という二重積分に表すこともできます。
　以下 $c$ を $c > - 1$ かつ $c \neq 0$ なる有理数とし、また $r, s (s > 0)$ を $c = r / s$ なる整数とします。

　非負整数 $k$ に対し
$\int_{0}^{1} \frac{x^{k}}{1 + c x} d x = \frac{1}{(- c)^{k}} (\frac{\log (1 + c)}{c} - \sum_{m = 1}^{k} \frac{(- c)^{m - 1}}{m})$
が成り立つ。特に
$I_{n} \in \frac{\log (1 + c)}{r^{n + 1}} Z + \frac{1}{r^{n} d_{n}} Z$
が成り立つ。

証明

$\begin{aligned} \int_{0}^{1} \frac{x^{k}}{1 + c x} d x & = \frac{1}{(- c)^{k}} \int_{0}^{1} \frac{1 - (1 - (- c x)^{k})}{1 + c x} d x \\ = \frac{1}{(- c)^{k}} \int_{0}^{1} (\frac{1}{1 + c x} - \sum_{m = 1}^{k} (- c x)^{m - 1}) d x \\ = \frac{1}{(- c)^{k}} (\frac{\log (1 + c)}{c} - \sum_{m = 1}^{k} \frac{(- c)^{m - 1}}{m}) \end{aligned}$
とわかる。

$0 < | I_{n} | < \frac{(\sqrt{1 + c} - 1)^{2 n}}{| c |^{n}} \frac{\log (1 + c)}{c}$
が成り立つ。

証明

$f (x) = \frac{x (1 - x)}{1 + c x}$
とおくとこれは $0 \leq x \leq 1$ において
$f (\frac{\sqrt{1 + c} - 1}{c}) = {(\frac{\sqrt{1 + c} - 1}{c})}^{2}$
を最大値に持つので
$\begin{aligned} 0 < | I_{n} | & < {(\frac{(\sqrt{1 + c} - 1)^{2}}{| c |})}^{n} \int_{0}^{1} \frac{d x}{1 + c x} \\ < {(\frac{(\sqrt{1 + c} - 1)^{2}}{| c |})}^{n} \frac{\log (1 + c)}{c} \end{aligned}$
を得る。

$(\sqrt{1 + c} - 1)^{2} s < e^{- 1}$
を満たすような有理数 $c$ に対し $\log (1 + c)$ は無理数となる。

証明

　補題16よりある整数 $p_{n}, q_{n}$ が存在して
$I_{n} = \frac{q_{n} \log (1 + c) - p_{n}}{r^{n + 1} d_{n}}$
と表せ、また補題17から十分大きい任意の $n$ に対し
$\begin{aligned} 0 < | q_{n} \log (1 + c) - p_{n} | & < (\sqrt{1 + c} - 1)^{2 n} s^{n} d_{n} \cdot s | \log (1 + c) | \\ < ((\sqrt{1 + c} - 1)^{2} s e^{1 + ε})^{n} \to 0 \end{aligned}$
が成り立つので補題2より $\log (1 + c)$ は無理数であることが示された。

　 $c = 1$ のとき $(\sqrt{2} - 1)^{2} \cdot 1 < e^{- 1}$ より $\log 2$ は無理数となる。

$I_{n}$ の満たす漸化式

　ちなみにこの $I_{n}$ は以下の漸化式を満たします。

　 $I_{n}$ は漸化式
$(n + 1) I_{n + 1} - (\frac{2}{c} + 1) (2 n + 1) I_{n} + n I_{n - 1} = 0$
を満たす。またその初項は
$I_{0} = \frac{1}{c} \log (1 + c), I_{1} = \frac{2 + c}{c^{2}} \log (1 + c) - \frac{2}{c}$
である。

証明

　 $P_{n} (x)$ は漸化式
$(n + 1) P_{n + 1} (x) + (2 x - 1) (2 n + 1) P_{n} (x) + n P_{n - 1} (x) = 0$
を満たすことが知られている(cf. ルジャンドル多項式)、特に
$(n + 1) I_{n + 1} + n I_{n - 1} = - \int_{0}^{1} \frac{(2 x - 1) (2 n + 1) P_{n} (x)}{1 + c x} d x$
が成り立つことに注意すると
$\begin{aligned} (n + 1) I_{n + 1} - (\frac{2}{c} + 1) (2 n + 1) I_{n} + n I_{n - 1} \\ = & - (2 n + 1) \int_{0}^{1} \frac{(\frac{2}{c} + 1) + (2 x - 1)}{1 + c x} P_{n} (x) d x \\ = & - (2 n + 1) \frac{2}{c} \int_{0}^{1} P_{n} (x) d x = 0 \end{aligned}$
を得る。

　また積分を直接計算することでその一般項は以下のように求まることがわかります。

$\begin{aligned} a_{n} & = \sum_{k = 0}^{n} (\binom{n}{k}) (\binom{n + k}{k}) \frac{1}{c^{k}} \sum_{m = 1}^{k} \frac{(- c)^{m - 1}}{m} \\ b_{n} & = \sum_{k = 0}^{n} (\binom{n}{k}) (\binom{n + k}{k}) \frac{1}{c^{k}} \end{aligned}$
とおくと
$b_{n} \frac{\log (1 + c)}{c} - a_{n} = \int_{0}^{1} \frac{P_{n} (x)}{1 + c x} d x$
が成り立つ。

　特に $c = 1$ としたもの
$\begin{aligned} a_{n} & = \sum_{k = 0}^{n} (\binom{n}{k}) (\binom{n + k}{k}) \sum_{m = 1}^{k} \frac{(- 1)^{m - 1}}{m} \\ b_{n} & = \sum_{k = 0}^{n} (\binom{n}{k}) (\binom{n + k}{k}) \end{aligned}$
はこの記事のおまけとして紹介した数列に一致します。

参考文献

[1]

F. Beukers, A note on the irrationality of ζ(2) and ζ(3), Bull. London Math. Soc., 1979, 268-272

[2]

K. Alladi, M. Robinson, Legendre polynomials and irrationality, J. Reine Angew. Math., 1980, 137-155

[3]

M. Apagodu, D. Zeilberger, Multi-variable Zeilberger and Almkvist–Zeilberger algorithms and the sharpening of Wilf–Zeilberger theory, Adv. in Appl. Math., 2006

[4]

R. Dougherty-Bliss, Integral Recurrences from A to Z, Amer. Math. Monthly, 2021, 805-815

[5]

R. Dougherty-Bliss, C. Koutschan, D. Zeilberger, Tweaking the Beukers Integrals In Search of More Miraculous Irrationality Proofs A La Apéry, Ramanujan Journal, 2022, 973-994

投稿日：2024年8月5日

更新日：2024年8月5日

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。

バッチを贈って投稿者を応援しよう

バッチを贈ると投稿者に現金やAmazonのギフトカードが還元されます。

投稿者

子葉

1066

259867

主に複素解析、代数学、数論を学んでおります。私の経験上、その証明が簡単に探しても見つからない、英語の文献を漁らないと載ってない、なんて定理の解説を主にやっていきます。同じ経験をしている人の助けになれば。最近は自分用のノートになっている節があります。

他の人のコメント

コメントはありません。

読み込み中

子葉

Apéryの定理とBeukersの手法

はじめに
あらすじ
Apéryの手法との比較
Beukersによる無理性証明
補題
証明
Apéryの数列との関係
余談
おわりに
おまけ： $e$ の無理性
$I_{n}$ の満たす漸化式
おまけ： $\log$ の無理性
$I_{n}$ の満たす漸化式
参考文献

Apéryの定理とBeukersの手法

はじめに

あらすじ

Apéryの手法との比較

Beukersによる無理性証明

補題

証明

Apéryの数列との関係

Almkvist-Zeilberger Algorithm

余談

おわりに

おまけ：eの無理性

Inの満たす漸化式

おまけ：logの無理性

Inの満たす漸化式

参考文献

この記事を高評価した人

この記事に送られたバッジ

投稿者

コメント

他の人のコメント

おまけ： $e$ の無理性

$I_{n}$ の満たす漸化式

おまけ： $\log$ の無理性

$I_{n}$ の満たす漸化式