按自己风格写的,以后就直接用了。
#include<iostream>
#include<cstdio>
using namespace std;
typedef long long ll;
const ll mod=1e9+7;
// [fn] = [1,1] * [fn-1]
// [fn-1] [1,0] [fn-2]
// 以0作为第一项
ll m[3][3],g[3][3];
ll fgfk(int n) {
m[1][1]=1,m[1][2]=1,m[2][1]=1,m[2][2]=0;
g[1][1]=1,g[1][2]=0,g[2][1]=0,g[2][2]=1;
if(n==1) {
return 0;
}else if(n==2) {
return 1;
}
n-=2;
while(n>0) {
if(n&1) {
ll a1=(g[1][1]*m[1][1]+g[1][2]*m[2][1])%mod;
ll a2=(g[1][1]*m[1][2]+g[1][2]*m[2][2])%mod;
ll a3=(g[2][1]*m[1][1]+g[2][2]*m[2][1])%mod;
ll a4=(g[2][1]*m[1][2]+g[2][2]*m[2][2])%mod;
g[1][1]=a1,g[1][2]=a2,g[2][1]=a3,g[2][2]=a4;
}
ll a1=(m[1][1]*m[1][1]+m[1][2]*m[2][1])%mod;
ll a2=(m[1][1]*m[1][2]+m[1][2]*m[2][2])%mod;
ll a3=(m[2][1]*m[1][1]+m[2][2]*m[2][1])%mod;
ll a4=(m[2][1]*m[1][2]+m[2][2]*m[2][2])%mod;
m[1][1]=a1,m[1][2]=a2,m[2][1]=a3,m[2][2]=a4;
n>>=1;
}
return g[1][1];
}
int main() {
int n;
while(scanf("%d",&n)==1) {
printf("第%d项 = %lld\n",n,fgfk(n));
}
}