当前位置: 代码迷 >> 综合 >> POJ 2391 floyd+二分+最大流
  详细解决方案

POJ 2391 floyd+二分+最大流

热度:24   发布时间:2024-01-20 20:27:31.0

思路构图还是十分的好想的... 比较的有经验了。对这题,没啥难度了。

题目大意:

FJ有很多田地,一些田地有雨棚,雨棚有容纳牛的容量,现在给你田中的牛和田里的雨棚的容纳量,还有田间路径走完花费的时间。

要求出使得所有牛都躲进雨棚里的最小时间。

构图还是很简单的:

将牛棚拆点为in和out点。

源点到in容量为牛棚中牛的数量,out到汇点为牛棚的容量。in--out为无穷大。

二分枚举时间,若两牛棚的时间小于枚举时间,则in1--out2有边容量为无穷大。

soso....

这么一写就过了....

这题恶心的地方就是时间枚举要long long 而坑爹的我找不到__int64 的最大值常量了= =

后来才发现有这样写的= = long long inf=1ll<<60;

#include<iostream>
#include<string>
#define MN 444
#define INF 0x3FFFFFFF
#define FF(i,N) for( int i=0;i<N;i++ )
#define CC(a) memset( a,0,sizeof(a) )
template<class T> void inline checkmin( T &a,T b ){ if( a>b||a==-1 ) a=b; }
using namespace std;const long long oo=1ll<<58;struct edge
{int u,v,len;
}E[2222];int F,P,s,t,tot;
__int64 map[MN][MN];
int maze[MN][MN],cap[MN][MN];bool cmp( edge a,edge b ){ return a.len<b.len; }void setG()
{memset( map,0x3F,sizeof(map) );CC(maze),CC(cap);s=0;t=2*F+1;tot=0;FF( i,F ){scanf( "%d%d",&maze[s][i+1],&maze[i+1+F][t] );maze[i+1][i+1+F]=INF;tot+=maze[s][i+1];}int u,v,len;FF( i,P ){scanf( "%d%d%d",&u,&v,&len );if( map[u][v]>len )map[u][v]=map[v][u]=len;}for( int i=1;i<=F;i++ )map[i][i]=0;for( int k=1;k<=F;k++ )for( int i=1;i<=F;i++ )for( int j=1;j<=F;j++ )if( map[i][j]>map[i][k]+map[k][j] )map[i][j]=map[i][k]+map[k][j];
}void initG( __int64 mid )
{for( int i=0;i<=t;i++ )for( int j=0;j<=t;j++ )cap[i][j]=maze[i][j];for( int i=1;i<=F;i++ )for( int j=1;j<=F;j++ )if( map[i][j]<=mid )cap[i][j+F]=INF;
}int cur[MN],pre[MN],dis[MN],gap[MN];
int sap()
{CC(cur),CC(pre),CC(dis),CC(gap);int u=pre[s]=s,maxflow=0,aug=-1;gap[0]=t+1;while( dis[s]<=t ){
loop:	   for( int v=cur[u];v<=t;v++ ){if( cap[u][v]&&dis[u]==dis[v]+1 ){pre[v]=u;cur[u]=v;checkmin( aug,cap[u][v] );u=v;if( v==t ){maxflow+=aug;for( u=pre[u];v!=s;v=u,u=pre[u] )cap[u][v]-=aug,cap[v][u]+=aug;aug=-1;}goto loop;}}int mind=t;for( int v=0;v<=t;v++ )if( cap[u][v]&&mind>dis[v] ){mind=dis[v];cur[u]=v;}if( --gap[dis[u]]==0 )break;gap[dis[u]=mind+1]++;u=pre[u];}return maxflow;
}int main()
{while( scanf("%d%d",&F,&P)!=EOF ){setG();__int64 l=1,r=(oo),mid;while( (mid=(l+r)/2)&&l<r ){initG(mid);if( tot==sap() ) r=mid;else l=mid+1;}if( mid==(oo) )printf( "-1\n" );elseprintf( "%I64d\n",mid );}return 0;
}