poj Secret Milking Machine

33 篇文章 0 订阅
7 篇文章 0 订阅

题意:有n个点,之间有p条路径(无向,可重复),要求求出要同时存在T条路径(每条边都不重复的条件下),路径中边的最短距离。

方法:二分搜索+最大流,二分枚举所有的距离,求出最短的,再构造网络(存在一条边两点之间的通路小于枚举值,两点间的流量自加1),最后的最大流如果是大于

#include<iostream>
#include<vector>
using namespace std;
const int inf = 0x7fffffff;
const int MAXN = 210;
struct edge
{
	int v,u,w;
}E[40001];
const int MAXM = MAXN*MAXN;
int map[MAXN][MAXN];
int dis[MAXN];
int queue[MAXM];
int bef[MAXN];
int num[MAXN];
void set_dis(int end)
{
	int top(0),base(0);
	memset(dis,-1,sizeof(dis));
	memset(num,0,sizeof(num));
	queue[top++] = end;
	dis[end]=0;num[0] = 1;
	while(top!=base)
	{
		int pre = queue[base++];
		for(int i=end-1;i!=0;i--)
		{
			if(dis[i] != -1 || map[i][pre] == 0)
				continue;
			dis[i] = dis[pre]+1;
			queue[top++] = i;
			num[dis[i]]++;
		}
	}
}
int cmp(int x,int y)
{
	if(x<y)
		return x;
	return y;
}
int remark(int n,int a)
{
	int tmp=inf;
    for(int i=1;i!=n+1;i++)
		if(map[a][i]!=0&&dis[i]>=0)
			tmp=cmp(tmp,dis[i]+1);
    if(tmp==inf)
		tmp=n;
    return tmp;
}
int SAP(int n,int s)
{
	int j,a = s;int flow(0);
	memset(bef,-1,sizeof(bef));
	while(dis[s]<n)
	{
		for(j=1;j!=n+1;j++)
		{
			if(map[a][j]>0&&dis[a]-1==dis[j])
			{
				break;
			}
		}
		if(j<=n)
		{
			bef[j] = a;a = j;
			if(a==n)
			{
				int i = n,max(inf);
				while(bef[i]!=-1)
				{
					max = cmp(max,map[bef[i]][i]);
					i = bef[i];
				}
				i = n;
				while(bef[i]!=-1)
				{
					map[bef[i]][i]-=max;
					map[i][bef[i]]+=max;
					i = bef[i];
				}
				flow+=max;a=s;
			}
		}
		else
		{
			int x = remark(n,a);
			num[x]++;num[dis[a]]--;
            if(num[dis[a]]==0)
				return flow;//间隙优化 gap
            dis[a]=x;
            if(a!=s)
				a=bef[a]; 
		}
	}
	return flow;	
}
void build(int n,int max,int p)
{
	for(int i=1;i!=p+1;i++)
	{
		if(E[i].w<=max)
		{
			map[E[i].u][E[i].v]++;
			map[E[i].v][E[i].u]++;
		}
	}
}
bool judge(int p,int n,int T,int max)
{
	memset(map,0,sizeof(map));
	build(n,max,p);
	set_dis(n);
	int f = SAP(n,1);
	//cout<<f<<endl;
	if(f >= T)
		return true;
	return false;
}
int main()
{
	int	N,P,T;
	while(cin>>N>>P>>T)
	{
		memset(E,0,sizeof(E));
		int max(-1);
		for(int i=1;i!=P+1;i++)
		{
			cin>>E[i].u>>E[i].v>>E[i].w;
			if(E[i].w>max)
				max = E[i].w+1;
		}
		int min = 0;int mid;
		while(min!=max)
		{
			mid = (max+min)>>1;
			if(judge(P,N,T,mid))
			{
				max = mid;
			}
			else
			{
				min = mid+1;
			}
		}
		cout<<max<<endl;
	}
	return 0;
}

或等于T的话就成立。

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值