Problem:
Problem G: The Morning after Halloween

You are working for an amusement park as an operator of an obakeyashiki, or a haunted house,
in which guests walk through narrow and dark corridors. The house is proud of their lively
ghosts, which are actually robots remotely controlled by the operator, hiding here and there in
the corridors. One morning, you found that the ghosts are not in the positions where they are
supposed to be. Ah, yesterday was Halloween. Believe or not, paranormal spirits have moved
them around the corridors in the night. You have to move them into their right positions before
guests come. Your manager is eager to know how long it takes to restore the ghosts.

In this problem, you are asked to write a program that, given a floor map of a house, finds the
smallest number of steps to move all ghosts to the positions where they are supposed to be.

A floor consists of a matrix of square cells. A cell is either a wall cell where ghosts cannot move
into or a corridor cell where they can.

At each step, you can move any number of ghosts simultaneously. Every ghost can either stay
in the current cell, or move to one of the corridor cells in its 4-neighborhood (i.e. immediately
left, right, up or down), if the ghosts satisfy the following conditions:

No more than one ghost occupies one position at the end of the step.

No pair of ghosts exchange their positions one another in the step.

For example, suppose ghosts are located as shown in the following (partial) map, where a sharp sign ('#) represents a wall cell and 'a', 'b', and 'c' ghosts.

####
ab#
#c##
####

The following four maps show the only possible positions of the ghosts after one step.

#### #### #### ####
ab# a b# acb# ab #
#c## #c## # ## #c##
#### #### #### ####

Input

The input consists of at most 10 datasets, each of which represents a floor map of a house. The format of a dataset is as follows.

w h n
c11c12...c1w
c21c22...c2w
. . .. .
.. .
.
.. .
ch1ch2...chw

w, h and n in the first line are integers, separated by a space. w and h are the floor width
and height of the house, respectively. n is the number of ghosts. They satisfy the following constraints.

4 ≤ w ≤ 16

4 ≤ h ≤ 16

1 ≤ n ≤ 3

Subsequent h lines of w characters are the floor map. Each of cij is either:

a '#' representing a wall cell,

a lowercase letter representing a corridor cell which is the initial position of a ghost,

an uppercase letter representing a corridor cell which is the position where the ghost corresponding to its lowercase letter is supposed to be, or

a space representing a corridor cell that is none of the above.

In each map, each of the first n letters from a and the first n letters from A appears once and
only once. Outermost cells of a map are walls; i.e. all characters of the first and last lines are
sharps; and the first and last characters on each line are also sharps. All corridor cells in a
map are connected; i.e. given a corridor cell, you can reach any other corridor cell by following
corridor cells in the 4-neighborhoods. Similarly, all wall cells are connected. Any 2 × 2 area on
any map has at least one sharp. You can assume that every map has a sequence of moves of ghosts that restores all ghosts to the positions where they are supposed to be.

The last dataset is followed by a line containing three zeros separated by a space.

Output

For each dataset in the input, one line containing the smallest number of steps to restore ghosts
into the positions where they are supposed to be should be output. An output line should not
contain extra characters such as spaces.

Sample Input

5 5 2
#####
#A#B#
# #
#b#a#
#####
16 4 3
################
## ########## ##
# ABCcba #
################
16 16 3
################
### ## # ##
## # ## # c#
# ## ########b#
# ## # # # #
# # ## # # ##
## a# # # # #
### ## #### ## #
## # # # #
# ##### # ## ##
#### #B# # #
## C# # ###
# # # ####### #
# ###### A## #
# # ##
################
0 0 0

Output for the Sample Input

7
36
77

Solution in C++11 (judge verdict: Time Limit Exceeded)

#include "bits/stdc++.h"
#include<unordered_map>
#include<unordered_set>
#pragma warning(disable:4996)
using namespace std;
using ld = long double;
template<class T>
using Table = vector<vector<T>>;
const ld eps=1e-9;

//// < "D:\D_Download\Visual Studio 2015\Projects\programing_contest_c++\Debug\a.txt"


int memo[15][15][15][15][15][15];
int dist[16][16][16][16];
const int dx[] = { -1,0,1,0,0 };
const int dy[] = { 0,1,0,-1,0 };
int W, H, N;

vector<int>gx;
vector<int>gy;
struct aa {
	vector<int>xs;
	vector<int>ys;
	int time;
	int geta()const {
		int a = time;
		int amin = 0;
		for (int i = 0; i < N; ++i) {
			amin=max(amin, dist[ys[i]][xs[i]][gy[i]][gx[i]]);
		}
		return a+amin;
	}
};
class Compare {
public:
	//aa?????????????????¶
	bool operator()(const aa&l, const aa&r) {
		return l.geta() > r.geta();
	}
};
int getmemo(vector<int>xs, vector<int>ys) {
	int an = xs.size();
	for (int i = 0; i < 3 - an; ++i) {
		xs.push_back(0);
		ys.push_back(0);
	}
	return memo[xs[0]][xs[1]][xs[2]][ys[0]][ys[1]][ys[2]];
}
void setmemo(vector<int>xs, vector<int>ys, const int value) {
	int an = xs.size();
	for (int i = 0; i < 3 - an; ++i) {
		xs.push_back(0);
		ys.push_back(0);
	}
	memo[xs[0]][xs[1]][xs[2]][ys[0]][ys[1]][ys[2]] = value;
}
bool check(const vector<int>&axs, const vector<int>&ays, vector<int>&bxs, vector<int>&bys) {
	int time;
	if (bxs.size() == 1)time = 0;
	else if (bxs.size() == 2)time = 1;
	else time = 3;
	for (int i = 0; i < time; ++i) {
		int a = i;
		int b = (i + 1) % 3;
		if (axs[a] == bxs[b] && ays[a] == bys[b]&&axs[b]==bxs[a]&&ays[b]==bys[a]) {
			return false;
		}
		if (bxs[a] == bxs[b] && bys[a] == bys[b]) {
			return false;
		}
	}
	return true;
}
int main() {
	while (1) {
		for (int i = 0; i < 15; ++i) {
			for (int j = 0; j < 15; ++j) {
				for (int k = 0; k < 15; ++k) {
					for (int l = 0; l < 15; ++l) {
						dist[i][j][k][l] = 1e9;
						for (int m = 0; m < 15; ++m) {
							for (int n = 0; n < 15; ++n) {
								memo[i][j][k][l][m][n] = 1e9;
							}
						}
					}
				}
			}
		}cin >> W >> H >> N;
		if (!W)break;
		vector<vector<int>>field(H, vector<int>(W));
		vector<int>sx(3);
		vector<int>sy(3);
		gx.clear();
		gy.clear();
		gx.resize(N);
		gy.resize(N);
		string st;
		getline(cin, st);
		for (int i = 0; i < H; ++i) {
			getline(cin, st);
			for (int j = 0; j < W; ++j) {
				if (st[j] == '#') {
					field[i][j] = true;
				}
				else if (st[j] >= 'a'&&st[j] <= 'c') {
					sx[st[j] - 'a'] = j;
					sy[st[j] - 'a'] = i;
				}
				else if (st[j] == ' ') {

				}
				else {
					gx[st[j] - 'A'] = j;
					gy[st[j] - 'A'] = i;
				}
			}
		}
		for (int i = 0; i < H - 1; ++i) {
			for (int j = 0; j < W; ++j) {
				if (!field[i][j] && !field[i + 1][j]) {
					dist[i][j][i + 1][j] = 1;
					dist[i+1][j][i][j] = 1;
				}
			}
		}
		for (int i = 0; i < H; ++i) {
			for (int j = 0; j < W-1; ++j) {
				if (!field[i][j] && !field[i][j+1]) {
					dist[i][j][i][j + 1] = 1;
					dist[i][j+1][i][j] = 1;
				}
			}
		}
		for (int i = 0; i < 256; ++i) {
			dist[i % 16][i / 16][i % 16][i / 16] = 0;
			
		}
		for (int k = 0; k < 256; ++k) {
			for (int i = 0; i < 256; ++i) {
				for (int j = 0; j < 256; ++j) {
					dist[i % 16][i / 16][j % 16][j / 16] = min(dist[i % 16][i / 16][j % 16][j / 16], dist[i % 16][i / 16][k % 16][k / 16] + dist[k % 16][k / 16][j % 16][j / 16]);
				}
			}
		}
		priority_queue<aa,vector<aa>,Compare>que;
		que.push(aa{ sx,sy,0 });
		setmemo(sx, sy, 0);
		int ans = -1;
		while (!que.empty()) {
			aa atop(que.top());
			//cout << atop.geta();
			vector<int>nowxs(atop.xs);
			vector<int>nowys(atop.ys);
			if (nowxs == gx&&nowys == gy) {
				while (!que.empty()) {
					aa a =que.top();
					que.pop();
					//cout << a.geta() << endl;
				}
				ans = atop.time;
				break;
			}
			que.pop();
			for (int i = 0; i < 125; ++i) {
				int n(i);
				vector<int>nexxs, nexys;
				bool aok = true;
				for (int j = 0; j < N; ++j) {
					const int way = n % 5;
					n /= 5;
					const int nexx = nowxs[j] + dx[way];
					const int nexy = nowys[j] + dy[way];
					if (field[nexy][nexx])aok = false;
					nexxs.push_back(nowxs[j] + dx[way]);
					nexys.push_back(nowys[j] + dy[way]);
				}
				if (!check(nowxs, nowys, nexxs, nexys))aok = false;
				if (aok) {
					if (getmemo(nexxs, nexys) > atop.time + 1) {
						setmemo(nexxs, nexys, atop.time + 1);
						que.push(aa{ nexxs,nexys,atop.time + 1 });
					}
				}
			}

		}
		cout << ans << endl;

	}
	return 0;
}