做上下界网络流的题姿势水平还有待提高啊~
把每个国家拆成两个点
Ai1,Ai2
,限制这两个点之间的流量,让其上下界都为
Vi
先建立超级源点
S
和超级汇点
Ai1−>T
,
S−>Ai2
流量都为
Vi
费用为零
在建立一个点
s−>Ai1
流量为无穷费用为零
S−>s
流量为
m
费用为零
对于原图中的边
然后跑最小费用最大流就是答案
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<vector>
#include<queue>
#include<cmath>
#include<set>
#include<map>
#define N 400006
#define mx 1011111111
using namespace std;
int sc()
{
int i=0,f=1;char c=getchar();
while(c>'9'||c<'0'){if(c=='-')f=-1;c=getchar();}
while(c>='0'&&c<='9')i=i*10+c-'0',c=getchar();
return i*f;
}
int dis[N],mark[N],from[N],num[N],q[N];
int head[N],lst[N],nxt[N],c[N],v[N];
int n,m,s,S,T,tot=1,ans;
void insert(int x,int y,int a,int b)
{
lst[++tot]=y;nxt[tot]=head[x];head[x]=tot;c[tot]=a;v[tot]=b;
lst[++tot]=x;nxt[tot]=head[y];head[y]=tot;c[tot]=0;v[tot]=-b;
}
bool spfa()
{
for(int i=1;i<=T;i++)dis[i]=mx;
dis[q[1]=S]=0;
for(int l=1,r=2;l!=r;l=(l+1)%N)
{
int x=q[l];mark[x]=0;
for(int i=head[x];i;i=nxt[i])
if(c[i]&&dis[lst[i]]>dis[x]+v[i])
{
dis[lst[i]]=dis[x]+v[i];
from[lst[i]]=x;
num[lst[i]]=i;
if(!mark[lst[i]])
mark[q[r]=lst[i]]=1,r=(r+1)%N;
}
}
//cout << dis[T]<<endl;
return dis[T]!=mx;
}
void mcf()
{
int w=mx;
for(int i=T;i!=S;i=from[i])
w=min(w,c[num[i]]);
for(int i=T;i!=S;i=from[i])
c[num[i]]-=w,c[num[i]^1]+=w;
ans+=w*dis[T];
}
int main()
{
n=sc(),m=sc();s=2*n+1,S=s+1,T=s+2;
insert(S,s,m,0);
for(int i=1;i<=n;i++)
{
int x=sc();
insert(s,i,mx,0);
insert(i,T,x,0);
insert(S,i+n,x,0);
}
for(int i=1;i<=n;i++)
for(int j=i+1;j<=n;j++)
{
int x=sc();
if(x==-1)continue;
insert(i+n,j,mx,x);
}
while(spfa())mcf();
cout<<ans;
return 0;
}