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
41
好久不见得中文题,可谓一股清流啊,由于刚接触舞蹈链说实话只是刚刚会套个模板而已,但看到这个题就蒙了,因为对于如何构造矩阵一头雾水,在翻阅了多个大神博客后才晓得以所给矩阵中的1为列(当然需要给1标号),然后以所给攻击范围的规模大小的矩阵为行,刚开始还不太明白,但后来明白了,其实如何确定矩阵很简单,首先咱们在模板的原型中要明白是找出几行然后使得每列中至少含有一个1(覆盖所有列),然后咱们再分析本题,本题是问最少攻击几次能把敌人全部消灭,每次攻击是一个给定规模的矩阵,所以简单点看就是让我们找最少的矩阵能覆盖所有1,那么矩阵的行列就好确定了,矩阵为行,1为列
ac代码:
#include <vector> #include <stdio.h> #include <string.h> #include <stdlib.h> #include <iostream> #include <algorithm> using namespace std; typedef long long ll; typedef unsigned long long ull; const int inf=0x3f3f3f3f; const ll INF=0x3f3f3f3f3f3f3f3fll; const int maxn=260; int L[maxn*maxn],R[maxn*maxn],U[maxn*maxn],D[maxn*maxn];//节点的上下左右四个方向的链表 int C[maxn*maxn],H[maxn],cnt[maxn],vis[maxn];//C列H行cnt列链表中元素个数 int n,m,id,fans; void init(int ll){ for(int i=0;i<=ll;i++){ cnt[i]=0;U[i]=D[i]=i; L[i+1]=i;R[i]=i+1; } R[ll]=0;id=ll+1; memset(H,-1,sizeof(H)); } void Link(int r,int c){ cnt[c]++;C[id]=c; U[id]=U[c];D[U[c]]=id; D[id]=c;U[c]=id; if(H[r]==-1) H[r]=L[id]=R[id]=id; else{ L[id]=L[H[r]];R[L[H[r]]]=id; R[id]=H[r];L[H[r]]=id; } id++; } void Remove(int Size){ for(int j=D[Size];j!=Size;j=D[j]) L[R[j]]=L[j],R[L[j]]=R[j]; } void Resume(int Size){ for(int j=D[Size];j!=Size;j=D[j]) L[R[j]]=R[L[j]]=j; } int h(){ int sum=0; memset(vis,0,sizeof(vis)); for(int i=R[0];i;i=R[i]){ if(vis[i]) continue; sum++; for(int j=D[i];j!=i;j=D[j]){ for(int k=R[j];k!=j;k=R[k]) vis[C[k]]=1; } } return sum; } void Dance(int k){ int mm=maxn,pos; if(k+h()>=fans) return; if(!R[0]){ if(k<fans) fans=k; return; } for(int i=R[0];i;i=R[i]) if(mm>cnt[i]) mm=cnt[i],pos=i; for(int i=D[pos];i!=pos;i=D[i]){ Remove(i); for(int j=R[i];j!=i;j=R[j]) Remove(j); Dance(k+1); for(int j=R[i];j!=i;j=R[j]) Resume(j); Resume(i); } } int A[20][20],B[20][20]; int main(){ int u,v,len,wid; while(scanf("%d%d",&n,&m)!=-1){ int kk=0,kkk=1; for(int i=1;i<=n;i++) for(int j=1;j<=m;j++){ scanf("%d",&A[i][j]); if(A[i][j]==1) B[i][j]=(++kk);//给1编号 } scanf("%d%d",&len,&wid); init(kk); for(int i=1;i<=n-len+1;i++){ for(int j=1;j<=m-wid+1;j++){ for(int k=0;k<len;k++){ for(int l=0;l<wid;l++){ if(A[i+k][j+l]==1) Link(kkk,B[i+k][j+l]); } } kkk++;//矩阵个数 } } fans=kkk-1; Dance(0); printf("%d\n",fans); } return 0; }