#include<iostream>
#include<cstring>
#include<algorithm>
#include<cstdio>
#include<cmath>
using namespace std;
struct mat
{ int date[2][2];
void init()//单位矩阵
{memset(date,0,sizeof(date));
date[0][0]=date[1][1]=1;
}
};
mat d,re;//d:幂矩阵
int a,b,n;
void init()//根据题目要求,将数据填入矩阵中
{ memset(d.date,0,sizeof(d.date));
d.date[0][0]=a;d.date[0][1]=b;
d.date[1][0]=1;
}
mat mul(mat u,mat v)//乘法
{ mat temp;
memset(temp.date,0,sizeof(d.date));
for(int i=0;i<2;i++)
for(int j=0;j<2;j++)
for(int k=0;k<2;k++)
temp.date[i][j]+=(u.date[i][k]*v.date[k][j])%7;//(可能是A^2*A^2或者是A^2*A)
return temp;
}
mat calc(mat op,int k)//矩阵快速幂
{
mat res;
res.init();//单位矩阵
while(k)
{
if(k&1)//判断是否是奇数,如果是奇数结果为1,偶数结果为0(用到二分快速幂的思想,如果奇数次方时,则乘以他本身一次)
res=mul(res,op);//(A^3时,结果是A^2*A)res:即为result结果
k>>=1;//(当A^4时,变为,A^2,A^2)
op=mul(op,op);//(A^4时,结果是A^2*A^2)
}
return res;
}
mat div(mat op,int k)//二分
{
if(k==1)
return op;
mat temp;
temp=calc(op,k/2);
temp=mul(temp,div(op,k-k/2));
return temp;
}
int main()
{
while(scanf("%d%d%d",&a,&b,&n)&&(a||b||n))
{
init();
if(n>2)
{
re=div(d,n-2);
cout<<(re.date[0][0]+re.date[0][1])%7<<endl;
}
else cout<<1<<endl;
}
return 0;
}
快速矩阵幂
最新推荐文章于 2022-09-06 10:58:10 发布