spoj 1693 Coconuts 最小割 二者取其一式

题目链接：http://www.spoj.com/problems/COCONUTS/

题意：N个城堡守卫就非洲的燕子能否搬运椰子而进行投票，每个人有自己的看法，但是为了避免跟自己的朋友持相反意见，他们时常会投相反的票。

现在给出每个人的初始看法以及朋友关系，求在某种投票方案下，违背自己意愿与持不同意见的朋友对数的总和最小。

思路：参考Edelweiss的《网络流建模汇总》。这是一个二者取其一式问题。

设置超级源点S和T。

每名守卫i如果赞成，连边s->i，容量为1;连边i->t，容量为0。

如果反对，连边i->t，容量为1;连边s-i,容量为0。

如果i和j是朋友，连边i->j和j->i，容量均为1。

然后求最小割就可以了。

#include <bits/stdc++.h>
using namespace std;
int n, m, s, t;
#define maxn 310
#define inf 0x3f3f3f3f
struct Edge
{
    int from, to, cap, flow;
    Edge(int f, int t, int c, int fl)
    {
        from = f; to = t; cap = c; flow = fl;
    }
};
vector <Edge> edges;
vector <int> G[maxn];
int cur[maxn], vis[maxn], d[maxn];
void AddEdge(int from, int to, int cap)
{
    edges.push_back(Edge(from, to, cap, 0));
    edges.push_back(Edge(to, from, 0, 0));
    m = edges.size();
    G[from].push_back(m-2);
    G[to].push_back(m-1);
}
bool bfs()
{
    memset(vis, 0, sizeof(vis));
    d[s] = 0; vis[s] = 1;
    queue <int> q;
    q.push(s);
    while(!q.empty())
    {
        int x = q.front(); q.pop();
        for(int i = 0; i < G[x].size(); i++)
        {
            Edge &e = edges[G[x][i]];
            if(!vis[e.to] && e.cap > e.flow)
            {
                d[e.to] = d[x]+1;
                vis[e.to] = 1;
                q.push(e.to);
            }
        }
    }
    return vis[t];
}
int dfs(int x, int a)
{
    if(x == t || a == 0) return a;
    int flow = 0, f;
    for(int &i = cur[x]; i < G[x].size(); i++)
    {
        Edge &e = edges[G[x][i]];
        if(d[x]+1 == d[e.to] && (f = dfs(e.to, min(e.cap - e.flow, a))) > 0)
        {
            e.flow += f;
            edges[G[x][i]^1].flow -= f;
            flow += f;
            a -= f;
            if(a == 0) break;
        }
    }
    return flow;
}
int maxflow()
{
    int flow = 0;
    while(bfs())
    {
        memset(cur, 0, sizeof(cur));
        flow += dfs(s, inf);
    }
    return flow;
}
int N, M;
int main() 
{
   // freopen("in.txt", "r", stdin);
  //  freopen("out.txt", "w", stdout);
    while(scanf("%d%d", &N, &M) && (N+M))
    {
        edges.clear();
        s = 0; t = N+1;
        for(int i = 0; i <= t; i++) G[i].clear();
        for(int i = 1; i <= N; i++)
        {
            int a; scanf("%d", &a);
            if(a)
            {
                AddEdge(s, i, 1); AddEdge(i, t, 0);
            }
            else
            {
                AddEdge(s, i, 0); AddEdge(i, t, 1);
            }
        }
        for(int i = 1; i <= M; i++)
        {
            int a, b; scanf("%d%d", &a, &b);
            AddEdge(a, b, 1); AddEdge(b, a, 1);
        }
        int flow = maxflow();
        printf("%d
", flow);
    }
    return 0;
}