POJ 1077 Eight(BFS + A* + 逆序对剪枝)

题意:

一个3*3的棋盘里,给8个格子填上1-8这8个数,另一个格子空着。每一步,可以将与空格相邻的格子与空格换位,换位方式可间接体现成空格进行上下左右的移动。先给定一个棋盘状态,问能否经过若干步后,将棋盘变成第一行123,第二行456,第三行78空的状态。若能,输出任意一种可行步骤。

思路:

1. 对于 123456789 的某个全排列,使用 hash 表来记录状态显得浪费。本题采用了康托展开,对于每一个排列生成唯一一个键值;

2. 普通 BFS 能过,但是时间复杂度太高。这里采用启发式的方法 f = h + g; g 表示状态 s 在搜索树中的深度,h 表示状态 s 距离目标状态的曼哈顿距离;

3. 通过 1, 2 时间空间上已经比较优化了,但是还可以再对其进行剪枝。对于 x,如果在所在行移动,逆序数对不变,如果在列上移动,逆序数对成偶数加减,根据这点来判断某个状态是否有可行解的潜质,最终代码优化到了 16ms. 不过还是在HDOJ 1043 WA;

http://zh.wikipedia.org/wiki/%E5%BA%B7%E6%89%98%E5%B1%95%E5%BC%80  康托展开

http://zh.wikipedia.org/wiki/%E6%9B%BC%E5%93%88%E9%A0%93%E8%B7%9D%E9%9B%A2 曼哈顿距离

 

#include <iostream>
#include <algorithm>
#include <queue>
#include <vector>
using namespace std;

const int MAXN = 362880 + 10;
const int fac[15] = {1, 1, 2, 6, 24, 120, 720, 5040, 40320, 362880};
const char dir[5] = {0, 'u', 'd', 'l', 'r'};
const int dis[9][2] = {{0,0}, {0,1}, {0,2}, {1,0}, {1,1}, {1,2}, {2,0}, {2,1}, {2,2}};

bool vis[MAXN];
int F[MAXN], ENDS;
char OP[MAXN];

struct ST {
    char e[10];
    int f, step;
    bool operator < (const ST& other) const {
        return f > other.f;
    }
};

int getvalue(const char e[]) {
    int ans = 0, count;
    for (int i = 0; i < 8; i++) {
        count = 0;
        for (int j = i + 1; j < 9; j++) {
            if (e[i] > e[j]) count += 1;
        }
        ans += fac[8-i] * count;
    }
    return ans;
}

bool challenge(ST& s, int pos, int op) {
    if (op == 1) {
        // up
        if (pos >= 3) {
            swap(s.e[pos-3], s.e[pos]);
            return true;
        }
    } else if (op == 2) {
        // down
        if (pos < 6) {
            swap(s.e[pos], s.e[pos+3]);
            return true;
        }
    } else if (op == 3) {
        // left
        if (pos % 3) {
            swap(s.e[pos-1], s.e[pos]);
            return true;
        }
    } else if (op == 4) {
        // right
        if ((pos + 1) % 3) {
            swap(s.e[pos], s.e[pos+1]);
            return true;
        }
    }
    return false;
}

int getdiff(const char e[]) {
    int ret = 0;
    for (int i = 0; i < 3; i++) {
        for (int j = 0; j < 3; j++) {
            int k = 3 * i + j;
            if (e[k] != 9) {
                ret += abs(i - dis[e[k]-1][0]) + abs(j - dis[e[k]-1][1]);
            }
        }
    }
    return ret;
}

bool judge(const char e[]) {
    int sum = 0;
    for (int i = 0; i < 9; i++) {
        if (e[i] != 9) {
            for (int j = i + 1; j < 9; j++)
                if (e[i] > e[j])  sum += 1;
        }
    }
    return !(sum & 0x01);
}

bool bfs(ST& s) {
    priority_queue<ST> Q;

    s.step = 0;
    s.f = s.step + getdiff(s.e);

    Q.push(s);
    
    int vals = getvalue(s.e);
    F[vals] = -1;
    OP[vals] = 0;
    vis[vals] = true;

    while (!Q.empty()) {
        ST u = Q.top();
        Q.pop();

        int pos;
        for (int i = 0; i < 9; i++)
            if (u.e[i] == 9) pos = i;

        for (int i = 1; i <= 4; i++) {
            ST v = u;
            if (challenge(v, pos, i)) {
                int state = getvalue(v.e);
                if (!vis[state] && judge(v.e)) {
                    vis[state] = true;
                    v.step = u.step + 1;
                    v.f = v.step + getdiff(v.e);
                    Q.push(v);

                    F[state] = getvalue(u.e);
                    OP[state] = dir[i];

                    if (state == ENDS)
                        return true;
                }
            }
        }
    }
    return false;
}

int main() {
    ST u, E;
    for (int i = 0; i < 9; i++)
        E.e[i] = i + 1;
    ENDS = getvalue(E.e);

    char ch[100];
    while (gets(ch)) {
        for (int i = 0, j = 0; ch[i]; i++) {
            if (ch[i] != ' ') {
                if (ch[i] == 'x')
                    u.e[j++] = 9;
                else
                    u.e[j++] = ch[i] - '0';
            }
        }

        if (getvalue(u.e) == ENDS) {
            //printf("lr\n");
            continue;
        }

        int j = 0;
        char opcode[1000];

        memset(vis, false, sizeof(vis));
        memset(F, -1, sizeof(F));
        memset(OP, 0, sizeof(OP));

        if (judge(u.e) && bfs(u)) {
            for (int i = ENDS; i != -1; i = F[i])
                opcode[j++] = OP[i];
            for (int i = j-2; i >= 0; i--)
                printf("%c", opcode[i]);
            printf("\n");
        } else {
            printf("unsolvable\n");
        }
    }
    return 0;
}

转载于:https://www.cnblogs.com/kedebug/archive/2013/03/23/2977232.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 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)$。 参考代码:
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值