POJ 2391 floyd+二分+最大流

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

题目大意:

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" );
  		   else
  		   	   printf( "%I64d\n",mid );
  	}
 	return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值