【bzoj2150】 部落战争 二分图匹配

裸的有向无环图最小不相交路径覆盖，直接拆点，然后原图点数-最大匹配就可以。数组开小了，WA了一发。

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<iostream>
#include<algorithm>
#define maxn 2610

using namespace std;

int a[maxn][maxn],lk[maxn],f[maxn][maxn],num[maxn][maxn];
bool vis[maxn];
int dx[5],dy[5];
char s[510];
int n,m,r,c,sum,ans,id;

bool find(int x)
{
	for (int i=1;i<=id;i++)
	  if (a[x][i] && !vis[i])
	  {
	  	vis[i]=1;
	  	if (!lk[i] || find(lk[i]))
	  	{
	  		lk[i]=x;
	  		return 1;
	  	}
	  }
	return 0;
}

int main()
{
	scanf("%d%d%d%d",&n,&m,&r,&c);
	id=0;
	for (int i=1;i<=n;i++)
	{
		scanf("%s",s);
		for (int j=0;j<m;j++)
		  if (s[j]=='.') f[i][j+1]=1,num[i][j+1]=++id;
		  else f[i][j+1]=0;
	}
	dx[0]=r;dx[1]=r;dx[2]=c;dx[3]=c;
	dy[0]=c;dy[1]=-c;dy[2]=r;dy[3]=-r;
	for (int i=1;i<=n;i++)
	  for (int j=1;j<=m;j++)
	    if (f[i][j])
	    for (int k=0;k<4;k++)
	    {
	    	int x=i+dx[k],y=j+dy[k];
	    	if (x<1 || x>n || y<1 || y>m) continue;
	    	if (f[x][y]) a[num[i][j]][num[x][y]]=1;
	    }
	for (int i=1;i<=id;i++)
	{
		memset(vis,0,sizeof(vis));
		if (find(i)) ans++;
	}
	printf("%d\n",id-ans);
	return 0;
}