結果

問題 No.1750 ラムドスウイルスの感染拡大-hard
ユーザー shinichishinichi
提出日時 2021-11-20 01:32:46
言語 PyPy3
(7.3.15)
結果
RE  
実行時間 -
コード長 715 bytes
コンパイル時間 724 ms
コンパイル使用メモリ 87,072 KB
実行使用メモリ 79,940 KB
最終ジャッジ日時 2023-08-30 14:29:18
合計ジャッジ時間 7,451 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 RE -
testcase_01 RE -
testcase_02 RE -
testcase_03 RE -
testcase_04 RE -
testcase_05 RE -
testcase_06 RE -
testcase_07 RE -
testcase_08 RE -
testcase_09 RE -
testcase_10 RE -
testcase_11 RE -
testcase_12 RE -
testcase_13 RE -
testcase_14 RE -
testcase_15 RE -
testcase_16 RE -
testcase_17 RE -
testcase_18 RE -
testcase_19 RE -
testcase_20 RE -
testcase_21 RE -
testcase_22 RE -
testcase_23 RE -
testcase_24 RE -
testcase_25 RE -
testcase_26 RE -
testcase_27 RE -
testcase_28 RE -
testcase_29 RE -
testcase_30 RE -
testcase_31 RE -
testcase_32 RE -
testcase_33 RE -
権限があれば一括ダウンロードができます

ソースコード

diff # 

from collections import defaultdict
import numpy as np

N, M, T = map(int, input().split())
A = np.zeros(shape=(N, N), dtype=np.int64)
for _ in range(M):
    s, t = map(int, input().split())
    A[s, t] = 1
    A[t, s] = 1


mod = 998244353

def mat_mul(A, B, MOD):
    A1, A2 = A >> 15, A & (1 << 15) - 1
    B1, B2 = B >> 15, B & (1 << 15) - 1
    AB1 = np.dot(A1, B1) % MOD
    AB2 = np.dot(A2, B2)
    AB3 = np.dot(A1+A2, B1+B2) - AB1 - AB2
    return ((AB1 << 30) + (AB3 << 15) + AB2) % MOD


def mat_pow(A, K):
    P = np.eye(len(A), dtype=np.int64)
    while K:
        if K & 1:
            P = mat_mul(A, P, mod)
        A = mat_mul(A, A, mod)
        K >>= 1
    return P

P = mat_pow(A, T)
print(P[0, 0])
0