题意:给出一个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;
}