HDOJ1005

#include<iostream>
#include<string.h>
#include<stdio.h>
using namespace std;
struct Matrix
{
int s[2][2];
}h;
int A,B,n;
Matrix Mul(Matrix a,Matrix b)
{
Matrix h;
int i,j,k;
memset(h.s,0,sizeof(h.s));
for (k=0;k<2;k++)
for (i=0;i<2;i++)
for (j=0;j<2;j++)
h.s[i][j]=(h.s[i][j]+a.s[i][k]*b.s[k][j])%7;
return h;
}
Matrix find(int n)
{
Matrix k;
if (n==1) return h;
k=find(n/2);
k=Mul(k,k);
if (n%2) k=Mul(k,h);
return k;
}
int main()
{
while (scanf("%d%d%d",&A,&B,&n))
{
if (!A && !B && !n) break;
h.s[0][0]=0; h.s[0][1]=1;
h.s[1][0]=B; h.s[1][1]=A;
if (n<=2)
{
printf("1\n");
continue;
}
h=find(n-2);
printf("%d\n",(h.s[1][0]+h.s[1][1])%7);
}
return 0;
}

转载于:https://www.cnblogs.com/ct0421/p/3606075.html