6
競技数学解説
文献あり

フィボナッチ数列の応用いろいろ

1229
0

フィボナッチ数は「指数-指数」と同じようなふるまいをすることに前回までの記事で触れてきました. ので, 「指数-指数」で有名な性質, LTEの補題とZsigmondyの定理を中心に, フィボナッチ数列の応用例をいくつか書いていきます!

証明中の記法など

  • Fnn番目のフィボナッチ数. F1=1,F2=1,F3=2,
  • a|babの約数である.
  • (ab)二項係数aCb
  • φ(n)Eulerのトーシェント関数, nと互いに素な1以上n以下の整数の個数
  • vp(n)整数nが素数pで割り切れる回数の最大値

前回の記事 およびそこから飛べる過去の記事の結果を色々使います.

まずはLTEから.

LTEの類似

前回の記事 のなかで扱ったフィボナッチ数列におけるLTEの補題の類似を示していきます. といっても, 証明は普通のLTEとほとんど同じです. (そして じゅんにーさんのこの記事 が以下の結果を包含しています. )

追記 OMC040Fの解説 にまったく同じ手法が用いられていることを発見しました. 議論の流れ上、該当部分の削除はいたしませんが、その点に十分ご留意ください

任意の正整数k, 整数nに対して以下の式が成立する.
2Fkn=(k1)FnLn(k1)5(k2)Fn2Fn(k2)+5(k3)Fn3Ln(k3)52(k4)Fn4Fn(k4)++5i1(k2i1)Fn2i1Ln(k(2i1))5i(k2i)Fn2iFn(k2i)+
ここで, L0=2,L1=1,Ln+2=Ln+1+Ln である(リュカ数列, ルカス数列と呼ばれるもの).
なお, 上の式は有限和である.

概略

α=1+52,β=152とおく. ((αkβk)+βk)nαknを二項定理で展開し, βn=Ln5Fn2,Fn=15(αnβn)を代入したのち5の係数を比較することにより所望の式を得る.

LTEの類似

pを奇素数, p|Fnとするとき
vp(Fkn)=vp(Fn)+vp(k)
が成立する

補題1の右辺第一項について考える.
p|Fnよりp|Fn(k1)であるが, もしp|Ln(k1)だと仮定するとp|5Fn(k1)2Ln(k1)2よりp|4であるが, これは不適. 従って,
vp(kFnLn(k1))=vp(Fn)+vp(k)
である. また, 右辺第二項以降について,
vp((ki)Fni)=ivp(Fn)+vp((ki))=ivp(Fn)+vp(ki(k1i1))i1+vp(Fn)+vp(k)vp(i)>vp(Fn)+vp(k)
であるので示された.

p=2は前回の方法で頑張りましょう.

使ってみよう!

LTEを使ってみます.

素数p, 正の整数nに対してFnpFn(N)の素因数をxp1x1の素因数と同じように考察してみます.

奇素数q|FnpFnを任意にとる.

  • qFnを割り切るとき
    定理2よりvq(Fpn)=vq(Fn)+vq(p)なので, p=qが従う.
  • qFnを割りらないとき
    Fnp,Fqorq±1のどちらもqの倍数なので, q|Fgcd(np,qorq±1)であり, 場合分けの仮定よりgcd(np,qorq±1)Fnである. これよりgcd(np,qorq±1)n, 従ってp|qorq±1である. つまりqmodp±1,0のいずれかに合同である.

また, 2|FnpFnとなるのは前回の補題よりp=2,3の場合のみなので結局FnpFnの素因数はmodp±1,0のいずれかに合同なものだけである. 特にすべてp1以上である.
多項式の場合と似た結果が得られました!

4で割って1余る素数の評価

整数コン 5についてです.

上の結果において, p13以上の素数, n=1とします. このとき, Fpp1以上の素因数しか持たず, pqの偶奇の考察およびp5なのでこれはp+1以上です. さらに,
Fp=F(p+1)/22+F(p1)/22
が成立し, F(p+1)/2F(p1)/2が互いに素であることを併せれば, Fpの素因数は全て4で割って1余ります(一般にa2+b2の素因数の議論). 従って, p=5の場合も考慮に入れれば, すべての4で割って1余る素数pに対してp+1以上Fp+8以下の, 4で割って1余る素数が存在することが言えました.

