Lucas数列の関係式

同義語:Lucas数列:関係式identities for Lucas sequences

概要

Lucas数列の関係式(identities for Lucas sequences)とは、特性多項式 $X^2-PX+Q$ から定まる第1種の数列 $U_n$ と第2種の数列 $V_n$ の添字・積・平方を結ぶ恒等式群である。$Q\ne0$、判別式 $D=P^2-4Q\ne0$ のもとではBinet表示を全整数添字へ拡張でき、負添字公式、加法・差の公式、$V_n^2-DU_n^2=4Q^n$、添字の合成公式、Cassini型・Catalan型恒等式が統一的に導かれる。Fibonacci数列は $(P,Q)=(1,-1)$ の特別な場合である。

$$\newcommand{C}[0]{\mathbb{C}} \newcommand{div}[0]{\mathbin{÷}} \newcommand{N}[0]{\mathbb{N}} \newcommand{Q}[0]{\mathbb{Q}} \newcommand{R}[0]{\mathbb{R}} \newcommand{Z}[0]{\mathbb{Z}} $$

前提知識: 複素数, 漸化式, 二項定理

規約と定義

Lucas数列には複数の添字・符号規約がある。本記事では $P,Q\in\mathbb{C}$ に対し
$$ D:=P^2-4Q $$
とおき、$Q\ne0$ かつ $D\ne0$ を仮定する。多項式
$$ X^2-PX+Q $$
の相異なる二根を $\alpha,\beta$ とする。このとき
$$ \alpha+\beta=P,\qquad \alpha\beta=Q,\qquad(\alpha-\beta)^2=D $$
である。

第1種・第2種のLucas数列

第1種と第2種のLucas数列を、すべての整数 $n$ に対して
$$ U_n:=\frac{\alpha^n-\beta^n}{\alpha-\beta}, \qquad V_n:=\alpha^n+\beta^n $$
と定める。
非負添字では、これは初期値
$$ U_0=0,\quad U_1=1, \qquad V_0=2,\quad V_1=P $$
と漸化式
$$ U_{n+2}=PU_{n+1}-QU_n, \qquad V_{n+2}=PV_{n+1}-QV_n $$
による定義と一致する。

名前の由来

2 次方程式 $x^2-Px+Q=0$ の根 $a,b$ から $u_n=(a^n-b^n)/(a-b)$、$v_n=a^n+b^n$ を作り、$u_{2n}=u_nv_n$、$2u_{m+n}=u_mv_n+u_nv_m$ などの関係式と整除の性質を体系的に調べたのは É. Lucas(1870 年代)であり、この名前はそれによる。Fibonacci 数列で、ある項の平方と両隣の項の積の差が $\pm1$ であることは、それより早く 18 世紀に R. Simson が示していた(Dic05 Chapter XVII、pp. 393–396)。

$Q\ne0$ により $\alpha,\beta$ はともに零でないので、負の指数も定義できる。$D\ne0$ は $\alpha\ne\beta$、従って $U_n$ の分母が零でないことを保証する。文献によって特性多項式を $X^2-PX-Q$ とする規約もあり、その場合は本記事の $Q$ を $-Q$ に置き換える必要がある。

基本恒等式

Lucas数列の基本恒等式

任意の整数 $m,n$ に対して次が成り立つ。

  1. 負添字公式
    $$ U_{-n}=-Q^{-n}U_n,\qquad V_{-n}=Q^{-n}V_n. $$
  2. ノルム恒等式
    $$ V_n^2-DU_n^2=4Q^n. $$
  3. 加法公式
    $$ 2U_{m+n}=U_mV_n+U_nV_m, $$
    $$ 2V_{m+n}=V_mV_n+DU_mU_n. $$
  4. 差の公式
    $$ 2Q^nU_{m-n}=U_mV_n-U_nV_m, $$
    $$ 2Q^nV_{m-n}=V_mV_n-DU_mU_n. $$

