結果

問題 No.659 徘徊迷路
ユーザー fine
提出日時 2018-03-02 23:45:38
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,756 bytes
コンパイル時間 1,551 ms
コンパイル使用メモリ 172,540 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-22 14:19:51
合計ジャッジ時間 3,076 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 5
other AC * 9 WA * 3
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;

using ll = long long;

const int dh[] = {1, -1, 0, 0};
const int dw[] = {0, 0, 1, -1};

double dp[2][10][10];
int cnt[10][10];

int main() {
    cin.tie(0);
    ios::sync_with_stdio(false);
    int r, c, t;
    cin >> r >> c >> t;
    int sy, sx, gy, gx;
    cin >> sy >> sx >> gy >> gx;
    vector<string> s(r);
    for (int i = 0; i < r; i++) cin >> s[i];
    for (int h = 0; h < r; h++) {
        for (int w = 0; w < c; w++) {
            for (int j = 0; j < 4; j++) {
                int nh = h + dh[j], nw = w + dw[j];
                if (nh < 0 || nh >= r || nw < 0 || nw >= c) continue;
                cnt[h][w] += (s[nh][nw] == '.');               
            }
        }
    }
    dp[0][sy][sx] = 1;
    if (cnt[sy][sx] == 0) {
        cout << setprecision(10) << dp[0][gy][gx] << endl;
        return 0;
    }
    const int MAX_T = 100000;
    for (int i = 0; i < min(t, MAX_T); i++) {
        for (int h = 0; h < r; h++) {
            for (int w = 0; w < c; w++) {
                if (s[h][w] == '#') continue;
                for (int j = 0; j < 4; j++) {
                    int nh = h + dh[j], nw = w + dw[j];
                    if (nh < 0 || nh >= r || nw < 0 || nw >= c) continue;
                    if (s[nh][nw] == '#') continue;
                    dp[(i + 1) % 2][nh][nw] += dp[i % 2][h][w] / cnt[h][w];
                }
                if (i + 1 < min(t, MAX_T)) dp[i % 2][h][w] = 0;
            }
        }
    }
    if (t <= MAX_T) cout << setprecision(10) << dp[t % 2][gy][gx] << endl;
    else {
        if ((t - MAX_T) % 2 == 1) cout << setprecision(10) << dp[1 - t % 2][gy][gx] << endl;
        else cout << setprecision(10) << dp[t % 2][gy][gx] << endl;
    }
    return 0;
}
0