Problem:
パチモンクリーチャー

某国で大人気のゲーム、パチモンクリーチャーが日本でリメイクされて発売されました。 ゲームが大好きなあなたは、 このゲームを何度もプレイするうちにどうしたら最速でクリアできるのか考えるようになりました。しかし、いくら考えても最速の攻略方法がわからなかったあなたは、どれだけ早くゲームをクリアできるかを求めるプログラムを作成することにしました。

ゲームの詳細は以下の通りです。

パチモンクリーチャー(以下、パチクリ)という生物が多く存在する世界がゲームの舞台です。各パチクリは、火属性、氷属性、木属性、土属性、水属性の 5 種類の属性のいずれか１つの属性を持ちます。ゲームの主人公は、ゲーム開始時に好きな属性のパチクリ一匹を冒険のパートナーとして選びます。そのパチクリと共にゴールを目指し、ゴールにいるライバルを倒してパチクリマスターになることがゲームの目的です。

しかし、ライバルを倒すためには全属性のパチクリがいないと勝てないので、途中で全属性のパチクリを捕まえなければなりません。パチクリを捕まえるには属性がカギとなります。火属性のパチクリは氷属性のパチクリを捕まえることができ、同様に、氷属性は木属性、木属性は土属性、土属性は水属性、水属性は火属性を捕まえることができます。属性の関連は以下の図のようになります。

以下の図はゲームが行われるマップの一例を表しています。

主人公はパチクリを一匹もってスタート地点である「S」から出発し、一マスずつ移動しながらゴール地点である「G」を目指します。その途中で、最初に持っているパチクリ以外の 4 つの属性のパチクリを捕まえ、ゴール地点であるマス目に移動するとゲーム終了となります。

主人公は、今いるマス目から、辺を共有する隣のマス目に移動することができ、それを一回の移動と数えます。主人公がパチクリのいるマスに移動した場合、そのパチクリを捕まえられる属性のパチクリを持っていればそのパチクリを捕まえたことになります。そのマスにいるパチクリを捕まえられるかの可否にかかわらず、すべてのマスに何度でも移動することができます。

マップの大きさ(横方向の列数、縦方向の行数)とマップの初期状態を入力とし、初めに選んだパチクリの属性と、それ以外の 4 つの属性のパチクリを捕まえるのにかかる、スタート地点からゴール地点に至る最小移動数を出力するプログラムを作成してください。

Input

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

W H
c11c12...c1W
c21c22...c2W
:
cH1cH2...cHW

1 行目にマップの横方向の列数 W と縦方向の行数 H (2 ≤ W, H ≤ 1000) が与えられます。続く H 行にマップの i 行目の情報が与えられます。入力されるマップには各マスの状態が与えられます。 「S」は主人公のスタート地点を、 「G」はゴール地点を、「1」「2」「3」「4」「5」はそこにいるパチクリの属性を( 1:火属性、 2:氷属性、 3:木属性、 4:土属性、 5:水属性 をそれぞれ表します)、 「.(ピリオド) 」は何もないマスをそれぞれ表します。
各属性のパチクリの数はそれぞれ 0 以上 1000 以下とします。

データセットの数は140 を超えません。また、データセットの 80 % について、W, H は100 を超えません。

Output

入力データセットごとに、最初に選んだパチクリの属性と最小移動数を１行に出力します。なお、どのように初めのパチクリを選んでも、どのような経路で移動してもそれら 4 つの属性のパチクリを捕まえることができない場合は NA と出力してください。 また、最小移動数が同じになるような最初のパチクリの選び方が複数ある場合は、属性の数字が小さいものを出力してください。

Sample Input

6 6
S.1..4
3..5..
..4.1.
4....5
.2.32.
5.1..G
3 2
...
S.G
0 0

Output for the Sample Input

2 10
NA

Solution in C (judge verdict: Wrong Answer)

// Aizu 0215: Pachimon Creature
// 2017.9.9 bal4u@uu

#include <stdio.h>
#include <string.h>

#define abs(a)  ((a)>=0?(a):-(a))
#define INF 0x7fffffff
#define MAX 50000
typedef struct { int t, r, c, k; } QUE;
QUE que[MAX]; int qsize;