$\alpha\beta=Q$ なので
$$ U_{-n} =\frac{\alpha^{-n}-\beta^{-n}}{\alpha-\beta} =\frac{\beta^n-\alpha^n}{Q^n(\alpha-\beta)} =-Q^{-n}U_n $$
であり、同様に $V_{-n}=Q^{-n}V_n$ である。
また $(\alpha-\beta)^2=D$ より
$$ V_n^2-DU_n^2 =(\alpha^n+\beta^n)^2-(\alpha^n-\beta^n)^2 =4(\alpha\beta)^n =4Q^n. $$
加法公式については、直接展開すると
$$ U_mV_n+U_nV_m =\frac{2\alpha^{m+n}-2\beta^{m+n}}{\alpha-\beta} =2U_{m+n} $$
を得る。また
$$ V_mV_n+DU_mU_n =(\alpha^m+\beta^m)(\alpha^n+\beta^n) +(\alpha^m-\beta^m)(\alpha^n-\beta^n) =2V_{m+n}. $$
加法公式で $n$ を $-n$ に置き、負添字公式を用いると
$$ 2U_{m-n} =Q^{-n}(U_mV_n-U_nV_m), $$
$$ 2V_{m-n} =Q^{-n}(V_mV_n-DU_mU_n) $$
となる。両辺に $Q^n$ を掛ければ差の公式を得る。

特に $m=n$ とおくと
$$ U_{2n}=U_nV_n, \qquad V_{2n}=V_n^2-2Q^n $$
を得る。さらに加法公式をもう一度用いると
$$ U_{3n}=U_n(V_n^2-Q^n)=U_n(DU_n^2+3Q^n), $$
$$ V_{3n}=V_n(V_n^2-3Q^n) $$
となる。

添字の合成

パラメータを明示するときは $U_n(P,Q)$、$V_n(P,Q)$ と書く。
この節では、内部に現れるLucas数列も非負添字については初期値と漸化式で定義する。この定義は内部パラメータ $(V_r(P,Q),Q^r)$ の判別式が $0$ の場合にも意味をもち、相異なる特性根を必要とするBinet表示には依存しない。

添字の積に関する合成公式

正の整数 $r$ と非負整数 $n$ に対して
$$ U_{rn}(P,Q)=U_r(P,Q)\,U_n(V_r(P,Q),Q^r), $$
$$ V_{rn}(P,Q)=V_n(V_r(P,Q),Q^r) $$
が成り立つ。

$A:=\alpha^r$、$B:=\beta^r$ とおくと
$$ A+B=V_r(P,Q),\qquad AB=Q^r. $$
従って $A,B$ は $X^2-V_r(P,Q)X+Q^r$ の根である。
$A\ne B$ の場合、Binet表示から
$$ U_n(V_r,Q^r)=\frac{A^n-B^n}{A-B}. $$
一方
$$ U_r(P,Q)=\frac{A-B}{\alpha-\beta} $$
だから、両式の積は
$$ \frac{A^n-B^n}{\alpha-\beta} =\frac{\alpha^{rn}-\beta^{rn}}{\alpha-\beta} =U_{rn}(P,Q) $$
である。また $V_n(V_r,Q^r)=A^n+B^n=V_{rn}(P,Q)$ である。
$A=B$ の場合は Binet 表示が使えないので、漸化式で示す。$S_n:=\sum_{j=0}^{n-1}A^{n-1-j}B^j$($S_0=0$)とおくと、$A,B$ の多項式としての恒等式
$$ A^n-B^n=(A-B)S_n,\qquad S_{n+2}=(A+B)S_{n+1}-AB\,S_n $$
が成り立つ(2 つ目は、両辺の $A^{n+1-j}B^j$ の係数を比べると確かめられる。$0< j< n+1$ では右辺の係数が $1+1-1=1$、$j=0,n+1$ では $1$ である)。$S_0=0$、$S_1=1$ なので、$S_n$ は $U_n(A+B,AB)=U_n(V_r,Q^r)$ と同じ初期値と漸化式をもち、両者は一致する。したがって $U_{rn}(P,Q)=\dfrac{A^n-B^n}{\alpha-\beta}=\dfrac{A-B}{\alpha-\beta}S_n=U_r(P,Q)\,U_n(V_r,Q^r)$ である(この計算は $A\ne B$ でも成り立つ)。第2種については、$A^n+B^n$ と $V_n(A+B,AB)$ がともに初期値 $2$、$A+B$ と同じ漸化式をもつので一致する。$\blacksquare$

