第 45 届ICPC亚洲区域赛（昆明）B Chessboard 上下界费用流 负环处理

Solution:

用的是题解的建图（详见出题人知乎回答）。然而上下界费用流为了平衡源汇点的流量，会有一条汇点到源点的边，而负权边的出现会导致原图存在负环。采用以下方法处理：

模板参考的是这位大佬。
存一下模板。

Code

#include<bits/stdc++.h>
using namespace std;
#define LL long long
#define pa pair<int,int>
const int Maxn=55;
const int N=160,M=2500000;
const int inf=1e9;
LL read()
{
	LL x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9')x=x*10ll+ch-'0',ch=getchar();
	return x*f;
}
struct Edge
{
	int x,y,d,c,next;
}e[M];
int last[N],len=1;
int n,m,sb[Maxn][Maxn],sw[Maxn][Maxn],l[Maxn],r[Maxn],L[Maxn],R[Maxn],ans=0;
int st,ed,S,T,SS,TT;
int f[N],from[N],D[N],E[N];bool in[N];
queue<int>q;
void ins(int x,int y,int d,int c)
{
	int t=++len;
	e[t].x=x;e[t].y=y;e[t].d=d;e[t].c=c;
	e[t].next=last[x];last[x]=t;
}
void addedge(int x,int y,int d,int c)
{
	if(c<0)
	{
		ans+=d*c;E[y]+=d,E[x]-=d; 
		ins(x,y,0,c),ins(y,x,d,-c);
	}
	else ins(x,y,d,c),ins(y,x,0,-c);
}
void Add(int x,int y,int lower,int upper,int c)
{
	ans+=c*lower;D[y]+=lower,D[x]-=lower;
	addedge(x,y,upper-lower,c);
}
bool spfa(int p1,int p2,bool flag)//true 最小费用最大流 false 最小费用可行流 
{
	for(int i=0;i<=TT;i++)f[i]=inf,in[i]=false,from[i]=-1;
	f[p1]=0;q.push(p1);
	while(!q.empty())
	{
		int x=q.front();q.pop();in[x]=false;
		for(int i=last[x];i;i=e[i].next)
		{
			int y=e[i].y;
			if(e[i].d&&f[x]+e[i].c<f[y])
			{
				f[y]=f[x]+e[i].c;
				from[y]=i;
				if(!in[y])in[y]=true,q.push(y);
			}
		}
	}
	if(flag)return(f[p2]!=inf);
	return(f[p2]<=0);
}
void MinCost(int p1,int p2,bool flag)
{
	while(spfa(p1,p2,flag))
	{
		int x=p2,mn=inf;
		while(x!=p1)
		{
			mn=min(mn,e[from[x]].d);
			x=e[from[x]].x;
		}
		x=p2;
		while(x!=p1)
		{
			ans+=mn*e[from[x]].c;
			e[from[x]].d-=mn,e[from[x]^1].d+=mn;
			x=e[from[x]].x;
		}
	}
}
int main()
{
	n=read(),m=read();
	for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)sb[i][j]=read(),ans+=sb[i][j];
	for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)sw[i][j]=read();
	for(int i=1;i<=n;i++)l[i]=read(),r[i]=read();
	for(int j=1;j<=m;j++)L[j]=read(),R[j]=read();
	st=n+m+1,ed=st+1,S=st+2,T=st+3,SS=st+4,TT=st+5;
	for(int i=1;i<=n;i++)Add(st,i,n-r[i],n-l[i],0);
	for(int j=1;j<=m;j++)Add(n+j,ed,m-R[j],m-L[j],0);
	for(int i=1;i<=n;i++)
	for(int j=1;j<=m;j++)
	{
		addedge(i,n+j,1,-sb[i][j]);
		addedge(i,n+j,1,sw[i][j]);
	}
	int l1,l2;
	l1=len;
	addedge(ed,st,inf,0);
	for(int i=1;i<=TT;i++)
	{
		if(D[i]>0)addedge(S,i,D[i],0);
		else addedge(i,T,-D[i],0);
	}
	l2=len;
	for(int i=1;i<=TT;i++)
	{
		if(E[i]>0)addedge(SS,i,E[i],0);
		else addedge(i,TT,-E[i],0);
	}
	MinCost(SS,TT,true);
	for(int x=1;x<=TT;x++)
	while(last[x]>l2)last[x]=e[last[x]].next;
	MinCost(S,T,true);
	for(int x=1;x<=TT;x++)
	while(last[x]>l1)last[x]=e[last[x]].next;
	MinCost(st,ed,false);
	printf("%d",ans);
}