結果

問題 No.659 徘徊迷路
コンテスト
ユーザー beet
提出日時 2018-03-02 22:48:41
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 58 ms / 2,000 ms
コード長 1,680 bytes
コンパイル時間 1,297 ms
コンパイル使用メモリ 169,196 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-22 04:37:09
合計ジャッジ時間 2,346 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 5
other AC * 12
権限があれば一括ダウンロードができます

ソースコード

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