FZU - 1686 神龙的难题 （舞蹈链，可重覆盖）

这是个剑与魔法的世界.英雄和魔物同在,动荡和安定并存.但总的来说,库尔特王国是个安宁的国家,人民安居乐业,魔物也比较少.但是.总有一些魔物不时会进入城市附近,干扰人民的生活.就要有一些人出来守护居民们不被魔物侵害.魔法使艾米莉就是这样的一个人.她骑着她的坐骑,神龙米格拉一起消灭干扰人类生存的魔物,维护王国的安定.艾米莉希望能够在损伤最小的前提下完成任务.每次战斗前,她都用时间停止魔法停住时间,然后米格拉他就可以发出火球烧死敌人.米格拉想知道,他如何以最快的速度消灭敌人,减轻艾米莉的负担.

Input

数据有多组,你要处理到EOF为止.每组数据第一行有两个数,n,m,(1<=n,m<=15)表示这次任务的地区范围. 然后接下来有n行,每行m个整数,如为1表示该点有怪物,为0表示该点无怪物.然后接下一行有两个整数,n1,m1 (n1<=n,m1<=m)分别表示米格拉一次能攻击的行,列数(行列不能互换，这里指的是可以攻击一个n1*m1的矩形范围),假设米格拉一单位时间能发出一个火球,所有怪物都可一击必杀.

 

Output

 

输出一行,一个整数,表示米格拉消灭所有魔物的最短时间.

 

Sample Input

4 4
1 0 0 1
0 1 1 0
0 1 1 0
1 0 0 1
2 2
4 4 
0 0 0 0
0 1 1 0
0 1 1 0
0 0 0 0
2 2

Sample Output

4
1

 

题解：在每一个格子内施法都可以杀掉  以该点为左上角的n1*m1的矩形，问题转换为：求最少的攻击次数，可以杀掉所有的怪物，很显然的可重覆盖问题了。行：格子，列：怪物，  如果在该格子内可以攻击到怪物，就Link。

#include<iostream>
#include<cstdio>
using namespace std;
#define rep(i,a,n) for (int i=a;i<n;i++)
#define per(i,a,n) for (int i=n-1;i>=a;i--)
#define oo cout<<"!!!"<<endl;
typedef long long ll;
typedef unsigned long long ull;
#define ms(s) memset(s, 0, sizeof(s))
const int inf = 0x3f3f3f3f;
const int mod = 1e9+7;
//head

const int maxn = 15*15+11;
const int maxm = 15*15+11;
const int maxnode = maxn*maxm;

struct DLX{

	int n,m,size;
	int U[maxnode],D[maxnode],R[maxnode],L[maxnode],Row[maxnode],Col[maxnode];
	int H[maxn],S[maxm];
	int ansd;
	void init(int _n,int _m)
	{
		n = _n,m = _m;
		for(int i = 0;i<=m;++i)
		{
			S[i] = 0;
			U[i] = D[i] = i;
			L[i] = i-1;
			R[i] = i+1; 
		}
		R[m] = 0,L[0] = m;
		size = m;
		for(int i = 1;i<=n;++i)H[i] = -1;

	}

	void link(int r,int c)
	{
		++S[ Col[++size] = c ];
		Row[size] = r;

		D[size] = D[c];
		U[ D[c] ] = size;
		U[size] = c;
		D[c] = size;

		if(H[r] < 0)H[r] = L[size] = R[size] = size;
		else
		{
			R[size] = R[H[r]];
			L[R[H[r]]] = size;
			L[size] = H[r];
			R[H[r]] = size;
		}
	}

	void remove(int c)
	{
		for(int i = D[c];  i != c;i = D[i])
			L[ R[i] ] = L[i],R[ L[i] ] = R[i];
	}

	void resume(int c)
	{
		for(int i = U[c];  i != c;i = U[i])
			L[R[i]] = R[ L[i] ] = i;
	}

	bool v[maxm];

	int f()
	{
		int res = 0;
		for(int c = R[0];c;c = R[c])v[c] = true;

		for(int c = R[0];c;c=R[c])
			if(v[c])
			{
				res++;
				v[c] = false;
				for(int i = D[c];i!=c;i = D[i])
					for(int j = R[i];j!=i;j = R[j])
						v[ Col[j] ] = false;
 			}

 		return res;
	}

	void dance(int d)
	{
		if(d+f() >= ansd)return;
		if(R[0] == 0)
		{
			if(d < ansd)ansd = d;
			return;
		}

		int c = R[0];
		for(int i = R[0];i;i = R[i])
			if(S[i] < S[c])
				c = i;

		for(int i = D[c];i != c;i = D[i])
		{
			remove(i);
			for(int j = R[i];j != i;j = R[j])remove(j);
			dance(d+1);
			for(int j = L[i];j != i; j = L[j])resume(j);
			resume(i);
		}	
	}	
}g;

int a[22][22];
int id[22][22];
int main() 
{
	int n,m;
	while(~scanf("%d%d",&n,&m))
	{
		int sz = 0;
		ms(id);
		for(int i = 0;i<n;++i)
			for(int j = 0;j<m;++j)
			{
				scanf("%d",&a[i][j]);
				if(a[i][j] == 1)id[i][j] = ++sz;
			}
		g.init(n*m,sz);
		sz = 1;
		int n1,m1;
		cin>>n1>>m1;
		for(int i = 0;i<n;++i)
			for(int j = 0;j<m;++j)
			{
				for(int x = 0;x<n1 && i+x<n;++x)
				for(int y = 0;y<m1 && j+y<m;++y)
					if(id[i+x][j+y])
						g.link(sz,id[i+x][j+y]);
				++sz;
			}
		g.ansd = inf;
		g.dance(0);
		cout << g.ansd<<endl;

	}
    return 0;
}