Problem:
Problem G: スプリング・タイル

ある朝、あなたが目覚めると、そこはバネだらけの迷宮の中だった。

なぜ自分がこんな見知らぬ場所にいるのかはわからない。この場でじっと助けを待つという選択肢も
あるにはあるが、こういった迷宮の中に長居していると、いつかは必ず突風が吹いてきて吹き飛ばさ
れてしまうことを、あなたは経験上知っていた。しかし、突風が吹いてくるまでの時間は判断のしよ
うがない。そこであなたは、迷宮から脱出するまでに必要な移動回数の期待値が最小になるように行
動するのが最善の戦略だと考えた。

足元に落ちていた正体不明の巻物を拾ってためしに読んでみたところ、偶然にも迷宮の全体マップを
知覚することができた。さらに、たまたま持っていた草を飲むことで、すべてのワナの位置も知るこ
とができた。どうやらこの迷宮には、危険なモンスターやバネ以外のワナは存在しないらしい。

この迷宮は、いくつかの種類のタイルが長方形に並んだ形をしており、たとえば以下のようにあらわ
される。

########
#..##g.#
#*#.*#.#
#......#
#*s#*.*#
########

「.」：床。この上なら自由に歩き回ることができる。

「#」：壁。あなたは壁の中に入ることはできない。

「s」：あなたの最初の位置。このタイルの下は床になっている。

「g」：階段。このタイルの上に乗れば、迷宮を脱出したことになる。

「*」：バネ。このタイルの上に乗ると、いずれかの床の上（階段、バネのタイルは除く）にランダムに飛ばされる。それぞれの床に飛ばされる確率はすべて等しい。

タイルの種類は、この５つのいずれかである。

あなたは、現在乗っているタイルから隣り合ったタイルへ、上下左右いずれかの方向に１マスずつ移
動することができる。ただし、ナナメに移動することはできない。

迷宮のマップが与えられたとき、最善の戦略をとった場合に迷宮から脱出するまでに必要な移動回数
の期待値を求めよ。バネで飛ばされるのは移動回数には含めない。

Input

W H

c11 c12 ... c1w

c21 c22 ... c2w

::    :

ch1 ch2 ... ccw

入力の１行目には、整数W（3 ≤ W ≤ 500）と整数H（3 ≤ H ≤ 500）が、この順に空白区切りで書かれている。整数W は迷宮の横幅を、整数H は迷宮の縦幅をあらわす。

続くH 行には、迷宮のマップをあらわすW 個の文字が書かれている（空白区切りではない）。この部分の形式は、先に挙げた例のとおりである。なお、データ中に「s」と「g」は、それぞれちょうど１回ずつ出現する。また、迷宮の周囲１マスは、必ず壁で囲まれている。

なお、与えられる迷宮では、どのように移動し、バネでどのように飛ばされたとしても、「バネで飛
ばされる床を任意に設定できたとしても脱出不可能な状態」に陥ることはないと仮定してよい。

Output

最善の戦略をとったときに、脱出までにかかる移動回数の期待値を出力せよ。出力は誤差を含んでい
てもよいが、真の値との相対誤差が10-9 未満でなければならない。

Sample Input 1

8 6
########
#..##g.#
#*#.*#.#
#......#
#*s#*.*#
########

Sample Output 1

5.857142857143

Sample Input 2

21 11
#####################
#*#*.*.*.*.*.*.*.*#*#
#*#...............#*#
#*#.#.#.#.#.#.#.#.#*#
##.#.#.#.#.#.#.#.#.##
#...................#
##.#.#.#.#.#.#.#.#.##
#s#.#.#.#.#.#.#.#.#g#
#...................#
#...................#
#####################

Sample Output 2

20.000000000000

Sample Input 3

9 8
#########
#...#*.*#
#.*.#...#
#...#*.*#
#########
#.s....g#
#.*#.*..#
#########

Sample Output 3

5.000000000000

Solution in C++11 (judge verdict: Wrong Answer)

#include <bits/stdc++.h>
using namespace std;
#define repl(i,a,b) for(int i=(int)(a);i<(int)(b);i++)
#define rep(i,n) repl(i,0,n)
#define mp(a,b) make_pair((a),(b))
#define pb(a) push_back((a))
#define all(x) (x).begin(),(x).end()
#define uniq(x) sort(all(x)),(x).erase(unique(all(x)),end(x))
#define fi first
#define se second
#define dbg(x) cout<<#x" = "<<((x))<<endl
template<class T,class U> ostream& operator<<(ostream& o, const pair<T,U> &p){o<<"("<<p.fi<<","<<p.se<<")";return o;}
template<class T> ostream& operator<<(ostream& o, const vector<T> &v){o<<"[";for(T t:v){o<<t<<",";}o<<"]";return o;}

#define INF 21474836001234567

int main(){
  int w,h;
  cin>>w>>h;
  vector<string> vec(h);
  rep(i,h) cin>>vec[i];

  int sx,sy,tx,ty;
  vector<int> spx,spy;
  rep(i,h)rep(j,w){
    if(vec[i][j]=='s') sx=i,sy=j;
    if(vec[i][j]=='g') tx=i,ty=j;
    if(vec[i][j]=='*') spx.pb(i), spy.pb(j);
  }

  vector<vector<long>> d(h, vector<long>(w, INF));
  // g?????§??????????????¢
  {
    queue<int> xs,ys;
    d[tx][ty] = 0;
    xs.push(tx); ys.push(ty);
    while(!xs.empty()){
      int x = xs.front(); xs.pop();
      int y = ys.front(); ys.pop();
      const int dx[] = {0,0,-1,1}, dy[] = {-1,1,0,0};
      rep(i,4){
        int nx = x+dx[i];
        int ny = y+dy[i];
        if(vec[nx][ny]!='#' && vec[nx][ny]!='*' && d[nx][ny]>d[x][y]+1){
          d[nx][ny] = d[x][y]+1;
          xs.push(nx);
          ys.push(ny);
        }
      }
    }
  }

  vector<vector<long>> ds(h,vector<long>(w,INF));
  // ????????????spring?????§??????????????¢
  {
    queue<int> xs,ys;
    rep(i,spx.size()){
      xs.push(spx[i]);
      ys.push(spy[i]);
      ds[spx[i]][spy[i]] = 0;
    }
    while(!xs.empty()){
      int x = xs.front(); xs.pop();
      int y = ys.front(); ys.pop();
      const int dx[] = {0,0,-1,1}, dy[] = {-1,1,0,0};
      rep(i,4){
        int nx = x+dx[i];
        int ny = y+dy[i];
        if(vec[nx][ny]!='#' && vec[nx][ny]!='*' && ds[nx][ny]>ds[x][y]+1){
          ds[nx][ny] = ds[x][y]+1;
          xs.push(nx);
          ys.push(ny);
        }
      }
    }
  }

  double p = INF; // ????????§?£???°??????????????¨?????????????§????????????????
  {
    int cnt=1;
    rep(i,h)rep(j,w)if(vec[i][j]=='.') cnt++;

    double l = 0, r = 500.0*500*500*500;
    rep(_,300){
      double m = (l+r)/2.0;
      double accm = min<double>(d[sx][sy], ds[sx][sy]+m);
      rep(i,h) rep(j,w)if(vec[i][j]=='.'){
        accm += min<double>(d[i][j], ds[i][j]+m);
      }
      p = accm/cnt;
      if(p<m) r = m;
      else l = m;
    }
    p = (r+l)/2.0;
  }

  printf("%.10f\n", min<double>(d[sx][sy], ds[sx][sy]+p));

  return 0;
}