結果
問題 | No.1750 ラムドスウイルスの感染拡大-hard |
ユーザー |
![]() |
提出日時 | 2025-03-31 17:40:51 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 907 ms / 2,000 ms |
コード長 | 1,401 bytes |
コンパイル時間 | 225 ms |
コンパイル使用メモリ | 82,432 KB |
実行使用メモリ | 77,064 KB |
最終ジャッジ日時 | 2025-03-31 17:42:21 |
合計ジャッジ時間 | 10,818 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 30 |
ソースコード
MOD = 998244353def main():import sysinput = sys.stdin.read().split()idx = 0N = int(input[idx]); idx +=1M = int(input[idx]); idx +=1T = int(input[idx]); idx +=1# Initialize adjacency matrixA = [[0]*N for _ in range(N)]for _ in range(M):s = int(input[idx]); idx +=1t = int(input[idx]); idx +=1A[s][t] += 1A[t][s] += 1# Matrix multiplication functiondef multiply(a, b):res = [[0]*N for _ in range(N)]for i in range(N):for k in range(N):if a[i][k] == 0:continuefor j in range(N):res[i][j] = (res[i][j] + a[i][k] * b[k][j]) % MODreturn res# Matrix exponentiation (binary method)def matrix_pow(mat, power):# Initialize result as identity matrixresult = [[0]*N for _ in range(N)]for i in range(N):result[i][i] = 1current = [row[:] for row in mat]while power > 0:if power % 2 == 1:result = multiply(result, current)current = multiply(current, current)power //= 2return result# Compute A^Tif T == 0:print(1 % MOD)returnA_power = matrix_pow(A, T)print(A_power[0][0] % MOD)if __name__ == "__main__":main()