BZOJ 2055:80人环游世界

无源汇有上下界最小费用流

（等我先再刷会儿网络流再去做鬼畜数据结构题）

比较果的吧，然而年代太久远几乎忘记了。。

大约是把每个点的度算出来

（度就是进来的下界的总减出去的下界的总流量）

然后各边的值就变为上界减下界

呣唔（本人专用叹词）

拆点，x1,x2之间流量费用为0，上界下界都为这个点的魅力值

其他边上界INF下界0，费用为这个城市到另一城市的费用

然后度大于0的由炒鸡源点S连一条值为du[i]的边，费用为0

度小于0的向炒鸡汇点T连一条流量为-du[i]的边，费用为0

因为总流量有限制所以加一条源点st向超级源点的边，流量为m，费用为0

对st与炒鸡汇点T跑一遍最小费用最大流

由于绝壁有解所以直接输出ans就好了呣唔。

代码：

#include<set>
#include<map>
#include<queue>
#include<stack>
#include<cmath>
#include<vector>
#include<cstdio>
#include<complex>
#include<cstring>
#include<iostream>
#include<algorithm>
#define g getchar()
#define S 201
#define T 202
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
inline ll read(){
    ll x=0,f=1;char ch=g;
    for(;ch<'0'||ch>'9';ch=g)if(ch=='-')f=-1;
    for(;ch>='0'&&ch<='9';ch=g)x=x*10+ch-'0';
    return x*f;
}
inline void out(ll x){
    int a[25],t=0;
    if(x<0)putchar('-'),x=-x;
    for(;x;x/=10)a[++t]=x%10;
    for(int i=t;i;--i)putchar('0'+a[i]);
    if(t==0)putchar('0');
    putchar('\n');
}
int n,m,ans,e=1;
struct re{int from,v,fl,w,next;}ed[100000];
int dis[600],dui[10000],head[205],from[205],du[205];
bool pd[600];
inline void ins(int u,int v,int fl,int w){
    ed[++e]=(re){u,v,fl,w,head[u]};head[u]=e;
    ed[++e]=(re){v,u,0,-w,head[v]};head[v]=e;
}
bool spfa(){
    int tou=1,wei=1;
    memset(dis,inf,sizeof(dis));
    memset(pd,0,sizeof(pd));
    dui[1]=0;pd[0]=1;dis[0]=0;
    for(;tou<=wei;pd[dui[tou++]]=0){
        int u=dui[tou];
        for(int i=head[u];i;i=ed[i].next)
        if(ed[i].fl&&ed[i].w+dis[u]<dis[ed[i].v]){
            dis[ed[i].v]=ed[i].w+dis[u];
            from[ed[i].v]=i;
            if(!pd[ed[i].v])pd[ed[i].v]=1,dui[++wei]=ed[i].v;
        }
    }
    return dis[T]!=inf;
}
int dfs(){
    int x=inf,tmp=0;
    for(int i=from[T];i;i=from[ed[i].from])x=min(x,ed[i].fl);
    for(int i=from[T];i;i=from[ed[i].from]){tmp+=x*ed[i].w;ed[i].fl-=x;ed[i^1].fl+=x;}
    return tmp;
}
int main(){
//  freopen("1.in","r",stdin);
    n=read();m=read();
    for(int i=1;i<=n;++i){
        int x=read();
        du[i]-=x;du[i+n]+=x;
    }
    ins(0,S,m,0);
    for(int i=1;i<=n;++i)ins(S,i,inf,0);
    for(int i=1;i<n;++i){
        for(int j=i+1;j<=n;++j){
            int x=read();
            if(x!=-1)ins(i+n,j,inf,x);
        }
    }
    for(int i=1;i<=2*n;++i){
        if(du[i]>0)ins(0,i,du[i],0);
        else ins(i,T,-du[i],0);
    }
    ans=0;
    while(spfa())ans+=dfs();
    out(ans);
    return 0;
}