#define PARENT(i) ((i)>>1)
#define LEFT(i)   ((i)<<1)
#define RIGHT(i)  (((i)<<1)+1)

void min_heapify(int i)
{
  int l, r;
  int smallest;

  l = LEFT(i), r = RIGHT(i);
  if (l < qsize && que[l].t < que[i].t) smallest = l; else smallest = i;
  if (r < qsize && que[r].t < que[smallest].t) smallest = r;
  if (smallest != i) {
    QUE t = que[i]; que[i] = que[smallest]; que[smallest] = t;
    min_heapify(smallest);
  }
}

int deq(int *r, int *c, int *t, int *k)
{
  if (qsize == 0) return 0;
  *r = que[0].r, *c = que[0].c, *t = que[0].t, *k = que[0].k;
  que[0] = que[--qsize];
  min_heapify(0);
  return 1;
}

void enq(int r, int c, int t, int k)
{
  int i, ii;

  i = qsize++;
  que[i].r = r, que[i].c = c, que[i].t = t, que[i].k = k;
  while (i > 0 && que[ii = PARENT(i)].t > que[i].t) {
    QUE tt = que[i]; que[i] = que[ii]; que[ii] = tt;
    i = ii;
  }
}

typedef struct { int r, c; } TBL;
TBL start, goal, tbl[5][1002]; int len[5];
int node[1002][1002][5]; char visited[1002][1002][5];
char map[1002][1002]; int w, h;

int search(int item)
{
	int r, c, r2, c2, i, k, fin, f, f2, t;
	TBL *tp;

	qsize = 0;
	memset(node, -1, sizeof(node)), memset(visited, 0, sizeof(visited));
	fin = (item + 4) % 5;
	node[start.r][start.c][item] = 0;
	enq(start.r, start.c, 0, item);
	while(deq(&r, &c, &t, &f)) {
		if (r == goal.r && c == goal.c) return t;
		if (visited[r][c][f]) continue;
		visited[r][c][f] = 1;
		f2 = f + 1; if (f2 >= 5) f2 = 0;
		if (f == fin && r != start.r) {
			r2 = goal.r, c2 = goal.c;
			if (visited[r2][c2][f]) continue;
			k = abs(r2 - r) + abs(c2 - c) + t;
			if (node[r2][c2][f] < 0 || k < node[r2][c2][f]) node[r2][c2][f] = k, enq(r2, c2, k, f);
		} else for (tp = tbl[f2], i = 0; i < len[f2]; i++, tp++) {
			r2 = tp->r, c2 = tp->c;
			if (visited[r2][c2][f2]) continue;
			k = abs(r2 - r) + abs(c2 - c) + t;
			if (node[r2][c2][f2] < 0 || k < node[r2][c2][f2]) node[r2][c2][f2] = k, enq(r2, c2, k, f2); 
		}
	}
	return -1;
}

int main()
{
	int r, c, s, g, i, ii, ans;
	char f[5];

	while (scanf("%d%d", &w, &h) && w) {
		for (r = 0; r < h; r++) scanf("%s", map[r]);
		memset(f, 0, sizeof(f)), memset(len, 0, sizeof(len));
		for (s = g = r = 0; r < h; r++) for (c = 0; c < w; c++) {
			char cc = map[r][c];
			if      (cc == 'S') start.r = r, start.c = c, s = 1;
			else if (cc == 'G') goal.r = r, goal.c = c, g = 1;
			else if (cc >= '1' && cc <= '5') {
				cc -= '1', f[cc] = 1;
				tbl[cc][len[cc]].r = r, tbl[cc][len[cc]].c = c, len[cc]++;
			}
		}
		for (i = 0; i < 5; i++) if (!f[i]) break;
		if (!s || !g || i < 5) { puts("NA"); continue; }
		for (ans = -1, i = 0; i < 5; i++) {
			int k = search(i);
			if (ans < 0 || k < ans) ans = k, ii = i;
		}
		printf("%d %d\n", ii+1, ans);
	}
	return 0;
}