题面是这样的，其实斐波那契我们之前也有接触过，并不是什么太陌生的玩意，第一个想到的方法其实是用递归来做，这样的话其实是非常轻松的，but同志们你们有没有关注过这样一个鬼东西
你以为蓝题是让你切着玩的吗？？？？？？
果不其然，递归写了一个，大红大紫啊喂
（雾
不过的确过不了就是了，直到我打开了题解，发现了一个叫矩阵快速幂的玩意
Fn表示数列的第n项
那么我们如果把Fn，Fn-1写成蒟阵的形式，可以按照如下推导过程对这个蒟阵进行拆分，从而写成便于计算的形式
其实我们就是把递归用矩阵的方式写了出来，然后想求第n项就直接输出矩阵的n次幂即可
快速幂在另一篇博客里看这里qaq
其实这道题的难点就是矩阵快速幂，既然会了这个的话就没什么大问题啦~
代码如下：
#include<iostream> #include<algorithm> #include<cstdio> #include<cstring> #include<cmath> #include<cstdlib> using namespace std; typedef long long ll; typedef long double ld; typedef pair<int,int> pr; const double pi=acos(-1); #define rep(i,a,n) for(ll i=a;i<=n;i++) #define per(i,n,a) for(ll i=n;i>=a;i--) #define Rep(i,u) for(int i=head[u];i;i=Next[i]) #define clr(a) memset(a,0,sizeof a) #define pb push_back #define mp make_pair #define fi first #define sc second ld eps=1e-9; ll pp=1000000007; ll mo(ll a,ll pp){if(a>=0 && a<pp)return a;a%=pp;if(a<0)a+=pp;return a;} ll powmod(ll a,ll b,ll pp){ll ans=1;for(;b;b>>=1,a=mo(a*a,pp))if(b&1)ans=mo(ans*a,pp);return ans;} ll read(){ ll ans=0; char last=' ',ch=getchar(); while(ch<'0' || ch>'9')last=ch,ch=getchar(); while(ch>='0' && ch<='9')ans=ans*10+ch-'0',ch=getchar(); if(last=='-')ans=-ans; return ans; } //head 从这里开始哦 struct matrix{ ll a[2][2]; };//注意这里要用ll保证不会爆 matrix operator *(matrix a, matrix b){//定义*运算 matrix c; rep(i,0,1)//简写的方式，for（int i =1;i<=n;++i)
rep(j,0,1){ c.a[i][j]=0; rep(k,0,1) c.a[i][j] = (c.a[i][j]+a.a[i][k]*b.a[k][j])%pp; } return c; } ll k; int main(){ cin>>k; matrix a; a.a[0][0]=0;a.a[0][1]=1; a.a[1][0]=1;a.a[1][1]=1; matrix ans; ans.a[0][0]=1;ans.a[0][1]=0; ans.a[1][0]=0;ans.a[1][1]=1;//把ans初始化为单位矩阵 ll b=k-1; while(b){ if(b&1)ans=ans*a; a=a*a; b/=2; }//一个快速幂 ll fk = (ans.a[0][0]+ ans.a[0][1])%pp; cout<<fk<<endl;//O(log B *2^3) }