結果
問題 | No.1706 Many Bus Stops (hard) |
ユーザー |
![]() |
提出日時 | 2021-10-10 13:46:43 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 48 ms / 2,000 ms |
コード長 | 1,033 bytes |
コンパイル時間 | 304 ms |
コンパイル使用メモリ | 82,092 KB |
実行使用メモリ | 62,576 KB |
最終ジャッジ日時 | 2024-09-14 08:15:09 |
合計ジャッジ時間 | 3,857 ms |
ジャッジサーバーID (参考情報) |
judge6 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 41 |
ソースコード
import sys#input = sys.stdin.readlineinput = sys.stdin.buffer.readlinedef mat_dot(A, B, MOD=pow(10,9)+7):#A:n*m行列, B:m*l行列n, m, l = len(A), len(A[0]), len(B[0])X = [[0]*l for i in range(n)]for i in range(n):for j in range(l):temp = 0for k in range(m):temp = (temp + A[i][k] * B[k][j]) % MODX[i][j] = tempreturn Xdef mat_pow(A, x, MOD):#A^xn = len(A)X = [[0] * n for i in range(n)]for i in range(n):X[i][i] = 1for i in range(x.bit_length()):if (x>>i)&1:X = mat_dot(X, A, MOD)A = mat_dot(A, A, MOD)return Xdef main():C,N,M = map(int,input().split()); MOD = pow(10,9) + 7invc = pow(C,MOD-2,MOD)X = [[invc,0,1,0],[0,invc,0,1],[0,invc,0,0],[1-invc,1-2*invc,0,0]]v = [[1],[0],[0],[0],[0],[0]]XN = mat_pow(X,N,MOD)g = mat_dot(XN,v,MOD)#print(XN)#1台の時、それがAにいない確率p = 1 - g[0][0]ans = 1 - pow(p,M,MOD)print(ans%MOD)if __name__ == '__main__':main()