Go: ///
// File:  0558.go
// Author: ymiyamoto
//
// Created on Sun Mar  4 20:58:45 2018
//
package main

import (
	"fmt"
	"strconv"
	"strings"
)

type pos struct {
	x, y int
}

var H, W, N int
var graph [][]string

func wfs(p pos, index int) (int, pos) {
	var q []pos
	dists := make([][]int, H)
	for i := range dists {
		dists[i] = make([]int, W)
	}

	q = append(q, p)
	for len(q) > 0 {
		p := q[0]
		dist := dists[p.x][p.y]
		q = q[1:]

		for _, d := range []pos{{-1, 0}, {1, 0}, {0, -1}, {0, 1}} {
			x := p.x + d.x
			y := p.y + d.y
			if 0 <= x && x < H && 0 <= y && y < W && graph[x][y] != "X" && dists[x][y] == 0 {
				dists[x][y] = dist + 1
				q = append(q, pos{x, y})
			}
		}
	}

	for i := range graph {
		for j := range graph[i] {
			if graph[i][j] == strconv.Itoa(index) {
				return dists[i][j], pos{i, j}
			}
		}
	}

	panic("not found")
}

func main() {
	fmt.Scan(&H, &W, &N)

	var start pos

	graph = make([][]string, H)
	for i := range graph {
		var s string
		fmt.Scan(&s)
		graph[i] = strings.Split(s, "")
		for j := range graph[i] {
			if graph[i][j] == "S" {
				start.x = i
				start.y = j
			}
		}
	}

	count := 0
	for i := 1; i <= N; i++ {
		c, s := wfs(start, i)
		start = s
		count += c
	}

	fmt.Println(count)
}

C++: #include <cstdio>
#include <cstdlib>
#include <cmath>
#include <climits>
#include <cfloat>
#include <map>
#include <utility>
#include <set>
#include <iostream>
#include <memory>
#include <string>
#include <vector>
#include <algorithm>
#include <functional>
#include <sstream>
#include <complex>
#include <stack>
#include <queue>
#include <cstring>
#include <sstream>
#include <cassert>
using namespace std;
static const double EPS = 1e-10;
typedef long long ll;
typedef pair<int,int> PI;
#define rep(i,n) for(int i=0;i<(int)n;++i)
#define FOR(i,c) for(__typeof((c).begin())i=(c).begin();i!=(c).end();++i)
#define ALL(c) (c).begin(), (c).end()
#define mp(a,b) make_pair(a,b)
#define pb(a) push_back(a)
int dx[]={0,1,0,-1},dy[]={1,0,-1,0};

typedef struct{
  short x,y;
  int cost;
  short ch;
}mouse;

short maxch[1000][1000];
char ku[1000][1000];

main(){
  int h,w,n;
  cin>>h>>w>>n;
  
  rep(i,h)scanf(" %s ",ku[i]);

  int sx,sy;
  rep(i,h)rep(j,w){
    if(ku[i][j]=='S')sx=i,sy=j;
    maxch[i][j]=-1;
  }

  queue<mouse> Q;

  Q.push((mouse){sx,sy,0,0});
  maxch[sx][sy]=0;


  int ans=0;
  while(!Q.empty()){
    //    cout<<Q.size()<<endl;
    mouse t=Q.front();Q.pop();
    int cx=t.x,cy=t.y,cost=t.cost;
    int ch=t.ch;
    if(ch==n){
      ans=cost;
      break;
    }



    rep(i,4){
      int nx=cx+dx[i],ny=cy+dy[i];
      if(nx<0 || nx>=h || ny<0 || ny>=w)continue;
	 if(ku[nx][ny]=='X')continue;
	 int nch=ch;
	 if(ku[nx][ny]<='9' && ku[nx][ny]>='1'){
	   int num=ku[nx][ny]-'0';
	   if(num==t.ch+1){
	     nch=ch+1;
	   }
	 }
	 if(maxch[nx][ny]>=nch)continue;
	 maxch[nx][ny]=nch;

	 Q.push((mouse){nx,ny,cost+1,nch});
    }
  }
  cout<<ans<<endl;
}

