POJ3469

题目链接：http://poj.org/problem?id=3469

题目大意：

要在核A和核B上运行N个模块，模块i在核A上执行的花费为Ai，在核B上执行的花费为Bi。有M个互相之间需要进行数据交换的模块组合(ai,bi),如果这两个模块在同一个核上执行则没有额外花费，否则会产生wi的花费，请计算执行所有模所需的最小花费

思路：

转化成最小割问题，最小花费就是固定花费加上额外花费。建立一个源点s（模块在A核上执行的集合）一个汇点t（模块在B核上执行的集合），将模块作为顶点，边的容量就是执行的花费，先将普通花费建立，之后再建立额外花费（双向边）

代码：

 

<span style="font-size:18px;">#include<iostream>
#include<cstdio>
#include<vector>
#include<queue>
#include<algorithm>
#include<cstring>
using namespace std;
const int maxn = 20010;
const int maxm = 200010;
const int inf = 0x3f3f3f3f;
struct edge
{
    int to, cap, rev;
};
vector<edge>G[maxm];
int level[maxm];
int iter[maxm];
void add_edge(int from, int to, int cap)
{
    G[from].push_back((edge){to, cap, G[to].size()});
    G[to].push_back((edge){from, 0, G[from].size() - 1});
}
void bfs(int s)
{
    memset(level, -1, sizeof(level));
    queue<int> que;
    level[s] = 0;
    que.push(s);
    while(!que.empty())
    {
        int v = que.front();
        que.pop();
        for(int i = 0; i < G[v].size(); i++)
        {
            edge &e = G[v][i];
            if(e.cap > 0 && level[e.to] < 0)
            {
                level[e.to] = level[v] + 1;
                que.push(e.to);
            }
        }
    }
}
int dfs(int u, int t, int f)
{
    if(u == t)
        return f;
    for(int &i = iter[u]; i < G[u].size(); i++)
    {
        edge &e = G[u][i];
        if( e.cap > 0 && level[u] <level[e.to])
        {
            int d = dfs(e.to, t, min(f, e.cap));
            if(d > 0)
            {
                e.cap -= d;
                G[e.to][e.rev].cap += d;
                return d;
            }
        }
    }
    return 0;
}
int max_flow(int s, int t)
{
    int flow = 0;
    for(;;)
    {
        bfs(s);
        if(level[t] < 0)
            return flow;
        memset(iter, 0, sizeof(iter));
        int f;
        while((f = dfs(s, t, inf))>0)
            flow += f;
    }
}

int main()
{
    int n,m;

    int a, b, ai, bi, wi;
    scanf("%d%d", &n, &m);
    int s = n, t = s + 1;
    for(int i = 0; i< n; i++)
    {
        scanf("%d%d",&a, &b);
        add_edge(s, i, a);
        add_edge(i, t, b);
    }
    for(int i = 0; i < m; i++)
    {
        scanf("%d%d%d", &ai, &bi, &wi);
        add_edge(ai-1, bi-1, wi);
        add_edge(bi-1, ai-1, wi);
    }
    printf("%d\n", max_flow(s, t));
    return 0;
}
</span>