Luogu 1962（矩阵快速幂）

传送门



图片转自：http://www.cnblogs.com/huxianglin/p/5995649.html

代码作为模板。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
const int MOD=1e9+7;
ll n;
struct Matrix {
int v[2][2];
Matrix() {memset(v,0,sizeof(v));}
}un,t;
Matrix mult(Matrix a,Matrix b) {
Matrix ret;
for (int k=0;k<2;++k)
for (int i=0;i<2;++i)
for (int j=0;j<2;++j)
ret.v[i][j]=(ret.v[i][j]+1ll*a.v[i][k]*b.v[k][j]%MOD)%MOD;
return ret;
}
int main() {
scanf("%lld",&n);
un.v[0][0]=un.v[1][1]=1;
t.v[0][0]=t.v[0][1]=t.v[1][0]=1;
while (n) {
if (n&1) un=mult(un,t);
t=mult(t,t),n>>=1;
}
printf("%lld\n",un.v[1][0]);
return 0;
}