BZOJ 3901 Magic 最小割

题目大意：给定一张有向图，每个点有一个权值，你可以付出b[i]的代价将第i个点的权值变成0，一个点的代价为所有出边指向的点的代价的最大值，要求代价之和最小

萌哒哒的zxr曾经给我们讲过这个建图方式~ 还记得真是太好了poi~

将每个点的所有出边指向的点按照权值从大到小建成从源点出发的一条链，一个点与之前的点的边流量为这个点的权值

然后把一个点所有分出去的点全都连到一起，然后向汇点连一条流量为这个点变成0的代价的边

跑最小割就是答案

这里一个点如果割掉和T的连边就相当于将这个点的权值清零

一个点的代价就是沿着它建出的链扫到的第一个没有清零的点

为了保证最小割因此要把这个点与之前的点的连边割掉 就会产生这个点权值的代价

图不是很好画请大家对着代码强行理解这段话吧- -

#include <vector>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define M 55000
#define S 0
#define T (M-1)
#define INF 0x3f3f3f3f
using namespace std;
int n,m,cnt,ans;
int a[M],b[M];
vector<int> to[M];
bool Compare(int x,int y)
{
	return a[x] < a[y];
}
namespace Max_Flow{
	struct abcd{
		int to,f,next;
	}table[1001001];
	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;
	}
	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);
		q[++r]=S;dpt[S]=1;
		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;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) );
				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,x,y,cnt;
	cin>>n>>m;cnt=n;
	for(i=1;i<=n;i++)
		scanf("%d",&a[i]);
	for(i=1;i<=n;i++)
		scanf("%d",&b[i]);
	for(i=1;i<=m;i++)
	{
		scanf("%d%d",&x,&y);
		to[x].push_back(y);
	}
	for(i=1;i<=n;i++)
	{
		vector<int>::iterator it;
		sort(to[i].rbegin(),to[i].rend(),Compare);
		int last=S;
		for(it=to[i].begin();it!=to[i].end();it++)
		{
			Link(last,++cnt,a[*it]);
			Link(cnt,*it,INF);
			last=cnt;
		}
		Link(i,T,b[i]);
	}
	while( BFS() )
		ans+=Dinic(S,INF);
	cout<<ans<<endl;
	return 0;
}