[插头DP] BZOJ2331 && SCOI2011 地板

最新推荐文章于 2018-09-15 11:21:58 发布

LowestJN

最新推荐文章于 2018-09-15 11:21:58 发布

阅读量398

点赞数

分类专栏：插头DP 文章标签：插头DP

本文链接：https://blog.csdn.net/Coldef/article/details/77056678

版权

插头DP 专栏收录该内容

2 篇文章 0 订阅

订阅专栏

插头DP裸题

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <string>

using namespace std;

const int N=60010,P=20110520;

int n,m;
char a[110][110],b[110][110];
int f[2][15][3][N];
int pw[110];

inline void turns(){
  swap(n,m);
  for(int i=1;i<=n;i++)
    for(int j=1;j<=m;j++)
      b[i][j]=a[j][i];
  memcpy(a,b,sizeof(b));
}

inline void add(int &x,int y){
  if((x+=y)>=P) x-=P;
}

inline int cg(int s,int p,int t){
  return ((s/pw[p])*3+t)*pw[p-1]+s%pw[p-1];
}

int main(){
  scanf("%d%d",&n,&m);
  for(int i=1;i<=n;i++)
    scanf("%s",a[i]+1);
  if(n<m) turns();
  pw[0]=1;
  for(int i=1;i<=m;i++) pw[i]=pw[i-1]*3;
  int st=0; f[st][0][0][0]=1;
  for(int i=1;i<=n;i++,st^=1){
    for(int j=0;j<pw[m];j++)
      add(f[st][0][0][j],f[st^1][m][0][j]);
    memset(f[st^1],0,sizeof(f[st^1]));
    for(int j=0;j<m;j++)
      for(int S=0;S<pw[m];S++)
    for(int t=0;t<3;t++){
      int s=(S/pw[j])%3,cur=f[st][j][t][S];
      if(cur==0) continue;
      int (*nxt)[N]=f[st][j+1];
      if(a[i][j+1]=='*'){
        if(s!=0 || t!=0) continue;
        add(nxt[0][cg(S,j+1,0)],cur);
      }
      else{
        if(s==0 && t==0)
          add(nxt[0][cg(S,j+1,1)],cur),add(nxt[1][cg(S,j+1,0)],cur),add(nxt[2][cg(S,j+1,2)],cur);
        if(t==1 && s==0)
          add(nxt[1][cg(S,j+1,0)],cur),add(nxt[0][cg(S,j+1,2)],cur);
        if(t==2 && s==0)
          add(nxt[2][cg(S,j+1,0)],cur),add(nxt[0][cg(S,j+1,0)],cur);
        if(s==1 && t==0)
          add(nxt[0][cg(S,j+1,1)],cur),add(nxt[2][cg(S,j+1,0)],cur);
        if(s==2 && t==0)
          add(nxt[0][cg(S,j+1,2)],cur),add(nxt[0][cg(S,j+1,0)],cur);
        if(s==1 && t==1)
          add(nxt[0][cg(S,j+1,0)],cur);
      }
      }
  }
  printf("%d\n",f[st^1][m][0][0]);
  return 0;
}