この形には $U_r$ による除算がない。$U_r\ne0$ の場合に限り
$$ \frac{U_{rn}(P,Q)}{U_r(P,Q)}=U_n(V_r(P,Q),Q^r) $$
と書けるが、商の形を無条件に使ってはいけない。

Cassini型・Catalan型恒等式

平方差の恒等式

任意の整数 $n,r$ に対して
$$ U_{n+r}^2-Q^rU_n^2=U_rU_{2n+r}, $$
$$ V_{n+r}^2-Q^rV_n^2=DU_rU_{2n+r} $$
が成り立つ。

Binet表示を用いると
$$ \begin{aligned} D(U_{n+r}^2-Q^rU_n^2) &=(\alpha^{n+r}-\beta^{n+r})^2 -(\alpha\beta)^r(\alpha^n-\beta^n)^2\\ &=(\alpha^r-\beta^r)(\alpha^{2n+r}-\beta^{2n+r})\\ &=D\,U_rU_{2n+r}. \end{aligned} $$
$D\ne0$ なので割ることができ、1 つ目の式を得る。
同様に
$$ \begin{aligned} V_{n+r}^2-Q^rV_n^2 &=(\alpha^{n+r}+\beta^{n+r})^2 -(\alpha\beta)^r(\alpha^n+\beta^n)^2\\ &=(\alpha^r-\beta^r)(\alpha^{2n+r}-\beta^{2n+r})\\ &=DU_rU_{2n+r}, \end{aligned} $$
となり 2 つ目の式を得る。

Catalan型恒等式

任意の整数 $n,r$ に対して
$$ U_{nr}^2-U_{(n-1)r}U_{(n+1)r} =Q^{(n-1)r}U_r^2, $$
$$ V_{nr}^2-V_{(n-1)r}V_{(n+1)r} =-DQ^{(n-1)r}U_r^2 $$
が成り立つ。

正負の符号を同時に扱うため
$$ (\alpha^{nr}\pm\beta^{nr})^2 -(\alpha^{(n-1)r}\pm\beta^{(n-1)r}) (\alpha^{(n+1)r}\pm\beta^{(n+1)r}) $$
を展開する。同符号の純粋な冪は相殺し、結果は
$$ \mp Q^{(n-1)r}(\alpha^r-\beta^r)^2 =\mp DQ^{(n-1)r}U_r^2 $$
となる。負号を選んだ式を $D=(\alpha-\beta)^2$ で割れば 1 つ目の式を得て、正号を選べば 2 つ目の式を得る。

対称な加減公式

加法公式と差の公式を組み合わせると、任意の整数 $m,n$ について
$$ U_{m+n}+Q^nU_{m-n}=U_mV_n, $$
$$ U_{m+n}-Q^nU_{m-n}=U_nV_m, $$
$$ V_{m+n}+Q^nV_{m-n}=V_mV_n, $$
$$ V_{m+n}-Q^nV_{m-n}=DU_mU_n $$
を得る。3 つ目の式の右辺に余分な係数 $2$ は付かない。

行列表示と整数性

Binet 表示は特性根を使うが、漸化式は $2\times2$ 行列の冪としても書ける。行列表示は $D=0$ の場合にもそのまま使え、隣り合う項の関係式を行列式として読める。

行列表示

$P,Q\in\mathbb{C}$ とし、$U_n,V_n$($n\geq0$)を初期値と漸化式で定める。
$$ M:=\begin{pmatrix}P&-Q\\1&0\end{pmatrix} $$
とおくと、すべての $n\geq1$ について
$$ M^n=\begin{pmatrix}U_{n+1}&-QU_n\\U_n&-QU_{n-1}\end{pmatrix} $$
が成り立つ。とくに行列式をとると、$n\geq1$ について
$$ U_{n+1}U_{n-1}-U_n^2=-Q^{n-1} $$
である(Cassini 型恒等式)。また $V_n=U_{n+1}-QU_{n-1}$($n\geq1$)であり、これは $M^n$ の跡に等しい。

