孤岛营救问题

【题目描述】

1944 年，特种兵麦克接到国防部的命令，要求立即赶赴太平洋上的一个孤岛，营救被敌军俘虏的大兵瑞恩。瑞恩被关押在一个迷宫里，迷宫地形复杂，但幸好麦克得到了迷宫的地形图。迷宫的外形是一个长方形，其南北方向被划分为 n𝑛 行，东西方向被划分为 m𝑚 列， 于是整个迷宫被划分为 n×m𝑛×𝑚 个单元。每一个单元的位置可用一个有序数对 (单元的行号, 单元的列号) 来表示。南北或东西方向相邻的 22 个单元之间可能互通，也可能有一扇锁着的门，或者是一堵不可逾越的墙。迷宫中有一些单元存放着钥匙，并且所有的门被分成 p𝑝 类， 打开同一类的门的钥匙相同，不同类门的钥匙不同。

大兵瑞恩被关押在迷宫的东南角，即 (n,m𝑛,𝑚) 单元里，并已经昏迷。迷宫只有一个入口， 在西北角。也就是说，麦克可以直接进入 (1,11,1) 单元。另外，麦克从一个单元移动到另一个相邻单元的时间为 11，拿取所在单元的钥匙的时间以及用钥匙开门的时间可忽略不计。

试设计一个算法，帮助麦克以最快的方式到达瑞恩所在单元，营救大兵瑞恩。

【输入】

第一行有三个整数,分别表示 n,m,p𝑛,𝑚,𝑝 的值。

第二行是一个整数k𝑘，表示迷宫中门和墙的总数。

第 i+2𝑖+2 行 (1≤i≤k1≤𝑖≤𝑘)，有 55 个整数，依次为 xi1,yi1,xi2,yi2,gi𝑥𝑖1,𝑦𝑖1,𝑥𝑖2,𝑦𝑖2,𝑔𝑖:当 gi≥1𝑔𝑖≥1 时，表示 (xi1,yi1𝑥𝑖1,𝑦𝑖1)单元与 (xi2,yi2𝑥𝑖2,𝑦𝑖2) 单元之间有一扇第 gi𝑔𝑖 类的门，当 gi=0𝑔𝑖=0 时， 表示 (xi1,yi1𝑥𝑖1,𝑦𝑖1) 单元与 (xi2,yi2𝑥𝑖2,𝑦𝑖2) 单元之间有一堵不可逾越的墙。

第 k+3𝑘+3 行是一个整数 s𝑠，表示迷宫中存放的钥匙总数。

第 k+3+j𝑘+3+𝑗 行 (1≤j≤s1≤𝑗≤𝑠) ，有 33 个整数，依次为 xi1,yi1,qi𝑥𝑖1,𝑦𝑖1,𝑞𝑖​​​​ ，表示第 j𝑗 把钥匙存放在 (xi1,yi1𝑥𝑖1,𝑦𝑖1) 单元里，并且第 j𝑗 把钥匙是用来开启第 qi𝑞𝑖 类门。

输入数据中同一行各相邻整数之间用一个空格分隔。

【输出】

输出麦克营救到大兵瑞恩的最短时间。如果问题无解，则输出 −1。

【输入样例】

4 4 9
9
1 2 1 3 2
1 2 2 2 0
2 1 2 2 0
2 1 3 1 0 
2 3 3 3 0
2 4 3 4 1
3 2 3 3 0
3 3 4 3 0
4 3 4 4 0
2
2 1 2 
4 2 1

【输出样例】

14

【提示】

|xi1−xi2|+|yi1−yi2|=1,0≤gi≤p|𝑥𝑖1−𝑥𝑖2|+|𝑦𝑖1−𝑦𝑖2|=1,0≤𝑔𝑖≤𝑝

1≤qi≤p1≤𝑞𝑖≤𝑝

n,m,p≤10,k<150

code:

#include<bits/stdc++.h>
using namespace std;
const int N = 12;
int dir[5][3]={{0,1},{1,0},{0,-1},{-1,0}};
int e[N][N][N][N],key[N][N][N],dkey[N][N];
int vis[N][N][1<<14];
int n,m;
struct node{
	int x,y,k,d;
	node(){}
	node(int x,int y,int k,int d):x(x),y(y),k(k),d(d){}
};
int getkey(int x,int y){
	int ans = 0;
	for(int i=1;i<=dkey[x][y];i++)
	ans|=(1<<(key[x][y][i]-1));
	return ans;
}
int bfs(){
	queue<node> q;
	int sk = getkey(1,1);
	q.push(node(1,1,sk,0)),vis[1][1][sk] = 1;
	while(q.size()){
		node u = q.front();q.pop();
		if(u.x==n && u.y==m) return u.d;
		int ux = u.x,uy = u.y;
		for(int i=0;i<4;i++){
			int nx = ux+dir[i][0];
			int ny = uy+dir[i][1];
			int opt = e[ux][uy][nx][ny];
			if(nx <1||ny<1||nx>n||ny>m||opt<0||(opt&&!(u.k&(1<<(opt-1))))) continue;
			int nkey = u.k|getkey(nx,ny);
			if(vis[nx][ny][nkey]) continue;
			q.push(node(nx,ny,nkey,u.d+1));
			vis[nx][ny][nkey] = 1;
		}
	}
	return -1;
}
void read(){
	int k,s,r;
	scanf("%d%d%d",&n,&m,&r);
	for(scanf("%d",&k);k--;){
		int x1,y1,x2,y2,g;
        scanf("%d%d%d%d%d",&x1,&y1,&x2,&y2,&g);
        if(g) e[x1][y1][x2][y2]=e[x2][y2][x1][y1]=g;
        else e[x1][y1][x2][y2]=e[x2][y2][x1][y1]=-1;
	}
   for(scanf("%d",&s);s--;) {
        int x,y,q;
        scanf("%d%d%d",&x,&y,&q);
        key[x][y][++dkey[x][y]]=q;
    }
} 
int main(){
	read();
    printf("%d\n",bfs());
	return 0;
}