POJ 3469 Dual Core CPU（最小割）

题意：

给你n个任务，每个任务要不在花费 $a_i$ 在a核上完成，要不花费 $b_i$ 在b核上完成。同时有m个约束条件，如果第i个任务和第j个任务在不同的核上完成需要额外的花费。

源点s，汇点t。
根据题意建图对于每个任务，s连i，容量为 $a_i$ 。
i连t，容量为 $b_i$ 。根据最小割的定义我们可以得到两个集合，一个以s为核心，一个以t为核心。现在考虑我们建的图，每个任务i的两条边中短的那条一定是最小割中的边。由此我们求出最小割即是答案。

然后我们考虑约束关系i，j，v。如果任务i和j在不同的集合中。那么建边 <i,j> <script type="math/tex" id="MathJax-Element-833"> </script> 和 <j,i> <script type="math/tex" id="MathJax-Element-834"> </script>，容量为v。
那么边 <i,j> <script type="math/tex" id="MathJax-Element-835"> </script> 一定是最小割中的边，所以一定会算进最小割中。至此本题考虑完全，直接上isap就可以了。

#include<iostream>
#include<cstdio>
#include<cstring>
#define inf 0x3f3f3f3f
typedef long long int lli;
using namespace std;
const lli maxn = 102000;
struct edge{
    lli to,v,next;
}ed[maxn+4*200000];
lli d[maxn],cur[maxn],pre[maxn],gap[maxn],q[maxn+2000],cnte,head[maxn];
void ae(lli x, lli y, lli v) {
    ed[cnte].to = y;
    ed[cnte].v = v;
    ed[cnte].next = head[x];
    head[x] = cnte++;
}
void rbfs (lli s,lli t) {
    lli fi,se;
    memset(gap,0,sizeof(gap));
    memset(d,-1,sizeof(d));
    d[t] = 0;
    gap[0] = 1;
    fi = se = 0;
    q[se++] = t;
    while (fi != se) {
        lli u = q[fi++];
        for (lli i=head[u];~i;i=ed[i].next) {
            lli v = ed[i].to;
            if (~d[v]) continue;
            d[v] = d[u] + 1;
            q[se++] = v;
            gap[d[v]]++;
        }
    }
}
lli isap(lli s,lli t){
    memcpy(cur,head,sizeof(head));
    rbfs (s,t);
    lli flow = 0, u = pre[s]=s,i;
    while(d[t] < t+1) {
        if(u==t) {
            lli f = inf,neck;
            for(i= s;i != t;i = ed[cur[i]].to){
                if(f > ed[cur[i]].v){
                    f = ed[cur[i]].v;
                    neck = i;
                }
            }
            for(i = s;i != t;i = ed[cur[i]].to){
                ed[cur[i]].v -= f;
                ed[cur[i]^1].v += f;
            }
            flow += f;
            u = neck;
        }
        for(i = cur[u];~i;i=ed[i].next) if(d[ed[i].to]+1 == d[u] && ed[i].v) break;
        if(~i) {
            cur[u] = i;
            pre[ed[i].to] = u;
            u = ed[i].to;
        }
        else{//否则回退，重新找增广路
            if(0 == (--gap[d[u]])) break;//GAP间隙优化，如果出现断层，可以知道一定不会再有增广路了
            lli mind = t+1;
            for(i = head[u];~i;i = ed[i].next){
                if(ed[i].v && mind > d[ed[i].to]){
                    cur[u] = i;
                    mind = d[ed[i].to];
                }
            }
            d[u] = mind + 1;
            gap[d[u]]++;
            u = pre[u];
        }
    }
    return flow;
}
void ini(){
    memset(head,-1,sizeof(head));cnte = 0;
}
int main(){
    int n,m,v,a,b;
    scanf("%d%d",&n,&m);
    ini();
    int s = 0,t = n+1;
    memset(head,-1,sizeof(head));cnte = 0;
    for(int i = 1;i <= n;i++){
        scanf("%d%d",&a,&b);
        ae(s,i,a);ae(i,s,0);
        ae(i,t,b);ae(t,i,0);
    }
    for(int i = 1;i <= m;i++){
        scanf("%d%d%d",&a,&b,&v);
        ae(a,b,v);ae(b,a,v);
    }
    printf("%lld\n",isap(s,t));
    return 0;
}