初学矩阵,建议看一下这个bloghttps://www.cnblogs.com/cmmdc/p/6936196.html讲得比较清楚简单。
codeforces450B
思路:很简单的水题,初学有几个注意点,矩阵幂一定要是方阵,填0凑成方阵;ans是单位矩阵,res要自己算,我根据上面博客最后求斐波那契数列的方式写的,一定要前后承接,多乘几个就能发现bug;负数的模运算,+mod凑成整数再模。
#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <stdio.h>
#include <iostream>
#include <cstdlib>
#include <cmath>
#include <cctype>
#include <string>
#include <cstring>
#include <algorithm>
#include <stack>
#include <queue>
#include <set>
#include <map>
#include <ctime>
#include <vector>
#include <fstream>
#include <list>
#include <iomanip>
#include <numeric>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
#define ms(s) memset(s, 0, sizeof(s))
const int inf = 0x3f3f3f3f;
const int mod =1000000007;
#define LOCAL
int x,y,n;
struct Maxtrix
{
int a[2][2];
}ans,res;
Maxtrix Mul(Maxtrix A,Maxtrix B)
{
Maxtrix temp;
for(int i=0;i<2;++i)
for(int j=0;j<2;++j)
temp.a[i][j]=0;
for(int i=0;i<2;++i)
for(int j=0;j<2;++j)
for(int k=0;k<2;++k)
temp.a[i][j]+=(A.a[i][k]%mod*B.a[k][j]%mod)%mod;
return temp;
}
void quick_power(int x,int y,int n)
{
ans.a[0][0]=ans.a[1][1]=1;
ans.a[0][1]=ans.a[1][0]=0;
res.a[0][0]=res.a[1][0]=1;
res.a[0][1]=-1;res.a[1][1]=0;
while(n)
{
if(n&1)ans=Mul(res,ans);
res=Mul(res,res);
/* for(int i=0;i<2;++i)
for(int j=0;j<2;++j)
cout<<ans.a[i][j]<<" ";*/
//cout<<endl;
n=n>>1;
}
}
int main()
{
// freopen("in.txt", "r", stdin);
while(cin>>x>>y>>n)
{
if(n==1){cout<<(x%mod+mod)%mod<<endl;}
else if(n==2){cout<<(y%mod+mod)%mod<<endl;}
else
{
quick_power(x,y,n-2);
cout<<((y*ans.a[0][0]%mod+x*ans.a[0][1]%mod)%mod+mod)%mod<<endl;
}
}
return 0;
}