[DLX重复覆盖] poj 1084 Square Destroyer

题意：

n*n的矩形阵(n<=5)，由2*n*(n+1)根火柴构成，那么其中会有很多诸如边长为1，，为2…为n的正方形，现在可以拿走一些火柴，那么就会有一些正方形被破坏掉。

求在已经拿走一些火柴的情况下，还需要拿走至少多少根火柴可以把所有的正方形都破坏掉。

思路：

对于每个位置遍历所有可能的边长，确定这个边长下的正方形的边对应的都是数字几，并且把正方形从1开始编号。

然后根据编号，把正方形和数字建边记录方便下面建图。

然后以火柴棍为行，正方形为列，建立dancing link

然后求解。

这里注意的是，需要强行插入某些行。

代码：

#include"stdio.h"#include"algorithm"#include"string.h"#include"iostream"#include"cmath"#include"queue"#include"map"#include"vector"#include"string"using namespace std;#define RN 70#define CN 70#define N 70*70int fuck[123];vector<int>edge[123];struct DLX{int n,m,C;int U[N],D[N],L[N],R[N],Row[N],Col[N];int H[RN],S[CN],cnt,ans[RN];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==0?m:i-1);R[i]=(i==m?0:i+1);}C=m;for(int i=1; i<=n; i++) H[i]=-1;}void link(int x,int y){C++;Row[C]=x;Col[C]=y;S[y]++;U[C]=U[y];D[C]=y;D[U[y]]=C;U[y]=C;if(H[x]==-1) H[x]=L[C]=R[C]=C;else{L[C]=L[H[x]];R[C]=H[x];R[L[H[x]]]=C;L[H[x]]=C;}}void del(int x){for(int i=D[x]; i!=x; i=D[i]){R[L[i]]=R[i];L[R[i]]=L[i];}}void rec(int x){for(int i=U[x]; i!=x; i=U[i]){R[L[i]]=i;L[R[i]]=i;}}int used[CN];int h(){int sum=0;for(int i=R[0]; i!=0; i=R[i]) used[i]=0;for(int i=R[0]; i!=0; i=R[i]){if(used[i]==0){sum++;used[i]=1;for(int j=D[i]; j!=i; j=D[j]) for(int k=R[j]; k!=j; k=R[k]) used[Col[k]]=1;}}return sum;}void dance(int x){if(x+h()>=cnt) return ;if(R[0]==0){cnt=min(cnt,x);return ;}int now=R[0];for(int i=R[0]; i!=0; i=R[i]){if(S[i]<S[now])now=i;}for(int i=D[now]; i!=now; i=D[i]){del(i);for(int j=R[i]; j!=i; j=R[j]) del(j);dance(x+1);for(int j=L[i]; j!=i; j=L[j]) rec(j);rec(i);}return ;}void DeleteTrick(int r) //强行先放入某行。{if(H[r] == -1) return ;for(int i = D[H[r]]; i != H[r]; i = D[i]){if(H[Row[i]] == i){if(R[i] == i){H[Row[i]] = -1;}else{H[Row[i]] = R[i];}}L[R[i]] = L[i];R[L[i]] = R[i];}for(int i = R[H[r]]; i != H[r]; i = R[i]){for(int j = D[i]; j != i; j = D[j]){if(H[Row[j]] == j){if(R[j] == j){H[Row[j]] = -1;}else{H[Row[j]] = R[j];}}L[R[j]] = L[j];R[L[j]] = R[j];}}}} dlx;int main(){int t;cin>>t;while(t–){int n,m;scanf("%d%d",&n,&m);for(int i=0; i<m; i++) scanf("%d",&fuck[i]);for(int i=0; i<=70; i++) edge[i].clear();int sum=0;for(int i=1; i<=n; i++){for(int j=1; j<=n; j++){int num=(i-1)*(n+n+1)+j;int lit=min(n-i+1,n-j+1);for(int k=1; k<=lit; k++){sum++;int up,down,left,right;up=num;left=num+n;right=left+k;down=up+k*(n+n+1);for(int o=0; o<k; o++){edge[up+o].push_back(sum);edge[left+o*(n+n+1)].push_back(sum);edge[right+o*(n+n+1)].push_back(sum);edge[down+o].push_back(sum);}}}}dlx.init(2*n*(n+1),sum);for(int i=1; i<=2*n*(n+1); i++){for(int j=0; j<(int)edge[i].size(); j++) dlx.link(i,edge[i][j]);}int fff[123];memset(fff,0,sizeof(fff));//for(int j=0; j<(int)edge[1].size(); j++) printf("%d ",edge[1][j]);//puts("");for(int i=0; i<m; i++){dlx.DeleteTrick(fuck[i]);}dlx.cnt=999;dlx.dance(0);printf("%d\n",dlx.cnt);}return 0;}/*44 101 2 3 4 5 6 7 8 9 10*/

版权声明：本文为博主原创文章，未经博主允许不得转载。

上帝助自助者。