HDU6214 2017青岛赛区网络赛

本弱并不会骚套路，当时并想了很久做不出。

将每个边的流量w变成w*(m+1)+1，因为边总数为m，所以后面那个1并不会影响最大流等于最小割的性质，再加上后面那个1体现了边数，跑出来的最大流为a*(m+1)+b，%掉m+1,就得到b，而每一种最大流等于最小割都可以得到相同的a，但b是不同的，由于最大流等于最小割，这个b最大就是最小割的边数。

#include<cstdio>
#include<cstring>
#define maxn 210
#define maxl 1010
#define inf 2000000001

int n,m,S,T,cnt,ans;
int cur[maxn],dis[maxn],a[maxn],ehead[maxn];
struct ed{int to,nxt,w;} e[maxl<<1];

void add(int u,int v,int w)
{
	e[++cnt].nxt=ehead[u];e[cnt].to=v;e[cnt].w=w;ehead[u]=cnt;
}

void prework()
{
	scanf("%d%d",&n,&m);
	scanf("%d%d",&S,&T);
	for(int i=1;i<=n;i++)
		ehead[i]=0;
	int u,v,w;
	cnt=0;
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d%d",&u,&v,&w);
		add(u,v,w*(m+1)+1);
		add(v,u,0);	
	}
}

bool bfs()
{
	int u,v,head=0,tail=1;
	for(int i=1;i<=n;i++)
		dis[i]=-1;
	a[1]=S;dis[S]=0;
	while(head<tail)
	{
		u=a[++head];
		for(int i=ehead[u];i;i=e[i].nxt)
		{
			v=e[i].to;
			if(dis[v]<0 && e[i].w)
			{
				dis[v]=dis[u]+1;
				a[++tail]=v;
			}
		}
	}
	if(dis[T]<0)
		return false;
	else
		return true;
}

int min(int a,int b)
{
	if(a<b)
		return a;
	else
		return b;
}

int find(int u,int low)
{
	int t,v;
	if(u==T)
		return low;
	for(int &i=cur[u];i;i=e[i].nxt)
	{
		v=e[i].to;
		if(dis[v]==dis[u]+1 && e[i].w && 
		(t=find(v,min(low,e[i].w))))
		{
			if(i&1)
				e[i+1].w+=t;
			else
				e[i-1].w+=t;
			e[i].w-=t;
			return t;
		}
	}
	return 0;
}

void mainwork()
{
	int t;
	ans=0;
	while(bfs())
	{
		for(int i=1;i<=n;i++)
			cur[i]=ehead[i];
		while((t=find(S,inf)))
			ans+=t;
	}
}

void print()
{
	ans%=m+1;
	printf("%d\n",ans);
}

int main()
{
	int t;
	scanf("%d",&t);
	for(int i=1;i<=t;i++)
	{
		prework();
		mainwork();
		print();
	}
	return 0;
}

据说以前有一道原题是跑完一遍以后吧满流设为1，否则为inf，再跑一遍，这样感觉没什么道理，最大流的时候最小割的数量是无法保证的，不同的跑法得到的满流边集都不一样，那道题是数据弱，这里就过不了了。然而又有人说反向边不改流量这样也可以过，这。。。就更没道理了，你不反悔那还叫网络流的算法？数据是会玩的。、、