“ 想必大家都看过成龙大哥的《80天环游世界》,里面的紧张刺激的打斗场面一定给你留下了深刻的印象。现在就有这么一个80人的团伙,也想来一次环游世界。
他们打算兵分多路,游遍每一个国家。因为他们主要分布在东方,所以他们只朝西方进军。设从东方到西方的每一个国家的编号依次为1...N。假若第i个人的游历路线为P1、P2......Pk(0≤k≤N),则P1 < P2 < ...... < Pk。
众所周知,中国相当美丽,这样在环游世界时就有很多人经过中国。我们用一个正整数Vi来描述一个国家的吸引程度,Vi值越大表示该国家越有吸引力,同时也表示有且仅有Vi个人会经过那一个国家。
为了节省时间,他们打算通过坐飞机来完成环游世界的任务。同时为了省钱,他们希望总的机票费最小。
明天就要出发了,可是有些人临阵脱逃,最终只剩下了M个人去环游世界。他们想知道最少的总费用,你能告诉他们吗?”
无源汇上下界费用流。
不存在源汇,自己先建一个超级源汇。
然后连上INF
构建一个发配源,连上超级源流m费0
构建一个新超级源汇。
”对于每个点,拆点:
超级源点向x连一条容量为v[x],费用为0的边;
x'向超级汇点连一条容量为v[x],费用为0的边;
源点向x连一条容量为inf,费用为0的边;
x'向汇点连一条容量为inf,费用为0的边。
对于原图一条u->v,费用为w的边:
u'向v连一条容量为inf,费用为w的边。“
这时便省略了发配源。原理是你不需要那个INF的超级汇到超级源的流量,它必然会被发配源限制。
#include<bits/stdc++.h>
using namespace std;
const int N=20000;
const int INF=1e9+7;
struct Front_satr{
int u,v,w,c,nxt;
}e[N*4];
int cnt=1;
int first[N]={0};
void addedge(int u,int v,int w,int c){
cnt++;
e[cnt].u=u;
e[cnt].v=v;
e[cnt].w=w;
e[cnt].c=c;
e[cnt].nxt=first[u];
first[u]=cnt;
}
void add(int u,int v,int w,int c){
addedge(u,v,w,c);
addedge(v,u,0,-c);
}
int n,m;
int dis[N]={0};
int pre[N]={0};
int inqueue[N]={0};
queue<int> q;
int S,SS,ST,T;
bool SPFA(){
for(int i=0;i<=T;i++)dis[i]=INF;
q.push(S);
dis[S]=0;
while(!q.empty()){
int u=q.front();
q.pop();
inqueue[u]=0;
for(int i=first[u];i;i=e[i].nxt){
int v=e[i].v;
if(e[i].w&&dis[v]>dis[u]+e[i].c){
dis[v]=dis[u]+e[i].c;
pre[v]=i;
if(!inqueue[v]){
q.push(v);
inqueue[v]=1;
}
}
}
}
return dis[T]!=INF;
}
void MCMF(){
int ret=0;
while(SPFA()){
int s=INF;
for(int i=pre[T];i;i=pre[e[i^1].v])
s=min(s,e[i].w);
for(int i=pre[T];i;i=pre[e[i^1].v]){
e[i].w-=s;
e[i^1].w+=s;
}
ret+=dis[T]*s;
}
cout<<ret;
}
int main(){
scanf("%d%d",&n,&m);
SS=n*2+1;
ST=n*2+2;
S=n*2+3;
T=n*2+4;
add(ST,SS,m,0);
for(int i=1;i<=n;i++){
int x;
scanf("%d",&x);
add(S,i+n,x,0);
add(i,T,x,0);
add(SS,i,INF,0);
add(i+n,ST,INF,0);
}
for(int i=1;i<=n;i++){
for(int j=i+1;j<=n;j++){
int x;
scanf("%d",&x);
if(x==-1)
continue;
add(i+n,j,INF,x);
}
}
MCMF();
}