找出递推关系式就好了
（fi+1）=（1 -1）（fi ）
（ fi）=（1 0）（fi-1）
不会打矩阵将就着看吧。。。
这是第一道矩阵快速幂。细节还是有很多没注意到的
本来想看挑战写的，结果上面的vector套vector看的我头都晕了。。
#include<map> #include<set> #include<cmath> #include<queue> #include<stack> #include<vector> #include<cstdio> #include<iomanip> #include<cstdlib> #include<cstring> #include<iostream> #include<algorithm> #define pi acos(-1) #define ll long long #define mod 1000000007 #define ls l,m,rt<<1 #define rs m+1,r,rt<<1|1 using namespace std; const double g=10.0,eps=1e-9; const int N=100+5,maxn=1<<10+5,inf=0x3f3f3f3f; struct Node{ ll row,col; ll a[N][N]; }; Node mul(Node x,Node y) { Node ans; memset(ans.a,0,sizeof ans.a); ans.row=x.row,ans.col=y.col; for(ll i=0;i<x.row;i++) for(ll j=0;j<y.row;j++) for(ll k=0;k<y.col;k++) ans.a[i][k]=(ans.a[i][k]+x.a[i][j]*y.a[j][k]+mod)%mod; return ans; } Node quick_mul(Node x,ll n) { Node ans; ans.row=x.row; ans.col=x.col; memset(ans.a,0,sizeof ans.a); for(int i=0;i<ans.row;i++)ans.a[i][i]=1; while(n){ if(n&1)ans=mul(ans,x); x=mul(x,x); n/=2; } return ans; } int main() { ios::sync_with_stdio(false); cin.tie(0); // cout<<setiosflags(ios::fixed)<<setprecision(2); int x,y,n; cin>>x>>y>>n; if(n==1) { cout<<(x+mod)%mod<<endl; return 0; } Node A,B; A.row=2,A.col=2; A.a[0][0]=1,A.a[0][1]=-1; A.a[1][0]=1,A.a[1][1]=0; /* for(int i=0;i<A.row;i++) { for(int j=0;j<A.col;j++) cout<<A.a[i][j]<<" "; cout<<endl; }*/ B.row=2,B.col=1; B.a[0][0]=y,B.a[1][0]=x; cout<<(mul(quick_mul(A,n-1),B).a[1][0]+mod)%mod<<endl; return 0; }