3

右辺が0でない三項間漸化式の問題(解説編)

3298
0

前回の記事 の解説をしていきたいと思います。
と言っても、これは私が最初に思いついた解法なので、無駄に遠回りをしている可能性があります。それと、前回の記事でいった(1)から(4)までのすべてに使える解法はまた別の記事に書きたいと思います。

さて、全ての問題において、結局はbn+25bn+1+6bn=0に帰着させています。やっぱり右辺が0のほうがいいんですね。

準備

これからbn+25bn+1+6bn=0をたくさん解くことになるので先に解けるところまで解いてしまします。
まず、いつもの特性方程式t25t+6=0を解いて上の漸化式を変形すると、
bn+22bn+1=3(bn+12bn)bn+23bn+1=2(bn+13bn)
となるのでb1,b2を用いて
bn+12bn=(b22b1)3n1bn+13bn=(b23b1)2n1
(上式)-(下式)よりbn=(b22b1)3n1(b23b1)2n1となります。

(1)an+25an+1+6an=12,a1=0,a2=1

方針は隣接二項間漸化式an+1=pan+qとまったく同様で、定数部分を特性方程式を用いて消すことを考えます。
つまりakに当たる部分をαに置き換えたα5α+6α=12を解いて出たα=6を用いて、与式を
(an+26)5(an+16)+6(an6)=0
と書きかえます。(これは特性方程式の原理が分かっていれば変形できます)
よって、an6={(a26)2(a16)}3n1{(a26)3(a16)}2n1
となって、a1=0,a2=1より
an=73n1132n1+6
が答えとなります。

(2)an+25an+1+6an=2n,a1=0,a2=1

これも、考え方は普通の(二項間)漸化式と同じで、都合の良い形を作ることを考えます。
つまり
(i){an+2A(n+2)B}5{an+1A(n+1)B}+6{anAnB}=0
という式を満たす定数A,Bを求めれば、右辺が0の隣接三項間漸化式として解くことができるということです。
早速展開してABを求めてみましょう。
(i)an+25an+1+6an=2An+(3A+2B)となるので、(2)の与式より
{2A=23A+2B=0
なので、A=1,B=32となります。
よって、(i){an+2(n+2)32}5{an+1(n+1)32}+6{ann32}=0となるので、ann32=bnと置けば、(i)bn+25bn+1+6bn=0となるので、これを解いて(anに戻して)、ann32=(52+252)3n1(52+352)2n1
つまり、
an=523n152n1+n+32
が答えとなります。

(3)an+25an+1+6an=5n,a1=0,a2=1

右辺が指数関数なので、その部分をまずは解消しようと思います。
an=5nbnと置くと、与式はbn+2bn+1+625bn=125,b1=0,b2=125となるので、右辺が定数、つまり(1)と同じ形にすることができました。こうなれば、あとは(1)と同様の特性方程式を解いて、α=16を得るので、
(bn+216)(bn+116)+625(bn16)=0となります。
また(隣接三項間漸化式の)特性方程式t2t+625=0を解くとt=25,35となるので、
(bn+216)25(bn+116)=35{(bn+116)25(bn16)}(bn+216)35(bn+116)=25{(bn+116)35(bn16)}
と変形して、
(bn+116)25(bn16)=350(35)n1(bn+116)35(bn16)=275(25)n1
(上式)-(下式)より
15bn=275(25)n1350(35)n1+130
よって、
an=232n1323n1+165n=132n123n+165n
が答えとなります。

(3)の別解

先ほどは指数関数の部分を無理やり定数に変えるという解法でした。こちらの別解では(2)の方法を用いて、(3)を解いてみようと思います。
まず、(2)の方法は、都合の良い形を作るというものでした。今回の都合の良い形というのは以下の形になります。
(an+2A5n+2)5(an+1A5n+1)+6(anA5n)=0
これと元の式を比較すると6A=1よりA=16となるので、
an165n={(11652)2(01651)}3n1{(11652)3(01651)2n1}
これを整理して、
an=132n123n+165n
が答えとなります。

(4)an+25an+1+6an=2n,a1=0,a2=1

これも右辺に指数関数がある形なので、(3)と同じ方針で解き進めてみましょう。
an=2nbnと置くと、与式は4bn+210bn+1+6bn=1,b1=0,b2=14となります。この漸化式の特性方程式は4α10α+6α=1となるのですが、αの係数を見ると、0になってしまいました。これはαが求まらない、つまり(1)のようなうまい変形ができない、ということを意味しています。なぜこのようなことになるかというと、an+25an+1+6an=0の特性方程式t25t+6=0の解の一つに2が含まれていることに起因します。つまり、上に出てきたbnについての漸化式を作った時、2nに次数を合わせるために係数に出てきた2bn+2なら22のように出てきたものが、t25t+6=0の解になったことでαの係数が0になってしまった((225×2+6)α=0)ということです。これは(3)でやった別解でも同様のことが起きます。
このままでは解けないので、何か別の方法を考えてみる必要があります。今、2が特性方程式t25t+6=0の解だったために、αの係数が0になってしまい、求めたいαの情報が消えてしまいました。しかし、逆に考えると、結局消えてくれるなら要らないものを増やしても問題がないのではないか?という考え方もできます。そう、nという"身代わり"を使ってαを守ろうということです。しかも、nならn+1n+2にも変化するので、αが残ってくれるだろうということも期待できます。ここまでくればあとは簡単です。an=2nbnでだめなら、an=n2nbnとしてみましょう。そうすると与式は4(n+2)bn+210(n+1)bn+1+6nbn=1,b1=0,b2=18となります。この式を展開してみるとn(4bn+210bn+1+6bn)+8bn+210bn+1=1となって、先ほどαが求められなかった原因である式の形がnにまとめられて狙い通りの形になっています。この式の特性方程式は8α10α=1となってα=12が分かります。これを使って式変形をすると、
4(n+2)(bn+2+12)10(n+1)(bn+1+12)+6n(bn+12)=0
となって、これを解くと(特性方程式うんぬんはもう省略します)
nbn=(32)n12n1
よって、
an=3nn2n12n
が答えとなります。

たった4問しかないのに、すごく長い解説になってしまいました…これは反省しないといけませんね。
こんなやり方も見つけた、という人はぜひ教えてください。では、今回はここまでにします。

投稿日:20201128
OptHub AI Competition

この記事を高評価した人

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

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

バッジはありません。
バッチを贈って投稿者を応援しよう

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

投稿者

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中