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