【GDSOI2018模拟4.19】排列

Description

有 n 个数 x1 ~xn 。你需要找出它们的一个排列，满足 m 个条件，每个条件形如 x_a 必须在x_b之前。在此基础上，你要最大化这个排列的最大子段和。
n<=500

Solution

神仙题
神奇的网络流建模，强行把两个最大权闭合子图搞在一起
每个点拆点，x和x’，对于每个正点，S->x和x’->T连边，边权为ai
对于每个负点，x->x’连边，边权为-ai
对于每条限制(u,v)，u->v，u’->v’连边，边权为inf
正确性可以感性理解一下

Code

#include <cstdio>
#include <cstring>
#include <algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
#define rep(i,a) for(int i=lst[a];i;i=nxt[i])
using namespace std;

const int N=1e3+5,M=1e4+5,inf=0x7fffffff;

int t[M],v[M],nxt[M],lst[N],l;
void add(int x,int y,int z) {
    t[++l]=y;v[l]=z;nxt[l]=lst[x];lst[x]=l;
    t[++l]=x;v[l]=0;nxt[l]=lst[y];lst[y]=l;
}

int n,m,S,T,x,y,ans;
int q[N],dis[N];

bool bfs() {
    fo(i,S,T) dis[i]=0;dis[S]=1;
    int i=0,j=1;q[1]=S;
    while (i<j)
        rep(k,q[++i])
            if (v[k]&&!dis[t[k]])
                dis[t[k]]=dis[q[i]]+1,q[++j]=t[k];
    return dis[T];
}

int dinic(int x,int y) {
    if (x==T) return y;
    int now=0;
    rep(i,x) 
        if (v[i]&&dis[t[i]]==dis[x]+1) {
            int k=dinic(t[i],min(y,v[i]));
            v[i]-=k;v[i^1]+=k;
            y-=k;now+=k;
            if (!y) break;
        }
    if (!now) dis[x]=-1;
    return now;
}

int main() {
    freopen("permutation.in","r",stdin);
    freopen("permutation.out","w",stdout);
    scanf("%d%d",&n,&m);l=1;S=0;T=2*n+1;
    fo(i,1,n) {
        scanf("%d",&x);
        if (x>0) add(S,i,x),add(i+n,T,x),ans+=x;
        else add(i,i+n,-x);
    }
    fo(i,1,m) {
        scanf("%d%d",&x,&y);
        add(x,y,inf);add(x+n,y+n,inf);
    }
    while (bfs()) ans-=dinic(S,inf);
    printf("%d\n",ans);
    return 0;
}