[bzoj2055]80人环游世界 上下界费用流

2055: 80人环游世界

Time Limit: 10 Sec   Memory Limit: 64 MB
[ Submit][ Status][ Discuss]

Description

    想必大家都看过成龙大哥的《80天环游世界》，里面的紧张刺激的打斗场面一定给你留下了深刻的印象。现在就有这么

    一个80人的团伙，也想来一次环游世界。

    他们打算兵分多路，游遍每一个国家。

    因为他们主要分布在东方，所以他们只朝西方进军。设从东方到西方的每一个国家的编号依次为1...N。假若第i个人的游历路线为P1、P2......Pk(0≤k≤N)，则P1<P2<......<Pk。

    众所周知，中国相当美丽，这样在环游世界时就有很多人经过中国。我们用一个正整数Vi来描述一个国家的吸引程度，Vi值越大表示该国家越有吸引力，同时也表示有且仅

有Vi个人会经过那一个国家。

    为了节省时间，他们打算通过坐飞机来完成环游世界的任务。同时为了省钱，他们希望总的机票费最小。

    明天就要出发了，可是有些人临阵脱逃，最终只剩下了M个人去环游世界。他们想知道最少的总费用，你能告诉他们吗？

Input

第一行两个正整数N，M。

第二行有N个不大于M正整数，分别表示V1，V2......VN。

接下来有N-1行。第i行有N-i个整数，该行的第j个数表示从第i个国家到第i+j个国家的机票费（如果该值等于-1则表示这两个国家间没有通航）。

Output

在第一行输出最少的总费用。

Sample Input

6 3
2 1 3 1 2 1
2 6 8 5 0
8 2 4 1
6 1 0
4 -1
4

Sample Output

27

HINT

1<= N < =100 1<= M <= 79

Source

上下界原图（很简单）

s->S m,m,0

S->xi 0,inf,0

yi->T 0,inf,0

xi->yi vi,vi,0

yi->xj 0,inf,ci

如何改造为上下界运行图

超超级源汇SS，TT

u->v  x,y,ci 改为u->v y-x,ci

bi = 流入下界和-流出下界和

对于bi>0的点,SS->i,bi,0

        bi<0的点,i->TT,bi,0

然后费用流

赶脚这个zkw费用流是错的，很多题都A不了

#include<iostream>
#include<cstring>
#include<cstdio>
#define inf 0x7fffffff
using namespace std;
const int N = 205;
int n,m,cnt=1,ans,S,T;
int last[N], b[N], q[N], d[N], inq[N], from[N];
struct Edge{
	int to,v,next,c,from;
}e[N*N*5];
void insert( int u, int v, int w, int c ){
	e[++cnt].to = v; e[cnt].c = c; e[cnt].v = w; e[cnt].next = last[u]; last[u] = cnt; e[cnt].from = u;
	e[++cnt].to = u; e[cnt].c = -c; e[cnt].v = 0; e[cnt].next = last[v]; last[v] = cnt; e[cnt].from = v;
}/*
bool spfa(){
	memset(inq,0,sizeof(inq));
	for( int i = 0; i <= T; i++ ) d[i] = inf;
	int head = 0, tail = 1;
	inq[T] = 1; d[T] = 0; q[0] = T;
	while( head != tail ){
		int now = q[head++]; if( head == T ) head = 0;
		for( int i = last[now]; i; i = e[i].next )
			if( e[i^1].v && d[now]+e[i^1].c < d[e[i].to] ){
				d[e[i].to] = d[now] + e[i^1].c;
				if( !inq[e[i].to] ){
					inq[e[i].to] = 1;
					q[tail++] = e[i].to; if( tail == T ) tail = 0;
				}
			}
	}
	return d[S] != inf;
}
int dfs( int x, int f ){
	inq[x] = 1;
	if( x == T ) return f;
	int w,used=0;
	for( int i = last[x]; i; i = e[i].next )
		if( d[e[i].to] == d[x]-e[i].c && e[i].v && !inq[e[i].to] ){
			w = dfs( e[i].to, min(f-used,e[i].v) );
			used += w; e[i].v -= w; e[i^1].v += w;
			ans += e[i].c*w; if( used == f ) return f;
		}
	return used;
}
void zkw(){
	while( spfa() ){
		inq[T] = 1;
		while( inq[T] ){
			memset(inq,0,sizeof(inq));
			dfs( S, inf );
		}
	}
}*/
bool spfa(){
	for( int i = 0; i <= T; i++ ) d[i] = inf;
	int head = 0, tail = 1;
	d[0] = 0; q[0] = 0; inq[0] = 1;
	while( head != tail ){
		int now = q[head++]; if( head == T ) head = 0;
		for( int i = last[now]; i; i = e[i].next )
			if( e[i].v && e[i].c+d[now] < d[e[i].to] ){
				d[e[i].to] = e[i].c+d[now]; from[e[i].to] = i;
				if( !inq[e[i].to] ){
					inq[e[i].to] = 1;
					q[tail++] = e[i].to; if( tail == T ) tail = 0;
				}
			}
		inq[now] = 0;
	}
	return d[T] != inf;
}
void mcf(){
	int x=inf;
	for( int i = from[T]; i; i = from[e[i].from] ) x=min(e[i].v,x);
    for( int i = from[T]; i; i = from[e[i].from] ){ 
		ans += x*e[i].c;
		e[i].v -= x; e[i^1].v += x;
	}
}
int main(){
	scanf("%d%d", &n, &m); S = n*2+1; T = n*2+2;
	for( int i = 1,x; i <= n; i++ ){
		scanf("%d", &x);
		insert( i, i+n, 0, 0 );
		b[i] -= x; b[i+n] += x;
	}
	insert( 0, S, m, 0 );
	for( int i = 1; i <= n; i++ ) insert( S, i, inf, 0 );
	for( int i = 1; i <= n; i++ )
		for( int j = i+1,x; j <= n; j++ ){
			scanf("%d", &x);
			if( x != -1 ) insert( i+n, j, inf, x );
		}
	for( int i = 1; i <= n*2; i++ ){
		if( b[i] > 0 ) insert( 0, i, b[i], 0 );
		if( b[i] < 0 ) insert( i, T, -b[i], 0 );
	}
	while(spfa()) mcf();
//	zkw();
	printf("%d", ans);
	return 0;
}