指数を含む不定方程式への利用

l(5)=5,F5=5であることと, LTEよりv5(Fn)=v5(n)が分かります. これより, Fn=5mなるn,mが存在したとき, v5(n)=m,n5mとなり, F5mFn=5mとなります. これは, m2で不成立なので, 5べきのフィボナッチ数は5だけです.
ところで, これにフィボナッチ数の判定法を用いると,
52m+1=k2±4
の解を求めたことになります. 52m=k2±4の解は, 平方数-平方数の形にできるので, 結局5m=k2±4の整数解を求められたことになります. 特に, 53=112+4という解はかなり非自明ですごいですよね. ということは, 逆に考えればこういった方程式は三項間漸化式に帰着させて解くことができそうな気がします. やってみましょう.

3m=n2+2を満たす正の整数n,mを全て求めよ.

これも33=52+2という解があるのでやばそうです.

mが偶数の場合は簡単. mが奇数の場合は3(m1)/2mで置きなおして, 3m2=n2+2になる.
ここで, an+2=can+1+an,a0=0,a1=1についても, フィボナッチ数列同様の乗法定理が成立する:
an+manm=an2(1)n+mam2
m=1として, 漸化式を代入すると
canan1+an12an2(1)n=0
an1についての判別式を見ると, (c2+4)an2+4(1)nが平方数になる. さて, これを利用したいのだが定数項の絶対値が4になってしまい困った. ので, 元の式から強引に4を作り出そう.
2nNとおけば, 6m2±4=N2となった. よって, c=2とすればよさそう. 数列を書いていくと,
0,1,2,3,42,11,152,41,562,153
となる. ここで注目すべきは, 奇数項目である. 偶数項目が2の整数倍であることから, 判別式(22+4)an24が平方数の2倍であり, N2の整数倍であることと整合性がとれた!

一見常にこのようにうまくいくように見えるが, 整合性がとれたのはただの数値設定上の偶然である.
具体的には,

  • 定数項が4の場合 (k+4)m2=kn2±4
  • 定数項が2の場合 (k+2)m2=kn2±2
  • 定数項が1の場合 (k+1)m2=kn2±1

の場合がこの方法で処理できる(多分). そのほかの場合にはより一般のペル方程式の議論を要する.

つまり, 3m2=n2+2の整数解は数列anの奇数項目で与えられる. 逆も, 詳細は省略するが, この記事の最後 っぽいことをすれば示せる.
さて, 上で考えた anについて, v3(an)=v3(n)が言える.

a3n=1an(a2n2(1)nan2)=an((an+1+an1)2(1)n)
あとはmod9頑張る(か, 冒頭のような方法を使う)

ここから, 先の例と同じように不等式評価をするとこの数列の奇数項目に現れる3べきは1,3だけだと分かった. 以上より, もとの方程式の解もしぼれ, (1,1),(3,5)のみだと分かる.


といった感じで, 一見難しそうな不定方程式を解くことができました. 常にというわけにはいかなさそうですが...

最後にZsigmodyの類似の紹介です.

Zsigmondyの定理の類似

INTEGERSのこの記事 の手法を応用し, フィボナッチ数列でも類似の結果を得られたので証明をかいてみます. なお, 実際はより一般の数列に対して成立することが知られています( この記事 で結果が紹介されています).

Zsigmondyの類似

任意の自然数n1,2,6,12に対し, Fnの素因数pであって, pF1,F2,Fn1を割り切らないものが存在する. このような素因数を原始的素因数と呼ぶ(上で紹介した記事における原始的約数とは若干異なるので注意).

この証明をしていきます.
証明の大まかな流れとしては, Fnの約数のうち, 原始的素因数でないものの積を上から評価するといった感じです. 実際には円分多項式にあたる概念, 円分フィボナッチ数を定義し, それについての評価をしていきます.


周期関数( この記事 で紹介)

