无源汇有上下界最小费用流
(等我先再刷会儿网络流再去做鬼畜数据结构题)
比较果的吧,然而年代太久远几乎忘记了。。
大约是把每个点的度算出来
(度就是进来的下界的总减出去的下界的总流量)
然后各边的值就变为上界减下界
呣唔(本人专用叹词)
拆点,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;
}