BZOJ_2039_[2009国家集训队]employ人员雇佣_ 最小割

BZOJ_2039_[2009国家集训队]employ人员雇佣_ 最小割

Description

作为一个富有经营头脑的富翁，小L决定从本国最优秀的经理中雇佣一些来经营自己的公司。这些经理相互之间合作有一个贡献指数,（我们用Ei,j表示i经理对j经理的了解程度），即当经理i和经理j同时被雇佣时，经理i会对经理j做出贡献，使得所赚得的利润增加Ei,j。当然，雇佣每一个经理都需要花费一定的金钱Ai，对于一些经理可能他做出的贡献不值得他的花费，那么作为一个聪明的人，小L当然不会雇佣他。 然而，那些没有被雇佣的人会被竞争对手所雇佣，这个时候那些人会对你雇佣的经理的工作造成影响，使得所赚得的利润减少Ei,j（注意：这里的Ei,j与上面的Ei,j 是同一个）。 作为一个效率优先的人，小L想雇佣一些人使得净利润最大。你可以帮助小L解决这个问题吗？

Input

第一行有一个整数N<=1000表示经理的个数 第二行有N个整数Ai表示雇佣每个经理需要花费的金钱 接下来的N行中一行包含N个数，表示Ei,j，即经理i对经理j的了解程度。（输入满足Ei,j=Ej,i）

Output

第一行包含一个整数，即所求出的最大值。

Sample Input

3
3 5 100
0 6 1
6 0 2
1 2 0

Sample Output

1
【数据规模和约定】
20%的数据中N<=10
50%的数据中N<=100
100%的数据中 N<=1000， Ei,j<=maxlongint, Ai<=maxlongint

---

转化为损失最少的最小割模型。

建立$S$，$T$分别表示选或不选，便于用$S$，$T$割集表示方案。

然后设两个经理$i$和$j$，连这样几条边$S->i(x1),S->j(x2),i->j(x3),i->T(x4),j->T(x5)$，括号内代表容量，$x$为设的未知数。

假设$i$和$j$都选，那么要割掉到$T$的两条边$x4+x5$。假设都不选，则需要割掉$S$连出的两条边$x1+x2$。

$i$选$j$不选，则需要割掉$x2+x3+x4$三条边。

然后分别和真正的收益对应一下，可以得到$x1=\sum\limits_{k=1}^{n}e[i][k]$ $x3=2*e[i][j]$ $x4=a[i]$。

于是建图跑最小割即可。

代码：

#include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std;
typedef long long ll;
inline char nc() {
    static char buf[100000],*p1,*p2;
    return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline int rd() {
    register int x=0;
    register char s=nc();
    while(s<'0'||s>'9') s=nc();
    while(s>='0'&&s<='9') x=(x<<3)+(x<<1)+s-'0',s=nc();
    return x;
}
#define N 1050
#define M 3000050
#define S (n+1)
#define T (n+2)
#define inf (1ll<<60)
int n,a[N],e[N][N];
struct Dinic {
    int head[N],to[M],nxt[M],cnt;
    int dep[N],Q[N],l,r;
    ll flow[M],sum;
    inline void add(int u,int v,ll f) {
        to[++cnt]=v; nxt[cnt]=head[u]; head[u]=cnt; flow[cnt]=f;
        to[++cnt]=u; nxt[cnt]=head[v]; head[v]=cnt; flow[cnt]=0;
    }
    bool bfs() {
        memset(dep,0,sizeof(dep));
        l=r=0; Q[r++]=S; dep[S]=1;
        while(l<r) {
            int x=Q[l++],i;
            for(i=head[x];i;i=nxt[i]) {
                if(!dep[to[i]]&&flow[i]) {
                    dep[to[i]]=dep[x]+1;
                    if(to[i]==T) return 1;
                    Q[r++]=to[i];
                }
            }
        }
        return 0;
    }
    ll dfs(int x,ll mf) {
        if(x==T) return mf;
        ll nf=0;
        int i;
        for(i=head[x];i;i=nxt[i]) {
            if(dep[to[i]]==dep[x]+1&&flow[i]) {
                ll tmp=dfs(to[i],min(mf-nf,flow[i]));
                if(!tmp) dep[to[i]]=0;
                nf+=tmp;
                flow[i]-=tmp;
                flow[i^1]+=tmp;
                if(nf==mf) break;
            }
        }
        return nf;
    }
    void dinic() {
        ll f;
        while(bfs()) while((f=dfs(S,inf))!=0) sum-=f;
        printf("%lld\n",sum);
    }
}G;
int main() {
    G.cnt=1;
    n=rd();
    int i,j;
    for(i=1;i<=n;i++) a[i]=rd();
    for(i=1;i<=n;i++) {
        ll now=0;
        for(j=1;j<=n;j++) {
            e[i][j]=rd();
            G.sum+=e[i][j];
            now+=e[i][j];
            G.add(i,j,2*e[i][j]);
        }
        G.add(S,i,now);
        G.add(i,T,a[i]);
    }
    G.dinic();
}