给定一个n*m的矩阵,使用1*2的小长方形覆盖矩阵,要求完全覆盖的同时不出现重合,也不允许超出边界,问有多少种可能的覆盖方法,方案数对1e9+7取模
2<=n<=1000
3<=m<=5
Input
整数n,m
Output
方案数
Sample input and output
Sample Input | Sample Output |
---|---|
2 4 | 5 |
Hint
注意取模
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <cmath>
#include <algorithm>
using namespace std;
typedef long long ll;
const int mod=1e9+7;
ll dp[1010][(1<<5)+10];//dp[i][j]=1表示第i行第j列放了骨牌
ll n,m;
bool init(ll s)
{
for(ll i=0;i<m;)
{
if(s&(1<<i))//提取十进制数s对应二进制的第i位,如果第1行第i列为1则一定要横放
{
if(i==m-1) return false;
if(s&(1<<(i+1))) i+=2;//提取十进制数s对应二进制的第i+1位,即第1行第i+1列要为1
else return false;
}
else i++;//如果第1行第i列为0
}
return true;
}
bool ok(ll s,ll ss)//ss表示上一行某一状态 s表示当前行某一状态
{
for(ll i=0;i<m;)
{
if(s&(1<<i))//提取十进制数s对应二进制的第i位
{
if(ss&(1<<i))//提取十进制数ss对应二进制的第i+1位
{
if(i==m-1||!(s&(1<<(i+1)))||!(ss&(1<<(i+1)))) return false;
i+=2;
// 第i行的第j列为1,第i-1行的第j列为1,这样的话,说明第i行的第j列一定
// 不是竖放而是横放否则会与第i-1行的第j列冲突,所以马上紧接着判断第i行
// 第j+1列,如果是1,那么满足横放的规则,同时也要第i-1行第j+1列也要为
// 1,否则的话这个格子没办法填充,成立后向左移动两格
}
else i++;
// 第i行第j列为1,第i-1行第j列为0,那么说明第i行第j列应该竖放并填充第
// i-1行第j列,成立后向左移动一格
}
else
{
if(ss&(1<<i)) i++;
else return false;
// 第i行第j列为0,说明不放方块,那么第i-1行第j列必须为1,否则没法填充
// 这个格子。若第i-1行第j列也为0,不兼容不合法(至于第i行第j列这个格子
// 空着干什么,其实就是留出来给第i+1行竖放的时候插进来的)
}
}
return true;
}
int main()
{
#ifdef LOCAL
freopen("input.in","r",stdin);
// freopen("output.in","w",stdout);
#endif // LOCAL
cin>>n>>m;
if(n%2==1&&m%2==1)
{
cout<<0<<endl;
return 0;
}
memset(dp,0,sizeof(dp));
ll maxn=(1<<m)-1;//2^m-1种可能,减一排除全为0的可能
for(ll i=0;i<=maxn;i++)//一行中所有放置可能对应的二进制的十进制数为i
{
if(init(i)) dp[1][i]=1;//
}
for(ll i=2;i<=n;i++)//i表示第i行
{
for(ll s=0;s<=maxn;s++)//枚举当前行所有可能状态,根据上一行状态判断
{
for(ll ss=0;ss<=maxn;ss++)//枚举上一行所有可能状态
{
if(ok(s,ss))
dp[i][s]=(dp[i][s]+dp[i-1][ss])%mod;
}
}
}
cout<<dp[n][maxn]<<endl;
return 0;
}