【bzoj2055】80人环游世界 有上下界的费用流

对有上下界的网络流理解的还不够。

有上下界的网络流
其实是无源无汇的有上下界的最小费用最大流(好像是吧?无源无汇是什么意思来着?)
每个人可以从任意一个点出发,每个点正好v个人
超级源点S,超级汇点T
源点s向每个地点连一条容量为inf费用为0的边
每个地点拆成两个点,连一条上下界为v费用为0的边
超级源点S向源点s连一条容量为m费用为0的边(表示s只能流出m的流量)
剩下的边随便连


#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#include<iostream>
#define maxn 510
#define maxm 500010 
#define inf 1000000000 

using namespace std;

int head[maxn],p[maxm],fr[maxn],to[maxm],c[maxm],len[maxm],next[maxm],q[maxn],dis[maxn];
int d[maxn];
bool vis[maxn];
int n,m,num,s,t,S,T,ans;

void addedge(int x,int y,int z,int w)
{
	num++;p[num]=x;to[num]=y;c[num]=z;len[num]=w;next[num]=head[x];head[x]=num;
	num++;p[num]=y;to[num]=x;c[num]=0;len[num]=-w;next[num]=head[y];head[y]=num;
}

bool spfa()
{
	for (int i=0;i<=T;i++) dis[i]=inf;
	int l=0,r=1;
	q[1]=S;dis[S]=0;vis[S]=1;
	while (l!=r)
	{
		l++;if (l==maxn) l=0;
		int x=q[l];
		for (int p=head[x];p;p=next[p])
		  if (c[p] && dis[x]+len[p]<dis[to[p]])
		  {
		  	dis[to[p]]=dis[x]+len[p];
		  	fr[to[p]]=p;
		  	if (!vis[to[p]])
		  	{
		  		r++;if (r==maxn) r=0;
		  		q[r]=to[p];vis[to[p]]=1;
		  	}
		  }
		vis[x]=0;
	}
	if (dis[T]==inf) return 0; else return 1;
}

void mcf()
{
	int x=inf;
	for (int i=fr[T];i;i=fr[p[i]]) x=min(x,c[i]);
	for (int i=fr[T];i;i=fr[p[i]]) ans+=x*len[i],c[i]-=x,c[i^1]+=x;
}

int main()
{
	scanf("%d%d",&n,&m);
	num=1;s=0;S=201;T=202;
	addedge(S,s,m,0);
	for (int i=1;i<=n;i++)
	{
		int x;
		scanf("%d",&x);
		addedge(s,i,inf,0);
		d[i]-=x;d[i+n]+=x;
	}
	for (int i=1;i<=n;i++)
	  for (int j=i+1;j<=n;j++)
	  {
	  	int x;
	  	scanf("%d",&x);
	  	if (x!=-1) addedge(i+n,j,inf,x);
	  }
	for (int i=1;i<=n*2;i++)
	  if (d[i]>0) addedge(S,i,d[i],0);
	  else if (d[i]<0) addedge(i,T,-d[i],0);
	while (spfa()) mcf();
	printf("%d\n",ans);
	return 0;
}


你可能感兴趣的:(【bzoj2055】80人环游世界 有上下界的费用流)