BOJ 1298 Hungary算法

//Result:1298	Accepted	3MS	308K	G++	 1669B	2012-06-13 23:39:57	Wizmann

#include <cstdio>

#include <cstring>

#include <cstdlib>

#include <iostream>

#include <algorithm>

#include <bitset>

#include <stack>

#include <queue>

#include <vector>



using namespace std;



#define print(x) cout<<x<<endl

#define input(x) cin>>x

#define N 10

#define SIZE 100



const int mx[]={1,0,-1,0};

const int my[]={0,1,0,-1};



vector<int> g[SIZE];

char maze[N][N];

int r,c;

bitset<SIZE> visit;

int match[SIZE];



inline int nr(int x,int y)

{

	return c*y+x;

}



void makeG()

{

	for(int i=0;i<SIZE;i++) g[i].clear();

	for(int i=0;i<r;i++)

	{

		for(int j=0;j<c;j++)

		{

			if(maze[i][j]=='*')

			{

				for(int k=0;k<4;k++)

				{

					int ny=i+my[k];

					int nx=j+mx[k];

					if(ny>=0 && ny<r && nx>=0 && nx<c && maze[ny][nx]=='*')

					{

						g[nr(j,i)].push_back(nr(nx,ny));

					}

				}

			}

		}

	}

}



bool dfs(int now)

{

	for(int i=0;i<(int)g[now].size();i++)

	{

		int next=g[now][i];

		if(!visit[next])

		{

			visit[next]=1;

			if(match[next]==-1 || dfs(match[next]))

			{

				match[next]=now;

				return true;

			}

		}

	}

	return false;

}



int hungary()

{

	int res=0;

	memset(match,-1,sizeof(match));

	for(int i=0;i<r*c;i++)

	{

		visit.reset();

		if(dfs(i)) res++;

	}

	return res;

}



int main()

{

	freopen("input.txt","r",stdin);

	int cas=1;

	while(input(r>>c))

	{

		printf("Case %d:\n",cas++);

		int block=0;

		for(int i=0;i<r;i++)

		{

			scanf("%s",maze[i]);

			for(int j=0;maze[i][j];j++)

			{

				if(maze[i][j]=='*') block++;

			}

		}

		if(block&1) print("no");

		else

		{

			makeG();

			int sum=hungary();

			if(sum==block) print("yes");

			else print("no");

		}

	}

	return 0;

}