插头dp初探:BZOJ1210&&Uva1519

这是一个坑……

几天前的考试一道插头dp我硬生生裸状压写了2个多小时然后考试挂了QAQ

想学插头dp的话还是陈丹琪的论文比较好

但是要明白最重要的一点：我们每个格子转移完毕后他的插头状态是要留给后面格子的

没什么好说的，注意细节

学完插头dp分类讨论什么鬼的再也不怕了QAQ

Uva1519(此题可在COGS上做)

链接：[URAL 1519]一级方程式赛车

这是CDQ的论问题，和HNOI2004邮递员一样是个入门题，这里只放前者

#include<queue>
#include<cstdio>
#include<vector>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#define LL long long
using namespace std;
const int maxs=20;
const int mod=9997;
const int maxl=3000000;
int code[maxs],x,y;
int ch[maxs],n,m;
bool a[15][15];
char s[15];

struct Hashmap{
	int next[maxl],h[mod],tot;
	LL val[maxl],As[maxl];
	void push(LL x,LL v){
		int now=x%mod; if (now<0) now+=mod;
		for (int i=h[now];i;i=next[i])
		    if (As[i]==x){val[i]+=v;return;}
		++tot; next[tot]=h[now]; h[now]=tot;
		As[tot]=x; val[tot]=v; return ;
	}
	void init(){
		memset(h,0,sizeof(h)); tot=0;
	}
}dp[2];

LL encode(){
	LL st=0; int cnt=0;
	memset(ch,-1,sizeof(ch)); ch[0]=0;
	for (int i=0;i<=m;++i){
		if (ch[code[i]]==-1) ch[code[i]]=++cnt;
		code[i]=ch[code[i]];
		st<<=3; st|=code[i];
	}return st;
}

void decode (LL st){
	for (int i=m;i>=0;--i) code[i]=st&7,st>>=3;
}

void shit(){
	for (int i=m;i>=1;--i)
	    code[i]=code[i-1];
	code[0]=0;
}

void dpblock(int i,int j,int cur){
	for(int k=1;k<=dp[cur].tot;++k){
		decode(dp[cur].As[k]);
		if(j==1){
			if(code[m])continue;
			shit();
		}
		int up=code[j],left=code[j-1];
		if (!up&&!left)
			dp[cur^1].push(encode(),dp[cur].val[k]);
	}
}

void dpright(int i,int j,int cur){
	for (int k=1;k<=dp[cur].tot;++k){
		decode(dp[cur].As[k]);
		if (j==1){
			if (code[m]) continue;
			shit();}
		int up=code[j],left=code[j-1];
		if (!up&&!left){
			if (!a[i+1][j]||!a[i][j+1]) continue;
			code[j-1]=code[j]=13;
			dp[cur^1].push(encode(),dp[cur].val[k]);
		}else if ((!up&&left)||(up&&!left)){
			int tmp=left?left:up;
			if (a[i+1][j]){
				code[j-1]=tmp; code[j]=0;
				dp[cur^1].push(encode(),dp[cur].val[k]);
			}if (a[i][j+1]){
				code[j]=tmp; code[j-1]=0;
				dp[cur^1].push(encode(),dp[cur].val[k]);
			}
		}else{
			if (up!=left){
				code[j]=code[j-1]=0;
				for (int u=0;u<=m;++u) if (code[u]==left) code[u]=up;
				dp[cur^1].push(encode(),dp[cur].val[k]);
			}else if (up==left){
				if (i==x&&j==y){
					code[j]=code[j-1]=0;
					dp[cur^1].push(encode(),dp[cur].val[k]);
				}
			}
		}
	}
}

void solve_DP(){
	int cur=0; dp[cur].push(0,1);
	for (int i=1;i<=n;++i){
		for (int j=1;j<=m;++j){
			dp[cur^1].init();
		    if (!a[i][j]) dpblock(i,j,cur);
		    else dpright(i,j,cur); cur^=1;
		}
	} printf("%lld",dp[cur].val[1]);
}

