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

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

版权声明：本文内容由网络用户投稿，版权归原作者所有，本站不拥有其著作权，亦不承担相应法律责任。如果您发现本站中有涉嫌抄袭或描述失实的内容，请联系我们jiasou666@gmail.com 处理，核实后本网站将在24小时内删除侵权内容。