![](http://acm.fzu.edu.cn/image/problem.gif)
Accept: 528 Submit: 1601
Time Limit: 1000 mSec Memory Limit : 32768 KB
Problem Description
Input
Output
n行m列 每个位置上都是0或1 现在攻击其中的n1行 m1列 问最少挑选多少行或列 使得每行或列至少有一个1
Dancing Links可重复覆盖模型:在0-1矩阵中选取最少的行,使得每一列最少有一个1
#include <cstdio>
#include <iostream>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <string.h>
#include <string>
#include <vector>
#include <queue>
#define MEM(a,x) memset(a,x,sizeof a)
#define eps 1e-8
#define MOD 10009
#define INF 99999999
#define ll __int64
#define bug cout<<"here"<<endl
#define fread freopen("ceshi.txt","r",stdin)
#define fwrite freopen("out.txt","w",stdout)
using namespace std;
const int N=15;
const int MAXM=N*N+10;
const int MAXN=N*N+10;
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 ret=0;
for(int c=R[0];c!=0;c=R[c]) v[c]=1;
for(int c=R[0];c!=0;c=R[c])
{
if(v[c])
{
ret++;
v[c]=0;
for(int i=D[c];i!=c;i=D[i])
for(int j=R[i];j!=i;j=R[j])
v[Col[j]]=0;
}
}
return ret;
}
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!=0;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);
}
}
};
DLX dlx;
int a[20][20];
int id[20][20];
//在0-1矩阵中选取最少的行 使得每一列最少有一个1
int main()
{
// fread;
int n,m;
while(scanf("%d%d",&n,&m)!=EOF)
{
int sz=0;
MEM(id,0);
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);
}
dlx.init(n*m,sz);
sz=1;
int n1,m1;
scanf("%d%d",&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])
dlx.Link(sz,id[i+x][j+y]);
sz++;
}
dlx.ansd=INF;
dlx.Dance(0);
printf("%d\n",dlx.ansd);
}
return 0;
}