【转】http://www.cnblogs.com/Knuth/archive/2009/09/04/1559951.html#2576160
[定理1] 标准Fibonacci序列（即第0项为0，第1项为1的序列）当N大于1时，一定有f(N)和f(N-1)互质
其实，结合“互质”的定义，和一个很经典的算法就可以轻松证明
对，就是辗转相除法
互质的定义就是最大公约数为1
数学归纳法是很有用的证明方法，我们接下来这个定理用数学归纳法就很好证明：
[定理2]若i为奇数， f(i)*f(i)=f(i-1)*f(i+1)+1，否则f(i)*f(i)=f(i-1)*f(i+1)-1
对，这个定理用数学归纳法可以轻松证明，大家有兴趣可以自己尝试
[定理3] f(n)=f(i)*f(n-i-1)+f(i+1)*f(n-i)
f(n)=f(1)*f(n-2)+ f(2)*f(n-1)
=f(2)*f(n-3)+ f(3)*f(n-2)
=f(3)*f(n-4)+ f(4)*f(n-3)
看来没有错，证明方法就是这样
这个公式也可以用来计算较大的fibonacci数除以某个数的余数
设i=n/2 不过，为了保证计算能延续下去 需要每次保留三个值
这样，下一次计算就可以利用这三个值来求出两个值，再相加就可以得到第三个值
譬如，计算出f(5),f(6),f(7)，可以计算出f(11)、f(12)，然后推出f(13)
就是刚才洛奇的悲鸣(364738334)所提到的矩阵方法
我们知道我们若要简单计算f(n)，有一种方法就是先保存
a=f(0),b=f(1)，然后每次设：
a'=b b'=a+b
并用新的a'和b'来继续这一运算
如果大家熟悉利用“矩阵”这一工具的话，就知道，如果把a、b写成一个向量[a,b]，完成上述操作相当于乘以矩阵
0 1
1 1
也就是说，如果我们要求第100个fibonacci数，只需要将矩阵
[0,1]乘上
0 1
1 1
的一百次方，再取出第一项
因为我们知道，矩阵运算满足结合律，一次次右乘那个矩阵完全可以用乘上那个矩阵的N次方代替，更进一步，那个矩阵的N次方就是这样的形式：
f(n-1) f(n)
f(n) f(n+1)
而求矩阵的N次方，由于矩阵乘法满足结合律，所以我们可以用log(N)的算法求出——这个算法大家都会么？
一个是二分，一个是基于二进制的求幂
二分的原理：要求矩阵的N次方A(N)，设i=N/2若N%2==1， 则 A(N)=A(i)*A(i)*A(1)若N%2==0， 则 A(N)=A(i)*A(i)
基于二进制的原理：将N拆为二进制数，譬如13=1101那么 A^13= A^8 * A^4 * A^1 （这里^表示幂运算）
也就是说，由A^1开始，自乘得到A^2，然后自乘得到A^4，如果N对应位为1，则将这个结果乘到目标上去
这样的话，将所有乘法改为模乘，就可以得到一个较大Fibonacci数除以M的余数
若不用递归，其实类似
http://acm.pku.edu.cn/JudgeOnline/problem?id=3070
这里用的fib矩阵略有不同，是
f(n+1) f(n)
f(n) f(n-1)
但实际上可以验证效果是一样的
这题是要求求F(n)的最后四位数，所有乘法过程增加一个模10000的步骤即可，大家可以收藏稍候AC
关于矩阵我们告一段落，等下会回来继续探讨利用矩阵来解决复杂些的Fibonacci问题
http://acm.hdu.edu.cn/showproblem.php?pid=1568
我们来看这题，这题要求求出Fibonacci某项的前四位
当然，用矩阵也可以解决这道题——只要将乘法改为乘并保留前四位
我们采用double 保留整数部分四位 这题最好还是double吧
不过显然有更好的解法——如果我们知道Fibonacci序列的通项公式
F(n) = (((1+Sqrt(5))/2)^n - ((1-Sqrt(5))/2)^n)*1/Sqrt(5)
不过组合数学里也有这一公式的推导方法 叫做“线性齐次递推式”
这个解法的核心是，通解是某个数的幂 将f(n)=x^n代入递推方程，可以解出三个通解 0和 (1+sqrt(5))/2
通常把“0”称作平凡解，那么特解就是通解的某个线性组合
再代入f(0)=0 f(1)=1，就可以得出我们刚才的公式
不过通常情况下，我们只需要记住那个公式就可以了
提醒大家，记忆公式的时候千万别忘记了系数1/sqrt(5)
因为(1-sqrt(5))/2的绝对值小于1
所以当i较大的时候，往往可以忽略掉这一项
f(i)≈((1+Sqrt(5))/2)^n/sqrt(5);
所以，刚才列举出的HDOJ的1568，可以很简单的30以内直接求解，30以上采用这个公式，还是用log(N)求幂的算法求解
恩，就是公式的前半部分
http://acm.hdu.edu.cn/showproblem.php?pid=1021
或 http://acm.zju.edu.cn/show_problem.php?pid=2060
Fibonacci某项是否被3整除
[定理5] 标准Fibonacci序列对任意大于2的正整数的余数序列，必然是以“0 1”为循环节开头的序列
显然0、1是序列开头，也就是说序列开头就是循环节开头
循环长度的计算貌似是个比较难的问题，我一时还没有想到有效解法，不过，要说明的是，计算复杂度时，这个循环节长度应该按复杂度O(N^2)计算
恩，证明方法是利用同余定理、反证法，还有我们之前证明过的相邻项一定互质的定理，留给大家家庭作业
http://acm.hdu.edu.cn/showproblem.php?pid=1588
这是前天比赛关于Fibonacci的一道题，大家先看看题。
Description看后半部分就行了
现在告诉大家一种正确解法，然后大家就可以去搞定这道题向别人炫耀了
首先，我们将问题整理一下，就是对等差数列 ai=k*i+b，求所有的f(ai)之和除以M的余数
当0<=i<N
大家有没有想到，因为ai是等差数列，倘若f(ai)也是个等什么序列，那说不定就有公式求了
f(ai)显然不是等差数列，直接看上去也不是等比数列
但是如果把f(ai)换成我们刚才所说的Fibonacci矩阵呢？
是的，可是我们对矩阵是直接求幂即可，由于矩阵加法的性质，我们要求A^ai的右上角元素之和，只要求A^ai之和的右上角元素
就矩阵这个东西来说，完全可以看作一个等比数列，
首项是：A^b
公比是：A^k
项数是：N
呵呵，我们可以把问题进一步简化
因为矩阵的加法对乘法也符合分配律，我们提出一个A^b来，形成这样的式子：
A^b*( I + A^k + (A^k)^2 + .... + (A^k)^(N-1) )
A^b 和 A^k 显然都可以用我们之前说过的方法计算出来，这剩下一部分累加怎么解决呢
简单起见，设A^k=B
要求 G(N)=I + ... + B^(N-1)，设i=N/2
若N为偶数，G(N)=G(i)+G(i)*B^i
若N为奇数，G(N)=I+ G(i)*B + G(i) * (B^(i+1))
呵呵，这个方法就是比赛当时ACRush用的
而农夫用的则是更精妙的方法，大家可想知道
我们来设置这样一个矩阵
B I
O I
其中O是零矩阵，I是单位矩阵
将它乘方，得到
B^2 I+B
O I
乘三方，得到
B^3 I+B+B^2
O I
乘四方，得到
B^4 I+B+B^2+B^3
O I
既然已经转换成矩阵的幂了，继续用我们的二分或者二进制法，直接求出幂就可以了
http://online-judge.uva.es/p/v110/11089.html
大家来读读这一题
Fibinary数是指没有相邻的两个1的二进制数。给N，求出第N大的Fibinary数
相对于二进制中每一位的值是2的幂，十进制中每一位的值是十的幂，
Fibonacci进制是每一位的值是对应Fibonacci数的一种计数系统。
8 5 3 2 1
1 1
2 1 0
3 1 0 0
4 1 0 1
5 1 0 0 0
6 1 0 0 1
7 1 0 1 0
8 1 0 0 0 0
9 1 0 0 0 1
10 1 0 0 1 0
11 1 0 1 0 0
12 1 0 1 0 1
以上是前12个数字对应的十进制到Fibonacci进制的表格
Fibonacci的运算方法很奇怪。首先，它每一位上非0即1，而且不同于二进制的逢二进一或者十进制的逢十进一，它的进位方法是逢连续两个1，则进1
譬如
1010110==> 1011000 ==> 1100000==>10000000
在最低位有个特殊情况，最低位既可以逢2进1，也可以和次低位一起逢相邻进1
这种奇怪的进位方法，换句话描述就是，不存在两个连续的1
因为Fibonacci数其实也增长很快，int范围内好像只有46个，本题只需要用最简单的办法转换成Fibonacii进制即可
其中一题是
http://www.mydrs.org/program/down/ ahoi2004da y1.pdf
中的第二题，叫做数字迷阵
还有一题是PKU上的很出名 的取石子问题
http://acm.pku.edu.cn/JudgeOnline/problem?id=1067
这题相当复杂，大家可以自己思考，往Fibonacci上想，也可以阅读这里的论文：
http://episte.math.ntu.edu.tw/articles/mm/mm_03_2_02/index.html
http://acm.pku.edu.cn/JudgeOnline/problem?id=2967
另外这题 可以利用Fibonacci判断数据范围进行优化设计。不过貌似可以水过去，仅仅给大家提供个思路罢
关于Fibonacci和黄金分割，还有很多更高明的结论和定理，希望大家也继续讨论，将自己的知识和他人共享
http://episte.math.ntu.edu.tw/articles/mm/mm_02_4_10/index.html
中例3详细讲述了用生成函数来计算Fibonacci数公式的运算过程。http://acm.hdu.edu.cn/showproblem.php?pid=1568
Fibonacci 求fibonacci前4位
http://acm.hdu.edu.cn/showproblem.php?pid=1588
Gauss Fibonacci
http://acm.pku.edu.cn/JudgeOnline/problem?id=1067
取石子问题
http://episte.math.ntu.edu.tw/articles/mm/mm_03_2_02/index.html 是报告
http://acm.pku.edu.cn/JudgeOnline/problem?id=3070
Fibonacci矩阵
http://acm.hdu.edu.cn/showproblem.php?pid=1021
或
http://acm.zju.edu.cn/show_problem.php?pid=2060
Fibonacci某项是否被3整除
http://acm.pku.edu.cn/JudgeOnline/problem?id=2116
Fibonacci进制计算
http://acm.pku.edu.cn/JudgeOnline/problem?id=2967
利用Fibonacci判断数据范围进行优化设计。
http://online-judge.uva.es/p/v110/11089.html
Fi-binary numbers，Fibonacci进制。
http://www.mydrs.org/program/down/ahoi2004day1.pdf
第二题 数字迷阵 这些，是今天涉及到的资料和网页