C++: #include <iostream>
#include <queue>

using namespace std;

/*
チーズ (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
*/

int H;
int W;
int N;
char** area;

void input(){
    cin>>H>>W>>N;
    area = new char*[H];
    for(int h=0; h<H; h++){
        area[h]=new char[W];
        for(int w=0; w<W; w++){
            cin>>area[h][w];
        }
    } 
}
void output(){
    cout<<"######output######"<<endl;
    cout<<"W:"<<W<<" H:"<<H<<" N:"<<N<<endl;
    for(int h=0; h<H; h++){
        for(int w=0; w<W; w++){
            cout<<area[h][w];
        }
        cout<<endl;
    }
    cout<<"######output######"<<endl;
}

typedef pair<int, int> P;
const int INF=1000000;
int dx[4]={1,0,-1,0}, dy[4]={0,1,0,-1};
void set2p(char start, char end, P& s_p, P& e_p){
    for(int h=0; h<H; h++){
        for(int w=0; w<W; w++){
            if(area[h][w]==start){
                s_p.first=h;
                s_p.second=w;
            }
            if(area[h][w]==end){
                e_p.first=h;
                e_p.second=w;
            }
        }
    }
    return;
}
int search(char start, char end){
    int** dis;
    dis = new int*[H];
    for(int h=0;h<H;h++){
        dis[h]=new int[W];
        for(int w=0;w<W;w++) dis[h][w]=INF;
    }

    P s_p, e_p;
    set2p(start, end, s_p, e_p);
    
    queue<P> que;
    que.push(s_p);
    dis[s_p.first][s_p.second]=0;
    while(que.size()){
        P p=que.front(); que.pop();
        if(p.first==e_p.first && p.second==e_p.second) break;

        for(int i=0; i<4; i++){
            int nx=p.second+dx[i], ny=p.first+dy[i];
            if(0<=nx&&nx<W&&0<=ny&&ny<H&&area[ny][nx]!='X'&&dis[ny][nx]==INF){
                que.push(P(ny, nx));
                dis[ny][nx]=dis[p.first][p.second]+1;
            }
        }        
    }
    int distance = dis[e_p.first][e_p.second];

    for(int h=0; h<H; h++){
        delete dis[h];
    }
    delete dis;
    return distance;
}

void solve(){
    int distance=0;
    for(int step=0; step<N; step++){
        char start='0'+step;
        char end='0'+(step+1);
        if(step==0) start='S';
        distance +=search(start, end);
    }
    cout << distance <<endl;
    return;
}

void run(){
    input();
    //output();
    solve();
    return;
}

int main(){
    run();
    return 0;
}

C++: #include <iostream>
#include <stdio.h>
#include <cstring>
#include <queue>
using namespace std;

char ma[1005][1005];int vis[1005][1005];
int w,h,n,xx,yy;int setp;

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

struct node
{
    int x,y,time;
    node(int a,int b,int c)
    {
        x=a;y=b;time=c;
    }
};
node start(0,0,0);

void bfs(node st,char ed)
{
    queue<node> que;
    que.push(st);
    while(!que.empty())
    {
        node no=que.front();
        que.pop();
        for(int i=0;i<4;i++)
        {
            int x=no.x+dx[i];
            int y=no.y+dy[i];
            if(vis[x][y]||x>=w||x<0||y>=h||y<0||ma[x][y]=='X')
                continue;
            if(ma[x][y]==ed)
            {
                start.x=x;
                start.y=y;
                start.time=no.time+1;
                return ;
            }
            vis[x][y]=1;
            que.push(node(x,y,no.time+1));
        }
    }
}

int main()
{
    while(cin>>w>>h>>n,w>0)
    {
        for(int i=0;i<w;i++)
        {
            cin>>ma[i];
            for(int j=0;j<h;j++)
                if(ma[i][j]=='S')
                {
                    start.x=i;
                    start.y=j;
                    start.time=0;
                    ma[i][j]='0';
                }

        }
            for(int i=0;i<n;i++)
            {
                memset(vis,0,sizeof(vis));
                bfs(start,ma[start.x][start.y]+1);
            }

            printf("%d\n",start.time);

    }
}