题意:给一张地图,1代表这个地方可以放,0代表不可以,可以放c-d个1*1的方块,其余的只能放1*2的,问有多少种方法。
插头DP,轮廓线DP,每一次都是从上一个状态转移到当前状态。。
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int MOD=1e9+7;
int dp[2][22][1050];
char mp[105][15];
int now,pre;
int main()
{
int i,j,k,d,c,s,n,m;
//freopen("in.txt","r",stdin);
while(scanf("%d%d%d%d",&n,&m,&c,&d)==4)
{
for(i=0;i<n;i++)
scanf("%s",mp[i]);
now=1,pre=0;
memset(dp[now],0,sizeof(dp[now]));
int all=(1<<m)-1;
dp[now][0][all]=1;
for(i=0;i<n;i++)
{
for(j=0;j<m;j++)
{
swap(now,pre);
memset(dp[now],0,sizeof(dp[now]));
if(mp[i][j]=='1')
{
for(k=0;k<=d;k++)
{
for(s=0;s<=all;s++)
{
if(k&&(s&(1<<j))) //放1*1的方块,上一个状态是s,当前状态和上一个状态一样
dp[now][k][s]=(dp[now][k][s]+dp[pre][k-1][s])%MOD;
if(!dp[pre][k][s])
continue;
if(j&&!(s&(1<<(j-1)))&&(s&(1<<j))) //横着放1*2的
dp[now][k][s|(1<<(j-1))]=(dp[now][k][s|(1<<(j-1))]+dp[pre][k][s])%MOD;
dp[now][k][s^(1<<j)]=(dp[now][k][s^(1<<j)]+dp[pre][k][s])%MOD;//竖着放1*2的,如果上一个状态这个格子有就放不了,没有就能放,所以刚好用异或表示
}
}
}
else
{
for(k=0;k<=d;k++)
{
for(s=0;s<=all;s++)
{
if((s&(1<<j)))
dp[now][k][s]=(dp[now][k][s]+dp[pre][k][s])%MOD;
}
}
}
}
}
int ans=0;
for(i=c;i<=d;i++)
{
ans=(ans+dp[now][i][all])%MOD;
}
printf("%d\n",ans);
}
return 0;
}