いわゆる分数型の漸化式の一般項

はじめに

分数型の漸化式というのは$\displaystyle{a_{n+1}=\dfrac{pa_n +q}{ra_n +s}}$のような形をしたものを指すらしいです．
解説サイトを見ていた時に唐突に$a_{n+1}=a_n=\alpha$とおいて方程式を解いていたことに腹が立ち，このような手法を用いずに一般項を求める方法を紹介するとともに公式化しようと思い，こうして記事にしようと思いました．

条件

今回考える漸化式は$a_{n+1}=\dfrac{pa_n +q}{ra_n +s}\ ,\ a_1=a\ ,\ r\neq 0\ ,\ ps -qr\neq 0$として，$ra_k +s=0$となる$k$が存在すれば数列$\lbrace a_n\rbrace$は第$k$項目までとします．
また，$p^2 +s^2 +4qr -2ps\neq 0$とします．$=0$の場合は公式の見やすさのため割愛しました．よければ計算してみてください．

解法

漸化式から，$a_n=\dfrac{p_ka_{n-k} +q_k}{r_ka_{n-k} +s_k}$となるような数列$\lbrace p_k\rbrace ,\lbrace q_k\rbrace ,\lbrace r_k\rbrace ,\lbrace s_k\rbrace$が存在します．この4つの数列の一般項を求めることができれば，数列$\lbrace a_n\rbrace$の一般項は$\displaystyle{a_1=a\ ,\ a_n=\dfrac{p_{n-1}a +q_{n-1}}{r_{n-1}a +s_{n-1}}\ (n\gt 1)}$となります．

$\dfrac{p_ka_{n-k} +q_k}{r_ka_{n-k} +s_k}=\dfrac{p_k\dfrac{pa_{n-k-1} +q}{ra_{n-k-1} +s} +q_k}{r_k\dfrac{pa_{n-k-1} +q}{ra_{n-k-1} +s} +s_k}=\dfrac{(pp_k +rq_k)a_{n-k-1} +(qp_k +sq_k)}{(pr_k +rs_k)a_{n-k-1} +(qr_k +ss_k)}$

