hdu 2842 Chinese Rings

 

http://acm.hdu.edu.cn/showproblem.php?pid=2842

 

题目大意：有n个环在一条线上，叫你把他们通通取下来最少要做几步操作，结果mod 200907.要拿下环来是有条件的，i-th环的前面的环如果通通拿下来了，那么(i+1)-th这个环，你可以放进去，也可以取下来.注意，第一个环，可以随意拿进取下，其他环就是要按照规则才能进行取下和拿进操作。

 

解题思路：构造矩阵乘法，首先要推出公式

 

f(n) = 2*f(n-2)+f(n-1)+1

 

怎么推得呢？拿5个环来推吧

 

设f(n)表示取下n个环的步数，那么f(5)就表示取下5个环要走几步啦。

首先我们想，要想取下最后一个环，我们就必须把前面3个环取下，第四个环留着，所以就需要f(3)步，取下最后一个环要一步，所以总共是f（3）+1步

之后我们把前面3个环放回去，就要f(3)步，之后我们只要再取走剩下的4个就好了，那不就是f(4)步吗。

 

所以  f(5) = 2*f(3)+f(4)+1,通过扩展就可得到上面的f(n)的公式了。

 

/* f(n) = 2*f(n-2)+f(n-1)+1 矩阵: Init: 1 2 1 1 0 0 0 0 1 e: 1 0 0 0 1 0 0 0 1 r: 2 0 0 1 0 0 1 0 0 结果为 Init^(n-2)*r */ #include <stdio.h> #define size 3 int n; struct Mat { __int64 mat[size][size]; }e,init,r; void InitMat() { int i,j; for (i=0;i<size;i++) { for (j=0;j<size;j++) { e.mat[i][j] = (i==j); r.mat[i][j] = init.mat[i][j] = 0; } } r.mat[1][0] = r.mat[2][0] = init.mat[0][0]=init.mat[0][2]=init.mat[1][0]=init.mat[2][2]=1; r.mat[0][0] = init.mat[0][1]= 2; } Mat Mul(Mat a,Mat b) { Mat c; int i,j,k; for (i=0;i<size;i++) { for (j=0;j<size;j++) { c.mat[i][j] = 0; for(k=0;k<size;k++) c.mat[i][j]+=(a.mat[i][k]*b.mat[k][j])%200907; c.mat[i][j]%=200907; } } return c; } __int64 Cal() { int t = n-2; Mat p,q; p = e,q = init; while (t) { if(t&1) p = Mul(p,q); t>>=1; q = Mul(q,q); } p = Mul(p,r); return p.mat[0][0]; } int main() { InitMat(); while (scanf("%d",&n),n) { if(n==1) printf("1/n"); else if(n==2) printf("2/n"); else printf("%I64d/n",Cal()); } return 0; }