2331: [SCOI2011]地板
Time Limit: 5 Sec Memory Limit: 128 MB
Submit: 909 Solved: 382
[Submit][Status][Discuss]
Description
lxhgww的小名叫“小L”,这是因为他总是很喜欢L型的东西。小L家的客厅是一个的矩形,现在他想用L型的地板来铺满整个客厅,客厅里有些位置有柱子,不能铺地板。现在小L想知道,用L型的地板铺满整个客厅有多少种不同的方案?
需要注意的是,如下图所示,L型地板的两端长度可以任意变化,但不能长度为0。铺设完成后,客厅里面所有没有柱子的地方都必须铺上地板,但同一个地方不能被铺多次。
Input
输入的第一行包含两个整数,R和C,表示客厅的大小。
接着是R行,每行C个字符。’_’表示对应的位置是空的,必须铺地板;’*’表示对应的位置有柱子,不能铺地板。
Output
输出一行,包含一个整数,表示铺满整个客厅的方案数。由于这个数可能很大,只需输出它除以20110520的余数。
Sample Input
2 2
*_
__
Sample Output
1
HINT
R*C<=100
f[i][j][s]
表示在格子
(i,j)
当前的状态是
s
。
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
#define Mod 20110520
const int N=110;
const int P=50000;
const int M=200000;
int n,m,map[N][N],now,f[2][M],hash[2][M],tot[2],point[P],next[M],bin[N],pre;
inline void add(int s,int value){
int i,x=s%P;
for(i=point[x];i;i=next[i])
if(hash[now][i]==s){
f[now][i]=(f[now][i]+value)%Mod;
return ;
}
f[now][++tot[now]]=value;
hash[now][tot[now]]=s;
next[tot[now]]=point[x];
point[x]=tot[now];
}
inline void dp(){
int i,j,k,p,q,s,value;
tot[now]=1;
hash[now][1]=0;
f[now][1]=1;
for(i=0;i<N;++i) bin[i]=i<<1;
for(i=1;i<=n;++i){
for(j=1;j<=tot[now];++j)
hash[now][j]=(hash[now][j]<<2)&((1<<bin[m+1])-1);
for(j=1;j<=m;++j){
pre=now;now^=1;
tot[now]=0;
memset(point,0,sizeof(point));
memset(f[now],0,sizeof(f[now]));
for(k=1;k<=tot[pre];++k)
if(f[pre][k]){
s=hash[pre][k];
value=f[pre][k];
p=(s>>bin[j-1])&3;
q=(s>>bin[j])&3;
if(!map[i][j]){
if(!p&&!q) add(s,value);
}
else if(!p&&!q){
if(map[i+1][j]==1) add(s+(1<<bin[j-1]),value);
if(map[i][j+1]==1) add(s+(1<<bin[j]),value);
if(map[i+1][j]==1&&map[i][j+1]==1) add(s+(2<<bin[j-1])+(2<<bin[j]),value);
}
else if(!p){
s-=(1<<bin[j])*q;
if(q==1){
if(map[i+1][j]==1) add(s+(1<<bin[j-1]),value);
if(map[i][j+1]==1) add(s+(2<<bin[j]),value);
}
else{
add(s,value);
if(map[i+1][j]==1) add(s|(2<<bin[j-1]),value);
}
}
else if(!q){
s-=(1<<bin[j-1])*p;
if(p==1){
if(map[i+1][j]==1) add(s+(2<<bin[j-1]),value);
if(map[i][j+1]==1) add(s+(1<<bin[j]),value);
}
else{
add(s,value);
if(map[i][j+1]==1) add(s+(2<<bin[j]),value);
}
}
else if(p==1&&q==1){
s-=((1<<bin[j-1])+(1<<bin[j]));
add(s,value);
}
}
}
}
}
int main(){
int i,j;
scanf("%d%d",&n,&m);
memset(map,127,sizeof(map));
if(n>m){
for(i=1;i<=n;++i)
for(j=1;j<=m;++j){
char ch=getchar();
while(ch!='*'&&ch!='_') ch=getchar();
map[i][j]=(ch=='_');
}
}
else{
for(i=1;i<=n;++i)
for(j=1;j<=m;++j){
char ch=getchar();
while(ch!='*'&&ch!='_') ch=getchar();
map[j][i]=(ch=='_');
}
swap(n,m);
}
dp();
printf("%d\n",f[now][1]);
}