Description

约翰的牛们非常害怕淋雨,那会使他们瑟瑟发抖.他们打算安装一个下雨报警器,并且安排了一个撤退计划.他们需要计算最少的让所有牛进入雨棚的时间.    牛们在农场的F(1≤F≤200)个田地上吃草.有P(1≤P≤1500)条双向路连接着这些田地.路很宽,无限量的牛可以通过.田地上有雨棚,雨棚有一定的容量,牛们可以瞬间从这块田地进入这块田地上的雨棚    请计算最少的时间,让每只牛都进入雨棚.

Input

第1行:两个整数F和P;

第2到F+1行:第i+l行有两个整数描述第i个田地,第一个表示田地上的牛数,第二个表示田地上的雨棚容量.两个整数都在0和1000之间.

第F+2到F+P+I行:每行三个整数描述一条路,分别是起点终点,及通过这条路所需的时间(在1和10^9之间).

Output

一个整数,表示最少的时间.如果无法使牛们全部进入雨棚,输出-1.

题解:

首先floyd处理出点到点之间的最短时间,然后二分一个时间,网络流判断是否全部进入雨棚。

代码:

#include <cstdio>
#include <cstring>
#include <algorithm>
#include<queue>
// by zrt
// problem:
// 无论你在什么时候开始,重要的是开始以后就不要停止。
using namespace std ;
typedef long long LL ;
const double eps(1e-10) ;
const int inf(0x3f3f3f3f) ;
int f,p;
LL INF=0x3f3f3f3f3f3f3f3f;
int ni[205],fi[205];
LL map[205][205];
int S,T;
int H[404],X[340000],P[340000],tot;
int flow[340000];
int d[406];
queue<int> q;
bool bfs(){
memset(d,0,sizeof d);
d[S]=1;q.push(S);
while(!q.empty()){
int x=q.front();q.pop();
for(int i=H[x];i;i=X[i]){
if(flow[i]>0&&!d[P[i]]){
d[P[i]]=d[x]+1;
q.push(P[i]);
}
}
}
return d[T];
}
int dfs(int x,int a){
int f=a,tmp;
if(x==T||a==0) return a;
for(int i=H[x];i;i=X[i]){
if(flow[i]>0&&d[P[i]]==d[x]+1){
tmp=dfs(P[i],min(a,flow[i]));
a-=tmp;
flow[i]-=tmp;
flow[i^1]+=tmp;
if(!a) break;
}
}
if(f==a) d[x]=-1;
// printf("%d\n",f-a);
return f-a;
}
int Dinic(){
int f=0;
while(bfs()) f+=dfs(S,inf);
return f;
}
inline void add(int x,int y,int z){
P[++tot]=y;X[tot]=H[x];H[x]=tot;flow[tot]=z;
}
int sum=0;
bool judge(LL lim){
tot=1;memset(H,0,sizeof H);
for(int i=1;i<=f;i++){
add(S,i<<1,ni[i]);
add(i<<1,S,0);
add(i<<1,i<<1|1,inf);
add(i<<1|1,i<<1,0);
add(i<<1|1,T,fi[i]);
add(T,i<<1|1,0);
}
for(int i=1;i<=f;i++){
for(int j=1;j<=f;j++){
if(i!=j) {
if(map[i][j]<=lim){
add(i<<1,j<<1|1,inf);
add(j<<1|1,i<<1,0);
}
}
}
}
if(Dinic()==sum) return true;
else return false;
}
int main(){
#ifdef LOCAL
freopen("in.txt","r",stdin) ;
freopen("out.txt","w",stdout) ;
#endif
scanf("%d%d",&f,&p);
S=0,T=1; for(int i=1;i<=f;i++) scanf("%d%d",&ni[i],&fi[i]),sum+=ni[i];
memset(map,0x3f,sizeof map);
for(int i=0,x,y,z;i<p;i++){
scanf("%d%d%d",&x,&y,&z);
map[x][y]=map[y][x]=min(map[x][y],(LL)z);
// map[x][y]=z;
}
for(int k=1;k<=f;k++){
for(int i=1;i<=f;i++){
if(i==k) continue;
for(int j=1;j<=f;j++){
if(i==j||k==j) continue;
map[i][j]=min(map[i][j],map[i][k]+map[k][j]);
}
}
}
if(!judge(INF-1)) {
puts("-1");
return 0;
}
LL l=-1,r=INF-1,m;
while(r-l>1){
m=(l+r)/2;
if(judge(m)){
r=m;
}else l=m;
}
printf("%lld\n",r);
return 0 ;
}
04-18 21:09