Problem:
氷の迷路

四角のマスを縦横にならべた長方形の迷路があります。この迷路では東西南北の隣接するマスへ移動しながら、スタートのマスSを出発し、ゴールのマスGを目指します。マスの種類には、平原、山、氷の３種類があります。SとGは、平原のマスに置かれています。平原のマスには移動できますが、山のマスには移動できません。氷のマスは移動できますが、条件によって氷が割れ動けなくなります。

氷のマスは東西南北に隣接するもの全体でひとつの塊として扱われます。

氷の塊のなかで、マスの個数の半分より多くのマスに移動すると、塊全体が割れます。

例えば、図１は与えられた迷路の最短経路のステップ数が 11 であることを示しています。

図1

しかし、図２のように氷のマスを通って近道をしようとすると、大きさが６の氷の塊に対する４回目の移動後に動けなくなってしまいGには到達できません。

図2

このような迷路の情報を入力とし、SからGまでの最短経路のステップ数を求めるプログラムを作成してください。マスの種類はそれぞれ以下の文字で表されます：

文字（半角）

マスの種類

. （ピリオド）

平原

#（シャープ）

山

X

氷

与えられる迷路は必ず解けるものとします。迷路はその東西方向のマスの数 X、南北方向のマスの数 Y と X × Y 個の文字として与えられます。

入力

複数のデータセットが与えられます。入力の終わりはゼロふたつの行で示されます。各データセットは以下の形式で与えられます。

X Y
line1
line2
:
lineY

1行目に迷路の大きさを表す整数 X, Y (2 ≤ X, Y ≤ 12) が与えられます。続く Y 行に迷路の i 行目の情報 linei (長さ X の半角英文字列) が与えられます。

データセットの数は 40 を超えない。

出力

データセットごとに、最小のステップ数を１行に出力します。

入力例

5 5
.X.S.
.X#..
.XX##
.#XG.
..X..
7 3
SXX.XXG
X.#.#X.
XXX.XX#
4 4
S...
X.X.
GX..
...X
10 10
..XXXXX.XX
.X.#.#X.XX
SX.#X.X..X
#X.##.X.XX
..XXXX#.XX
##.##.##XX
....X.XX#X
.##X..#X#X
....XX#..X
...#XXG..X
0 0

出力例

11
10
10
33

Solution in C++ (judge verdict: Runtime Error)

#include <algorithm>
#include <cmath>
#include <climits>
#include <cstdio>
#include <cstring>
#include <fstream>
#include <iostream>
#include <list>
#include <map>
#include <queue>
#include <set>
#include <sstream>
#include <stack>
#include <string>
#include <vector>
 
using namespace std;
 
//#define LOG(...) printf(__VA_ARGS__)
#define LOG(...)
#define FOR(i,a,b) for(int i=(a);i<=(b);++i)
#define REP(i,n) for(int i=0;i<(int)(n);++i)
#define ALL(a) (a).begin(),(a).end()
#define EXIST(s,e) ((s).find(e)!=(s).end())
#define SORT(c) sort((c).begin(),(c).end())
 
typedef long long ll;
typedef unsigned long long ull;
typedef vector<bool> vb;
typedef vector<int> vi;
typedef vector<ll> vll;
typedef vector<vb> vvb;
typedef vector<vi> vvi;
typedef vector<vll> vvll;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;

int w,h;
vector<string> field;

int dx[]={-1,0,1,0};
int dy[]={0,-1,0,1};

int dfs(int x, int y, int n){
	int cnt=1;
	field[y][x]=-n;
	REP(i, 4){
		int sx=x+dx[i];
		int sy=y+dy[i];
		if(0<=sx&&sx<w&&0<=sy&&sy<h&&field[sy][sx]=='X')
			cnt+=dfs(sx, sy, n);
	}
	return cnt;
}

struct P
{
	int x,y,cost;
	vi iceList,l;
};
 
int main() {
	while(cin >> w >> h, w|h){
		field.clear();
		REP(y, h){
			string line;
			cin >> line;
			field.push_back(line);
		}

		P s,g;
		vi l;
		REP(y, h){REP(x, w){
			switch(field[y][x]){
				case 'X':
					l.push_back(dfs(x,y,l.size()+1)/2);
					break;
				case 'S':
					s=(P){x,y,0};
					break;
				case 'G':
					g=(P){x,y};
					break;
			}
		}}

		REP(y, h){REP(x, w){
			LOG("%3d ",field[y][x]);
		}
		LOG("\n");
		}
		s.iceList=l;

		REP(i, l.size()){
			LOG("ice life %d\n",l[i]);
		}

		LOG("gool %d %d \n",g.x,g.y);

		queue<P> que;
		que.push(s);
		while(!que.empty()){
			LOG("---\n");
			P p=que.front();que.pop();
			LOG("now %d %d %d\n",p.x,p.y,p.cost);
			if(p.y==g.y&&p.x==g.x){
				cout << p.cost << endl;
				break;
			}
			REP(i, 4){
				vi l=p.iceList;
				vi ll=p.l;
				int sx=p.x+dx[i];
				int sy=p.y+dy[i];
				if(!(0<=sx&&sx<w&&0<=sy&&sy<h))
					continue;
				int d=field[sy][sx];
				bool t=false;
				REP(j, p.l.size()){
					if(p.l[j]==sy*w+sx){
						t=true;
						break;
					}
				}
				if(t)continue;
				ll.push_back(sy*w+sx);
				if(d<0){
					d=-d-1;
					LOG("ice %d = %d\n",d,l[d]);
					if(l[d]>0){
						l[d]--;
						que.push((P){sx,sy,p.cost+1,l,ll});
					}
				}else if(d=='.'||d=='G')
					que.push((P){sx,sy,p.cost+1,l,ll});
			}
		}
	}

	return 0;
}