2019.01.23 hdu1964 Pipes(轮廓线dp)

时间:2022-03-24 04:03:27

传送门

题意简述:给一个没有障碍的网格图,任意两个格子连通需要花费一定代价,现在求一条覆盖所有格子的哈密顿回路的总权值的最小值。


思路:

这道题一毛一样,除了把求和变成求最小值以外。

代码:

#include<bits/stdc++.h>
#define ri register int
using namespace std;
inline int read(){
	int ans=0;
	char ch=getchar();
	while(!isdigit(ch))ch=getchar();
	while(isdigit(ch))ans=(ans<<3)+(ans<<1)+(ch^48),ch=getchar();
	return ans;
}
typedef long long ll;
const int mod=1e6+7;
int mp[15][15],a[15][15],b[15][15],n,m,ans;
bool cur;
char s[30][30];
struct Statement{
	int mn[mod],sta[mod],idx[mod],tot;
	inline void clear(){tot=0,memset(idx,-1,sizeof(idx));}
	inline void insert(int stat,int mnn){
		int pos=stat%mod;
		if(!pos)++pos;
		while(~idx[pos]&&sta[idx[pos]]!=stat)pos=pos==mod-1?1:pos+1;
		if(~idx[pos])mn[idx[pos]]=min(mn[idx[pos]],mnn);
		else mn[idx[pos]=++tot]=mnn,sta[tot]=stat;
	}
}f[2];
inline int getbit(int x,int p){return (x>>((p-1)<<1))&3;}
inline void update(int&x,int p,int v){x^=((getbit(x,p)^v)<<((p-1)<<1));}
inline void solve(){
	ans=1e9,f[cur=0].clear(),f[cur].insert(0,0);
	for(ri i=1;i<=n;++i){
		for(ri j=1;j<=m;++j){
			f[(cur^=1)].clear();
			for(ri tt=1;tt<=f[cur^1].tot;++tt){
				int stat=f[cur^1].sta[tt],p=getbit(stat,j),q=getbit(stat,j+1),mn=f[cur^1].mn[tt];
				if(mn>=ans)continue;
				if(!(p+q)){
					if(mp[i+1][j]&&mp[i][j+1]){
						update(stat,j,1),update(stat,j+1,2);
						f[cur].insert(stat,mn+a[i][j]+b[i][j]);
					}
					continue;
				}
				if(!p){
					if(mp[i][j+1])f[cur].insert(stat,mn+a[i][j]);
					if(mp[i+1][j])update(stat,j,q),update(stat,j+1,0),f[cur].insert(stat,mn+b[i][j]);
					continue;
				}
				if(!q){
					if(mp[i+1][j])f[cur].insert(stat,mn+b[i][j]);
					if(mp[i][j+1])update(stat,j,0),update(stat,j+1,p),f[cur].insert(stat,mn+a[i][j]);
					continue;
				}
				if(p==1&&q==2){
					if(i==n&&j==m)ans=mn;
					continue;
				}
				update(stat,j,0),update(stat,j+1,0);
				if(p==2&&q==1){f[cur].insert(stat,mn);continue;}
				if(p==1){
					for(ri cnt=1,k=j+1,bit;k<=m;++k){
						bit=getbit(stat,k);
						if(bit==1)++cnt;
						if(bit==2)--cnt;
						if(!cnt){update(stat,k,1),f[cur].insert(stat,mn);break;}
					}
					continue;
				}
				for(ri cnt=-1,k=j-1,bit;k>1;--k){
					bit=getbit(stat,k);
					if(bit==1)++cnt;
					if(bit==2)--cnt;
					if(!cnt){update(stat,k,2),f[cur].insert(stat,mn);break;}
				}
			}
		}
		for(ri j=1;j<=f[cur].tot;++j)f[cur].sta[j]<<=2;
	}
}
int main(){
	for(ri tt=read();tt;--tt){
		n=read(),m=read(),memset(mp,0,sizeof(mp)),memset(a,0,sizeof(a)),memset(b,0,sizeof(b));
		for(ri i=1;i<=n*2+1;++i)gets(s[i]+1);
		for(ri i=1;i<=n;++i){
			for(ri j=1;j<=m;++j){
				int x=i*2,y=j*2;
				a[i][j]=s[x][y+1]=='#'?0:s[x][y+1]-'0',b[i][j]=s[x+1][y]=='#'?0:s[x+1][y]-'0',mp[i][j]=1;
			}
		}
		solve(),cout<<ans<<'\n';
	}
	return 0;
}