[BZOJ2095]Bridges

最大值最小,是二分

转化为判定问题:给定一个混合图,问是否存在欧拉回路

首先,有向图存在欧拉回路的充要条件是每个点的入度等于出度,现在我们有一个混合图,我们要做的就是给其中的无向边定向,使得它变成有向图之后存在欧拉回路

记点$x$的入度为$in_x$,出度为$out_x$,我们的目标是使得每个点$x$满足$in_x-out_x=0$

先随便给每条无向边定向,这样不一定满足要求,所以我们必须让某些边反向,反向$a\rightarrow b$会让$in_a-out_a$增加$2$,让$in_b-out_b$减少$2$,所以如果存在$in_x-out_x\equiv1(\text{mod }2)$那么无解

为了决定是否反向某些无向边,我们这样建图:

对于$in_x\lt out_x$的$x$,连边$x\rightarrow T$,容量为$\dfrac{out_x-in_x}2$

对于$in_x\gt out_x$的$x$,连边$S\rightarrow x$,容量为$\dfrac{in_x-out_x}2$

对于一条无向边,如果一开始硬点它的方向为$x\rightarrow y$,那么连边$y\rightarrow x$,权值为$1$

这样建图跑最大流,每流过一条原图中的无向边就相当于将它反向(这条无向边的两个端点的流量之和就是$\left|in_x-out_x\right|$的改变量),跑最大流是因为我们想要尽可能地缩小$in_x$和$out_x$的差距,如果满流,自然就存在欧拉回路了

#include<stdio.h>
#include<string.h>
const int inf=2147483647;
int abs(int x){return x>0?x:-x;}
int min(int a,int b){return a<b?a:b;}
int h[1010],cur[1010],nex[10010],to[10010],cap[10010],dis[1010],q[10010],M,S,T;
void add(int a,int b,int c){
	M++;
	to[M]=b;
	cap[M]=c;
	nex[M]=h[a];
	h[a]=M;
	M++;
	to[M]=a;
	cap[M]=0;
	nex[M]=h[b];
	h[b]=M;
}
bool bfs(){
	int head,tail,x,i;
	memset(dis,-1,sizeof(dis));
	head=tail=1;
	q[1]=S;
	dis[S]=0;
	while(head<=tail){
		x=q[head];
		head++;
		for(i=h[x];i;i=nex[i]){
			if(cap[i]&&dis[to[i]]==-1){
				dis[to[i]]=dis[x]+1;
				if(to[i]==T)return 1;
				tail++;
				q[tail]=to[i];
			}
		}
	}
	return 0;
}
int dfs(int x,int flow){
	if(x==T)return flow;
	int i,f;
	for(i=cur[x];i;i=nex[i]){
		if(cap[i]&&dis[to[i]]==dis[x]+1){
			f=dfs(to[i],min(flow,cap[i]));
			if(f){
				cap[i]-=f;
				cap[i^1]+=f;
				if(cap[i])cur[x]=i;
				return f;
			}
		}
	}
	dis[x]=-1;
	return 0;
}
int dicnic(){
	int ans=0,tmp;
	while(bfs()){
		memcpy(cur,h,sizeof(h));
		while(tmp=dfs(S,inf))ans+=tmp;
	}
	return ans;
}
int n,m,in[1010],ou[1010];
struct edge{
	int x,y,a,b;
}e[2010];
bool check(int lim){
	int i,s;
	memset(h,0,sizeof(h));
	memset(in,0,sizeof(in));
	memset(ou,0,sizeof(ou));
	M=1;
	for(i=1;i<=m;i++){
		if(e[i].a<=lim&&e[i].b<=lim){
			add(e[i].y,e[i].x,1);
			ou[e[i].x]++;
			in[e[i].y]++;
		}else if(e[i].a<=lim){
			ou[e[i].x]++;
			in[e[i].y]++;
		}else if(e[i].b<=lim){
			ou[e[i].y]++;
			in[e[i].x]++;
		}else
			return 0;
	}
	s=0;
	for(i=1;i<=n;i++){
		if(abs(in[i]-ou[i])&1)return 0;
		if(in[i]<ou[i])add(i,T,(ou[i]-in[i])>>1);
		if(in[i]>ou[i]){
			add(S,i,(in[i]-ou[i])>>1);
			s+=(in[i]-ou[i])>>1;
		}
	}
	return s==dicnic();
}
int main(){
	int i,l,r,mid,ans;
	scanf("%d%d",&n,&m);
	for(i=1;i<=m;i++)scanf("%d%d%d%d",&e[i].x,&e[i].y,&e[i].a,&e[i].b);
	S=n+1;
	T=n+2;
	ans=-1;
	l=0;
	r=1000;
	while(l<=r){
		mid=(l+r)>>1;
		if(check(mid)){
			ans=mid;
			r=mid-1;
		}else
			l=mid+1;
	}
	if(ans==-1)
		puts("NIE");
	else
		printf("%d",ans);
}

猜你喜欢

转载自www.cnblogs.com/jefflyy/p/8905848.html