poj2286 唯美的IDA*(80行)

ida*果然比a*强太多了,少了hash判重和堆维护,是估价函数和DFS的完美结合啊!

嗯,主要是DFS函数中第二行的那个强剪枝太给力了!IDA*的灵魂啊!

基本的架构就是如此了,可以任意套。

//poj2286
#include <iostream>
#include <cmath>
using namespace std;
const int b[5][8] = { 7,8,9,12,13,16,17,18,
                      1,3,7,12,16,21,23,0,
                      2,4,9,13,18,22,24,0,
                      11,10,9,8,7,6,5,0,
                      20,19,18,17,16,15,14,0};
const int c[9] = {0,1,2,3,4,6,5,8,7};
int a[25],stk[1010],bound,x,final;
bool ans;
int h(int *a) {
    int tmp[3]={0,0,0};
    for (int i=0;i<8;i++)
        tmp[a[b[0][i]]-1]++;
    return 8-max(max(tmp[0],tmp[1]),tmp[2]);
}
void trans(int *a,int p) {
     if (p<=4) {
               int tmp = a[b[p][0]];
               for (int i=0;i<6;i++)
                   a[b[p][i]]=a[b[p][i+1]];
               a[b[p][6]]=tmp;
     } else {
               p-=4;
               int tmp = a[b[p][6]];
               for (int i=6;i>0;i--)
                   a[b[p][i]]=a[b[p][i-1]];
               a[b[p][0]]=tmp;
     }
}
int dfs(int *a,int dep,int pre_move) {
    int hv = h(a);
    if (hv + dep > bound) return hv+dep;
    if (hv==0) {
               ans = true;
               final = a[7];
               return dep;
    }
    int next_bound = 10000;
    for (int i=1;i<=8;i++) {
        if (abs(c[i]-pre_move)==4) continue;   //与上一步相反的移动
        int tmp[25];
        for (int j=1;j<=24;j++) tmp[j]=a[j];
        trans(tmp,c[i]);
        stk[dep]=i;
        int new_bound = dfs(tmp,dep+1,c[i]);
        if (ans) return new_bound;
        next_bound = min(next_bound,new_bound);
    }
    return next_bound;
}
void ida_star() {
     memset(stk,0,sizeof(stk));
     bound = h(a);
     ans = false;
     while (!ans && bound<=100)
           bound = dfs(a,0,-100);
     for (int i=0;i<bound;i++) printf("%c",stk[i]+64);
     printf("\n%d\n",final);
}
int main() {
    while (~scanf("%d",&x)) {
          if (x==0) break;
          a[1]=x;
          for (int i=2;i<=24;i++)
              scanf("%d",&a[i]);
                  
          if (h(a)==0) {
                    printf("No moves needed\n");
                    printf("%d\n",a[7]);
                    continue;
          }
          ida_star();
    }
    return 0;
}


这是一道比较经典的计数问题。题目描述如下: 给定一个 $n \times n$ 的网格图,其中一些格子被标记为障碍。一个连通块是指一些被标记为障碍的格子的集合,满足这些格子在网格图中连通。一个格子是连通的当且仅当它与另一个被标记为障碍的格子在网格图中有公共边。 现在,你需要计算在这个网格图中,有多少个不同的连通块,满足这个连通块的大小(即包含的格子数)恰好为 $k$。 这是一道比较经典的计数问题,一般可以通过计算生成函数的方法来解决。具体来说,我们可以定义一个生成函数 $F(x)$,其中 $[x^k]F(x)$ 表示大小为 $k$ 的连通块的个数。那么,我们可以考虑如何计算这个生成函数。 对于一个大小为 $k$ 的连通块,我们可以考虑它的形状。具体来说,我们可以考虑以该连通块的最左边、最上边的格子为起点,从上到下、从左到右遍历该连通块,把每个格子在该连通块中的相对位置记录下来。由于该连通块的大小为 $k$,因此这些相对位置一定是 $(x,y) \in [0,n-1]^2$ 中的 $k$ 个不同点。 现在,我们需要考虑如何计算这些点对应的连通块是否合法。具体来说,我们可以考虑从左到右、从上到下依次处理这些点,对于每个点 $(x,y)$,我们需要考虑它是否能够与左边的点和上边的点连通。具体来说,如果 $(x-1,y)$ 和 $(x,y)$ 都在该连通块中且它们在网格图中有公共边,那么它们就是连通的;同样,如果 $(x,y-1)$ 和 $(x,y)$ 都在该连通块中且它们在网格图中有公共边,那么它们也是连通的。如果 $(x,y)$ 与左边和上边的点都不连通,那么说明这个点不属于该连通块。 考虑到每个点最多只有两个方向需要检查,因此时间复杂度为 $O(n^2 k)$。不过,我们可以使用类似于矩阵乘法的思想,将这个过程优化到 $O(k^3)$ 的时间复杂度。 具体来说,我们可以设 $f_{i,j,k}$ 表示状态 $(i,j)$ 所代表的点在连通块中,且连通块的大小为 $k$ 的方案数。显然,对于一个合法的 $(i,j,k)$,我们可以考虑 $(i-1,j,k-1)$ 和 $(i,j-1,k-1)$ 这两个状态,然后把点 $(i,j)$ 加入到它们所代表的连通块中。因此,我们可以设计一个 $O(k^3)$ 的 DP 状态转移,计算 $f_{i,j,k}$。 具体来说,我们可以考虑枚举连通块所包含的最右边和最下边的格子的坐标 $(x,y)$,然后计算 $f_{x,y,k}$。对于一个合法的 $(x,y,k)$,我们可以考虑将 $(x,y)$ 所代表的点加入到 $(x-1,y,k-1)$ 和 $(x,y-1,k-1)$ 所代表的连通块中。不过,这里需要注意一个细节:如果 $(x-1,y)$ 和 $(x,y)$ 在网格图中没有相邻边,那么它们不能算作连通的。因此,我们需要特判这个情况。 最终,$f_{n,n,k}$ 就是大小为 $k$ 的连通块的个数,时间复杂度为 $O(n^2 k + k^3)$。 参考代码:
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值