hdu 1946 Pipes

本文介绍了一种求解所有路径中权值最小的最短路径算法,并通过具体代码实现展示了算法的工作原理。该算法适用于动态规划问题,通过状态转移方程进行求解,确保每个可行状态取到最小值。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

这题基本上和上两题是一样的,不过这题是求所有路径权值最小,转移就直接带权转,每个可行状态取最小即可。

Run IDSubmit TimeJudge StatusPro.IDExe.TimeExe.MemoryCode Len.LanguageAuthor
68496812012-10-01 21:43:55Accepted196431MS616K3026 BC++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;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值