ssoj2279: 磁力阵（IDA*）

题意：给一个由小正方形组成的大正方形，其中有k条边被删除，问最少再删除几条边图中没有正方形

IDA*：

        1、总边数（2*n*(n+1)<=60）可以状压。预处理出各个边长的正方形的状态。

        2、设置估价函数h*（到目标还要删除的最小正方形数）<=h（到目标还要删除的最小边数）。若（h*+stp>ans）剪枝。

        3、枚举第一个找到的正方形的边，删除，继续dfs。

贴代码：

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
using namespace std;
typedef long long ll;
const int maxn=66;
int n,k,cnt,sum,dep;
ll tar,a[6][6],sta[maxn];
inline int get(){
    char c;while(!isdigit(c=getchar()));
    int v=c-48;while(isdigit(c=getchar()))v=v*10+c-48;
    return v;
}
inline ll shl(int x){return 1ll<<(x-1);}
inline int heng(int x,int y){return (x-1)*(2*n+1)+y;}
inline int shu(int x,int y){return x*(2*n+1)-n-1+y;}
inline void predeal(){
    for(int i=1;i<=n;++i)
        for(int j=1;j<=n;++j){
		    a[i][j]=shl(heng(i,j))|shl(heng(i+1,j))|shl(shu(i,j))|shl(shu(i,j+1));
		    sta[++cnt]=a[i][j];
		}
	for(int l=2;l<=n;++l)
	    for(int i=1;i<=n+1-l;++i)
	        for(int j=1;j<=n+1-l;++j){
			    sta[++cnt]=0;
			    for(int x=i;x<i+l;++x)
			        for(int y=j;y<j+l;++y)
			            sta[cnt]^=a[x][y];
			}
}
inline bool dfs(ll tst,int stp){
    ll tmp=tst;int h=0;ll fir=0;
    for(int i=1;i<=cnt;++i)if((tmp&sta[i])==sta[i]){//<span style="font-family: Arial, Helvetica, sans-serif;">tmp&sta[i]   ！！！记得加括号！！！</span>
	    ++h;tmp^=sta[i];
	    if(!fir)fir=sta[i];
	}
	if(h==0)return 1;             
	else if(h+stp>dep)return 0;  //估价函数
	for(int i=1;i<=sum;++i)if((fir&shl(i)) && dfs(tst^shl(i),stp+1))return 1;
	return 0;
}
int main(){
	int t=get();
	while(t--){
	    n=get();k=get();sum=2*n*(n+1);
	    tar=shl(sum+1)-1;cnt=0;dep=0;
	    predeal();
	    while(k--){
		    int x=get();
		    tar^=shl(x);
		}
		while(!dfs(tar,0))++dep;
		printf("%d\n",dep);
	}
	return 0;
}