【POI 2010】Bridges

最新推荐文章于 2019-07-26 22:01:00 发布

forever_dreams

最新推荐文章于 2019-07-26 22:01:00 发布

阅读量214

点赞数

分类专栏： --------基础算法-------- # 二分 --------图论-------- # 欧拉回路 # 网络流文章标签：混合图欧拉回路二分网络流

本文链接：https://blog.csdn.net/forever_dreams/article/details/83540308

版权

--------基础算法-------- 同时被 3 个专栏收录

77 篇文章 4 订阅

订阅专栏

--------图论--------

75 篇文章 1 订阅

订阅专栏

网络流

27 篇文章 0 订阅

订阅专栏

【题目】

题目描述：

YYD 为了减肥，他来到了瘦海，这是一个巨大的海，海中有 $n$ 个小岛，小岛之间有 $m$ 座桥连接，两个小岛之间不会有两座桥，并且从一个小岛可以到另外任意一个小岛。现在 YYD 想骑单车从小岛 $1$ 出发，骑过每一座桥，到达每一个小岛，然后回到小岛 $1$ 。霸中同学为了让 YYD 减肥成功，召唤了大风，由于是海上，风变得十分大，经过每一座桥都有不可避免的风阻碍 YYD，YYD 十分 ddt，于是用泡芙贿赂了你，希望你能帮他找出一条承受的最大风力最小的路线。

输入格式：

输入：第一行为两个用空格隔开的整数 $n$ （ $2 \leq n \leq 1000$ ）， $m$ （ $1 \leq m \leq 2000$ ），接下来读入 $m$ 行由空格隔开的 $4$ 个整数 $a$ ， $b$ （ $1 \leq a, b \leq n$ ， $a\ne b$ ）， $c$ ， $d$ （ $1 \leq c, d \leq 1000$ ），表示第 $i + 1$ 行第 $i$ 座桥连接小岛 $a$ 和 $b$ ，从 $a$ 到 $b$ 承受的风力为 $c$ ，从 $b$ 到 $a$ 承受的风力为 $d$ 。

输出格式：

输出：如果无法完成减肥计划，则输出 NIE，否则第一行输出承受风力的最大值（要使它最小)

样例数据：

输入
4 4
1 2 2 4
2 3 3 4
3 4 4 4
4 1 5 4

输出
4

提示：
在这里插入图片描述
注意：通过桥为欧拉回路

【分析】

看到最后的"使承受风力的最大值最小"，就知道多半要二分

二分之后，对于每个中间值 $m i d$ ，小于 $m i d$ 的边就丢掉，然后相当于在剩下的图中判断有没有欧拉回路

但是这个图不是无向图，也不是有向图，是一个混合图，该怎么办呢？

先对每条无向边任意定向，统计出入度和出读，如果有一个点的入度和出度的奇偶性不同，那么就肯定无解（因为最后的目标是入度和出度相等，而改变无向边方向的话，会让它们的差值变化 $2$ ，无法使它们相等）

然后现在就想如何将无向边反向，使得最后的如何和出度相等，这是我们的网络流就闪亮登场了

建立一个源点 $s$ ，一个汇点 $t$ ，假设入度 $-$ 出度 $= x$ ，若 $x > 0$ ，就向 $t$ 连一条 $\frac{x}{2}$ 的边（代表这个点的入度要减 $\frac x2$ 才能和出度相等，即要使 $\frac x2$ 条边反向），否则就从源点连一条 $-\frac{x}2$ 的边（同理），对于原来定向的无向边 $(a, b)$ ，建立一条从 $b$ 到 $a$ ，容量为 $1$ 的边

其实说的简单一点，在新建的图中，流一条流就代表将一条无向边反向

那么跑一个最大流，如果最后从 $s$ 出去的边都能流满，则说明可以找到一种方案使得原图是欧拉回路

然后就华丽丽（并不）的结束啦

【代码】

#include<queue>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 1005
#define M 200005
#define inf (1ll<<31ll)-1
using namespace std;
int n,m,s,t,num,tot;
int a[M],b[M],c[M],d[M];
int f[N],du[N],dis[N],first[N],v[M],w[M],nxt[M];
void add(int x,int y,int z)
{
	num++;
	nxt[num]=first[x];
	first[x]=num;
	v[num]=y;
	w[num]=z;
}
void create(int x)
{
	int i;tot=0,num=1;
	memset(du,0,sizeof(du));
	memset(first,0,sizeof(first));
	for(i=1;i<=m;++i)
	{
		if(c[i]<=x)  du[a[i]]--,du[b[i]]++;
		if(d[i]<=x)  add(b[i],a[i],1),add(a[i],b[i],0);
	}
	for(i=1;i<=n;++i)
	{
		if(du[i]<0)  add(i,t,-du[i]/2),add(t,i,0);
		else  tot+=du[i]/2,add(s,i,du[i]/2),add(i,s,0);
	}
}
bool bfs()
{
	int x,y,i;
	memcpy(f,first,sizeof(f));
	memset(dis,-1,sizeof(dis));
	queue<int>q;q.push(s);dis[s]=0;
	while(!q.empty())
	{
		x=q.front();q.pop();
		for(i=first[x];i;i=nxt[i])
		{
			y=v[i];
			if(w[i]&&dis[y]==-1)
			{
				dis[y]=dis[x]+1;
				if(y==t)
				  return true;
				q.push(y);
			}
		}
	}
	return false;
}
int dinic(int now,int flow)
{
	if(now==t)  return flow;
	int x,delta,ans=0;
	for(int &i=f[now];i;i=nxt[i])
	{
		x=v[i];
		if(w[i]&&dis[x]==dis[now]+1)
		{
			delta=dinic(x,min(flow,w[i]));
			w[i]-=delta,w[i^1]+=delta;
			flow-=delta,ans+=delta;
			if(!flow)  return ans;
		}
	}
	return ans;
}
bool check(int mid)
{
	int i,ans=0;
	create(mid);
	for(i=1;i<=n;++i)
	  if(du[i]&1)
	    return false;
	while(bfs())
	  ans+=dinic(s,inf);
	return ans==tot;
}
int main()
{
	int i;
	scanf("%d%d",&n,&m);
	s=0,t=n+1;
	int Min=1001,Max=0;
	for(i=1;i<=m;++i)
	{
		scanf("%d%d%d%d",&a[i],&b[i],&c[i],&d[i]);
		if(c[i]>d[i])  swap(a[i],b[i]),swap(c[i],d[i]);
		Min=min(Min,c[i]),Max=max(Max,d[i]);
	}
	int l=Min,r=Max;
	while(l<r)
	{
		int mid=(l+r)>>1;
		if(check(mid))  r=mid;
		else  l=mid+1;
	}
	if(!check(l))  printf("NIE");
	else  printf("%d",l);
	return 0;
}