BZOJ 1266 AHOI2006 上学路线route Floyd+最小割

题目大意:给定一张图,每条边有一个长度和一个花费,要求删掉一些边使1到n的最短路变长,求最小花销

首先求出最短路(用什么求随便,反正数据范围小),然后将所有在最短路上的边连到新图中,求最小割就是答案

图没有重边- - 数组开小WA了半篇- -

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define M 510
#define S 1
#define T n
#define INF 0x7fffffff
using namespace std;
struct edge{
	int x,y,f,c;
}edges[124800];
int n,m,ans;
int map[M][M];
namespace Max_Flow{
	struct abcd{
		int to,next,f;
	}table[124800<<2];
	int head[M],tot=1;
	int dpt[M];
	void Add(int x,int y,int z)
	{
		table[++tot].to=y;
		table[tot].f=z;
		table[tot].next=head[x];
		head[x]=tot;
	}
	inline void Link(int x,int y,int z)
	{
		Add(x,y,z);
		Add(y,x,0);
	}
	bool BFS()
	{
		static int q[M];
		int i,r=0,h=0;
		memset(dpt,-1,sizeof dpt);
		dpt[S]=1;q[++r]=S;
		while(r!=h)
		{
			int x=q[++h];
			for(i=head[x];i;i=table[i].next)
				if(table[i].f&&!~dpt[table[i].to])
				{
					dpt[table[i].to]=dpt[x]+1;
					q[++r]=table[i].to;
					if(table[i].to==T)
						return true;
				}
		}
		return false;
	}
	int Dinic(int x,int flow)
	{
		int i,left=flow;
		if(x==T) return flow;
		for(i=head[x];i&&left;i=table[i].next)
			if(table[i].f&&dpt[table[i].to]==dpt[x]+1)
			{
				int temp=Dinic(table[i].to,min(left,table[i].f) );
				//if(!temp) dpt[table[i].to]=-1;
				left-=temp;
				table[i].f-=temp;
				table[i^1].f+=temp;
			}
		if(left) dpt[x]=-1;
		return flow-left;
	}
}
int main()
{
	using namespace Max_Flow;
	int i,j,k,x,y,z;
	cin>>n>>m;
	memset(map,0x3f,sizeof map);
	for(i=1;i<=n;i++)
		map[i][i]=0;
	for(i=1;i<=m;i++)
	{
		scanf("%d%d%d%d",&x,&y,&z,&edges[i].c);
		edges[i].x=x;edges[i].y=y;edges[i].f=z;
		map[x][y]=map[y][x]=z;
	}
	for(k=1;k<=n;k++)
		for(i=1;i<=n;i++)
			for(j=1;j<=n;j++)
				map[i][j]=min(map[i][j],map[i][k]+map[k][j]);
	cout<<map[1][n]<<endl;
	for(i=1;i<=m;i++)
	{
		x=edges[i].x;y=edges[i].y;z=edges[i].f;
		if( map[1][x]+map[y][n]+z==map[1][n] )
			Link(x,y,edges[i].c);
		if( map[1][y]+map[x][n]+z==map[1][n] )
			Link(y,x,edges[i].c);
	}
	while( BFS() )
		ans+=Dinic(S,INF);
	cout<<ans<<endl;
	return 0;
}


你可能感兴趣的:(最小割,floyd,dinic,bzoj,BZOJ1266)