省选专练【80人环游世界】

“　想必大家都看过成龙大哥的《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();
}

转载于:https://www.cnblogs.com/Leo-JAM/p/10079283.html