hdu 1946 Pipes

這題基本上和上兩題是一樣的,不過這題是求所有路徑權值最小,轉移就直接帶權轉,每個可行狀態取最小即可。

Run ID Submit Time Judge Status Pro.ID Exe.Time Exe.Memory Code Len. Language Author
6849681 2012-10-01 21:43:55 Accepted 1964 31MS 616K 3026 B C++ xym2010
#include<cstdio>
#include<cstring>
#include<string>
#include<iostream>
#include<algorithm>
#include<cmath>
#define LL long long
using namespace std;
const int maxn=30001;
int n,m;
int gp1[20][20],gp2[20][20];
int mov[20]={0,2,4,6,8,10,12,14,16,18,20};
char s[100];
struct node
{
	int size,head[maxn],next[maxn],sum[maxn];
	LL sta[maxn];
	void clear()
	{
		memset(head,-1,sizeof(head));
		size=0;
	}
	void push(LL st,const int v)
	{
		LL hash=st%maxn;
		for(int i=head[hash];i>=0;i=next[i])
		{
			if(sta[i]==st)
			{
				sum[i]=min(sum[i],v);
				return ;
			}
		}
		sta[size]=st,sum[size]=v;
		next[size]=head[hash],head[hash]=size++;
	}
}dp[2];
inline LL getbit(LL st,int k){return 3&(st>>mov[k]);}
inline LL pybit(LL st,int k){return st<<mov[k];}
inline LL clrbit(LL st,int a,int b){return st&(~(3<<mov[a]))&(~(3<<mov[b]));}
int fl(int st,int k)
{
	int cnt=1;
	for(int i=k+1;i<=m;i++)
	{
		int e=getbit(st,i);
		if(e==2) cnt--;
		else if(e==1) cnt++;
		if(!cnt) return i;
	}
}
int fr(int st,int k)
{
	int cnt=1;
	for(int i=k-1;i>=0;i--)
	{
		int e=getbit(st,i);
		if(e==2) cnt++;
		else if(e==1) cnt--;
		if(!cnt) return i;
	}
}
int DP()
{
	dp[0].clear();dp[0].push(0,0);
	int now=0,pre=1;
	for(int i=1;i<=n;i++)
	{
		pre=now,now^=1;dp[now].clear();
		for(int j=0;j<dp[pre].size;j++)
			dp[now].push(dp[pre].sta[j]<<2,dp[pre].sum[j]);
		for(int j=1;j<=m;j++)
		{
			pre=now,now^=1;dp[now].clear();
			for(int k=0;k<dp[pre].size;k++)
			{
				LL l=getbit(dp[pre].sta[k],j-1);
				LL up=getbit(dp[pre].sta[k],j);
				LL st=clrbit(dp[pre].sta[k],j-1,j);
				if(!l&&!up)
				{
					if(i<n&&j<m)
						dp[now].push(st|pybit(1,j-1)|pybit(2,j),dp[pre].sum[k]+gp1[i][j]+gp2[i][j]);
				}
				else if(!l||!up)
				{
					int e=l==0?up:l;
					if(i<n)
						dp[now].push(st|pybit(e,j-1),dp[pre].sum[k]+gp2[i][j]);
					if(j<m)
						dp[now].push(st|pybit(e,j),dp[pre].sum[k]+gp1[i][j]);
				}
				else if(l==1&&up==1)
					dp[now].push(st^pybit(3,fl(st,j)),dp[pre].sum[k]);
				else if(l==2&&up==2)
					dp[now].push(st^pybit(3,fr(st,j-1)),dp[pre].sum[k]);
				else if(l==2&&up==1)
					dp[now].push(st,dp[pre].sum[k]);
				else if(i==n&&j==m)
					dp[now].push(st,dp[pre].sum[k]);
			}
		}
	}
	int ans=0;
	for(int i=0;i<dp[now].size;i++)
		if(dp[now].sta[i]==0)
			ans+=dp[now].sum[i];
	return ans;
}
int main()
{
	int T;
	scanf("%d",&T);
	while(T--)
	{
		scanf("%d%d",&n,&m);
		getchar();
		gets(s);
		for(int i=1;i<n;i++)
		{
			gets(s);
			for(int j=1,k=2;j<m;j++,k+=2)
				gp1[i][j]=s[k]-'0';
			gets(s);
			for(int j=1,k=1;j<=m;j++,k+=2)
				gp2[i][j]=s[k]-'0';
		}
		gets(s);
		for(int j=1,k=2;j<m;j++,k+=2)
			gp1[n][j]=s[k]-'0';
		gets(s);
		cout<<DP()<<endl;
	}
	return 0;
}


發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章