快速矩阵幂

#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;
}