题意:有n个牛棚,现在每个牛棚都有ai牛,下雨的时候每一个牛棚可以放bi只牛才不会淋湿。把牛放入另一个牛棚需要一些时间。。问最短要多少时间。能把牛放到其他的牛盆里,牛才不会被淋湿。。
拆点,2*a个牛棚 , 源点连a个牛棚,流量为ai , (s,i,ai) ,每个牛棚连汇点t , (i+a,t,bi) ; 弗洛伊德算出两两牛棚之间的最小花费时间,枚举花费时间;
对于枚举到的花费时间mid ,如果两个牛棚间花费时间小于等于mid ;两个牛棚就连边 , (i,j+a,inf) ; 跑一边最大流,最大流=所有的牛 ,即有解 ;
说明mid答案合理 ,继续找更小的合理答案 ;
注意很多细节,数据会超int,两点间可能有很多边 ;
#include<cstdio>
#include<cstring>
#include<vector>
#include<cmath>
#include<cstdlib>
#include<stack>
#include<queue>
#include <iomanip>
#include<iostream>
#include<algorithm>
using namespace std ;
#define MAXN 405
#define INF 1e8
struct edge
{
int u,v,w,next;
}E[200000];
int head[MAXN],ecnt;
int gap[MAXN],cur[MAXN],pre[MAXN],dis[MAXN];
long long map[205][205];
int N,M,scr,sink,vn,num;
int A[205],B[205];
void Insert(int u,int v,int w)
{
E[ecnt].u=u;
E[ecnt].v=v;
E[ecnt].w=w;
E[ecnt].next=head[u];
head[u]=ecnt++;
E[ecnt].u=v;
E[ecnt].v=u;
E[ecnt].w=0;
E[ecnt].next=head[v];
head[v]=ecnt++;
}
void Init()
{
int i,u,v;
long long c;
memset(map,-1,sizeof(map));
scr=0;sink=2*N+1;vn=sink+1;num=0;
for(i=1;i<=N;i++)
{
scanf("%d%d",&A[i],&B[i]);
num+=A[i];
}
for(i=1;i<=M;i++)
{
scanf("%d%d%I64d",&u,&v,&c);
if(map[u][v]==-1||map[u][v]>c)
{map[u][v]=map[v][u]=c;}
}
}
void Floyd()
{
int i,j,k;
for(k=1;k<=N;k++)
{
for(i=1;i<=N;i++)
{
if(i==k||map[i][k]==-1) continue;
for(j=1;j<=N;j++)
{
if(j==k||i==j||map[j][k]==-1) continue;
if(map[i][j]==-1||map[i][j]>map[i][k]+map[k][j])
map[i][j]=map[i][k]+map[k][j];
}
}
}
}
int Sap(int s,int t,int n)//核心代码(模版)
{
int ans=0,aug=INF;//aug表示增广路的流量
int i,v,u=pre[s]=s;
for(i=0;i<=n;i++)
{
cur[i]=head[i];
dis[i]=gap[i]=0;
}
gap[s]=n;
bool flag;
while(dis[s]<n)
{
flag=false;
for(int &j=cur[u];j!=-1;j=E[j].next)//一定要定义成int &j,why
{
v=E[j].v;
if(E[j].w>0&&dis[u]==dis[v]+1)
{
flag=true;//找到容许边
aug=min(aug,E[j].w);
pre[v]=u;
u=v;
if(u==t)
{
ans+=aug;
while(u!=s)
{
u=pre[u];
E[cur[u]].w-=aug;
E[cur[u]^1].w+=aug;//注意
}
aug=INF;
}
break;//找到一条就退出
}
}
if(flag) continue;
int mindis=n;
for(i=head[u];i!=-1;i=E[i].next)
{
v=E[i].v;
if(E[i].w>0&&dis[v]<mindis)
{
mindis=dis[v];
cur[u]=i;
}
}
if((--gap[dis[u]])==0) break;
gap[dis[u]=mindis+1]++;
u=pre[u];
}
return ans;
}
int Judge(long long mid)//构图+判定
{
int i,j;
memset(head,-1,sizeof(head));ecnt=0;
for(i=1;i<=N;i++)
{
Insert(scr,i,A[i]);
Insert(i+N,sink,B[i]);
Insert(i,i+N,INF);
}
for(i=1;i<=N;i++)
{
for(j=i+1;j<=N;j++)
{
if(map[i][j]==-1) continue;
if(map[i][j]<=mid)//注意要连入点-->出点(INF)
{
Insert(i,j+N,INF);
Insert(j,i+N,INF);
}
}
}
return Sap(scr,sink,vn)==num;
}
void Solve()
{
int i,j;
long long l,r,mid,Max=0;
Floyd();
for(i=1;i<=N;i++)
{
for(j=i+1;j<=N;j++)
Max=max(Max,map[i][j]);
}
l=0;r=Max;
while(l!=r)
{
mid=(l+r)/2;
if(Judge(mid))
r=mid;
else
l=mid+1;
}
if(Judge(l))
printf("%I64d\n",l);
else
printf("-1\n");
}
int main()
{
while(scanf("%d%d",&N,&M)!=EOF)
{
Init();
Solve();
}
return 0;
}
,