結果
問題 | No.1750 ラムドスウイルスの感染拡大-hard |
ユーザー |
![]() |
提出日時 | 2021-11-21 12:46:44 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,837 ms / 2,000 ms |
コード長 | 1,010 bytes |
コンパイル時間 | 358 ms |
コンパイル使用メモリ | 82,512 KB |
実行使用メモリ | 76,988 KB |
最終ジャッジ日時 | 2024-06-22 18:54:36 |
合計ジャッジ時間 | 21,943 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 30 |
ソースコード
import sys #input = sys.stdin.readline #input = sys.stdin.buffer.readline def 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 = 0 for k in range(m): temp = (temp + A[i][k] * B[k][j]) % MOD X[i][j] = temp return X def mat_pow(A, x, MOD): #A^x n = len(A) X = [[0] * n for i in range(n)] for i in range(n): X[i][i] = 1 for i in range(x.bit_length()): if (x>>i)&1: X = mat_dot(X, A, MOD) A = mat_dot(A, A, MOD) return X def main(): N,M,T = map(int,input().split()); MOD = 998244353 X = [[0]*N for _ in range(N)] for _ in range(M): a,b = map(int,input().split()) X[a][b] = 1 X[b][a] = 1 v = [[0] for _ in range(N)] v[0][0] = 1 XT = mat_pow(X,T,MOD) ans = mat_dot(XT,v,MOD) #print(ans) print(ans[0][0]) if __name__ == '__main__': main()