$n$ についての帰納法で示す。$n=1$ では $U_2=P$、$U_1=1$、$U_0=0$ なので右辺は $M$ に等しい。$n$ で成り立つとすると
$$ M^{n+1}=M^nM=\begin{pmatrix}U_{n+1}&-QU_n\\U_n&-QU_{n-1}\end{pmatrix}\begin{pmatrix}P&-Q\\1&0\end{pmatrix}=\begin{pmatrix}PU_{n+1}-QU_n&-QU_{n+1}\\PU_n-QU_{n-1}&-QU_n\end{pmatrix} $$
であり、漸化式 $U_{n+2}=PU_{n+1}-QU_n$、$U_{n+1}=PU_n-QU_{n-1}$ により $n+1$ の場合の形になる。
行列式は乗法的なので $\det M^n=(\det M)^n=Q^n$ である。一方、右辺の行列式は $-QU_{n+1}U_{n-1}+QU_n^2$ である。$Q\ne0$ なら両辺を $-Q$ で割って $U_{n+1}U_{n-1}-U_n^2=-Q^{n-1}$ を得る。$Q=0$ でも、両辺は $P,Q$ の多項式であり($U_n$ は漸化式から $P,Q$ の整数係数多項式)、$Q\ne0$ のすべての点で一致するので、多項式として等しく、$Q=0$ でも成り立つ。
最後に $W_n:=U_{n+1}-QU_{n-1}$ とおくと、$W_1=P-0=P=V_1$、$W_2=U_3-QU_1=(P^2-Q)-Q=P^2-2Q=V_2$ であり、$W_n$ は $U$ の同じ漸化式を満たす 2 つの列の一次結合なので、同じ漸化式を満たす。よって $W_n=V_n$ である。跡は $U_{n+1}-QU_{n-1}$ そのものである。$\blacksquare$

prop-lucas-sequence-identities-catalan で $r=1$ とした 1 つ目の式 $U_n^2-U_{n-1}U_{n+1}=Q^{n-1}$ と比べると、行列式から得た式と一致する。行列表示の証明は判別式 $D$ を使わないので、$D=0$ でもこの式は成り立つ。たとえば $(P,Q)=(2,1)$ で $U_n=n$ なら $(n+1)(n-1)-n^2=-1=-Q^{n-1}$ である。

整数性と $\alpha^n$ の表示

$P,Q$ が整数なら、すべての $n\geq0$ で $U_n,V_n$ は整数であり、$V_n\equiv PU_n\pmod 2$ である。さらに $Q\ne0$、$D\ne0$ のとき、すべての整数 $n$ について
$$ \alpha^n=\frac{V_n+U_n(\alpha-\beta)}{2},\qquad \beta^n=\frac{V_n-U_n(\alpha-\beta)}{2} $$
が成り立つ。$\alpha-\beta$ は $D$ の平方根の一つである。

整数性は初期値 $0,1,2,P$ が整数で、漸化式の係数 $P,-Q$ が整数であることから、$n$ についての帰納法で従う。
$V_n\equiv PU_n\pmod2$ も帰納法で示す。$n=0$ では $V_0=2\equiv0=PU_0$、$n=1$ では $V_1=P=PU_1$ である。$n,n+1$ で成り立てば
$$ V_{n+2}=PV_{n+1}-QV_n\equiv P\cdot PU_{n+1}-Q\cdot PU_n=PU_{n+2}\pmod 2 $$
である。
後半は定義から直接計算できる。$V_n+U_n(\alpha-\beta)=(\alpha^n+\beta^n)+(\alpha^n-\beta^n)=2\alpha^n$ であり、符号を変えれば $\beta^n$ の式を得る。$(\alpha-\beta)^2=D$ は冒頭で見た。$\blacksquare$

例

Pell数列

$(P,Q)=(2,-1)$ では $D=8$、$\alpha=1+\sqrt2$、$\beta=1-\sqrt2$ であり、$U_n$ は Pell 数 $0,1,2,5,12,29,\dots$、$V_n$ は $2,2,6,14,34,82,\dots$ である。整数性の命題の表示は
$$ (1+\sqrt2)^n=\frac{V_n}{2}+U_n\sqrt2 $$
となる($\alpha-\beta=2\sqrt2$)。たとえば $(1+\sqrt2)^3=7+5\sqrt2$ で、$V_3/2=7$、$U_3=5$ である。ノルム恒等式は $V_n^2-8U_n^2=4(-1)^n$ となり、$n=3$ では $196-200=-4$ である。$V_n$ はすべて偶数で、$P=2$ なので $V_n\equiv PU_n\equiv0\pmod2$ と合う。

