BZOJ3511 土地划分-最小割

各种retle也真是不想说什么了

这种类型的最小割模板题大全集

题目链接：右转进入题目

题目大意：请自行参考原题

解：

对第i个城市建点C=i，对第i条边两端点X和Y建点A=i+n和B=i+n+m，建源点S=0，汇点T=n+2*m+1。

建边：(S,C,VA),(C,T,VB),(S,A,EA),(B,T,EB),(A,X,INF),(A,Y,INF),(X,B,INF),(Y,B,INF),(X,Y,EC),(Y,X,EC)。

ans=sigma{VA}+sigma{VB}+sigma{EA}+sigma{EB}-最小割

特别的，对于城市1，只建立(S,1,INF)。

对城市n，只建立(n,T,INF)。

没了

#include<queue>
#include<cstdio>
#include<climits>
#include<cstring>
#include<iostream>
#include<algorithm>
#define MAXN 100010
#define MAXM 1000000
#define INF INT_MAX
#define build_edge(u,v,c) add_edge(u,v,c),add_edge(v,u,0)
using namespace std;
struct edges{
	int to,pre,resf;
}e[MAXM];
int etop,h[MAXN],lev[MAXN],cur[MAXN];
bool vis[MAXN];queue<int> q;
int add_edge(int u,int v,int c)
{
	etop++;
	e[etop].to=v;
	e[etop].pre=h[u];
	h[u]=etop;
	e[etop].resf=c;
	return 0;
}
bool bfs(int s,int t)
{
	memset(vis,false,sizeof(vis));
	memset(lev,0,sizeof(lev));
	while(!q.empty()) q.pop();
	q.push(s);lev[s]=0;vis[s]=true;
	while(!q.empty())
	{
		int x=q.front();q.pop();
		for(int i=h[x];i;i=e[i].pre)
			if(!vis[e[i].to]&&e[i].resf)
			{
				vis[e[i].to]=true;
				lev[e[i].to]=lev[x]+1;
				q.push(e[i].to);
			}
	}
	return vis[t];
}
int dfs(int s,int t,int a)
{
	if(s==t||!a) return a;
	int flow=0,f;
	for(int &i=cur[s];i;i=e[i].pre)
		if(lev[e[i].to]==lev[s]+1&&(f=dfs(e[i].to,t,min(a,e[i].resf)))>0)
		{
			a-=f;flow+=f;e[i].resf-=f;
			e[((i-1)^1)+1].resf+=f;
			if(!a) break;
		}
	return flow;
}
int main()
{
	int n,m;scanf("%d%d",&n,&m);
	int s=0,t=n+2*m+1,totans=0;
	build_edge(s,1,INF);
	build_edge(n,t,INF);
	for(int i=2;i<n;i++)
	{
		int va;scanf("%d",&va);
		build_edge(s,i,va);totans+=va;
	}
	for(int i=2;i<n;i++)
	{
		int vb;scanf("%d",&vb);
		build_edge(i,t,vb);totans+=vb;
	}
	for(int i=1;i<=m;i++)
	{
		int x,y,a,b,c;scanf("%d%d%d%d%d",&x,&y,&a,&b,&c);
		int A=i+n,B=i+n+m;totans+=a+b;
		add_edge(x,y,c);add_edge(y,x,c);
		build_edge(s,A,a);build_edge(B,t,b);
		build_edge(A,x,INF);build_edge(A,y,INF);
		build_edge(x,B,INF);build_edge(y,B,INF);
	}
	int flow=0;
	while(bfs(s,t))
	{
		for(int i=s;i<=t;i++)
			cur[i]=h[i];
		flow+=dfs(s,t,INF);
	}
	printf("%d\n",totans-flow);
	return 0;
}