goj 1460 dag的最小路径覆盖

思路:直接把每一个相邻的植物都用边连起来，然后求一次最小路径覆盖就可以了。

注意:匈牙利算法返回1代表找到一条增广路，对应一个匹配。

最小路径覆盖=最大独立集=节点数-最大匹配数。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
using namespace std;
const int maxn=100;

char s[maxn][maxn];
int h,w,cnt,dirx[]={1,-1,0,0},diry[]={0,0,1,-1},g[maxn][maxn],num[maxn][maxn],ny[maxn*maxn],T;
bool done[maxn*maxn];

struct Flower{
	int x,y;
}f[maxn*maxn];

bool judge(int x,int y){
	if(x>=0&&x<h&&y>=0&&y<w) return true;
	return false;
}

int hungary(int u){
	for(int i=0;i<cnt;i++){
		if(g[u][i]&&!done[i]){
			done[i]=true;
			if(!ny[i]||hungary(ny[i])){
				ny[i]=u;
				return 1;
			}
		}
	}
	return 0;
}

int main(){
	while(scanf("%d",&T)!=EOF){
		for(int cas=1;cas<=T;cas++){
			memset(f,0,sizeof(f));
			memset(num,0,sizeof(num));
			memset(g,0,sizeof(g));
			memset(ny,0,sizeof(ny));
			cnt=0;
			scanf("%d%d",&h,&w);
			for(int i=0;i<h;i++)
			for(int j=0;j<w;j++){
			scanf(" %c",&s[i][j]);///这里前面留一个空格吸收前面的回车 
			if(s[i][j]=='*')
				f[cnt].x=i,f[cnt].y=j,num[i][j]=cnt++;
			}
			for(int i=0;i<cnt;i++)
				for(int j=0;j<4;j++){
					int xx=f[i].x+dirx[j];
					int yy=f[i].y+diry[j];
					if(judge(xx,yy)&&s[xx][yy]=='*'){
						int id=num[xx][yy];
						if(g[i][id]) continue;
						g[i][id]=g[id][i]=1;
					}
				}
			int ans=0;
			for(int i=0;i<cnt;i++){
				memset(done,0,sizeof(done));
				ans+=hungary(i);
			}
			printf("%d\n",cnt-(ans>>1));
		}
	}
	return 0;
}