2331: [SCOI2011]地板

Time Limit: 5 Sec  Memory Limit: 128 MB
Submit: 541  Solved: 239
[Submit][Status]

Description

lxhgww的小名叫“小L”,这是因为他总是很喜欢L型的东西。小L家的客厅是一个bzoj 2331: [SCOI2011]地板 插头DP-LMLPHP的矩形,现在他想用L型的地板来铺满整个客厅,客厅里有些位置有柱子,不能铺地板。现在小L想知道,用L型的地板铺满整个客厅有多少种不同的方案?

需要注意的是,如下图所示,L型地板的两端长度可以任意变化,但不能长度为0。铺设完成后,客厅里面所有没有柱子的地方都必须铺上地板,但同一个地方不能被铺多次。

bzoj 2331: [SCOI2011]地板 插头DP-LMLPHP

Input

输入的第一行包含两个整数,R和C,表示客厅的大小。

接着是R行,每行C个字符。’_’表示对应的位置是空的,必须铺地板;’*’表示对应的位置有柱子,不能铺地板。

Output

输出一行,包含一个整数,表示铺满整个客厅的方案数。由于这个数可能很大,只需输出它除以20110520的余数。

Sample Input

2 2

*_

__

Sample Output

1

HINT

R*C<=100

  太久没做插头DP,一来就把R*C<=100看成R,C<=100,这道题主要问题还是在数据范围里面,交wa了一次,因为我默认R最大为10,而实际上R可以取到100.

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define MOD 20110520
#define MAXN 12
#define MAXV 531442
#define gv(x,ps) ((x)/p[(ps)]%3)
#define cd(x,ps) ((x)-((x)/p[ps]%3*p[ps]))
#define sd(x,ps,v) ((x)+v*p[ps])
int n,m;
void pm(int x)
{
int i;
for (i=;i<=m;i++)
{
printf("%d",x%);
x/=;
}
printf("\n");
}
void deal(int &x,int y)
{
x=(x+y)%MOD;
}
char mp[][];
int dp[][][MAXV];
int p[];
int main()
{
freopen("input.txt","r",stdin);
int i,j,k,x,y,z;
p[]=;
scanf("%d%d\n",&n,&m);
for (i=;i<;i++)
p[i]=p[i-]*;
for (i=;i<n;i++)
{
for (j=;j<m;j++)
{
scanf("%c",&mp[i][j]);
}
scanf("\n");
}
int l;
if (n<m)
{
for (i=;i<max(n,m);i++)
{
for (j=i+;j<max(n,m);j++)
{
swap(mp[i][j],mp[j][i]);
}
}
swap(n,m);
}
l=p[m+];
dp[][][]=;
int v;
for (i=;i<n;i++)
{
for (j=;j<m;j++)
{
if (j==m-)//{{{
{
if (mp[i][j]=='*')
{
for (k=;k<l;k++)
{
if (!dp[i&][j][k])continue;
x=gv(k,j);
y=gv(k,j+);
if (x||y)continue;
deal(dp[(i+)&][][k*],dp[i&][j][k]);
}
memset(dp[i&],,sizeof(dp[i&]));
}else
{
for (k=;k<l;k++)
{
if (!dp[i&][j][k])continue;
x=gv(k,j);
y=gv(k,j+);
v=cd(k,j);
v=cd(v,j+);
if (x== && y==)
{
z=sd(v,j,);
deal(dp[(i+)&][][z*],dp[i&][j][k]);
}else if (x== && y==)
{
z=sd(v,j,);
deal(dp[(i+)&][][z*],dp[i&][j][k]);
}else if (x== && y==)
{
z=sd(v,j,);
deal(dp[(i+)&][][z*],dp[i&][j][k]);
z=v;
deal(dp[(i+)&][][z*],dp[i&][j][k]);
}else if (x== && y==)
{
z=sd(v,j,);
deal(dp[(i+)&][][z*],dp[i&][j][k]);
}else if (x== && y==)
{
z=v;
deal(dp[(i+)&][][z*],dp[i&][j][k]);
}else if (x== && y==)
{
//DO nothing
}else if (x== && y==)
{
z=v;
deal(dp[(i+)&][][z*],dp[i&][j][k]);
}else if (x== && y==)
{
//Do nothing
}else if (x== && y==)
{
//Do nothing
}
}
memset(dp[i&],,sizeof(dp[i&]));
}
}else//}}}
{
if (mp[i][j]=='*')
{
for (k=;k<l;k++)
{
if (!dp[i&][j][k])continue;
x=gv(k,j);
y=gv(k,j+);
if (x||y)continue;
deal(dp[i&][j+][k],dp[i&][j][k]);
}
}else
{
for (k=;k<l;k++)
{
if (!dp[i&][j][k])continue;
x=gv(k,j);
y=gv(k,j+);
v=cd(k,j);
v=cd(v,j+);
if (x== && y==)
{
z=sd(v,j+,);
deal(dp[i&][j+][z],dp[i&][j][k]);
z=sd(v,j,);
deal(dp[i&][j+][z],dp[i&][j][k]);
z=sd(v,j+,);
z=sd(z,j,);
deal(dp[i&][j+][z],dp[i&][j][k]);
}else if (x== && y==)
{
z=sd(v,j+,);
deal(dp[i&][j+][z],dp[i&][j][k]);
z=sd(v,j,);
deal(dp[i&][j+][z],dp[i&][j][k]);
}else if (x== && y==)
{
z=v;
deal(dp[i&][j+][z],dp[i&][j][k]);
z=sd(v,j,);
deal(dp[i&][j+][z],dp[i&][j][k]);
}else if (x== && y==)
{
z=sd(v,j+,);
deal(dp[i&][j+][z],dp[i&][j][k]);
z=sd(v,j,);
deal(dp[i&][j+][z],dp[i&][j][k]);
}else if (x== && y==)
{
z=v;
deal(dp[i&][j+][z],dp[i&][j][k]);
}else if (x== && y==)
{
//Do nothing
}else if (x== && y==)
{
z=sd(v,j+,);
deal(dp[i&][j+][z],dp[i&][j][k]);
z=v;
deal(dp[i&][j+][z],dp[i&][j][k]);
}else if (x== && y==)
{
//Do nothing
}else if (x== && y==)
{
//Do nothing
}
}
}
}
}
}
cout<<dp[n&][][]<<endl;
return ;
}
04-28 02:23