很久在学校的时候看一个人微软的讲面试讲到斐波拉契数列。当初也没做太深的研究只记得可以求出一个通项公式。最近无聊翻出来研究下。
第一种:递归(最简单,最low)
时间复杂度O(n2)
这里就不做大数处理了。超过50跑起来时间就超过了我的忍耐了。
private int fib(int i) {
if(i==0)
return 0;
else if(i==1)
return 1;
return fib(i-1)+fib(i-2);
}
第二种:循环
时间复杂度 O(n)
private long fib(int n) {
long x;
long fi1 = 1;
long fi0 = 0;
for(int i = 2;i<=n;i++){
x = fi1+fi0;
fi0 = fi1%1000000007;
fi1 = x%1000000007;
}
return fi1;
}
第三种:矩阵乘法规律
时间复杂度 O(log n)
矩阵可以得到下面公式
上面公式可以推出下面公式,
“`
由上面公式直接求矩阵 矩阵幂就好了。
因为f(1)=1
f(0)=0
即: f(n) = a[0][0]
代码如下,注意快速幂那块,加分项。
private int fib(int i) {
if(i==0)
return 0;
else if(i==1)
return 1;
long[][] fs = pow(f,i-1);
return (int) fs[0][0];
}
private long[][] pow(long[][] f2, int i) {
if(i==1)
return f;
if(i==2)
return fun(f,f);
if( i % 2 == 0){//偶数
return fun(pow(f,i/2),pow(f,i/2));
}else{
return fun(pow(f,i/2),pow(f,i/2 + 1));
}
}
private long[][] fun(long[][] f2, long[][] f3) {
long[][] temp = new long[2][2];
temp[0][0] = (f2[0][0]*f3[0][0] + f2[0][1]*f3[1][0])%1000000007;
temp[0][1] = (f2[0][0]*f3[0][1] + f2[0][1]*f3[1][1])%1000000007;
temp[1][0] = (f2[1][0]*f3[0][0] + f2[1][1]*f3[1][0])%1000000007;
temp[1][1] = (f2[1][0]*f3[0][1] + f2[1][1]*f3[1][1])%1000000007;
return temp;
}
第四种,通项公式求解。
其实fibonacci数列是有一个通项公式的如下:
因此就可以直接根据公式返回其求出值double fi;
然后返回return (int)fi+0.5;
代码就不贴了。时间复杂度O(1).