結果
問題 | No.1547 [Cherry 2nd Tune *] 偶然の勝利の確率 |
ユーザー |
|
提出日時 | 2022-03-15 21:57:04 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 550 ms / 2,000 ms |
コード長 | 1,577 bytes |
コンパイル時間 | 151 ms |
コンパイル使用メモリ | 82,452 KB |
実行使用メモリ | 77,148 KB |
最終ジャッジ日時 | 2024-09-22 14:19:10 |
合計ジャッジ時間 | 9,427 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 36 |
ソースコード
P = 998244353m,n,S = map(int,input().split())pa = m * pow(n,P-2,P) % Pqa = (1 - pa) % Pm,n,T = map(int,input().split())pb = m * pow(n,P-2,P) % Pqb = (1 - pb) % PK = int(input())C = 205SA = [0] * CSB = [0] * Cnow = qaSA[0] = qafor i in range(1,C):now = now * pa % PSA[i] = nownow = qbSB[0] = nowfor i in range(1,C):now = now * pb % PSB[i] = nown = S + T + 1def seki(x,y):l = [[0] * n for _ in range(n)]for i in range(n):for j in range(n):for k in range(n):l[i][j] = (l[i][j] + x[i][k] * y[k][j]) % Preturn le = [[0] * n for _ in range(n)]for i in range(n):e[i][i] = 1R = [[0] * n for _ in range(n)]for j in range(-T,S+1):for i in range(-T,S+1):if i == S:if j != S:continueR[S-j][S-i] = 1continueif i == -T:if j != -T:continueR[S-j][S-i] = 1continueif j == S:tmp = pow(pa,S - i,P)R[S-j][S-i] = tmpcontinueif j == -T:tmp = 0for k in range(i,S):tmp += SA[k-i] * pow(pb,k-j,P) % Ptmp %= PR[S-j][S-i] = tmpcontinueelse:delta = max(i,j)tmp = 0for k in range(delta,S):tmp += SA[k-i] * SB[k - j] % Ptmp %= PR[S-j][S-i] = tmpwhile K:if K & 1:e = seki(e,R)R = seki(R,R)K >>= 1a = e[0][S]b = e[-1][S]print(a)print(b)