初探插头dp

开学那个月学了点新东西，不知道还记不记得了，mark一下

感觉cdq的论文讲的很详细

题主要跟着kuangbin巨做了几道基础的

http://www.cnblogs.com/kuangbin/archive/2012/10/02/2710343.html

还有几道没做，留着坑

感觉广义括号表示法虽然神奇，但一般最小表示法就够用了吧，感觉最小表示法更直观一点

hdu1693

#include<cstdio>
#include<iostream>
#include<cstring>
typedef long long ll;
using namespace std;
ll f[13][13][1<<13];
int a[13][13];
int n,m,cas;

int main()
{
  scanf("%d",&cas);
  for (int tt=1; tt<=cas; tt++)
  {
    scanf("%d%d",&n,&m);
    for (int i=1; i<=n; i++)
      for (int j=1; j<=m; j++)
        scanf("%d",&a[i][j]);
    printf("Case %d: ",tt);
    memset(f,0,sizeof(f));
    f[0][m][0]=1;
    for (int i=1; i<=n; i++)
    {
      for (int j=0; j<(1<<m); j++)
        f[i][0][j<<1]=f[i-1][m][j];
      for (int j=1; j<=m; j++)
      {
        for (int st=0; st<(1<<(m+1)); st++)
        {
          int x=1<<j,y=1<<(j-1);
          if (a[i][j])
          {
            if ((st&x)&&(st&y)) f[i][j][st]=f[i][j-1][st^x^y];
            else if (!(st&x)&&!(st&y)) f[i][j][st]=f[i][j-1][st|x|y];
            else f[i][j][st]=f[i][j-1][st]+f[i][j-1][st^x^y];
          }
          else {
            if (!(st&x)&&!(st&y)) f[i][j][st]=f[i][j-1][st];
            else f[i][j][st]=0;
          }
        }
      }
    }
    printf("There are %lld ways to eat the trees.
",f[n][m][0]);
  }
  return 0;
}

poj1793（男人八题!）

#include<iostream>
#include<cstdio>
#include<cstring>
typedef long long ll;
const int mo=30007;
const int maxl=1000010;
using namespace std;
struct node{
  int p[mo],nex[maxl],len;
  ll st[maxl],f[maxl];
  void clr()
  {
    len=0; memset(p,255,sizeof(p));
  }
  void push(ll nw,ll s)
  {
    int x=nw%mo;
    for (int i=p[x];i!=-1; i=nex[i])
      if (st[i]==nw)
      {
        f[i]+=s;
        return;
      }
    st[++len]=nw; f[len]=s;
    nex[len]=p[x]; p[x]=len;
  }
} h[2];
int n,m,cas,p,ex,ey;
int a[15][15],b[15],v[15];

void get(ll st)
{
  for (int i=m; i>=0; i--)
  {
    b[i]=st&7;
    st>>=3;
  }
}

ll put()
{
  memset(v,255,sizeof(v)); v[0]=0;
  ll st=0; int t=0;
  for (int i=0; i<=m; i++)
  {
    if (v[b[i]]==-1) v[b[i]]=++t;
    b[i]=v[b[i]];
    st<<=3; st|=b[i];
  }
  return st;
}

void shift()
{
  for (int i=m;i; i--) b[i]=b[i-1];
  b[0]=0;
}

void work(int j,int k)
{
  if (j==m) shift();
  h[p].push(put(),h[1-p].f[k]);
}

void dp(int i,int j)
{
  for (int k=1; k<=h[1-p].len; k++)
  {
    get(h[1-p].st[k]);
    if (!a[i][j])
    {
      b[j]=b[j-1]=0; work(j,k);
      continue;
    }
    int x=b[j],y=b[j-1];
    if (x&&y)
    {
      if ((x==y&&i==ex&&j==ey)||(x!=y))
      {
        b[j]=b[j-1]=0;
        if (x!=y)
          for (int r=0; r<=m; r++) if (b[r]==x) b[r]=y;
        work(j,k);
      }
    }
    else if ((x&&!y)||(!x&&y))
    {
      int r=x+y;
      if (a[i][j+1]) {b[j]=r; b[j-1]=0; work(j,k);}
      if (a[i+1][j]) {b[j-1]=r; b[j]=0; work(j,k);}
    }
    else if (a[i][j+1]&&a[i+1][j])
    {
      b[j]=b[j-1]=m+1;
      work(j,k);
    }
  }
}

int main()
{
  char s[12];
  scanf("%d%d",&n,&m);
  while (n&&m)
  {
    memset(a,0,sizeof(a));
    for (int i=1; i<=n; i++)
    {
      scanf("%s",s+1);
      for (int j=1; j<=m; j++)
        if (s[j]=='.') a[i][j]=1;
    }
    for (int i=1; i<=m; i++)
    {
      a[n+2][i]=1;
      a[n+1][i]=(i>1&&i<m)?0:1;
    }
    n+=2;  ex=n; ey=m;
    h[0].clr();
    h[0].push(0,1);  p=0;
    for (int i=1; i<=n; i++)
      for (int j=1; j<=m; j++)
      {
        p^=1; h[p].clr();
        dp(i,j);
      }
    ll ans=0;
    for (int i=1; i<=h[p].len; i++)
      ans+=h[p].f[i];
    printf("%lld
",ans);
    scanf("%d%d",&n,&m);
  }
  return 0;
}

fzu1977

#include<iostream>
#include<cstdio>
#include<cstring>
typedef long long ll;
const int mo=30007;
const int maxl=300010;
using namespace std;
struct node{
  int p[mo],nex[maxl],len;
  ll st[maxl],f[maxl];
  void clr()
  {
    len=0; memset(p,255,sizeof(p));
  }
  void push(ll nw,ll s)
  {
    int x=nw%mo;
    for (int i=p[x];i!=-1; i=nex[i])
      if (st[i]==nw)
      {
        f[i]+=s;
        return;
      }
    st[++len]=nw; f[len]=s;
    nex[len]=p[x]; p[x]=len;
  }
} h[2];
int n,m,cas,p,ex,ey;
int a[15][15],b[15],v[15];

void get(ll st)
{
  b[m+1]=st&1; st>>=1;
  for (int i=m; i>=0; i--)
  {
    b[i]=st&7;
    st>>=3;
  }
}

ll put()
{
  memset(v,255,sizeof(v)); v[0]=0;
  ll st=0; int t=0;
  for (int i=0; i<=m; i++)
  {
    if (v[b[i]]==-1) v[b[i]]=++t;
    b[i]=v[b[i]];
    st<<=3; st|=b[i];
  }
  st<<=1; st|=b[m+1];
  return st;
}

void shift()
{
  for (int i=m;i; i--) b[i]=b[i-1];
  b[0]=0;
}

void dp(int i,int j)
{
  for (int k=1; k<=h[p^1].len; k++)
  {
    get(h[p^1].st[k]);
    if (!a[i][j])
    {
      b[j]=b[j-1]=0;
      h[p].push(put(),h[p^1].f[k]);
      continue;
    }
    int x=b[j],y=b[j-1];
    if (b[m+1])
    {
      if (!x&&!y&&a[i][j]!=1)
      {
        b[j]=b[j-1]=0;
        h[p].push(put(),h[p^1].f[k]);
      }
      continue;
    }
    if (x&&y)
    {
      b[j]=b[j-1]=0;
      if (x!=y){
        for (int r=0; r<=m; r++) if (b[r]==x) b[r]=y;
      }
      else b[m+1]=1;
      h[p].push(put(),h[p^1].f[k]);
    }
    else if ((x&&!y)||(!x&&y))
    {
      int r=x+y;
      if (a[i][j+1])
      {
        b[j]=r; b[j-1]=0;
        h[p].push(put(),h[p^1].f[k]);
      }
      if (a[i+1][j])
      {
        b[j-1]=r; b[j]=0;
        h[p].push(put(),h[p^1].f[k]);
      }
    }
    else {
      if (a[i][j+1]&&a[i+1][j])
      {
        b[j]=b[j-1]=m+1;
        h[p].push(put(),h[p^1].f[k]);
      }
      if (a[i][j]==2)
      {
        b[j]=b[j-1]=0;
        h[p].push(put(),h[p^1].f[k]);
      }
    }
  }
}

int main()
{
  char s[15]; int cas;
  scanf("%d",&cas);
  for (int tt=1; tt<=cas; tt++)
  {
    scanf("%d%d",&n,&m);
    memset(a,0,sizeof(a));
    for (int i=1; i<=n; i++)
    {
      scanf("%s",s+1);
      for (int j=1; j<=m; j++)
        if (s[j]=='O') a[i][j]=1;
        else if (s[j]=='*') a[i][j]=2;
    }
    h[0].clr();
    h[0].push(0,1);  p=0;
    for (int i=1; i<=n; i++)
    {
      p^=1; h[p].clr();
      for (int k=1; k<=h[p^1].len; k++)
      {
        get(h[p^1].st[k]); shift();
        h[p].push(put(),h[p^1].f[k]);
      }
      for (int j=1; j<=m; j++)
      {
        p^=1; h[p].clr();
        dp(i,j);
      }
    }
    ll ans=0;
    for (int i=1; i<=h[p].len; i++)
      ans+=h[p].f[i];
    printf("Case %d: %I64d
",tt,ans);
  }
  return 0;
}

hdu3377

#include<iostream>
#include<cstdio>
#include<cstring>
typedef long long ll;
const int mo=30007;
const int maxl=1000010;
const int inf=-100000007;
using namespace std;
struct node{
  int p[mo],nex[maxl],len,f[maxl];
  ll st[maxl];
  void clr()
  {
    len=0; memset(p,255,sizeof(p));
  }
  void push(ll nw,int s)
  {
    int x=nw%mo;
    for (int i=p[x];i!=-1; i=nex[i])
      if (st[i]==nw)
      {
        f[i]=max(f[i],s);
        return;
      }
    st[++len]=nw; f[len]=s;
    nex[len]=p[x]; p[x]=len;
  }
} h[2];
int n,m,cas,p,ex,ey;
int a[15][15],b[15],v[15];

void get(ll st)
{
  for (int i=m; i>=0; i--)
  {
    b[i]=st&7;
    st>>=3;
  }
}

ll put()
{
  memset(v,255,sizeof(v)); v[0]=0;
  ll st=0; int t=0;
  for (int i=0; i<=m; i++)
  {
    if (v[b[i]]==-1) v[b[i]]=++t;
    b[i]=v[b[i]];
    st<<=3; st|=b[i];
  }
  return st;
}

void shift()
{
  for (int i=m;i; i--) b[i]=b[i-1];
  b[0]=0;
}

void dp(int i,int j)
{
  for (int k=1; k<=h[1-p].len; k++)
  {
    get(h[1-p].st[k]);
    int x=b[j],y=b[j-1];
    if (i==1&&j==1)
    {
      if (i<n) {b[j-1]=m+1; b[j]=0; h[p].push(put(),h[1-p].f[k]+a[i][j]);}
      if (j<m) {b[j-1]=0; b[j]=m+1; h[p].push(put(),h[1-p].f[k]+a[i][j]);}
      continue;
    }
    if (i==n&&j==m)
    {
      if ((!x&&y)||(!y&&x))
      {
        b[j]=b[j-1]=0;
        h[p].push(put(),h[1-p].f[k]+a[i][j]);
      }
      continue;
    }
    if (x&&y)
    {
      if (x!=y)
      {
        b[j]=b[j-1]=0;
        if (x!=y)
          for (int r=0; r<=m; r++) if (b[r]==x) b[r]=y;
        if (j==m) shift();
        h[p].push(put(),h[1-p].f[k]+a[i][j]);
      }
    }
    else if ((x&&!y)||(!x&&y))
    {
      int r=x+y;
      if (j<m) {
        b[j]=r; b[j-1]=0;
        h[p].push(put(),h[1-p].f[k]+a[i][j]);
      }
      if (i<n)
      {
        b[j-1]=r; b[j]=0;
        if (j==m) shift();
        h[p].push(put(),h[1-p].f[k]+a[i][j]);
      }
    }
    else {
      if (i<n&&j<m)
      {
        b[j]=b[j-1]=m+1;
        h[p].push(put(),h[1-p].f[k]+a[i][j]);
      }
      b[j]=b[j-1]=0;
      if (j==m) shift();
      h[p].push(put(),h[1-p].f[k]);
    }
  }
}

int main()
{
  int tt=0;
  while (scanf("%d%d",&n,&m)!=EOF)
  {
    for (int i=1; i<=n; i++)
      for (int j=1; j<=m; j++)
        scanf("%d",&a[i][j]);
    if (n==1&&m==1)
    {
      printf("Case %d: %d
",++tt,a[1][1]);
      continue;
    }
    h[0].clr();
    h[0].push(0,0);  p=0;
    for (int i=1; i<=n; i++)
      for (int j=1; j<=m; j++)
      {
        p^=1; h[p].clr();
        dp(i,j);
      }
    int ans=inf;
    for (int i=1; i<=h[p].len; i++)
      ans=max(ans,h[p].f[i]);
    printf("Case %d: %d
",++tt,ans);
  }
  return 0;
}

poj3133

#include<iostream>
#include<cstdio>
#include<cstring>

const int maxl=1100010;
using namespace std;
struct node{
  int len,f[maxl],st[maxl];
  void clr()
  {
    for (int i=1; i<=len; i++) f[st[i]]=-1;
    len=0;
  }
  void push(int nw,int s)
  {
    if (f[nw]==-1)
    {
      st[++len]=nw;
      f[nw]=s;
    }
    else f[nw]=min(f[nw],s);
  }
} h[2];
int n,m,cas,p;
int a[15][15],b[15];

void get(int st)
{
  for (int i=m; i>=0; i--)
  {
    b[i]=st&3;
    st>>=2;
  }
}

int put()
{
  int st=0;
  for (int i=0; i<=m; i++) {st<<=2; st|=b[i];}
  return st;
}

void shift()
{
  for (int i=m;i; i--) b[i]=b[i-1];
  b[0]=0;
}

void dp(int i,int j)
{
  for (int k=1; k<=h[1-p].len; k++)
  {
    get(h[1-p].st[k]);
    int x=b[j],y=b[j-1],st=h[1-p].st[k];
    if (a[i][j]!=1)
    {
      if (!a[i][j]&&!x&&!y) h[p].push(st,h[1-p].f[st]);
      else if (a[i][j]>1)
      {
        if (((x&&!y)||(!x&&y))&&(x+y==a[i][j])) {b[j]=b[j-1]=0; h[p].push(put(),h[1-p].f[st]+1);}
        else if (!x&&!y)
        {
          if (a[i][j+1]==a[i][j]||a[i][j+1]==1) {b[j]=a[i][j]; b[j-1]=0; h[p].push(put(),h[1-p].f[st]+1);}
          if (a[i+1][j]==a[i][j]||a[i+1][j]==1) {b[j-1]=a[i][j]; b[j]=0; h[p].push(put(),h[1-p].f[st]+1);}
        }
      }
      continue;
    }
    if (x&&y&&x==y) {b[j]=b[j-1]=0; h[p].push(put(),h[1-p].f[st]+1);}
    else if ((x&&!y)||(!x&&y))
    {
      int r=x+y;
      if (a[i][j+1]==1||a[i][j+1]==r) {b[j]=r; b[j-1]=0; h[p].push(put(),h[1-p].f[st]+1);}
      if (a[i+1][j]==1||a[i+1][j]==r) {b[j-1]=r; b[j]=0; h[p].push(put(),h[1-p].f[st]+1);}
    }
    else if (!x&&!y)
    {
      if (a[i][j+1]&&a[i+1][j])
      {
        if (a[i][j+1]==1&&a[i+1][j]==1)
        {
          b[j]=b[j-1]=2;  h[p].push(put(),h[1-p].f[st]+1);
          b[j]=b[j-1]=3;  h[p].push(put(),h[1-p].f[st]+1);
        }
        else if (a[i][j+1]==2&&a[i+1][j]==1||a[i][j+1]==1&&a[i+1][j]==2||a[i+1][j]==2&&a[i][j+1]==2) {b[j]=b[j-1]=2; h[p].push(put(),h[1-p].f[st]+1);}
        else if (a[i][j+1]==3&&a[i+1][j]==1||a[i][j+1]==1&&a[i+1][j]==3||a[i+1][j]==3&&a[i][j+1]==3) {b[j]=b[j-1]=3; h[p].push(put(),h[1-p].f[st]+1);}
      }
      h[p].push(st,h[1-p].f[st]);
    }
  }
}

int main()
{
  h[1].len=h[0].len=0;
  memset(h[0].f,255,sizeof(h[0].f));
  memset(h[1].f,255,sizeof(h[1].f));
  scanf("%d%d",&n,&m);
  while (n&&m)
  {
    memset(a,0,sizeof(a));
    for (int i=1; i<=n; i++)
    {
      for (int j=1; j<=m; j++)
      {
        scanf("%d",&a[i][j]);
        if (a[i][j]==0||a[i][j]==1) a[i][j]^=1;
      }
    }
    h[0].push(0,0);  p=0;
    for (int i=1; i<=n; i++)
    {
      p^=1;  h[p].clr();
      for (int k=1; k<=h[1-p].len; k++)
      {
        get(h[1-p].st[k]); int st=h[1-p].st[k]; shift();
        h[p].push(put(),h[1-p].f[st]);
      }
      for (int j=1; j<=m; j++)
      {
        p^=1; h[p].clr();
        dp(i,j);
      }
    }
    int ans=0;
    for (int i=1; i<=h[p].len; i++)
    {
      int st=h[p].st[i];
      ans+=h[p].f[st];
    }
    if (ans>0) ans-=2;
    printf("%d
",ans);
    h[1].clr(); h[0].clr();
    scanf("%d%d",&n,&m);
  }
  return 0;
}

hdu4285

#include<iostream>
#include<cstdio>
#include<cstring>
typedef long long ll;
const int ha=30007;
const int maxl=300010;
const int mo=1000000007;
using namespace std;
struct node{
  int p[ha],nex[maxl],len,f[maxl];
  ll st[maxl];
  void clr()
  {
    len=0; memset(p,255,sizeof(p));
  }
  void push(ll nw,int s)
  {
    int x=nw%ha;
    for (int i=p[x];i!=-1; i=nex[i])
      if (st[i]==nw)
      {
        f[i]=(f[i]+s)%mo;
        return;
      }
    st[++len]=nw; f[len]=s;
    nex[len]=p[x]; p[x]=len;
  }
} h[2];
int n,m,cas,p,t;
int a[15][15],b[15],v[15];

void get(ll st)
{
  b[m+1]=st&63; st>>=6;
  for (int i=m; i>=0; i--)
  {
    b[i]=st&7;
    st>>=3;
  }
}

ll put()
{
  memset(v,255,sizeof(v)); v[0]=0;
  ll st=0; int t=0;
  for (int i=0; i<=m; i++)
  {
    if (v[b[i]]==-1) v[b[i]]=++t;
    b[i]=v[b[i]];
    st<<=3; st|=b[i];
  }
  st<<=6; st|=b[m+1];
  return st;
}

void shift()
{
  for (int i=m;i; i--) b[i]=b[i-1];
  b[0]=0;
}

bool check(int j)
{
  bool ff=1;
  for (int i=0; i<j-1; i++)
    if (b[i]) ff^=1;
  if (!ff) return 0;
  for (int i=j+1; i<=m; i++)
    if (b[i]) ff^=1;
  return ff;
}

void dp(int i,int j)
{
  for (int k=1; k<=h[p^1].len; k++)
  {
    get(h[p^1].st[k]);
    if (!a[i][j])
    {
      b[j]=b[j-1]=0;
      h[p].push(put(),h[p^1].f[k]);
      continue;
    }
    int x=b[j],y=b[j-1];
    if (b[m+1]==t) continue;
    if (x&&y)
    {
      b[j]=b[j-1]=0;
      if (x!=y){
        for (int r=0; r<=m; r++) if (b[r]==x) b[r]=y;
      }
      else {
        if (check(j)) ++b[m+1]; else continue;
      }
      h[p].push(put(),h[p^1].f[k]);
    }
    else if ((x&&!y)||(!x&&y))
    {
      int r=x+y;
      if (a[i][j+1])
      {
        b[j]=r; b[j-1]=0;
        h[p].push(put(),h[p^1].f[k]);
      }
      if (a[i+1][j])
      {
        b[j-1]=r; b[j]=0;
        h[p].push(put(),h[p^1].f[k]);
      }
    }
    else {
      if (a[i][j+1]&&a[i+1][j])
      {
        b[j]=b[j-1]=m+1;
        h[p].push(put(),h[p^1].f[k]);
      }
    }
  }
}

int main()
{
  char s[15]; int cas;
  scanf("%d",&cas);
  while (cas--)
  {
    scanf("%d%d%d",&n,&m,&t);
    memset(a,0,sizeof(a));
    for (int i=1; i<=n; i++)
    {
      scanf("%s",s+1);
      for (int j=1; j<=m; j++)
        if (s[j]=='.') a[i][j]=1;
    }
    h[0].clr();
    h[0].push(0,1);  p=0;
    for (int i=1; i<=n; i++)
    {
      p^=1; h[p].clr();
      for (int k=1; k<=h[p^1].len; k++)
      {
        get(h[p^1].st[k]); shift();
        h[p].push(put(),h[p^1].f[k]);
      }
      for (int j=1; j<=m; j++)
      {
        p^=1; h[p].clr();
        dp(i,j);
      }
    }
    int ans=0;
    for (int i=1; i<=h[p].len; i++)
    {
      get(h[p].st[i]);
      if (b[m+1]==t) ans=(ans+h[p].f[i])%mo;
    }
    printf("%d
",ans);
  }
  return 0;
}

矩阵乘法加速插头dp的几题一直想做，抽时间要补一下