POJ 3131 Cubic Eight-Puzzle 双向BFS + HASH

超赞的一道搜索题，麻烦到没朋友，不过思路还是很清楚的。

双向搜索的时候需要注意一个地方就是起始状态只有一个，目标状态最多的时候感觉会有512个，所以要控制两个方向搜索的层数。第一次知道双向搜索还能控制层数，简直点赞。

#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <queue>
#include <cmath>
#include <stack>
#include <map>

#pragma comment(linker, "/STACK:1024000000");
#define EPS (1e-8)
#define LL long long
#define ULL unsigned long long
#define _LL __int64
#define _INF 0x3f3f3f3f
#define Mod 6000007

using namespace std;

char Map[5][5];

int Judge(char a,char b,char c)
{
    if(a == 'W')
    {
        if(b == 'R')
            return 0;
        return 5;
    }
    else if(a == 'R')
    {
        if(b == 'W')
            return 1;
        return 3;
    }
    else if(a == 'B')
    {
        if(b == 'R')
            return 2;
        return 4;
    }
    else return 6;
}

struct Q
{
    int sta,x,y,fl,ans;
};

int mark[3000010][2];

int jx[] = { 0, 0, 1,-1};
int jy[] = {-1, 1, 0, 0};

int re[10][10];

int bit[10];

struct H
{
    int site,data,next;
}st[Mod];

int head[Mod];

int Top_H;

int Query(int data,int Top)
{
    int temp = data%Mod;

    int p = head[temp];
    for(;p != -1; p = st[p].next)
    {
        if(st[p].data == data)
            return st[p].site;
    }

    p = Top_H++;
    st[p].data = data;
    st[p].site = Top;
    st[p].next = head[temp];
    head[temp] = p;
    return -1;
}

int bfs(int x,int y)
{
    Top_H = 0;
    memset(mark,-1,sizeof(mark));

    int Top = 0;

    int i,j;

    for(i = 0;i < Mod; ++i)
        head[i] = -1;

    int sta = 0;

    for(i = 0; i < 3; ++i)
    {
        for(j = 0; j < 3; ++j)
            if(i == x && j == y)
                sta += Judge('E','E','E')*bit[i*3+j+1];
            else
                sta += Judge('W','B','R')*bit[i*3+j+1];
    }

    Q f,t;
    queue<Q> q;

    f.sta = sta;
    f.ans = 0,f.x = x,f.y = y,f.fl = 0;
    q.push(f);

    Query(f.sta,Top);

    mark[Top][0] = 0;
    Top++;

    f.ans = 0;
    f.fl = 1;

    for(i = 0; i < 512; ++i)
    {
        sta = 0,x = 0,y = 0;
        for(j = 0; j < 9; ++j)
        {
            if(Map[x][y] == 'E')
            {
                f.x = x;
                f.y = y;
                sta += Judge('E','E','E')*bit[x*3+y+1];
            }
            else if(i&(1<<j))
            {
                if(Map[x][y] == 'W')
                    sta += Judge('W','R','B')*bit[x*3+y+1];
                else if(Map[x][y] == 'R')
                    sta += Judge('R','B','W')*bit[x*3+y+1];
                else
                    sta += Judge('B','R','W')*bit[x*3+y+1];
            }
            else
            {
                if(Map[x][y] == 'W')
                    sta += Judge('W','B','R')*bit[x*3+y+1];
                else if(Map[x][y] == 'R')
                    sta += Judge('R','W','B')*bit[x*3+y+1];
                else
                    sta += Judge('B','W','R')*bit[x*3+y+1];
            }
            y++;
            if(y == 3)
                y = 0,x++;
        }
        f.sta = sta;
        q.push(f);
        if(Query(f.sta,Top) == -1)
        {
            mark[Top++][1] = 0;
        }
        else
        {
            mark[Query(f.sta,Top)][1] = 0;
        }
    }

    int site,tsite,tdata;

    while(q.empty() == false)
    {
        f = q.front();
        q.pop();

        if((f.fl == 1 && f.ans > 9) || (f.fl == 0 && f.ans > 20))
            continue;

        if(mark[Query(f.sta,Top)][f.fl^1] != -1)
        {
            return mark[Query(f.sta,Top)][f.fl^1] + f.ans;
        }
        site = f.x*3 + f.y;

        t.fl = f.fl;
        t.ans = f.ans+1;

        for(i = 0; i < 4; ++i)
        {
            t.x = f.x + jx[i];
            t.y = f.y + jy[i];

            if(0 <= t.x && t.x <= 2 && 0 <= t.y && t.y <= 2)
            {
                tsite = t.x*3+t.y;
                tdata = f.sta/bit[tsite+1]%7;

                t.sta = f.sta/bit[tsite]*bit[tsite] + 6*bit[tsite+1] + f.sta%bit[tsite+1];
                t.sta = t.sta/bit[site]*bit[site] + re[tdata][i]*bit[site+1] + t.sta%bit[site+1];

                if(Query(t.sta,Top) == -1)
                {
                    mark[Top++][t.fl] = t.ans;
                    q.push(t);
                }
                else if(mark[Query(t.sta,Top)][t.fl] == -1)
                {
                    mark[Query(t.sta,Top)][t.fl] = t.ans;
                    q.push(t);
                }
            }
        }
    }
    return -1;
}

int main()
{
    re[0][0] = 1;
    re[0][1] = 1;
    re[0][2] = 2;
    re[0][3] = 2;

    re[1][0] = 0;
    re[1][1] = 0;
    re[1][2] = 4;
    re[1][3] = 4;

    re[2][0] = 3;
    re[2][1] = 3;
    re[2][2] = 0;
    re[2][3] = 0;

    re[3][0] = 2;
    re[3][1] = 2;
    re[3][2] = 5;
    re[3][3] = 5;

    re[4][0] = 5;
    re[4][1] = 5;
    re[4][2] = 1;
    re[4][3] = 1;

    re[5][0] = 4;
    re[5][1] = 4;
    re[5][2] = 3;
    re[5][3] = 3;

    int i,j,x,y;

    bit[9] = 1;

    for(i = 8; i >= 0; --i)
        bit[i] = bit[i+1]*7;

    while(scanf("%d %d",&y,&x) && (x||y))
    {

        for(i = 0; i < 3; ++i)
        {
            for(j = 0; j < 3; ++j)
            {
                scanf("%*c%c",&Map[i][j]);
            }
        }

        printf("%d\n",bfs(x-1,y-1));
    }

    return 0;
}