結果

問題 No.659 徘徊迷路
ユーザー sekiya9311sekiya9311
提出日時 2018-03-04 20:54:15
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 40 ms / 2,000 ms
コード長 2,130 bytes
コンパイル時間 1,746 ms
コンパイル使用メモリ 175,532 KB
実行使用メモリ 4,380 KB
最終ジャッジ日時 2023-09-22 23:31:31
合計ジャッジ時間 2,883 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 7 ms
4,380 KB
testcase_01 AC 13 ms
4,376 KB
testcase_02 AC 34 ms
4,380 KB
testcase_03 AC 35 ms
4,376 KB
testcase_04 AC 34 ms
4,376 KB
testcase_05 AC 8 ms
4,376 KB
testcase_06 AC 10 ms
4,376 KB
testcase_07 AC 10 ms
4,376 KB
testcase_08 AC 15 ms
4,376 KB
testcase_09 AC 34 ms
4,380 KB
testcase_10 AC 39 ms
4,376 KB
testcase_11 AC 39 ms
4,376 KB
testcase_12 AC 9 ms
4,376 KB
testcase_13 AC 35 ms
4,376 KB
testcase_14 AC 35 ms
4,376 KB
testcase_15 AC 40 ms
4,376 KB
testcase_16 AC 40 ms
4,380 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;

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

using row = std::vector<double>;
using matrix = std::vector<row>;
//(A * B) mod
matrix matrixMul(matrix &A, matrix &B, const int mod) {
    const int asz = A.size();
    const int bsz = B.size();
    const int bcolsz = B.front().size();
    matrix res(asz, row(bcolsz, 0));
    for (int i = 0; i < asz; i++) {
        for (int k = 0; k < bsz; k++) {
            for (int j = 0; j < bcolsz; j++) {
                (res[i][j] += A[i][k] * B[k][j]);
            }
        }
    }
    return res;
}
//(mat^p) mod
matrix matPowMod(matrix mat, long long p, const int mod) {
    const int matsz = mat.size();
    assert((int)mat.front().size() == matsz);//must square matrix!!
    matrix res(matsz, row(matsz, 0));
    for (int i = 0; i < matsz; i++) {
        res[i][i] = 1;
    }
    while (p) {
        if (p & 1) res = matrixMul(res, mat, mod);
        mat = matrixMul(mat, mat, mod);
        p >>= 1;
    }
    return res;
}

matrix make_matrix(int r, int c) {
    return matrix(r, row(c));
}

int main() {
    int r, c, t;
    cin >> r >> c >> t;
    int sy, sx, gy, gx;
    cin >> sy >> sx >> gy >> gx;
    auto mat = make_matrix(100, 100);
    vector<string> B(r);
    for_each(begin(B), end(B), [](auto &e){cin >> e;});
    for (int i = 0; i < r; i++) {
        for (int j = 0; j < c; j++) {
            const int rIdx = i * 10 + j;
            if (B[i][j] == '.') {
                int cnt = 0;
                for (int k = 0; k < 4; k++) {
                    cnt += B[i + dy[k]][j + dx[k]] == '.';
                }
                if (cnt) {
                    for (int k = 0; k < 4; k++) {
                        const int cIdx = (i + dy[k]) * 10 + j + dx[k];
                        mat[rIdx][cIdx] = 1. * (B[i + dy[k]][j + dx[k]] == '.') / cnt;
                    }
                } else {
                    mat[rIdx][rIdx] = 1;
                }
            }
        }
    }
    auto mat2 = matPowMod(mat, t, -1);
    cout << fixed << setprecision(10) << mat2[sy * 10 + sx][gy * 10 + gx] << endl;;
}
0