HDOJ 1005 - 上周新队员做的题～我也来水水～矩阵乘法的说..

    其实这道题开始那个谁选的时候是用的找循环节的方法～～～反正是mod 7嘛～～并且每一项只与前两项有关～～所以最多最多7*7=49次就会进行循环了～～所以找到循环节再mod一下循环节长度对应找到就可以了～

    而我是用矩阵乘法来解决的～～反正这是一类问题～～用矩阵乘法解递推数列～～递推关系都给出来了～～那么就可以构造矩阵～～

     0  1

     B  A

     当n<=2时～～还没超出初值～～直接输出1～～当n>2时～～计算这个矩阵的n-2次方～～得到的矩阵的( [ 1 ] [ 0 ] + [ 1 ] [ 1 ] )%7 即使答案～～

   PS : 我这题是再UBUNTU下裸敲的～～～MOhahahaha~~~

Program:

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