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

你可能感兴趣的:(HDOJ 1005 - 上周新队员做的题~我也来水水~矩阵乘法的说..)