正整数nに対し, Fmnの倍数になる正整数mが存在するので, このうち最小のものをl(n)とかく.

Fmnの倍数のとき, ml(n)の倍数である.

証明は演習とする.

任意のkvp(Fl(p)), 奇素数pに対し, l(pk)=pkvp(Fl(p))l(p), 任意のk3に対してl(2k)=2k23が成立する.

補題2より, Fmpkの倍数のときml(p)の倍数なので, LTEの類似よりkvp(Fm/l(p))=vp(m)+vp(Fl(p))が成立. これより, ml(p)pkvp(Fl(p))となるので示された(逆は同様にLTEを使おう). p=2の場合も同様.


メビウス関数

μをメビウス関数, すなわち正整数nに対しnが平方因子を持つ場合μ(n)=0, それ以外の場合nの素因数の個数tについてμ(n)=(1)tとする.

メビウス関数には以下の重要な性質がある.

メビウスの反転公式

正の整数に定義される関数(厳密には終域は任意のアーベル群)f,g
g(n)=d|nf(n),f(n)=d|ng(n)μ(n/d)
のそれぞれを満たすことは同値である.

円分フィボナッチ数

正の整数nに対し, 円分フィボナッチ数列F~nを以下で定める:
F~n=d|nFnμ(n/d)

F~1,F~2,...1,1,2,3,5,4,13,7,17,11,89,6,となる.

円分フィボナッチ数の性質

2以上の整数nに対して,
F~n=p|Fn:素数pmax(0,vp(Fn)max1kn1vp(Fk))=gcd(m,n)=1,1mn(1+52ζnm152)
が成立する(n=1は左の=のみ成立). ここで, ζn=e2πi/n1の原始n乗根の一つである. 特に, F~nは整数である.

(1つ目の=)
素数qp|Fn:素数pmax(0,vp(Fn)max1kn1vp(Fk))を割り切るのは, n=l(q),l(q2),l(q3),...の場合のみであることに注意する.
補題6より, d|np|Fd:素数pmax(0,vp(Fd)max1kd1vp(Fk))
を示せばよいが, 補題4より, l(q),l(q2),,l(qvp(Fn))はいずれもnを割り切るため,
vq(d|np|Fd:素数pmax(0,vp(Fd)max1kd1vp(Fk)))=vq(qvq(Fn))=vq(Fq)
である. よって示された.

(2つ目の=)
α=(1+5)/2,β=(15)/2とおく. Binetの公式より,
F~n=d|nFnμ(n/d)=d|n15(αnβn)μ(n/d)=Φ(α,β)15d|nμ(n/d)=Φ(α,β)=gcd(m,n)=1,1mn(αζnmβ)
ここで, Φは斉次の円分多項式であり, 最後から三つ目の=には一般の円分多項式の性質, 最後から二つ目の=にはd|nμ(d)=0を用いた.

円分フィボナッチ数自体を操作することが円分多項式に比べて難しいため, あらかじめ素因数で記述した. そして, さらに以下の補題によりより明示的に記述できる.

円分フィボナッチ数の性質

