結果

問題 No.659 徘徊迷路
ユーザー beetbeet
提出日時 2018-03-02 22:43:05
言語 C++11
(gcc 11.4.0)
結果
WA  
実行時間 -
コード長 1,682 bytes
コンパイル時間 1,481 ms
コンパイル使用メモリ 169,140 KB
実行使用メモリ 6,948 KB
最終ジャッジ日時 2024-06-22 04:13:21
合計ジャッジ時間 2,469 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,816 KB
testcase_01 WA -
testcase_02 AC 2 ms
6,940 KB
testcase_03 AC 2 ms
6,944 KB
testcase_04 AC 19 ms
6,944 KB
testcase_05 AC 2 ms
6,944 KB
testcase_06 WA -
testcase_07 AC 1 ms
6,940 KB
testcase_08 WA -
testcase_09 WA -
testcase_10 WA -
testcase_11 AC 56 ms
6,940 KB
testcase_12 AC 13 ms
6,944 KB
testcase_13 AC 50 ms
6,944 KB
testcase_14 AC 51 ms
6,940 KB
testcase_15 AC 58 ms
6,944 KB
testcase_16 AC 59 ms
6,940 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(i,j)][idx(ny,nx)]=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