poj3420

题意：用1 X 2的矩形填充4 X n的矩形，共有多少种不同方法。
//思路：网上的代码大部由状态压缩与矩阵乘法的结合，这里给大家说一由递推+矩阵乘的方法。
// 我们把4当作列数，n当作行数。当第n行填满时，第(n+1)行会出现以下几种情况：
//a[n] = a[n-1] + b[n-1] + c[n-1] + dx[n-1] + dy[n-1];
//b[n] = a[n-1];
//c[n] = a[n-1] + e[n-1];
//dx[n]= a[n-1] + dy[n-1];
//dy[n]= a[n-1] + dx[n-1];
//e[n] = c[n-1].
// 令d[n] = dx[n] + dy[n],则
// a[n] = a[n-1] + b[n-1] + c[n-1] + d[n-1];
//b[n] = a[n-1];
//c[n] = a[n-1] + e[n-1];
//d[n] = a[n-1] * 2 + d[n-1];
//e[n] = c[n-1].
={ a[n-1] b[n-1] c[n-1] d[n-1] e[n-1] } X
///  1  1  1  2  0 \
// |  1  0  0  0  0 |
// |  1  0  0  0  1 |
// |  1  0  0  1  0 |
// \  0  0  1  0  0 /
//A[0] = 1;

//172K 0MS

#include <iostream>
#include <vector>
using namespace std;

static int N,M;
typedef vector<int>vec;
typedef vector<vec>mat;
typedef long long ll;

static mat mul(mat&A,mat&B){
	mat C(A.size(),vec(B[0].size()));
	for (int i=0;i<A.size();++i)
	{
		for (int k=0;k<B.size();++k)
		{
			for (int j=0;j<B[0].size();++j)
			{
				C[i][j] = (C[i][j] + A[i][k]*B[k][j])%M;
			}
		}
	}

	return C;
}

static mat pow(mat A,int n){
	mat B(A.size(),vec(A.size()));
	for (int i=0;i<A.size();++i)
		B[i][i] = 1;

	while(n>0){
		if (n&1)
			B = mul(B,A);
		A = mul(A,A);
		n>>=1;
	}

	return B;
}

int main(){
	while(scanf("%d %d",&N,&M)){
		if (!N&&!M)
		 break;
		mat A(5,vec(5));
		A[0][0] = 1;A[0][1] = 1;A[0][2] = 1;A[0][3] = 2;A[0][4] = 0;
		A[1][0] = 1;A[1][1] = 0;A[1][2] = 0;A[1][3] = 0;A[1][4] = 0;
		A[2][0] = 1;A[2][1] = 0;A[2][2] = 0;A[2][3] = 0;A[2][4] = 1;
		A[3][0] = 1;A[3][1] = 0;A[3][2] = 0;A[3][3] = 1;A[3][4] = 0;
		A[4][0] = 0;A[4][1] = 0;A[4][2] = 1;A[4][3] = 0;A[4][4] = 0;
		
		A = pow(A,N);

		printf("%d\n",A[0][0]);
	}

	return 0;
}