【集训队互测2013】供电网络

最新推荐文章于 2021-01-18 22:10:06 发布

alan_cty

最新推荐文章于 2021-01-18 22:10:06 发布

阅读量1.3k

点赞数

分类专栏：上下界网络流费用流文章标签：集训队互测 2013 供电网络上下界网络流费用流

本文链接：https://blog.csdn.net/alan_cty/article/details/50849394

版权

费用流同时被 2 个专栏收录

6 篇文章 0 订阅

订阅专栏

上下界网络流

3 篇文章 0 订阅

订阅专栏

Description

给出一个n个点，m条边的有向图，每个点初始的值有正有负。你可以在每一个点加上一个数或是减去一个数，代价为你所选择的数*你所选择的点的加或减的费用（加减费用可以不同）。你也可以通过一条边把一定大小的数从起点“运”到终点，设为所运输数为x，则代价为 $a*x^2+b*x$ 。a，b为这条边的基础属性。且每条边有运输上界r和下界l，求把所有点的数都变为0的最小代价。
n <= 200, m <= 600, 1 <= a, b <=3, 1 <= l<= 10, 1 <= r<=100。

Solution

很明显的网络流模型。从原点向每个点连上界为∞，下界为0，费用为它加一个数的代价的边。表示你可以任意的加上一个数。从每个点向汇点连上界为∞，下界为0，费用为它减一个数的代价的边。表示你可以任意的减少一个数。然后对于每一个正权点，从原点向它连上界和下界都为它的权值，费用为0的边，表示你这个点之前就有这么大的数。对于每一个负权点，从它向汇点连上界和下界都为它的权值的相反数，费用为0的边，表示你还需要这么大的数。对于每一条边，由于它是二次的，于是我们可以把它拆成一堆边，每条边的上界为1，下界为0。即 $a+b,3*a+b,5*a+b...$ 这样就能表示所有情况了。最后再跑一遍上下界最小费用可行流就行了。注意，由于这样拆边会导致边数过大，于是我们可以使用动态加边大法。

Code

#include<cstdio>
#include<cstring>
#include<algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define rep(i,a) for(int i=last[a];i;i=next[i])
#define inf 1000000
#define ll long long
#define N 205
#define M 150500
using namespace std;
ll ans;
int n,m,x,y,z,j,k,l,S,T,SS,TT,used[N],flag[N][N];
int t[M],f[M],v[M],next[M],last[N],dis[N],bz[N],d[M],p[N];
void add(int x,int y,int z,int k) {
    t[++l]=y;f[l]=z;v[l]=k;next[l]=last[x];last[x]=l;
    t[++l]=x;v[l]=-k;next[l]=last[y];last[y]=l;
}
void link(int x,int y,int low,int up,int cost) {
    add(x,y,up-low,cost);
    if (low) used[y]+=low,used[x]-=low,ans+=cost*low;
}
struct note{
    int x,y,a,b,l,r;
}a[605];
bool spfa() {
    fo(i,1,m)
        if (flag[a[i].x][a[i].y]==a[i].l&&a[i].l<a[i].r) {
            a[i].l++;add(a[i].x,a[i].y,1,
            a[i].a*(2*a[i].l-1)+a[i].b);
        }
    memset(dis,127,sizeof(dis));int mx=dis[SS];dis[SS]=0;
    memset(bz,0,sizeof(bz));bz[SS]=1;
    int i=0,j=1;d[1]=SS;
    while (i<j) {
        rep(k,d[++i])
            if (f[k]&&dis[t[k]]>dis[d[i]]+v[k]) {
                dis[t[k]]=dis[d[i]]+v[k];p[t[k]]=k;
                if (!bz[t[k]]) bz[t[k]]=1,d[++j]=t[k];
            }
        bz[d[i]]=0;
    }
    return dis[TT]!=mx;
}
void find() {
    int i=TT;
    while (i!=SS) {
        x=t[p[i]^1];y=t[p[i]];
        f[p[i]]--;f[p[i]^1]++;
        flag[x][y]++;flag[y][x]--;
        ans+=v[p[i]];i=x;
    }
}
int main() {
    scanf("%d%d",&n,&m);S=0,T=n+1;SS=n+2;TT=n+3;l=1;
    fo(i,1,n) {
        scanf("%d%d%d",&x,&y,&z);
        if (x>0) link(S,i,x,x,0);else link(i,T,-x,-x,0);
        link(S,i,0,inf,y);link(i,T,0,inf,z);
    }
    fo(i,1,m) {
        scanf("%d%d%d%d%d%d",&a[i].x,&a[i].y,
        &a[i].a,&a[i].b,&a[i].l,&a[i].r);
        link(a[i].x,a[i].y,a[i].l,a[i].l,0);
        ans+=a[i].l*a[i].l*a[i].a+a[i].l*a[i].b;
        flag[a[i].x][a[i].y]+=a[i].l;
        flag[a[i].y][a[i].x]-=a[i].l;
    }
    link(T,S,0,inf,0);
    fo(i,S,T) if (used[i]>0) link(SS,i,0,used[i],0);
    else link(i,TT,0,-used[i],0);
    while (spfa()) find();
    printf("%lld",ans);
}