より，ある0でない実数$m$が存在して
\begin{eqnarray} \left\{ \begin{array}{l} p_{k+1}=m(pp_k +rq_k)\ ,\ p_1=p\\ q_{k+1}=m(qp_k +sq_k)\ ,\ q_1=q \end{array} \right. \end{eqnarray}

\begin{eqnarray} \left\{ \begin{array}{l} r_{k+1}=m(pr_k +rs_k)\ ,\ r_1=r\\ s_{k+1}=m(qr_k +ss_k)\ ,\ s_1=s \end{array} \right. \end{eqnarray}

となります．少し計算するとわかるように$m$は最終的に約分されて消えるので$m=1$で良いでしょう．また，2つの連立漸化式を見比べたときに初項が異なるだけだとわかるので，第1式のみを考えます．

$p_{k+2}=(p +s)p_{k+1} +(qr -ps)p_k$

となります．

$p +s=u +v\ ,\ qr -ps=-uv$となるように複素数$u,v$をとると次のように式変形できます．
\begin{align} p_{k+2} +up_{k+1}&=v(p_{k+1} +up_k) \\ &=v^k(p_2 +up_1) \\ p_{k+1} +up_k&=v^{k-1}(p_2 +up_1)\ \ (k\gt 1) \end{align}
同様にして，
$$ p_{k+1} +vp_k=u^{k-1}(p_2 +vp_1)\ \ (k\gt 1) $$
を得ます．これらは$k=1$のときでも成り立ちます．$u,v$の条件から
$$ u=\dfrac{p+s +\sqrt{(p +s)^2 -4(ps -qr)}}{2}\ ,\ v=\dfrac{p+s -\sqrt{(p +s)^2 -4(ps -qr)}}{2} $$
とわかるので，$p^2 +s^2 +4qr -2ps\neq 0$より$u\neq v$です．
さて，
\begin{eqnarray} \left\{ \begin{array}{l} p_{k+1} +up_k=v^{k-1}(p_2 +up_1) \\ p_{k+1} +vp_k=u^{k-1}(p_2 +vp_1) \end{array} \right. \end{eqnarray}
とわかりました．第1式から第2式を引いて両辺を$u-v$で割ることで，以下の式を得ます．
$$ p_k=\dfrac{1}{u-v}\Big(v^{k-1}(p_2 +up_1) -u^{k-1}(p_2 +vp_1)\Big) $$
また漸化式から，
\begin{align} q_k&=\dfrac{1}{r} (p_{k+1} -pp_k) \\ &=\dfrac{1}{u-v}\dfrac{1}{r}\Big(v^k(p_2 +up_1) -u^k(p_2 +vp_1) -pv^{k-1}(p_2 +up_1) +pu^{k-1}(p_2 +vp_1)\Big) \\ &=\dfrac{1}{u-v}\dfrac{1}{r}\Big((v-p)v^{k-1}(p_2 +up_1) -(u-p)u^{k-1}(p_2 +vp_1)\Big) \end{align}
とわかります．ここで$p_1=p\ ,\ p_2=p^2 +qr\ ,\ r_1=r\ ,\ r_2=(p+s)r$なので，それぞれ代入して
\begin{eqnarray} \left\{ \begin{array}{l} p_k=\dfrac{1}{u-v}\Big(v^{k-1}(p^2 +qr -up) -u^{k-1}(p^2 +qr -vp)\Big) \\ q_k=\dfrac{1}{u-v}\dfrac{1}{r}\Big((v-p)v^{k-1}(p^2 +qr -up) -(u-p)u^{k-1}(p^2 +qr -up)\Big) \\ r_k=\dfrac{1}{u-v}\Big(v^{k-1}(p +s -u)r -u^{k-1}(p +s -v)r\Big) \\ s_k=\dfrac{1}{u-v}\Big((v -p)v^{k-1}(p +s -u) -(u -p)u^{k-1}(p +s -v)\Big) \end{array} \right. \end{eqnarray}
となります．これで各数列の一般項を求めることができました．最後に仕上げです．
\begin{align} a_n&=\dfrac{p_{n-1}a +q_{n-1}}{r_{n-1}a +s_{n-1}} \\ &=\dfrac{\Big(v^{n-2}(p^2 +qr -up) -u^{n-2}(p^2 +qr -vp)\Big)a +\dfrac{1}{r}\Big((v-p)v^{n-2}(p^2 +qr -up) -(u-p)u^{n-2}(p^2 +qr -vp)\Big)}{\Big(v^{n-2}(p +s -u)r -u^{n-2}(p +s -v)r\Big)a +\Big((v-p)v^{n-2}(p +s -u) -(u-p)u^{n-2}(p +s -v)\Big)} \\ &=\dfrac{1}{r}\dfrac{u^{n-2}(p^2 +qr -vp)(ar +u -p) -v^{n-2}(p^2 +qr -up)(ar +v -p)}{u^{n-2}(p +s -v)(ar +u -p) -v^{n-2}(p +s -v)(ar +v -p)} \\ &=\dfrac{1}{r}\dfrac{u^{n-2}(p(p -v) +qr)(ar +u -p) -v^{n-2}(p(p -u) +qr)(ar +v -p)}{u^{n-1}(ar +u -p) -v^{n-1}(ar +v -p)} \\ &=\dfrac{1}{r}\dfrac{u^{n-2}(p(u -s) +qr)(ar +u -p) -v^{n-2}(p(v -s) +qr)(ar +v -p)}{u^{n-1}(ar +u -p) -v^{n-1}(ar +v -p)} \\ \therefore a_n&=\dfrac{1}{r}\dfrac{u^{n-1}(p -v)(ar +u -p) -v^{n-1}(p -u)(ar +v -p)}{u^{n-1}(ar +u -p) -v^{n-1}(ar +v -p)}\ \ \ (n\gt 1) \end{align}
$n=1$のときを考えると，
\begin{align} a_1&=\dfrac{1}{r}\dfrac{(p -v)(ar +u -p) -(p -u)(ar +v -p)}{(ar +u -p) -(ar +v -p)} \\ &=\dfrac{1}{r}\dfrac{p(u -v) +(u -v)(ar -p)}{u -v} \\ &=a \end{align}
となり$n=1$でも成り立ちます．

以上のような一般項の求め方はとても自然で難しくありません．まとめると次のようになります．

最後に

一般項は$u,v$の満たす性質を用いてもう少し変形することができます．
\begin{align} a_n&=\dfrac{1}{r}\dfrac{u^{n-1}(p -v)(ar +u -p) -v^{n-1}(p -u)(ar +v -p)}{u^{n-1}(ar +u -p) -v^{n-1}(ar +v -p)} \\ &=\dfrac{1}{r}\dfrac{u^{n-1}\big(ar(p -v) +(p -v)(u -p)\big) -v^{n-1}\big(ar(p -u) +(p -u)(v -p)\big)}{u^{n-1}(ar +u -p) -v^{n-1}(ar +v -p)} \\ &=\dfrac{1}{r}\dfrac{u^{n-1}\big(ar(p -v) +qr\big) -v^{n-1}\big(ar(p -u) +qr\big)}{u^{n-1}(ar +u -p) -v^{n-1}(ar +v -p)} \\ &=\dfrac{u^{n-1}(ap -av +q) -v^{n-1}(ap -au +q)}{u^{n-1}(ar +u -p) -v^{n-1}(ar +v -p)} \end{align}

となります．この式は$r=0$の場合にも対応しており(実際，$r=0\ ,\ s=1$として$u=p\ ,\ v=1$とおくと$a_n=p^{n-1}\left(a +\dfrac{q}{p -1}\right) -\dfrac{q}{p -1}$となり，これは漸化式$a_{n+1}=pa_n +q\ ,\ a_1=a$の一般項になっている)とても便利なのですが，式の見やすさの観点から敢えてそのままにしました．

ここまで読んでいただきありがとうございます．