奇素数p, 正整数nに対して, 以下が成立する.
max(0,vp(Fn)max1kn1vp(Fk))={vp(Fl(p))(n=l(p))1(n=pl(p),p2l(p),...)0(Otherwise),

max(0,v2(Fn)max1kn1v2(Fk))={1(n=3,223,233,243,...)2(n=6)0(Otherwise)

左辺が正になるのはn=l(p),l(p2),の場合のみなので, LTE, 補題5より直ちに従う.

さて, d|nF~n=Fnより, F~nFnを割り切るので, F~nFnの原始的素因数があることを言えばよい. 参照元の記事の記号を用い, 以下のように定義する. 同様にλnを評価していく.

F~n=λnPnとおく. ここで, PnFnの原始的素因数のみを素因数に持ち, λnは持たない.

λnの素因数の考察

n=6,12の場合を除き, λn1または素数である.

素数pに対して, pFl(p)の原始的素因数でないので, 補題7, 8よりpλnの素因数のとき, n=psl(p)なる正整数sをとれる.
λnが相異なる素因数p,qを持つとする.
このとき, ある正整数s,tが存在し,
n=psl(p)=qtl(q)
が成立する. p,qは相異なるので, l(p),l(q)はそれぞれqt,psの倍数である. 前回の記事 の定理3より, l(p)p+1,l(q)q+1なのでこれをあわせると,
qqtl(p)p+1,ppsl(q)q+1
となる. 従って, 偶奇を考えればp,qのいずれかは2であり, このときもう一方はl(2)=3となる. このとき, n=12である.
それ以外の場合, λn=psと記述できるが, 補題8よりn=6の場合を除きs=1である.

これより, F~nが十分大きいことを言えれば, Pn>1となりZsigmondyの類似が証明される.

F~nの評価

n=1,2,6,12の場合を除き, F~n>λnである.

α=(1+5)/2,β=(15)/2とおけば, 補題7より|F~n|=gcd(m,n)=1,1mn|αζnmβ|である. 三角不等式より右辺の各因数は1以上であり, n3の場合1より大きいものが存在するので, λn=1の場合は示された.

λn=p:素数とし, np2の倍数とする. このとき, n=p2Nとおけば, 1mNpgcd(m,n)=1を満たすとき, gcd(Npm,n)=1である. これらをペアにして考え, nNpmnについても同様のペアを考えれば, 合計でφ(pN)(p1)組のペアができる. また, 各ペアについて|αζnmβ||αζn(Npm)β|は, 「中心α半径βの円の原点における方べきの値(=5)」以上であるため(図を描いて位置関係の議論による), |αζnmβ||αζn(Npm)β|5が示された.
よって, |F~n|=gcd(m,n)=1,1mn|αζnmβ|5(p1)/2が言え, これはpよりも大きい.

そうでない場合, n=pNとおく. (以下の変形は参照元で利用されているものである.)
|F~n|=gcd(m,pN)=1,1mpN|αζpNmβ|=gcd(m,N)=1,1mpN|αζpNmβ|gcd(m,N)=1,p|m,1mpN|αζpNmβ|=gcd(m,N)=1,1mN|αpζNmβp|gcd(m,N)=1,1mN|αζNmβ|gcd(m,N)=1,1mN|αp(β)p|gcd(m,N)=1,1mN|αβ||αp(β)p||αβ|αp15
この最右辺はp7pより大きい. p=2,3,5のとき, λn=pおよびvp(n)=1より, それぞれn=6,12,(存在しない)となり, いずれも除外済みである.

以上より, n1,2,6,12のときPn>1となり, すなわちFnには原始的素因数が存在する.


これでZsigmodyの類似も示されました!
これを利用すると, 例えばJMO2011本選2のフィボナッチ数列版のさらに一般化を考察できたりします(あまりいい感じの応用が思いつかなかった...).

JMO2011本選2, 強化版

Fm=Ft1×Ft2××Ftn
を満たす2以上の整数nおよび3以上の整数m,t1,t2,,tnの組を全て求めよ.

とか. Fmの原始的素因数を取ればm=6,12と分かります.


さて, 今までフィボナッチ数列の性質を色々扱ってきましたが, 実はほとんどの命題においてフィボナッチ数列であることは本質的に使っていません. 一般のa0=0,a1=1なる三項間漸化式によって定義される数列(特にan+2=can+1+anという形のもの)において多少形は違えどほとんど同じ結果が得られます. ただ, その中でもフィボナッチ数列というのは一番シンプルなだけあり, 得られる結果も綺麗になりやすいです. 今後はもっと二次体とかと絡めた議論をしてみて, より一般的な結果を得たいですね.

最後まで読んでいただきありがとうございました!また次回(フィボナッチ数列関連は今回で最後かも?)

参考文献

投稿日:20231015
更新日:202433
OptHub AI Competition

この記事を高評価した人

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

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

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

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

投稿者

りぼーす
りぼーす
139
29952

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中
  1. LTEの類似
  2. 使ってみよう!
  3. Zsigmondyの定理の類似
  4. 参考文献