ZOJ1654:Place the Robots

http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=654

这个题目看题后会很直接以空白处作为顶点建图，这样问题转化为求最大独立集问题，当然我们可以进一步转化为二分图匹配问题，利用二分图的性质：最大顶点独立集个数= 顶点个数 —最小顶点覆盖，而最小顶点覆盖=最大匹配，这一点在前面的Boys And Girls的解题报告中已经使用了。这里最麻烦的就是建图，一般的建图时间代价是O(n^3)，多次提交的结果是一只超时，优化后效果不大，我暂未想到更好的建图方法来降低代价。不过我们可以从另一种方式来考虑，即：将水平方向与垂直方向分别编号作为图的顶点，这样就形成两个独立的集合，这样就直接转化为二分图的最大匹配，只需求出最大匹配的边即可。源代码如下：

#include<iostream>
#include<string.h>
#include<stdio.h>
using namespace std ;

const int maxn = 51 ;

int m ;
int n ;

char map[maxn][maxn] ;

int x[maxn*maxn] ;
int y[maxn*maxn] ;
int xs[maxn][maxn] ;  
int ys[maxn][maxn] ;

int xn  ;
int yn1 ;
bool g[maxn * maxn][maxn * maxn] ;

int t[maxn * maxn] ;

bool path(int u)
{
	int v ;
	for(v = 1 ; v <= yn1 ; v ++)
	{
		if(g[u][v] && !t[v])
		{
			t[v] = 1 ;
			
			if(!y[v] || path(y[v]))
			{
				x[u] = v ;
				y[v] = u ;
				return 1 ;
			}	
		}
	}
	return 0 ;
} 

void MaxMatch()
{
	int i;
	int ans = 0 ;
	memset(x , 0 , sizeof(x)) ;
	memset(y , 0 , sizeof(y)) ;
	
	for(i = 1 ; i <= xn ; i ++)
	{
		if(!x[i])
		{
			memset(t , 0 , sizeof(t)) ;
			if(path(i))
				ans ++ ;
		}
	}
	printf("%d\n" , ans) ;
}

int main()
{
	int k , kcase ;
	int i ;
	int j ;
	int number ;
	int flag ;
	scanf("%d" , &kcase) ;
	for(k = 0 ; k < kcase ; k ++)
	{
		printf("Case :%d\n" , k+1) ;
		scanf("%d%d", &m ,&n) ;
		memset(xs , 0, sizeof(xs)) ;
		
		for(i = 0 ; i < m ; i ++)
		{
			scanf("%s" , map[i]) ;
		}
		number = 0 ;
		for(i = 0 ; i < m  ; i ++)
		{
			flag = 0 ;
			for(j = 0 ; j < n ; j ++)
			{
				if(map[i][j]=='o')
				{
					if(flag ==0)
						number++ ;
					xs[i][j] = number ;
					flag = 1 ;	
				}
				else if(map[i][j]=='#')
					flag = 0 ;
			}
		}
		xn = number ;
		number = 0 ;
		for(j = 0 ; j < n ; j ++)
		{
			flag = 0 ;
			for(i = 0 ; i < m ; i ++)
			{
				if(map[i][j] == 'o')
				{
					if(flag ==0)
						number ++ ;
					ys[i][j] = number ;
					flag = 1 ;	
				}
				else if(map[i][j]=='#')
					flag = 0 ;
			}
		}
		
		yn1 = number ;
		memset(g , 0, sizeof(g)) ;
		
		for(i = 0 ; i < m ; i ++)
		{
			for(j = 0 ; j < n ; j ++)
			{
				if(xs[i][j])
					g[xs[i][j]][ys[i][j]] = 1 ;
			}
		}
		
		MaxMatch() ;
	}
	return 0 ;
}