当前位置: 代码迷 >> 综合 >> FZU - 1686 神龙的难题 (舞蹈链,可重覆盖)
  详细解决方案

FZU - 1686 神龙的难题 (舞蹈链,可重覆盖)

热度:85   发布时间:2024-01-22 01:31:52.0

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

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;
//headconst 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;
}