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