ssoj2279: 磁力阵（IDA*）

来源：互联网发布：ipad视频软件编辑：程序博客网时间：2024/05/17 09:41

题意：给一个由小正方形组成的大正方形，其中有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;}

0 0