Fibonacci数列と通常のLucas数列

$(P,Q)=(1,-1)$ とすると、$U_n$ はFibonacci数列 $F_n$、$V_n$ は通常のLucas数列 $L_n$ になる。$D=5$ なのでノルム恒等式は
$$ L_n^2-5F_n^2=4(-1)^n $$
となる。また平方差の式で $r=1$ とすれば
$$ F_{n+1}^2+F_n^2=F_{2n+1} $$
を得る。

根が2と1の場合

$(P,Q)=(3,2)$ では $\alpha=2$、$\beta=1$ と取れるので
$$ U_n=2^n-1,\qquad V_n=2^n+1. $$
例えば加法公式
$$ 2U_{m+n}=U_mV_n+U_nV_m $$
は
$$ 2(2^{m+n}-1) =(2^m-1)(2^n+1)+(2^n-1)(2^m+1) $$
という直接確認できる恒等式になる。

反例と適用条件

記事の恒等式が使う仮定と、それを外したときに崩れることを表にまとめる。

外す条件反例成り立たなくなること
$Q\ne0$$(P,Q)=(1,0)$負添字の定義と負添字公式($Q^{-n}$ が定まらない)
$D\ne0$$(P,Q)=(2,1)$相異なる根による Binet 表示($0/0$ になる)
$U_r\ne0$$(P,Q)=(0,-1)$、$r=2$商の形の合成公式 $U_{rn}/U_r=U_n(V_r,Q^r)$
特性多項式の符号の規約$X^2-X-1$ を $X^2-PX-Q$ の形で $(P,Q)=(1,1)$ と書く流儀この記事の式に $(P,Q)=(1,1)$ を代入してよいこと(正しくは $(1,-1)$)

表の 4 行目は数学的な誤りではなく規約の違いだが、ノルム恒等式 $V_n^2-DU_n^2=4Q^n$ の右辺の符号や $D=P^2-4Q$ の値が変わるので、文献の式を使うときにはまず規約をそろえる必要がある。$(P,Q)=(1,1)$ をこの記事の規約でそのまま使うと、特性多項式は $X^2-X+1$ で $D=-3$ になり、$U_n$ は $0,1,1,0,-1,-1,0,\dots$ と周期 $6$ の列になって Fibonacci 数列にならない。

$Q=0$ では負添字公式を使えない

$(P,Q)=(1,0)$ では特性根は $1$ と $0$ であり、非負添字の漸化式は定義できる。しかし $0^{-n}$ が現れるため、Binet表示を負添字へ拡張できず、$Q^{-n}$ を含む負添字公式も定義できない。これは「非負添字のLucas数列が定義できれば、負添字公式も無条件に使える」という含意を破り、仮定 $Q\ne0$ の必要性を示す。

$D=0$ では相異なる根によるBinet表示を使えない

$(P,Q)=(2,1)$ では $D=0$、特性多項式は $(X-1)^2$ である。漸化式から $U_n=n$、$V_n=2$ が得られるが
$$ \frac{\alpha^n-\beta^n}{\alpha-\beta} $$
は $0/0$ となり定義できない。これは本記事のBinet表示に $D\ne0$ が必要であることを示す。重根の場合は極限または漸化式から別に扱う必要がある。

$U_r$ による除算は常にはできない

$(P,Q)=(0,-1)$ では特性根を $1,-1$ と取れ、$D=4$ かつ $Q\ne0$ であるが
$$ U_2=\frac{1^2-(-1)^2}{1-(-1)}=0 $$
である。従って $U_{2n}/U_2$ は定義できない。一方、除算を避けた合成公式
$$ U_{2n}=U_2\,U_n(V_2,Q^2) $$
は両辺が $0$ で正しく成立する。これは商表示には $U_r\ne0$ という追加条件が必要であることを示す。

関連項目

参考文献

Mathpediaは寄付と、参考文献の書籍リンク(Amazonアソシエイト)の紹介料で運営されています。 支援について / 寄付する