POJ-3041-Asteroids

最新推荐文章于 2019-07-12 19:50:31 发布

南宮逸辰

最新推荐文章于 2019-07-12 19:50:31 发布

阅读量478

点赞数

分类专栏： ACM POJ 文章标签：二分图匹配 Dinic

本文链接：https://blog.csdn.net/z309241990/article/details/22940329

版权

ACM 同时被 2 个专栏收录

716 篇文章 0 订阅

订阅专栏

POJ

378 篇文章 0 订阅

订阅专栏

二分图匹配一道~用的网络流解。

代码：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<queue>
using namespace std;
const int inf=1<<29;
const int maxn=40010;
int n,k,e,st,des,pnt[maxn],nxt[maxn],head[maxn],flow[maxn],level[maxn];
queue<int> q;
void addEdge(int u,int v,int f)
{
    pnt[e]=v;nxt[e]=head[u];flow[e]=f;head[u]=e++;
    pnt[e]=u;nxt[e]=head[v];flow[e]=0;head[v]=e++;
}
bool BFS()
{
    memset(level,0,sizeof(level));
    while(!q.empty())
	q.pop();
    q.push(st);
    level[st]=1;
    while(!q.empty())
    {
	int u=q.front();
	q.pop();
	if(u==des)
	    return true;
	for(int i=head[u];i!=-1;i=nxt[i])
	    if(flow[i]&&!level[pnt[i]])
	    {
		level[pnt[i]]=level[u]+1;
		q.push(pnt[i]);
	    }
    }
    return level[des];
}
int DFS(int u,int sum)
{
    if(u==des)
	return sum;
    for(int i=head[u],t;i!=-1;i=nxt[i])
    {
	if(flow[i]&&level[pnt[i]]==level[u]+1&&(t=DFS(pnt[i],min(sum,flow[i]))))
	{
	    flow[i]-=t;
	    flow[i^1]+=t;
	    return t;
	}
    }
    return level[u]=0;
}
int maxflow()
{
    int ans=0;
    while(BFS())
	while(1)
	{
	    int val=DFS(st,inf);
	    if(!val)
		break;
	    ans+=val;
	}
    return ans;
}
int main()
{
    while(scanf("%d%d",&n,&k)!=EOF)
    {
	e=st=0;
	des=n*2+1;
	memset(head,-1,sizeof(head));
	for(int i=0;i<k;i++)
	{
	    int u,v;
	    scanf("%d%d",&u,&v);
	    addEdge(u,n+v,1);
	}
	for(int i=1;i<=n;i++)
	    addEdge(0,i,1);
	for(int i=1;i<=n;i++)
	    addEdge(n+i,des,1);
	printf("%d\n",maxflow());
    }
    return 0;
}