結果

問題 No.659 徘徊迷路
コンテスト
ユーザー chocorusk
提出日時 2018-12-17 00:40:57
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 149 ms / 2,000 ms
コード長 2,039 bytes
コンパイル時間 985 ms
コンパイル使用メモリ 106,236 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-09-25 06:41:50
合計ジャッジ時間 3,002 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 5
other AC * 12
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cstdio>
#include <cstring>
#include <iostream>
#include <string>
#include <cmath>
#include <bitset>
#include <vector>
#include <map>
#include <set>
#include <queue>
#include <deque>
#include <algorithm>
#include <complex>
#include <unordered_map>
#include <unordered_set>
#include <random>
#include <cassert>
using namespace std;
typedef long long int ll;
typedef pair<int, int> P;
vector<vector<double>> matrixmul(int l, int m, int n, vector<vector<double>> a, vector<vector<double>> b){
	vector<vector<double>> c;
	for(int i=0; i<l; i++){
		vector<double> v;
		for(int j=0; j<n; j++){
			double x=0;
			for(int k=0; k<m; k++){
				x+=(a[i][k]*b[k][j]);
			}
			v.push_back(x);
		}
		c.push_back(v);
	}
	return c;
}
vector<vector<double>> matrixpow(int n, vector<vector<double>> a, ll k){
	vector<vector<double>> ap=a, ans;
	for(int i=0; i<n; i++){
		vector<double> v;
		for(int j=0; j<n; j++){
			if(i==j){
				v.push_back(1);
			}else{
				v.push_back(0);
			}
		}
		ans.push_back(v);
	}
	while(k){
		if(k&1) ans=matrixmul(n, n, n, ap, ans);
		ap=matrixmul(n, n, n, ap, ap);
		k>>=1;
	}
	return ans;
}
int main()
{
	int r, c; ll t;
	cin>>r>>c>>t;
	int sy, sx, gy, gx;
	cin>>sy>>sx>>gy>>gx;
	string b[10];
	for(int i=0; i<r; i++) cin>>b[i];
	int dx[4]={1, 0, -1, 0}, dy[4]={0, 1, 0, -1};
	vector<vector<double>> a;
	for(int i=0; i<r; i++){
		for(int j=0; j<c; j++){
			vector<double> v;
			if(b[i][j]=='#'){
				for(int k=0; k<r*c; k++){
					v.push_back(0);
				}
				a.push_back(v);
				continue;
			}
			double ct=0;
			for(int k=0; k<4; k++){
				if(b[i+dx[k]][j+dy[k]]!='#') ct+=1;
			}
			if(ct==0){
				for(int k=0; k<r*c; k++){
					if(k==i*c+j) v.push_back(1);
					else v.push_back(0);
				}
				a.push_back(v);
				continue;
			}
			for(int k=0; k<r*c; k++){
				if(((abs(k/c-i)==1 && k%c==j) || (abs(k%c-j)==1 && k/c==i)) && b[k/c][k%c]!='#'){
					v.push_back(1/ct);
				}else{
					v.push_back(0);
				}
			}
			a.push_back(v);
		}
	}
	a=matrixpow(r*c, a, t);
	printf("%.7lf\n", a[sy*c+sx][gy*c+gx]);
	return 0;
}
0