int main(){
	scanf("%d%d",&n,&m);
	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,x=i,y=j;
	        else a[i][j]=0;
	}
	solve_DP();
}

Bzoj1187 [HNOI2007]神奇游乐园

没有了坏格子，但是我们可以在任意一个地方结束，所以在每个地方只要上插头左插头都存在就可以判断并更新答案

判断的依据是当前轮廓线上没有多余插头，否则不是标准回路

同时因为我们并没有要求每个格子都经过，所以我们可以空一个格子不选直接进行下一个，这种情况当且仅当左插头和上插头都没有

还有就是要一个回路QAQ，所以只要形成了一个符合要求的标准回路他是不会加入下一个格子的状态集合进行转移的

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#define LL long long
using namespace std;
const int mod=1331;
const int maxl=1000000+10;
const LL inf=1e15;
LL ans=-inf;
int code[15],ch[15],n,m;
int a[110][110];

struct Hashmap{
    int next[maxl],h[mod],tot;
    LL val[maxl],As[maxl];
    void push(LL x,LL v){
        int now=x%mod; if (now<0) now+=mod;
        for (int i=h[now];i;i=next[i])
            if (As[i]==x){val[i]=max(val[i],v);return;}
        ++tot; next[tot]=h[now]; h[now]=tot;
        As[tot]=x; val[tot]=v; return ;
    }
    void init(){
        memset(h,0,sizeof(h)); tot=0;
    }
}dp[2];

LL encode(){
    LL st=0; int cnt=0;
    memset(ch,-1,sizeof(ch)); ch[0]=0;
    for (int i=0;i<=m;++i){
        if (ch[code[i]]==-1) ch[code[i]]=++cnt;
        code[i]=ch[code[i]];
        st<<=3; st|=code[i];
    }return st;
}

void decode (LL st){
    for (int i=m;i>=0;--i) code[i]=st&7,st>>=3;
}

void shit(){
    for (int i=m;i>=1;--i)
        code[i]=code[i-1];
    code[0]=0;
}

void dpright(int i,int j,int cur){
    for (int k=1;k<=dp[cur].tot;++k){
        decode(dp[cur].As[k]);
        if (j==1){
            if (code[m]) continue;
            shit();}
        int up=code[j],left=code[j-1];
        if (!up&&!left){
            code[j]=code[j-1]=0;
            dp[cur^1].push(encode(),dp[cur].val[k]);
            if (i+1>n||j+1>m) continue;
            code[j]=code[j-1]=11;
            dp[cur^1].push(encode(),dp[cur].val[k]+a[i][j]);
        }else if ((!up&&left)||(up&&!left)){
            int tmp=left?left:up;
            if (i+1<=n){
                code[j-1]=tmp; code[j]=0;
                dp[cur^1].push(encode(),dp[cur].val[k]+a[i][j]);
            }if (j+1<=m){
                code[j]=tmp; code[j-1]=0;
                dp[cur^1].push(encode(),dp[cur].val[k]+a[i][j]);
            }
        }else{
            if (up!=left){
                code[j]=code[j-1]=0;
                for (int u=0;u<=m;++u) if (code[u]==left) code[u]=up;
                dp[cur^1].push(encode(),dp[cur].val[k]+a[i][j]);
            }else if (up==left){//!!!!!!!!!!!!!!!!!!!!!!!!!!!
                code[j]=code[j-1]=0;
                if (!encode()) ans=max(ans,dp[cur].val[k]+a[i][j]);
            }
        }
    }
}

void solve_DP(){
    int cur=0; dp[cur].push(0,0);
    for (int i=1;i<=n;++i){
        for (int j=1;j<=m;++j){
            dp[cur^1].init();
            dpright(i,j,cur); cur^=1;
        }
    }printf("%lld",ans);
}

int main(){
    scanf("%d%d",&n,&m);
    for (int i=1;i<=n;++i)
        for (int j=1;j<=m;++j)
            scanf("%d",&a[i][j]);
    solve_DP();
}