結果
問題 | No.659 徘徊迷路 |
ユーザー |
![]() |
提出日時 | 2024-07-09 04:13:21 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,100 ms / 2,000 ms |
コード長 | 1,498 bytes |
コンパイル時間 | 536 ms |
コンパイル使用メモリ | 81,904 KB |
実行使用メモリ | 81,008 KB |
最終ジャッジ日時 | 2024-07-09 04:13:36 |
合計ジャッジ時間 | 12,934 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 5 |
other | AC * 12 |
ソースコード
import sysinput = sys.stdin.readlineR,C,T=map(int,input().split())x,y=map(int,input().split())z0,w0=map(int,input().split())def prod(A,B,k,l,m):# A:k*l,B:l*mC=[[0 for i in range(m)] for j in range(k)]for i in range(k):for j in range(m):ANS=0for pl in range(l):ANS=(ANS+A[i][pl]*B[pl][j])C[i][j]=ANSreturn Cdef plus(A,B,k,l):# a,B:k*lC=[[0 for i in range(l)] for j in range(k)]for i in range(k):for j in range(l):C[i][j]=(A[i][j]+B[i][j])return CB=[input().strip() for i in range(R)]A=[[0]*(R*C) for i in range(R*C)]for i in range(R):for j in range(C):fr=C*i+jif B[i][j]==".":pos=0for z,w in [(i+1,j),(i-1,j),(i,j+1),(i,j-1)]:if 0<=z<R and 0<=w<C:if B[z][w]==".":pos+=1if pos>0:for z,w in [(i+1,j),(i-1,j),(i,j+1),(i,j-1)]:if 0<=z<R and 0<=w<C:if B[z][w]==".":to=z*C+wA[fr][to]=1/poselse:A[i*C+j][i*C+j]=1POWA=[A]for i in range(31):POWA.append(prod(POWA[-1],POWA[-1],R*C,R*C,R*C)) # ベキを求めてn=TX=[0]*(R*C)X[x*C+y]=1X=[X]while n:X=prod(X,POWA[n.bit_length()-1],1,R*C,R*C) # n乗の場合n-=1<<(n.bit_length()-1)print(X[0][z0*C+w0])