Problem:
チーズ (Cheese)

問題

今年も JOI 町のチーズ工場がチーズの生産を始め，ねずみが巣から顔を出した．JOI 町は東西南北に区画整理されていて，各区画は巣，チーズ工場，障害物，空き地のいずれかである．ねずみは巣から出発して全てのチーズ工場を訪れチーズを 1 個ずつ食べる．

この町には，N 個のチーズ工場があり，どの工場も１種類のチーズだけを生産している．チーズの硬さは工場によって異なっており，硬さ 1 から N までのチーズを生産するチーズ工場がちょうど 1 つずつある．

ねずみの最初の体力は 1 であり，チーズを 1 個食べるごとに体力が 1 増える．ただし，ねずみは自分の体力よりも硬いチーズを食べることはできない．

ねずみは，東西南北に隣り合う区画に 1 分で移動することができるが，障害物の区画には入ることができない．チーズ工場をチーズを食べずに通り過ぎることもできる．すべてのチーズを食べ終えるまでにかかる最短時間を求めるプログラムを書け．ただし，ねずみがチーズを食べるのにかかる時間は無視できる．

入力

入力は H+1 行ある．1 行目には 3 つの整数 H，W，N (1 ≤ H ≤ 1000，1 ≤ W ≤ 1000，1 ≤ N ≤ 9) がこの順に空白で区切られて書かれている．2 行目から H+1 行目までの各行には，'S'，'1', '2', ..., '9'，'X'，'.' からなる W 文字の文字列が書かれており，各々が各区画の状態を表している．北から i 番目，西から j 番目の区画を (i,j) と記述することにすると (1 ≤ i ≤ H, 1 ≤ j ≤ W)，第 i+1 行目の j 番目の文字は，区画 (i,j) が巣である場合は 'S' となり，障害物である場合は 'X' となり，空き地である場合は '.' となり，硬さ 1, 2, ..., 9 のチーズを生産する工場である場合はそれぞれ '1', '2', ..., '9' となる．入力には巣と硬さ 1, 2, ..., N のチーズを生産する工場がそれぞれ 1 つずつある．他のマスは障害物または空き地であることが保証されている．ねずみは全てのチーズを食べられることが保証されている．

出力

すべてのチーズを食べ終えるまでにかかる最短時間（分）を表す整数を 1 行で出力せよ．

入出力例

入力例 1

3 3 1
S..
...
..1

出力例 1

4

入力例 2

4 5 2
.X..1
....X
.XX.S
.2.X.

出力例 2

12

入力例 3

10 10 9
.X...X.S.X
6..5X..X1X
...XXXX..X
X..9X...X.
8.X2X..X3X
...XX.X4..
XX....7X..
X..X..XX..
X...X.XX..
..X.......

出力例 3

91

問題文と自動審判に使われるデータは、情報オリンピック日本委員会が作成し公開している問題文と採点用テストデータです。

Solution in C (judge verdict: Accepted)

#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#define recp(I,Y,X) {q[I].x=X;q[I].y=Y;m[Y][X]='X';I++;}
char m[1005][1005],m2[1005][1005];
typedef struct point{
	int x,y;
}POINT;
POINT sq[500],q[500],g[11];
int h,w,n,i,j,ans=0,gx,gy,cx,cy,sqn,qn;
int main(){
	char s[1005];
	scanf("%d%d%d",&h,&w,&n);
	for(i=0;i<=w+1;i++)m[0][i]=m[h+1][i]='X';
	for(i=1;i<=h;i++){
		scanf("%s",s);
		strcat(m[i]+1,s);
		for(j=1;j<=w;j++){
			if(m[i][j]=='S')g[0].y=i,g[0].x=j;
			else if(m[i][j]!='.' && m[i][j]!='X')g[m[i][j]-'0'].y=i,g[m[i][j]-'0'].x=j;
		}
		m[i][0]=m[i][w+1]='X';
	}
	memcpy(m2,m,sizeof(m));
	for(i=1;i<=n;i++){
		sq[0].x=g[i-1].x,sq[0].y=g[i-1].y;
		m[g[i-1].y][g[i-1].x]='X',gx=g[i].x,gy=g[i].y;sqn=1;
		//printf("%d %d\n",sq[0].y,sq[0].x);
		//printf("%d %d\n",gy,gx);
		while(sqn){
			qn=0;
			//printf("\n");
			for(j=0;j<sqn;j++){
				cx=sq[j].x;cy=sq[j].y;
				//printf("%d %d\n",cy,cx);
				if(cx==gx && cy==gy)goto nex;
				if(m[cy-1][cx]!='X')recp(qn,cy-1,cx);
				if(m[cy+1][cx]!='X')recp(qn,cy+1,cx);
				if(m[cy][cx-1]!='X')recp(qn,cy,cx-1);
				if(m[cy][cx+1]!='X')recp(qn,cy,cx+1);
			}
			memcpy(sq,q,sizeof(POINT)*qn);
			//printf("%d\n",qn);
			sqn=qn;
			ans++;
		}
		nex:
		memcpy(m,m2,sizeof(m2));
	}
	printf("%d\n",ans);
	return 0;
}