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