还需改进

#include<bits/stdc++.h>
using namespace std;
//第一步：用结构体来存一个矩阵
struct matrix{
int a[2][2];

};

matrix mul(matrix a,matrix b)
{
matrix ans;
/*需要用三重循环来运算矩阵乘法
乘完还是一个矩阵 */
for(int i=0;i<2;i++){
for(int j=0;j<2;j++){
ans.a[i][j]=0;
for(int k=0;k<2;k++){
ans.a[i][j]=(ans.a[i][j]+(a.a[i][k]*b.a[k][j])%7)%7;
}
}
}
return ans;

}
int main()
{
int A,B,n;
while(~scanf("%d%d%d",&A,&B,&n)){
if(A0&&B0&&n0)
break;
//n<=2时需要特殊说明
if(n<=2){
printf(“1\n”);
continue;
}
matrix a,ans;//定义一个单位矩阵ans,和关系矩阵a;
//关系矩阵 对应系数
a.a[0][0]=A;
a.a[0][1]=B;
a.a[1][0]=1;
a.a[1][1]=0;
//初始化为单位矩阵
ans.a[0][0]=ans.a[1][1]=1;
ans.a[0][1]=ans.a[1][0]=0;
n=n-2;
//快速幂
while(n){
if(n%21){
ans=mul(ans,a);
}
n/=2;
a=mul(a,a);
}
printf("%d\n",(ans.a[0][0]+ans.a[0][1])%7);
}

return 0;

}