結果

問題 No.659 徘徊迷路
コンテスト
ユーザー どらら
提出日時 2018-03-02 23:24:38
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 1,282 ms / 2,000 ms
コード長 1,353 bytes
コンパイル時間 1,887 ms
コンパイル使用メモリ 178,296 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-22 11:07:35
合計ジャッジ時間 11,612 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 5
other AC * 12
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
#define REP(i,a,n) for(int i=(a); i<(int)(n); i++)
#define rep(i,n) REP(i,0,n)
#define FOR(it,c) for(__typeof((c).begin()) it=(c).begin(); it!=(c).end(); ++it)
#define ALLOF(c) (c).begin(), (c).end()
typedef long long ll;
typedef unsigned long long ull;

int vy[4] = {-1,1,0,0};
int vx[4] = {0,0,-1,1};

int main(){
  int R, C, T;
  cin >> R >> C >> T;
  int Sy, Sx, Gy, Gx;
  cin >> Sy >> Sx;
  cin >> Gy >> Gx;
  vector<string> B(R);
  rep(i,R){
    cin >> B[i];
  }

  vector<vector<double>> prev(R,vector<double>(C,0)), next(R,vector<double>(C,0));
  prev[Sy][Sx] = 1.0;
  rep(t,min(T,1000000 + T%2)){
    /*
    cout << t << endl;
    rep(i,R){
      rep(j,C){
        cout << prev[i][j] << "\t";
      }
      cout << endl;
    }
    cout << endl;
     */
    rep(i,R) rep(j,C) next[i][j] = 0;
    rep(i,R) rep(j,C){
      if(B[i][j] == '#') continue;
      int cnt = 0;
      rep(k,4){
        int ny = i + vy[k];
        int nx = j + vx[k];
        if(B[ny][nx] == '.') cnt++;
      }
      if(cnt == 0){
        next[i][j] = prev[i][j];
      }else{
        rep(k,4){
          int ny = i + vy[k];
          int nx = j + vx[k];
          if(B[ny][nx] == '.') next[ny][nx] += 1.0 / cnt * prev[i][j];
        }
      }
    }
    prev = next;
  }

  printf("%.12lf\n", prev[Gy][Gx]);
  
  return 0;
}
0