[插头DP] BZOJ2331 && SCOI2011 地板

插头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;}