【bzoj2055】80人環遊世界 有上下界的費用流
阿新 • • 發佈:2019-01-23
對有上下界的網路流理解的還不夠。
有上下界的網路流
其實是無源無匯的有上下界的最小費用最大流(好像是吧?無源無匯是什麼意思來著?)
每個人可以從任意一個點出發,每個點正好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; }