BZOJ 1934 善意的投票（最小割）

最新推荐文章于 2019-08-09 15:14:00 发布

蘑菇小哲

最新推荐文章于 2019-08-09 15:14:00 发布

阅读量589

点赞数

分类专栏：网络流文章标签： ACM ACM 图论

本文链接：https://blog.csdn.net/moguxiaozhe/article/details/38276603

版权

网络流专栏收录该内容

27 篇文章 0 订阅

订阅专栏

建立一个源点S和一个汇点T，同意的向S连边，不同意的向T连边，每对朋友之间互相连边，然后跑一遍最大流算法。

#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <algorithm>
#include <queue>
#include <vector>
#include <cmath>
#define LL long long
using namespace std;
const int maxn = 1000 + 10;
const int INF = 10000000;
struct Edge
{
	int from,to,cap,flow;
	Edge(int u,int v,int c,int f) : from(u),to(v),cap(c),flow(f) { }
};
int N , M;
vector<Edge>edges;
vector<int>G[maxn];
int a[maxn];
int p[maxn];
void init(int n)
{
	for(int i=0;i<=n;i++)
		G[i].clear();
	edges.clear();
}
void AddEdge(int from,int to,int cap)
{
	edges.push_back(Edge(from,to,cap,0));
	edges.push_back(Edge(to,from,0,0));
	int m = edges.size();
	G[from].push_back(m-2);
	G[to].push_back(m-1);
}
int Maxflow(int s,int t)
{
	int flow = 0;
	for(;;)
	{
		memset(a,0,sizeof(a));
		queue<int>Q;
		Q.push(s);
		a[s] = INF;
		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(!a[e.to] && e.cap > e.flow)
				{
					p[e.to] = G[x][i];
					a[e.to] = min(a[x],e.cap-e.flow);
					Q.push(e.to);
				}
			}
			if(a[t]) break;
		}
		if(!a[t]) break;
		for(int u=t;u!=s;u = edges[p[u]].from)
		{
			edges[p[u]].flow += a[t];
			edges[p[u]^1].flow -= a[t];
		}
		flow += a[t];
	}
	return flow;
}
int main()
{
	while(scanf("%d%d",&N,&M)!=EOF)
	{
		int tmp;
		for(int i=1;i<=N;i++)
		{
			scanf("%d",&tmp);
			if(tmp == 1) AddEdge(0,i,1);
			else AddEdge(i,N+1,1);
		}
		int u , v;
		while(M--)
		{
			scanf("%d%d",&u,&v);
			AddEdge(u,v,1);
			AddEdge(v,u,1);
		}
		int ans = Maxflow(0,N+1);
		printf("%d\n",ans);
	}
	return 0;
}