【图论最短路】【CQBZOJ 2432】Maze

Maze问题 H(2432): Maze
时间限制: 1 Sec 内存限制: 256 MB

题目描述
众所周知(怎么又是众所周知)，仙剑的迷宫是难走的要命，某人就在仙四的女罗岩困了很长时间。我们可以把女罗岩的地图抽象成n*n的地图，我们现在在((1,1)处，而出口在(n, n)。每次行动只能向上下左右移动一格。图中有M个机关，只有打开某个机关之后，与该机关相对应的地方才可以行走。当然，地图中还会有小怪兽，他们能够监视着他所在地区以及上下左右共五个方格，我们不愿意与他们战斗(为了节省时间)，所以他们能监视到的地方不可以行走。同样的，地图上的障碍物处也不能经过。我们需要求出从起点到终点的最少步数。

输入
第1行，两个整数N,M。表示地图大小为M∗ N，机关的数量为M。
第2~N + 1行，每行N个整数，每个整数 i 可能是0，−1，−2或者一个正整数。i =0表示该位置为一块空地，i = −1表示该位置为一个障碍物，i = −2表示该位置为一个小怪兽。如果i是一个属于[1,M]的正整数，则表示该位置为一个机关，其编号为i。如果i是一个大于M的正整数，则表示该位置为一个机关控制区域，它由编号为i ~ M的机关控制。

输出
一个整数，为走出迷宫所需的最少的步数。

样例输入
6 2
0 0 0 -2 -1 2
-1 0 0 0 -1 0
-2 0 0 0 3 3
-2 0 0 -1 -1 4
0 -1 0 0 -1 0
1 0 0 0 -1 0

样例输出
24

提示
地图如下图，S为入口，T为目标，黑色的单元格为障碍物。每个E表示一个小怪兽，(E)为小怪兽的监视范围。K1表示机关1，K2表示机关2。D1表示机关为1的控制区域，D2表示机关为2的控制区域。
最优的路线为(1,1) → (1,2) → (2,2) → (2,3) → (3,3) → (4,3) → (5,3) → (6,3) →(6,2) → (6,1) → (6,2) → (6,3) → (5,3) → (4,3) → (3,3) → (3,4) → (3,5) →(3,6) → (2,6) → (1,6) → (2,6) → (3,6) → (4.6) → (5,6) → (6,6)

【数据规模】
1 ≤ N ≤ 50
0 ≤ M ≤ 16

简单的拆点Dijsktra。我们将最短路状态拆成$dist[pos][keys]$表示当前在pos节点，拥有钥匙的状态为keys的最短路径。跑最短路的时候注意不能去的节点的预处理和开门的判定。

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<queue>
using namespace std;

#define MAXN 50
#define MAXM 65536
#define INF 0x3f3f3f3f
typedef long long int LL;

const int dir[4][2]={{-1,0},{1,0},{0,-1},{0,1}};

int N,M,S;
int G[MAXN+5][MAXN+5];

struct point
{
    int x,y,s;
    int w;
    point(){}
    point(int a,int b,int c,int d):x(a),y(b),s(c),w(d){}
};

bool vis[MAXN+1][MAXN+1][MAXM+10];
queue<point>que;
int bfs(point S)
{
    memset(vis,0,sizeof(vis));
    while(!que.empty())que.pop();

    vis[S.x][S.y][S.s]=1;
    que.push(S);

    point now;
    int nx,ny,ns,nw;
    while(!que.empty())
    {
        now=que.front();que.pop();

        if(now.x==N&&now.y==N)return now.w;

        for(int i=0;i<4;++i)
        {
            nx=now.x+dir[i][0],ny=now.y+dir[i][1];

            if(nx<1||ny<1||nx>N||ny>N||G[nx][ny]==-1)continue;
            if(M<G[nx][ny]&&((now.s&(1<<(G[nx][ny]-M-1)))==0))continue;

            ns=now.s|(0<G[nx][ny]&&G[nx][ny]<=M?(1<<(G[nx][ny]-1)):0);
            nw=now.w+1;

            if(!vis[nx][ny][ns])
            {
                vis[nx][ny][ns]=1;
                que.push(point(nx,ny,ns,nw));
            }
        }
    }

    return -1;
}

int main()
{
    scanf("%d%d",&N,&M);
    int i,j;
    int tmp;

    for(i=1;i<=N;++i)
        for(j=1;j<=N;++j)
        {
            scanf("%d",&tmp);
            if(tmp==-2)
            {
                for(int k=0;k<4;++k)
                    G[i+dir[k][0]][j+dir[k][1]]=-1;
                G[i][j]=-1;
            }
            else if(G[i][j]!=-1)G[i][j]=tmp;
        }

    printf("%d\n",bfs(point(1,1,(0<G[1][1]&&G[1][1]<=M?(1<<(G[1][1]-1)):0),0)));
}