【HDU】2444 The Accomodation of Students 二分匹配

传送门:【HDU】2444 The Accomodation of Students


题目分析:能不能分成两个集合主要看这个图是否是一个二分图,我们可以用黑白染色法判定二分图。如果一个点的后继是已经染色的且该后继的颜色和该点相同,说明这个图一定不是二分图,否则判断到所有点都成功染色后,那么该图就是二分图了。由于不知道是否是连通图,所以用并查集得到各个连通块,分别判断是否是二分图。接下来就是二分匹配的事了。


代码如下:


#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std ;

#define REP( i , a , b ) for ( int i = a ; i < b ; ++ i )
#define FOR( i , a , b ) for ( int i = a ; i <= b ; ++ i )
#define REV( i , a , b ) for ( int i = a ; i >= b ; -- i )
#define CLR( a , x ) memset ( a , x , sizeof a )

const int MAXN = 205 ;
const int MAXQ = 205 ;
const int MAXE = 80005 ;
const int INF = 0x3f3f3f3f ;

struct Edge {
	int v , n ;
	Edge () {}
	Edge ( int v , int n ) : v ( v ) , n ( n ) {}
} E[MAXE] ;

int H[MAXN] , cntE ;
int dx[MAXN] , dy[MAXN] ;
int Lx[MAXN] , Ly[MAXN] ;
int Q[MAXQ] , head , tail ;
int vis[MAXN] ;
int dis ;
int n , m ;

int p[MAXN] ;
int color[MAXN] ;

void init () {
	cntE = 0 ;
	CLR ( H , -1 ) ;
}

void addedge ( int u , int v ) {
	E[cntE] = Edge ( v , H[u] ) ;
	H[u] = cntE ++ ;
}

int Hopcroft_Karp () {
	dis = INF ;
	head = tail = 0 ;
	CLR ( dx , -1 ) ;
	CLR ( dy , -1 ) ;
	FOR ( i , 1 , n )
		if ( Lx[i] == -1 ) {
			Q[tail ++] = i ;
			dx[i] = 0 ;
		}
	while ( head != tail ) {
		int u = Q[head ++] ;
		if ( dx[u] >= dis )
			continue ;
		for ( int i = H[u] ; ~i ; i = E[i].n ) {
			int v = E[i].v ;
			if ( dy[v] == -1 ) {
				dy[v] = dx[u] + 1 ;
				if ( Ly[v] == -1 )
					dis = dy[v] ;
				else {
					dx[Ly[v]] = dy[v] + 1 ;
					Q[tail ++] = Ly[v] ;
				}
			}
		}
	}
	return dis != INF ;
}

int hungary ( int u ) {
	for ( int i = H[u] ; ~i ; i = E[i].n ) {
		int v = E[i].v ;
		if ( !vis[v] && dy[v] == dx[u] + 1 ) {
			vis[v] = 1 ;
			if ( ~Ly[v] && dy[v] == dis )
				continue ;
			if ( Ly[v] == -1 || hungary ( Ly[v] ) ) {
				Lx[u] = v ;
				Ly[v] = u ;
				return 1 ;
			}
		}
	}
	return 0 ;
}

int match () {
	int res = 0 ;
	CLR ( Lx , -1 ) ;
	CLR ( Ly , -1 ) ;
	while ( Hopcroft_Karp () ) {
		CLR ( vis , 0 ) ;
		FOR ( i , 1 , n )
			if ( Lx[i] == -1 )
				res += hungary ( i ) ;
	}
	return res ;
}

int find ( int x ) {
	return p[x] == x ? x : ( p[x] = find ( p[x] ) ) ;
}

int draw () {
	head = tail = 0 ;
	CLR ( color , 0 ) ;
	FOR ( i , 1 , n )
		if ( find ( i ) == i ) {
			Q[tail ++] = i ;
			color[i] = 1 ;
		}
	while ( head != tail ) {
		int u = Q[head ++] ;
		for ( int i = H[u] ; ~i ; i = E[i].n ) {
			int v = E[i].v ;
			if ( !color[v] ) {
				color[v] = 3 - color[u] ;
				Q[tail ++] = v ;
			}
			if ( color[v] + color[u] != 3 )
				return 0 ;
		}
	}
	return 1 ;
}

void solve () {
	int u , v ;
	init () ;
	FOR ( i , 1 , n )
		p[i] = i ;
	while ( m -- ) {
		scanf ( "%d%d" , &u , &v ) ;
		addedge ( u , v ) ;
		addedge ( v , u ) ;
		u = find ( u ) ;
		v = find ( v ) ;
		p[u] = v ;
	}
	if ( !draw () )
		printf ( "No\n" ) ;
	else
		printf ( "%d\n" , match () >> 1 ) ;
}

int main () {
	while ( ~scanf ( "%d%d" , &n , &m ) )
		solve () ;
	return 0 ;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值