[CF311E]生物学家

题目描述

网络流

先想到可以把式子改写，用总资助金额减去拒绝的资助、请喝茶的费用、变性的费用。
然后变成了最小化，我们可以考虑最小割。
解决变性问题，规定靠S的这边是雄性，i是雄性s向其连vi的边表示变形就要割掉，否则i向t连vi的边。
解决条件问题，如果需求为雄性，s向i连wi(+g)，i向需要节点全部连正无穷的边，因为正无穷的边不会被割，而雌性靠T，所以要么切靠T边表示变性，否则切靠S边表示条件不满足。
需求为雌性同理。
然后跑最大流。

#include<cstdio>#include<iostream>#include<algorithm>#define fo(i,a,b) for(i=a;i<=b;i++)using namespace std;const int maxn=20000+10,maxm=200000*2+10,inf=1000000000;int now[maxn],d[maxn],h[maxn],go[maxm],dis[maxm],fx[maxm],next[maxm];int v[maxn],w[maxn],len[maxn];bool sex[maxn],nsex[maxn],bz[maxn],czy;int i,j,k,l,r,s,t,n,m,q,tot,ans;void add(int x,int y,int z,int d){    go[++tot]=y;    dis[tot]=z;    fx[tot]=tot+d;    next[tot]=h[x];    h[x]=tot;}int dfs(int x,int flow){    if (x==t){        ans+=flow;        return flow;    }    bz[x]=1;    int r=now[x],k;    while (r){        if (!bz[go[r]]&&dis[r]&&d[x]==d[go[r]]+1){            k=dfs(go[r],min(flow,dis[r]));            if (k){                dis[r]-=k;                dis[fx[r]]+=k;                now[x]=r;                return k;            }        }        r=next[r];    }    now[x]=0;    return 0;}bool change(){    int i,r,tmp=inf;    fo(i,s,t)        if (bz[i]){            r=h[i];            while (r){                if (!bz[go[r]]&&dis[r]&&d[go[r]]+1-d[i]<tmp) tmp=d[go[r]]+1-d[i];                r=next[r];            }        }    if (tmp==inf) return 0;    fo(i,s,t)        if (bz[i]) d[i]+=tmp;    return 1;}int main(){    scanf("%d%d%d",&n,&m,&q);    s=1;t=n+m+2;    fo(i,1,n) scanf("%d",&sex[i]);    fo(i,1,n) scanf("%d",&v[i]);    fo(i,1,n)        if (sex[i]){            add(s,i+1,v[i],1);            add(i+1,s,0,-1);        }        else{            add(i+1,t,v[i],1);            add(t,i+1,0,-1);        }    fo(i,1,m){        scanf("%d",&k);        nsex[i]=k;        scanf("%d%d",&w[i],&len[i]);        fo(j,1,len[i]){            scanf("%d",&k);            if (nsex[i]){                add(i+n+1,k+1,inf,1);                add(k+1,i+n+1,0,-1);            }            else{                add(k+1,i+n+1,inf,1);                add(i+n+1,k+1,0,-1);            }        }        scanf("%d",&k);        czy=k;        if (nsex[i]){            add(s,i+n+1,w[i]+czy*q,1);            add(i+n+1,s,0,-1);        }        else{            add(i+n+1,t,w[i]+czy*q,1);            add(t,i+n+1,0,-1);        }    }    ans=0;    do{        fo(i,s,t) now[i]=h[i];        fill(bz+s,bz+t+1,0);        while (dfs(s,inf)) fill(bz+s,bz+t+1,0);    }while (change());    ans=-ans;    fo(i,1,m) ans+=w[i];    printf("%d\n",ans);}