POJ 3420 Quad Tiling 线性递推 矩阵快速幂

题意：给出一个4 * N的方格。用1 * 2的格子去填充，问有方案数。

思路：以前做过填充格子的轮廓线DP。但是因为N<=10^9,根本无法进行DP。

          但是注意到所有格子都是可以放1*2的格子的，而不是像其他的题有某些障碍，是没有办法放置。

          类比轮廓线DP的思路，我们可以以一行为转移单位。枚举前一行的所有可能的情况，对于每种情况，将没有填满的空填满。这样就转移到了下一行，继续进行上面的过程。

          因为没有障碍，上面的转移过程就是一个固定的模式。

          但是怎样求出这样的转移方式呢？和轮廓线DP的思路继续类似，我们填充一行就能的到下一行的结果，填充的方式也是横着和竖着两种。具体细节，看dfs函数。

          得到对应的转移方式，我们就能用矩阵快速幂来进行加速。

注意：得到的转移方式，我们也可以看做构造出了一张有向图，而矩阵乘法的过程就是求从某点到另外一点的所有路径计数。而最终答案可以这样描述；从111状态回到111状态且长度为N的路径有多少种。

代码如下：

#include <cstdio>
#include <vector>

using namespace std;

typedef long long ll;
typedef vector<ll> vec;
typedef vector<vec> mat;

ll N,M;
mat A(16,vec(16));

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;
}
mat pow(mat A, ll n)
{
	mat B(A.size(), vec(A.size()));
	for (int i = 0; i < A.size(); ++i)
		B[i][i] = 1;
	while (n){
		if (n & 1) B = mul(B, A);
		A = mul(A, A);
		n >>= 1;
	}
	return B;
}

void dfs(int l , int now, int pre)//pre为当前行的状态,now是将pre行填完后，下一行的状态
{
    if(l > 4) return;
    if(l == 4){
        A[now][pre] = 1;
        return;
    }
    dfs(l+1,(now<<1)|1,pre<<1);//竖着放
    dfs(l+1,(now<<1),(pre<<1)|1);//当前列已经被放置，没法在放
    dfs(l+2,(now<<2)|3,(pre<<2)|3);//下两列已被放置，横着放
    return;
}


int main(void)
{
    //freopen("input.txt","r",stdin);
    dfs(0,0,0);
    while(scanf("%lld%lld",&N,&M),N||M){
        mat B(16,vec(16));
        B = pow(A,N);
        printf("%lld\n",B[15][15]);
    }
    return 0;
}