【BZOJ4405】挑战NPC 带花树模板 一般图最大匹配

原创 2016年05月31日 21:10:48

这段时间被主教练找的各种poj英文题虐的头昏脑涨也没时间来更新博客QAQ

好不容易找了时间来学带花树QAQ

#include<cstdlib>
#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<queue>
#include<vector>
using namespace std;
#define MAXN 2005
#define MAXM 200005
int N,M,Q,head[MAXN],edge_ct,f[MAXN],mark[MAXN],next[MAXN],match[MAXN];
struct edge{int to,next;}e[MAXM];
void add(int x,int y){e[++edge_ct]=(edge){y,head[x]}; head[x]=edge_ct; e[++edge_ct]=(edge){x,head[y]};head[y]=edge_ct;}
int Find(int x){return f[x]==x ? x : f[x]=Find(f[x]);};
int q[MAXN],front,rear;
void _read(int &x)
{
	x=0; char ch=getchar(); bool flag=false;
	while(ch<'0' || ch>'9'){if(ch=='-')flag=true;ch=getchar(); }
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0'; ch=getchar();} if(flag)x=-x; return ;
}
int tmp[MAXN],_T;
int LCA(int u,int v) //把奇环给找出来 
{
	_T++;
	while(1)
	{
		if(u)
		{
			u=Find(u);
			if(tmp[u]==_T)return u;
			tmp[u]=_T;
			if(match[u])u=next[match[u]];
			else u=0;
		}
		swap(u,v);
	}
}
void Link(int a,int p)
{
	int b,c;
	while(a!=p)
	{
		b=match[a]; c=next[b];
		if(Find(c)!=p)next[c]=b;
		if(mark[b]==2)mark[q[rear++]=b]=1;
		if(mark[c]==2)mark[q[rear++]=c]=1;
		f[Find(b)]=Find(a);f[Find(b)]=Find(c);
		a=c;
	}
	return ;
}
void _Find(int now)
{
	for(int i=1;i<=N+3*M;i++)mark[i]=next[i]=0,f[i]=i;
	front=1;rear=1;q[rear++]=now; mark[now]=1; int u,v,p,t;
	while(front<rear)
	{
		u=q[front++];
		for(int id=head[u];id;id=e[id].next)
		{
			v=e[id].to;
			if(v==match[u] || mark[v]==2 || Find(u)==Find(v))continue;
			if(mark[v]==1)
			{
				p=LCA(u,v);
				if(Find(u)!=p)next[u]=v;
				if(Find(v)!=p)next[v]=u;
				Link(u,p); Link(v,p);
			}
			else if(!match[v])
			{
				next[v]=u;
				while(v)
				{
					u=next[v];
					t=match[u];
					match[v]=u; match[u]=v;
					v=t;
				}
				return ;
			}
			else
			{
				next[v]=u; mark[v]=2; mark[match[v]]=1; q[rear++]=match[v];
			}
		}
	}
}
void Init()
{
	_read(N); _read(M); _read(Q);
	_T=0; edge_ct=0; for(int i=1;i<=N+3*M;i++)match[i]=head[i]=tmp[i]=0;
	int x,y;
	for(int i=1;i<=Q;i++)
	{
		_read(x); _read(y); add(x,N+3*y-2);add(x,N+3*y-1);add(x,N+3*y-0);
	}
	for(int i=1;i<=M;i++)add(N+3*i-2,N+3*i-1),add(N+3*i-1,N+3*i-0),add(N+3*i-2,N+3*i-0);
	return ;
}
void work()
{
	int ans=0;
	for(int i=1;i<=N+3*M;i++)if(!match[i])
	_Find(i);
	for(int i=1;i<=N+3*M;i++)if(match[i])ans++;
	ans=ans/2; ans=ans-N; printf("%d\n",ans); 
//	for(int i=1;i<=N;i++)printf("%d ",(match[i]-N-1)/3+1);
//	putchar('\n');
	return ;
}
int main()
{
//	freopen("in.txt","r",stdin);
	int T; _read(T);
	while(T--)
	{
		Init();
		work();
	}
	return 0;
}


版权声明:本文为博主原创文章,未经博主允许不得转载。

相关文章推荐

UOJ171. 挑战NPC【WC2016】【带花树】【一般图匹配】

UOJ#171. 挑战NPC【WC2016】 一般图最大匹配 题目飞机票 带花树(一般最大匹配)讲解 没想到建图居然建一般图,考场上想了很久的二分图最大匹配,大样例都过不了 (๑•́ ...

uoj171 bzoj 4405: [wc2016]挑战NPC 一般图最大匹配

不得不说建图太巧妙了,我在知道是一般图最大匹配的情况下都不会建图。。。        首先每个筐放三个,这个“三”是很关键的。        把一个筐拆成三个点,然后每一个球向可以放的三个筐都连边;同...

一般图最大匹配-带花树-带注释模板

一般图最大匹配 带花树算法真·存模板系列

[模板]带花树算法(一般图最大匹配)

->什么是增广路增广路:长度为奇数,第一条和最后一条边都是未匹配的边,增光路径是一条匹配边,一条未匹配边这样交错的路径,也叫做交错路 ->带花树的做法像匈牙利算法那样不断枚举点寻找增广路, 当找到环时...
  • Frods
  • Frods
  • 2017年01月15日 20:27
  • 255

uoj79 一般图最大匹配&【带花树】初学

上次听毛主力讲完课忘记写笔记了o(╯□╰)o。。发现现在忘得差不多了。于是补一发。        对于一般图 ,最大匹配,采用和二分图一样的思路。就是从每一个点出发,看能不能找到一条新的增广路,如果有...

一般图最大匹配问题-带花树开花算法

以前用这个算法写过一两个水题,当时纯粹是套用模板,对算法本身是一知半解。然后Watashi的多校题中有个带花树模板题,现成的模板都套出了各种死循环,RE问题,弱爆了。这两天重新看了看论文和博客,重新理...

【UOJ#97】一般图最大匹配【带花树】

模板题咯。 match[]记录匹配结果。fri[]记录增广时的匹配对象。fa[]为并查集的父亲数组。top[]记录路径,为了方便找lca。ring[]记录当前点所在的环的lca。 odd[...

带花树(一般图最大匹配)详解 ZOJ 3316

{CSDN:CODE:3423010}

一般图最大匹配问题-带花树开花算法

以前用这个算法写过一两个水题,当时纯粹是套用模板,对算法本身是一知半解。然后Watashi的多校题中有个带花树模板题,现成的模板都套出了各种死循环,RE问题,弱爆了。这两天重新看了看论文和博客,重新理...

【UOJ#171】【WC2016】挑战NPC【带花树】

vfk出的(好)题。赛时无人AC,一大堆人60分... 前两个点,暴搜即可。 第三个点,贪心。 第四到六个点,最大流。 std:带花树 把一个篮子拆为3个槽,每两个槽之间互相...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:【BZOJ4405】挑战NPC 带花树模板 一般图最大匹配
举报原因:
原因补充:

(最多只允许输入30个字)