結果

問題 No.659 徘徊迷路
ユーザー beetbeet
提出日時 2018-03-02 22:48:41
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 66 ms / 2,000 ms
コード長 1,680 bytes
コンパイル時間 1,519 ms
コンパイル使用メモリ 155,460 KB
実行使用メモリ 4,380 KB
最終ジャッジ日時 2023-09-04 05:09:26
合計ジャッジ時間 2,930 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,380 KB
testcase_01 AC 14 ms
4,380 KB
testcase_02 AC 2 ms
4,376 KB
testcase_03 AC 2 ms
4,376 KB
testcase_04 AC 21 ms
4,380 KB
testcase_05 AC 3 ms
4,380 KB
testcase_06 AC 2 ms
4,380 KB
testcase_07 AC 2 ms
4,380 KB
testcase_08 AC 24 ms
4,380 KB
testcase_09 AC 56 ms
4,376 KB
testcase_10 AC 64 ms
4,376 KB
testcase_11 AC 65 ms
4,376 KB
testcase_12 AC 13 ms
4,376 KB
testcase_13 AC 56 ms
4,380 KB
testcase_14 AC 57 ms
4,380 KB
testcase_15 AC 66 ms
4,376 KB
testcase_16 AC 66 ms
4,376 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;
using Int = long long;

typedef vector<double> arr;
typedef vector<arr> mat;
inline arr mul(const mat &a,arr &b){
  arr res(b.size(),0);
  for(Int i=0;i<(Int)b.size();i++)
    for(Int j=0;j<(Int)a[i].size();j++)
      res[i]+=a[i][j]*b[j];
  return res;
}
inline mat mul(const mat &a,const mat &b){
  mat res(a.size(),arr(b[0].size(),0));
  for(Int i=0;i<(Int)a.size();i++)
    for(Int j=0;j<(Int)b[0].size();j++)
      for(Int k=0;k<(Int)b.size();k++)
	res[i][j]+=a[i][k]*b[k][j];
  return res;
}
inline mat mat_pow(mat a,Int n){
  mat res(a);
  for(Int i=0;i<(Int)a.size();i++)
    for(Int j=0;j<(Int)a[i].size();j++)
      res[i][j]=(i==j);
  while(n){
    if(n&1) res=mul(a,res);
    a=mul(a,a);
    n>>=1;
  }
  return res;
}


struct Precision{
  Precision(){
    cout<<fixed<<setprecision(12);
  }
}precision_beet;

//INSERT ABOVE HERE
signed main(){
  Int h,w,t;
  cin>>h>>w>>t;
  Int sy,sx,gy,gx;
  cin>>sy>>sx>>gy>>gx;
  vector<string> s(h);
  for(Int i=0;i<h;i++) cin>>s[i];

  Int N=h*w;
  auto idx=[&](Int y,Int x){return y*w+x;};
  mat m(N,arr(N,0));
  Int dy[]={1,-1,0,0};
  Int dx[]={0,0,1,-1};
  
  for(Int i=0;i<h;i++){
    for(Int j=0;j<w;j++){
      if(s[i][j]=='#') continue;
      Int cnt=0;
      for(Int k=0;k<4;k++){
	Int ny=i+dy[k],nx=j+dx[k];
	if(s[ny][nx]=='#') continue;
	cnt++;
      }
      for(Int k=0;k<4;k++){
	Int ny=i+dy[k],nx=j+dx[k];
	if(s[ny][nx]=='#') continue;
	m[idx(ny,nx)][idx(i,j)]=1.0/cnt;
      }
      if(!cnt) m[idx(i,j)][idx(i,j)]=1.0;
    }
  }
  
  {
    auto a=mat_pow(m,t);
    arr b(N,0);
    b[idx(sy,sx)]=1;
    auto c=mul(a,b);
    cout<<c[idx(gy,gx)]<<endl;
  }
  return 0;
}
0