SPOJ_1693 COCONUTS
对于每个人的选择一共有两种，根据他们最后的选择就可以将这些人分成两个集合，显然集合间如果有朋友形成的边的话是要删掉的，这样实际上就形成了一个割，对应的持不同意见的朋友关系的数量。为了体现出改变自己的意见产生的代价，可以将S和所有起始选择1的人连一条容量为1的边，再将起始选择0的人和T连一条容量为1的边，朋友关系的边都是容量为1的无向边。这样对这个图求最小割就可以得到最后的结果。
#include#include #include #define MAXD 310#define MAXM 90610#define INF 0x3f3f3f3fint N, M, first[MAXD], e, next[MAXM], v[MAXM], flow[MAXM];int S, T, d[MAXD], q[MAXD], work[MAXD];void add(int x, int y, int z){ v[e] = y, flow[e] = z; next[e] = first[x], first[x] = e ++; }void init(){ int i, x, y; S = 0, T = N + 1; memset(first, -1, sizeof(first[0]) * (T + 1)); e = 0; for(i = 1; i <= N; i ++) { scanf("%d", &x); if(x) add(S, i, 1), add(i, S, 0); else add(i, T, 1), add(T, i, 0); } for(i = 0; i < M; i ++) { scanf("%d%d", &x, &y); add(x, y, 1), add(y, x, 1); }}int bfs(){ int i, j, rear = 0; memset(d, -1, sizeof(d[0]) * (T + 1)); d[S] = 0, q[rear ++] = S; for(i = 0; i < rear; i ++) for(j = first[q[i]]; j != -1; j = next[j]) if(flow[j] && d[v[j]] == -1) { d[v[j]] = d[q[i]] + 1, q[rear ++] = v[j]; if(v[j] == T) return 1; } return 0;}int dfs(int cur, int a){ if(cur == T) return a; for(int &i = work[cur]; i != -1; i = next[i]) if(flow[i] && d[v[i]] == d[cur] + 1) if(int t = dfs(v[i], std::min(a, flow[i]))) { flow[i] -= t, flow[i ^ 1] += t; return t; } return 0;}int dinic(){ int ans = 0, t; while(bfs()) { memcpy(work, first, sizeof(first[0]) * (N + 1)); while(t = dfs(S, INF)) ans += t; } return ans;}void solve(){ printf("%d\n", dinic()); }int main(){ while(scanf("%d%d", &N, &M), N || M) { init(); solve(); } return 0; }