HDU 3081 Marriage Match II

2011-05-14

18:08：47

第一次用SAP记录备忘。

http://acm.hdu.edu.cn/showproblem.php?pid=3081

Marriage Match II

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 540    Accepted Submission(s): 152

Problem Description

Presumably, you all have known the question of stable marriage match. A girl will choose a boy; it is similar as the game of playing house we used to play when we are kids. What a happy time as so many friends playing together. And it is normal that a fight or a quarrel breaks out, but we will still play together after that, because we are kids.
Now, there are 2n kids, n boys numbered from 1 to n, and n girls numbered from 1 to n. you know, ladies first. So, every girl can choose a boy first, with whom she has not quarreled, to make up a family. Besides, the girl X can also choose boy Z to be her boyfriend when her friend, girl Y has not quarreled with him. Furthermore, the friendship is mutual, which means a and c are friends provided that a and b are friends and b and c are friend.
Once every girl finds their boyfriends they will start a new round of this game—marriage match. At the end of each round, every girl will start to find a new boyfriend, who she has not chosen before. So the game goes on and on.
Now, here is the question for you, how many rounds can these 2n kids totally play this game?

 

 

Input

There are several test cases. First is a integer T, means the number of test cases.
Each test case starts with three integer n, m and f in a line (3<=n<=100,0<m<n*n,0<=f<n). n means there are 2*n children, n girls(number from 1 to n) and n boys(number from 1 to n).
Then m lines follow. Each line contains two numbers a and b, means girl a and boy b had never quarreled with each other.
Then f lines follow. Each line contains two numbers c and d, means girl c and girl d are good friends.

 

 

Output

For each case, output a number in one line. The maximal number of Marriage Match the children can play.

 

 

Sample Input

  
  

   
   1
4 5 2
1 1
2 3
3 2
4 2
4 4
1 4
2 3
  
  

 

 

Sample Output

  
  

   
   2
  
  
  
  

   
   题目大意：
  
  
  
  

   
   n个男孩n个女孩，女孩选男孩，每个女孩都要选到不同的人
  
  
  
  

   
   k对女孩有相同选择标准，
  
  
  
  

   
   女孩每轮都选择没选过的男孩，
  
  
  
  

   
   问总共能选几轮。
  
  
  
  

   
   题解：
  
  
  
  

   
   女孩1..n，男孩n+1..2*n编号
  
  
  
  

   
   由女孩到男孩建容量为1的边
  
  
  
  

   
   起点st=2*n+1，到1..n建边；
  
  
  
  

   
   n+1..2*n到终点ed=2*n+2建边
  
  
  
  

   
   二分搜索最大容量即为答案；
  
  
  
  

   
   #include <iostream>
using namespace std;
const int inf=0x3f3f3f3f;
struct rec
{
	int u,v,w,link;
}edge[40005],E[40005];
struct re
{
	int u,v;
}p[20005];
int head[405],cur[405],f[405],dis[405],gap[405],pre[405];
int n,m,NE,st,ed;
bool map[105][105];
inline void addedge(int u,int v,int w)
{
	edge[NE].u=u;
	edge[NE].v=v;
	edge[NE].w=w;
	edge[NE].link=head[u];
	head[u]=NE++;
	edge[NE].u=v;
	edge[NE].v=u;
	edge[NE].w=0;
	edge[NE].link=head[v];
	head[v]=NE++;
}
int max_flow()
{
	int ans=0,i,k,aug=inf,u,v;
	memset(dis,0,sizeof(dis));
	memset(gap,0,sizeof(gap));
	gap[0]=ed;
	u=pre[st]=st;
	for(i=1;i<ed;i++)
		cur[i]=head[i];
	while(dis[st]<ed)
	{
loop:	for(i=cur[u];i;i=cur[u]=edge[i].link)
		{
			v=edge[i].v;
			if(edge[i].w && dis[u]==dis[v]+1)
			{
				if(aug>edge[i].w)
					aug=edge[i].w;
				pre[v]=u;
				u=v;
				if(v==ed)
				{
					ans+=aug;
					for(u=pre[v];v!=st;v=u,u=pre[u])
					{
						edge[cur[u]].w-=aug;
						edge[cur[u]^1].w+=aug;
					}
					aug=inf;
				}
				goto loop;
			}
		}
		k=ed;
		for(i=head[u];i;i=edge[i].link)
		{
			v=edge[i].v;
			if(edge[i].w && k>dis[v])
			{
				k=dis[v];
				cur[u]=i;
			}
		}
		if(--gap[dis[u]]==0)
			break;
		gap[dis[u]=k+1]++;
		u=pre[u];
	}
	return ans;
}
int find(int x)
{
	if(f[x]==x)
		return x;
	return f[x]=find(f[x]);
}

int main()
{
	int T,i,j,k,x,y,fx,fy,ans;
	scanf("%d",&T);
	while(T--)
	{
		scanf("%d%d%d",&n,&m,&k);
		for(i=1;i<=m;i++)
		{
			scanf("%d%d",&p[i].u,&p[i].v);
		}
		for(i=1;i<=n;i++)
			f[i]=i;
		while(k--)
		{
			scanf("%d%d",&x,&y);
			fx=find(x);
			fy=find(y);
			if(fx==fy)
				continue;
			f[fx]=fy;
		}
		for(i=1;i<=n;i++)
			find(i);
		st=2*n+1;
		ed=2*n+2;
		NE=2;
		memset(head,0,sizeof(head));
		for(i=1;i<=n;i++)
		{
			addedge(st,i,0);
		}
		for(i=1;i<=n;i++)
		{
			addedge(i+n,ed,0);
		}
		memset(map,false,sizeof(map));
		for(i=1;i<=m;i++)
		{
			x=p[i].u;
			y=p[i].v;
			for(j=1;j<=n;j++)
			{
				if(f[x]==f[j] && !map[j][y])
				{
					addedge(j,y+n,1);
					map[j][y]=true;
				}
			}
		}
		for(i=2;i<NE;i++)
			E[i]=edge[i];
		int l=0,r=n,mid,tmp;
		while(l<=r)
		{
			mid=(l+r)>>1;
			for(i=2;i<=4*n;i+=2)
			{ 
	     		edge[i].w=mid;
				edge[i+1].w=0;
			} 
			for(i=4*n+2;i<NE;i++)
				edge[i]=E[i];
			/*for(i=2;i<NE;i++)
				cout<<i<<' '<<edge[i].u<<' '<<edge[i].v<<' '<<edge[i].w<<endl;*/
			tmp=max_flow();
		//	cout<<tmp<<endl;
			if(tmp==n*mid)
			{
				ans=mid;
				l=mid+1;
			}
			else
				r=mid-1;
		}
		cout<<ans<<endl;
	}
	return 0;
}