POJ 3469 最小割

题意:

有一些模块（modules）和一个双核处理器，一个模块可以在任意一个核上处理，每个核对应每个模块有个开销。现在有一些模块间需要数据交换，如果需要数据交换的模块在一个核上处理，则不需要额外开销，否则需要加上一个开销。现在需要完成所有模块，问最小需要多少开销。

如果没有这个额外的开销，那么每个模块只要选择开销小的那个核就行了。额外的开销给选择加上了限制。

题解：

最小割模型。

S和T分别表示两个核，每个模块和S，T连边求最小割就是了~

 

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>

#define N 201000
#define M  2001000
#define INF 100000000

using namespace std;

int head[N],next[M],to[M],len[M],cnt,S,T,a[N],b[N],aa[M],bb[M],n,m,w[N],layer[N],q[M<<4]; 

inline void add(int u,int v,int w)
{
    to[cnt]=v; len[cnt]=w; next[cnt]=head[u]; head[u]=cnt++;
    to[cnt]=u; len[cnt]=0; next[cnt]=head[v]; head[v]=cnt++;
}

void read()
{
    memset(head,-1,sizeof head);
    for(int i=1;i<=n;i++) scanf("%d%d",&a[i],&b[i]);
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d%d",&aa[i],&bb[i],&w[i]);
        add(aa[i],bb[i],w[i]); add(bb[i],aa[i],w[i]);
    }
    S=0; T=n+1;
    for(int i=1;i<=n;i++) add(S,i,a[i]),add(i,T,b[i]);
}

bool bfs()
{
    memset(layer,-1,sizeof layer);
    int h=1,t=2,sta;
    q[1]=S; layer[S]=0;
    while(h<t)
    {
        sta=q[h++];
        for(int i=head[sta];~i;i=next[i])
            if(len[i]>0&&layer[to[i]]<0)
            {
                layer[to[i]]=layer[sta]+1;
                q[t++]=to[i];
            }
    }
    return layer[T]!=-1;
}

int find(int u,int cur_flow)
{
    if(u==T) return cur_flow;
    int result=0,tmp=0;
    for(int i=head[u];~i&&result<cur_flow;i=next[i])
        if(len[i]>0&&layer[to[i]]==layer[u]+1)
        {
            tmp=find(to[i],min(len[i],cur_flow-result));
            len[i]-=tmp; len[i^1]+=tmp; result+=tmp;
        }
    if(!result) layer[u]=-1;
    return result; 
}

int dinic()
{
    int ans=0;
    while(bfs()) ans+=find(S,INF);
    return ans;
}

int main()
{
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        read();
        printf("%d",dinic());
    }
    return 0;
}