[bzoj2331][SCOI2011]地板

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$。
$0$表示没有，$1$表示没有折过的，$2$表示折过的。

#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]);
}