bzoj 1814 Ural 1519 Formula 1 ——插头DP

题目：https://www.lydsy.com/JudgeOnline/problem.php?id=1814

普通的插头 DP 。但是调了很久。注意如果合并两个 1 的话，不是 “把向右第一个 2 该成 1 ”，而是 “把向右第一个没有与 1 匹配的 2 改成 1 ”。

原来获取哈希值是用字符串哈希的方法，遍历12个位置；太慢。直接对某数取模作为哈希值，手写哈希表保证不会找错状态。大概 1e5 个状态？

在转移的时候看一下下方和右边有没有障碍、只做合法的转移的话，取答案的时候就不用在判断 “当前两个位置是 1 和 2 ”之外再判断 “其余位置都是 0 ” 了。

自己写的是 “如果当前位置是结束位置且当前两个插口是 1 和 2 就输出答案并 break ” 。这样无解的时候要手动输出 0 。

滚动数组。并且只遍历有效的状态。方法是用手写队列存 “得到转移” 的状态，用 vis 判断该状态是否已经在队列里；把下一层的队列做好之后，遍历下一层（不是这一层）的队列把 vis 清空即可。

注意 long long 。

哈希表的映射方式是：原状态 -> 哈希值 -> ++tot 地分配空间 -> 原状态 。“ -> ” 就是要记下的对应关系。因为哈希值要用数组记它对应哪个位置，所以哈希值不宜太大。

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const int N=,M=1e5+;//mention
int n,m,bin[N],q[][M],t[]; ll dp[][M];
bool b[N][N],vis[M],en[N][N];
namespace H{
  const int md=1e6+;//
  int hd[md+],xnt,nxt[M],vl[M];
  int Ps(int s)
  {
    int h=s%md;
    for(int i=hd[h];i;i=nxt[i])
      if(vl[i]==s){return i;}
    vl[++xnt]=s; nxt[xnt]=hd[h]; hd[h]=xnt;
    return xnt;
  }
}
int get(int cr,int j){ return (cr&(bin[j+]-bin[j]))>>(*j);}
void cz(int d,bool v){ if(!vis[d]) vis[d]=,q[v][++t[v]]=d;}
void solve()
{
  bin[]=; int lm=m+;
  for(int i=;i<=lm;i++)bin[i]=bin[i-]<<;
  bool flag=;
  for(int i=n;i&&!flag;i--)
    for(int j=m;j;j--)
      if(!b[i][j]){en[i][j]=flag=;break;}
  bool u=,v=; flag=;
  dp[u][H::Ps()]=; q[u][++t[u]]=;
  for(int i=;i<=n;i++)
    {
      for(int j=;j<=m;j++)
    {
      for(int c=;c<=t[u];c++)
        {
          int cr=H::vl[q[u][c]];ll tp=dp[u][q[u][c]];//ll!!!
          int d0=get(cr,j-), d1=get(cr,j), t0,t1;
          if(b[i][j])
        {
          if(!d0&&!d1)
            { t0=H::Ps(cr); cz(t0,v); dp[v][t0]+=tp;}
          continue;
        }
          if(en[i][j])
        {
          if(d0==&&d1==)
            { printf("%lld\n",tp);flag=;break;}
          continue;
        }
          if((!d0&&d1)||(!d1&&d0))
        {
          t0=cr; t1=cr+d1*bin[j-]-d0*bin[j-]+d0*bin[j]-d1*bin[j];
          if(i<n&&!b[i+][j])
            {if(d0){t0=H::Ps(t0); cz(t0,v); dp[v][t0]+=tp;}
              else{t1=H::Ps(t1); cz(t1,v); dp[v][t1]+=tp;}}
          if(j<m&&!b[i][j+])
            {if(d0){t1=H::Ps(t1); cz(t1,v); dp[v][t1]+=tp;}
              else{t0=H::Ps(t0); cz(t0,v); dp[v][t0]+=tp;}}
          continue;
        }
          if(!d0&&!d1&&i<n&&j<m&&!b[i+][j]&&!b[i][j+])
        {
          t0=cr+bin[j-]+*bin[j];
          t0=H::Ps(t0); cz(t0,v);
          dp[v][t0]+=tp; continue;
        }
          if(d0==&&d1==)
        {
          t0=cr-bin[j-]-bin[j];
          for(int p=j+,top=,d;p<=m;p++)
            {
              d=get(t0,p);
              if(d==)top++; else if(d==&&top)top--;
              else if(d==){t0-=bin[p];break;}
            }
          t0=H::Ps(t0); cz(t0,v); dp[v][t0]+=tp; continue;
        }
          if(d0==&&d1==)
        {
          t0=cr-*bin[j-]-*bin[j];//2*!!!
          for(int p=j-,top=,d;p>=;p--)
            {
              d=get(t0,p);
              if(d==)top++; else if(d==&&top)top--;
              else if(d==){t0+=bin[p];break;}
            }
          t0=H::Ps(t0); cz(t0,v); dp[v][t0]+=tp;
          continue;
        }
          if(d0==&&d1==)
        {
          t0=cr-*bin[j-]-bin[j];
          t0=H::Ps(t0); cz(t0,v); dp[v][t0]+=tp; continue;
        }
        }
      if(flag)break;
      for(int c=;c<=t[u];c++)dp[u][q[u][c]]=;
      for(int c=;c<=t[v];c++)vis[q[v][c]]=;//v not u!!!
      t[u]=; swap(u,v);
    }
      if(flag)break;
      for(int c=;c<=t[u];c++)
    {
      int cr=H::vl[q[u][c]];ll tp=dp[u][q[u][c]];//ll!!!
      if(get(cr,m))continue;
      int d=cr<<;
      d=H::Ps(d); cz(d,v); dp[v][d]=tp;
    }
      for(int c=;c<=t[u];c++)dp[u][q[u][c]]=;
      for(int c=;c<=t[v];c++)vis[q[v][c]]=;
      t[u]=; swap(u,v);
    }
  if(!flag)puts("");///
}
int main()
{
  scanf("%d%d",&n,&m); char ch[N];
  for(int i=;i<=n;i++)
    {
      scanf("%s",ch+);
      for(int j=;j<=m;j++)
    b[i][j]=(ch[j]=='*');
    }
  solve();
  return ;
}