poj1469(最大二分匹配)

题目大意:

一共有N个学生跟P门课程,一个学生可以任意选一

门或多门课,问是否达成:

  1.每个学生选的都是不同的课(即不能有两个学生选同一门课)

  2.每门课都有一个代表(即P门课都被成功选过)

输入为:

P N(课程数跟学生数)

接着有P行,格式为Count studenti studenti+1 ……studentcount

(Count表示对课程1感兴趣的学生数,接着有Count个学生)

如第一行3 1 2 3表示学生1跟学生2跟学生3对课程1感兴趣

输出为:

若能满足上面两个要求这输出”YES”,否则为”NO”

思路:其实只是一道二分最大匹配的裸题,关键在于算法的选择上,一开始选了最大流法,超时。后来选了匈牙利算法,不仅代码简单,而且速度快。碰上二分最大匹配还是用匈牙利算法吧,心累。

//用dinic算法,当最大流来做,超时!!!!所以可以看看后面的匈牙利算法

匈牙利算法:

#include<cstdio>
#include<iostream>
#include<cstring>
#include<queue>
#include<vector>
#define inf 0x3f3f3f3f
using namespace std;
int n,p;
int map[300][300],chk[300],used[300];
int find(int u)
{
	for(int i=1;i<=n;i++)
	{
		if(map[u][i]&&!chk[i])
		{
			chk[i]=1;
			if(!used[i]||find(used[i]))
			{
				used[i]=u;
				return 1;
			}
		}
	}
	return 0;
}
void solve()
{
	int res=0;
	memset(used,0,sizeof(used));
	for(int i=1;i<=p;i++)
	{
		memset(chk,0,sizeof(chk));
		res+=find(i);
	}
	if(res==p) printf("YES\n");
	else printf("NO\n");
}
int main()
{
	//freopen("t.txt","r",stdin);
	int w,t,a;
	scanf("%d",&w);
	while(w--)
	{
		memset(map,0,sizeof(map));
		scanf("%d%d",&p,&n);
		for(int i=1;i<=p;i++)
		{
			scanf("%d",&t);
			while(t--)
			{
				scanf("%d",&a);
				map[i][a]=1;
			}
		}
		solve();
	}
	return 0;
 } 

dinic法:

#include<cstdio>
#include<iostream>
#include<cstring>
#include<queue>
#include<vector>
#define inf 0x3f3f3f3f
using namespace std;
int p,n;
int dis[809],M;
struct node
{
	int v,cap,id;
	node(int v1,int cap1,int id1)
	{
		v=v1,cap=cap1,id=id1;
	}
};
vector<node> map[405];
void add(int u,int v,int cap)
{
	int l1=map[u].size();
	int l2=map[v].size();
	map[u].push_back(node(v,cap,l2));
	map[v].push_back(node(u,0,l1));
}
bool bfs()
{
	queue<int> que;
	memset(dis,0,sizeof(dis));
	dis[0]=1;
	que.push(0);
	while(!que.empty())
	{
		int t=que.front();que.pop();
		for(int i=0;i<map[t].size();i++)
		{
			if(!map[t][i].cap) continue;
			int v=map[t][i].v;
			if(!dis[v])
			{
				dis[v]=dis[t]+1;
				if(v==M) return true;
				que.push(v);
			}
		}
	}
	return false;
}
int dfs(int u,int cap)
{
	if(u==M) return cap;
	int tp=cap,t;
	for(int i=0;i<map[u].size();i++) 
	{
		if(!map[u][i].cap) continue;
		int v=map[u][i].v;
		if(dis[v]==dis[u]+1)
		{
			t=dfs(v,min(tp,map[u][i].cap));
			tp-=t;
			map[u][i].cap-=t;
			map[v][map[u][i].id].cap+=t;
		}
	}
	return cap-tp;
}
void dinic()
{
	int ans=0;
	while(bfs())
	{
		ans+=dfs(0,inf);
	}
	if(ans==p) printf("YES\n");
	else printf("NO\n");
}
int main()
{
	//freopen("t.txt","r",stdin);
	int w,t,a;
	scanf("%d",&w);
	while(w--)
	{
		for(int i=0;i<405;i++) map[i].clear();
		scanf("%d%d",&p,&n);
		M=p+n+1;
		for(int i=1;i<=p;i++)	add(0,i,1);
		for(int i=1;i<=n;i++) add(p+i,p+n+1,1);
		for(int i=1;i<=p;i++)
		{
			scanf("%d",&t);
			while(t--)
			{
				scanf("%d",&a);
				add(i,a+p,1);
			}
		}
		dinic();
	}
	return 0;
 } 

 

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值