結果

問題 No.659 徘徊迷路
ユーザー ryoissyryoissy
提出日時 2018-03-02 22:53:43
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 36 ms / 2,000 ms
コード長 1,673 bytes
コンパイル時間 2,554 ms
コンパイル使用メモリ 146,956 KB
実行使用メモリ 7,028 KB
最終ジャッジ日時 2023-09-04 05:29:40
合計ジャッジ時間 2,601 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 5 ms
5,424 KB
testcase_01 AC 26 ms
6,868 KB
testcase_02 AC 2 ms
4,376 KB
testcase_03 AC 2 ms
4,376 KB
testcase_04 AC 16 ms
6,452 KB
testcase_05 AC 5 ms
4,956 KB
testcase_06 AC 6 ms
5,280 KB
testcase_07 AC 3 ms
4,492 KB
testcase_08 AC 35 ms
6,800 KB
testcase_09 AC 36 ms
6,812 KB
testcase_10 AC 36 ms
6,828 KB
testcase_11 AC 36 ms
6,788 KB
testcase_12 AC 35 ms
6,908 KB
testcase_13 AC 35 ms
7,028 KB
testcase_14 AC 35 ms
6,828 KB
testcase_15 AC 34 ms
6,892 KB
testcase_16 AC 35 ms
6,996 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#define MOD 1000000007LL
using namespace std;
typedef long long ll;
typedef pair<int,int> P;

int r,c,t;
int fie[11][11];
int sx,sy;
int gx,gy;
double per[35][11][11][11][11];
int dx[4]={1,0,-1,0};
int dy[4]={0,1,0,-1};
double dp[2][11][11];

int main(void){
	scanf("%d%d%d",&r,&c,&t);
	scanf("%d%d",&sy,&sx);
	scanf("%d%d",&gy,&gx);
	for(int i=0;i<r;i++){
		string str;
		cin >> str;
		for(int j=0;j<c;j++){
			if(str[j]=='#'){
				fie[i][j]=-1;
			}
		}
	}
	for(int i=0;i<r;i++){
		for(int j=0;j<c;j++){
			if(fie[i][j]==0){
				int cnt=0;
				for(int k=0;k<4;k++){
					int nx=j+dx[k];
					int ny=i+dy[k];
					if(nx>=0 && nx<c && ny>=0 && ny<r){
						if(fie[ny][nx]==0)cnt++;
					}
				}
				if(cnt==0){
					per[0][i][j][i][j]=1.0;
				}else{
					for(int k=0;k<4;k++){
						int nx=j+dx[k];
						int ny=i+dy[k];
						if(nx>=0 && nx<c && ny>=0 && ny<r){
							if(fie[ny][nx]==0)per[0][i][j][ny][nx]=(double)1.0/cnt;
						}
					}
				}
			}
		}
	}
	for(int i=0;i<30;i++){
		for(int j=0;j<r;j++){
			for(int k=0;k<c;k++){
				for(int a=0;a<r;a++){
					for(int b=0;b<c;b++){
						for(int a2=0;a2<r;a2++){
							for(int b2=0;b2<c;b2++){
								per[i+1][j][k][a2][b2]+=per[i][j][k][a][b]*per[i][a][b][a2][b2];
							}
						}
					}
				}
			}
		}
	}
	int cur=0,next=1;
	dp[cur][sy][sx]=1.0;
	for(int i=0;i<30;i++){
		if(t>>i & 1){
			for(int j=0;j<r;j++){
				for(int k=0;k<c;k++){
					for(int a=0;a<r;a++){
						for(int b=0;b<c;b++){
							dp[next][a][b]+=dp[cur][j][k]*per[i][j][k][a][b];
						}
					}
				}
			}
			swap(cur,next);
			memset(dp[next],0,sizeof(dp[next]));
		}
	}
	printf("%.10f\n",dp[cur][gy][gx]);
	return 0;
}
0