1287. 躲雨

Description

　　FJ的奶牛真的很怕雨，所以一下雨他们就急着要找地方躲雨。为了做好防雨准备，他们装了警报器，一有雨就拉响警报器，他们利用天气预报来预测雨的到来，但有时候不准导致拉错警报，为了减少拉错警报的可能，他们想在保证奶牛有足够时间躲起来的前提下尽可能晚地拉响警报。
　　农场有F(1<=F<=200)个奶牛吃草的区域，有P(1<=P<=1500)条路连接一些区域，路是双向的。有些区域有遮雨棚，每个遮雨棚有自己的容量限制。
　　现在要你计算所有牛都能找到地方躲雨最少需要多少时间。

Input

　　第1行：两个空格隔开的整数F和P
　　第2到F+1行：每行两个空格隔开的整数X,Y，X表示当前这个区域的牛的数量，Y表示该区域的遮雨棚最多只能容纳Y头牛。
　　第F+2到F+P+1行：3个空格隔开的整数a,b,c,表示a和b之间有一条长度为c的路直接相连。

Output

　　输出一个整数，表示所有牛都能在这个时间内到达遮雨棚。

Sample Input

3 4
7 2
0 4
2 6
1 2 40
3 2 70
2 3 90
1 3 120

Sample Output

110

Data Constraint

Hint

【样例说明】
　　1号区域里保留2头牛，4头牛去2号区域，还有1个去3号区域，3号区域的2头牛保持不动，可能还有其他方案，但110是最少的时间。

Solution

首先预处理出两两之间的最短路，然后二分答案。对于那些距离小于mid的边，我们才用它。然后我们将一个点拆成两个点，每个点向S连一条容量为当前有的牛的个数，再在两个相同的点于对应的附加点之间连一条INF的边，再将距离小于mid的边（x,y）,用点x连向点y的附加点，最后将所有附加的点连线汇点T，这样我们做一遍网络流之后看一下汇点是否满流即可（所有的牛都能走到汇点）。

Code

#include<cstdio>
#include<cstring>
#include<algorithm>
#define inf 1e18
#define N 1000
#define M 81200
#define ll long long
using namespace std;
int n,m,x,y,S,T,a[N],b[N],d[N*100],q[N*N],h[N],nx[M*2],t[M*2],tot;
ll l,r,mid,ans=-1,dis[N][N],z,s[M*2],f[N],sum,sa=0;
void add(int x,int y,ll z){
	t[++tot]=y;nx[tot]=h[x];s[tot]=z;h[x]=tot;
	t[++tot]=x;nx[tot]=h[y];s[tot]=0;h[y]=tot;
}
int bfs(){
	memset(d,0,sizeof(d));
	d[S]=1;
	q[1]=S;int i=1,j=1;
	while(i<=j){
		int x=q[i++];
		for(int k=h[x];k;k=nx[k]) if(s[k]&&!d[t[k]]){
			d[t[k]]=d[x]+1;
			q[++j]=t[k];
			if(t[k]==T) return 1;
		}
	}
	return 0;
}
ll dinic(int x,ll flow){
	if(x==T) return flow;
	ll rest=flow,p;
	for(int k=h[x];k&&rest;k=nx[k]) if(s[k]&&d[t[k]]==d[x]+1){
		p=dinic(t[k],min(rest,s[k]));
		if(!p) d[t[k]]=0;
		s[k]-=p;s[k^1]+=p;
		rest-=p;
		if(!flow) break;
	}
	return flow-rest;
}
int check(ll x){
	memset(h,0,sizeof(h));
	tot=1;
	for(int i=1;i<=n;i++){
		if(a[i]) add(S,i,a[i]);if(b[i]) add(i+n,T,b[i]);
		for(int j=1;j<=n;j++) if(dis[i][j]<=x) add(i,j+n,inf);
	}
	sum=0;
	while(bfs()) sum+=dinic(S,inf);
	return  sum==sa;
}
int main(){
	scanf("%d%d",&n,&m);S=0;T=(n<<1)+1;
	memset(dis,0x3f,sizeof(dis));
	for(int i=1;i<=n;i++){
		dis[i][i]=0;
		scanf("%d%d",&a[i],&b[i]);
		sa+=a[i];
	}
	for(int i=1;i<=m;i++){
		scanf("%d%d%lld",&x,&y,&z);
		r+=z;
		dis[x][y]=dis[y][x]=min(dis[x][y],z);
	}
	for(int k=1;k<=n;k++)
		for(int i=1;i<=n;i++)
			for(int j=1;j<=n;j++)
				dis[i][j]=min(dis[i][j],dis[i][k]+dis[k][j]);
	while(l<=r){
		mid=(l+r)>>1;
		if(check(mid)){ans=mid;r=mid-1;}
		else l=mid+1;
	}
	printf("%lld\n",ans);
	return 0;
}

---

---

作者：zsjzliziyang 
QQ:1634151125 
转载及修改请注明 
本文地址:https://blog.csdn.net/zsjzliziyang/article/details/88764304