結果

問題 No.2435 Order All Company
ユーザー yuyu_5510yuyu_5510
提出日時 2023-08-15 00:30:04
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
AC  
実行時間 802 ms / 2,000 ms
コード長 1,782 bytes
コンパイル時間 98 ms
コンパイル使用メモリ 10,780 KB
実行使用メモリ 8,864 KB
最終ジャッジ日時 2023-08-15 15:11:15
合計ジャッジ時間 12,412 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 17 ms
8,120 KB
testcase_01 AC 16 ms
8,160 KB
testcase_02 AC 16 ms
8,208 KB
testcase_03 AC 16 ms
8,036 KB
testcase_04 AC 21 ms
8,148 KB
testcase_05 AC 802 ms
8,724 KB
testcase_06 AC 787 ms
8,860 KB
testcase_07 AC 585 ms
8,756 KB
testcase_08 AC 791 ms
8,732 KB
testcase_09 AC 796 ms
8,748 KB
testcase_10 AC 45 ms
8,136 KB
testcase_11 AC 24 ms
8,096 KB
testcase_12 AC 23 ms
8,020 KB
testcase_13 AC 44 ms
8,560 KB
testcase_14 AC 393 ms
8,044 KB
testcase_15 AC 484 ms
8,632 KB
testcase_16 AC 339 ms
8,128 KB
testcase_17 AC 760 ms
8,860 KB
testcase_18 AC 309 ms
8,656 KB
testcase_19 AC 300 ms
8,580 KB
testcase_20 AC 330 ms
8,080 KB
testcase_21 AC 351 ms
8,048 KB
testcase_22 AC 393 ms
8,028 KB
testcase_23 AC 233 ms
8,644 KB
testcase_24 AC 383 ms
8,644 KB
testcase_25 AC 368 ms
8,196 KB
testcase_26 AC 369 ms
8,608 KB
testcase_27 AC 382 ms
8,076 KB
testcase_28 AC 350 ms
8,052 KB
testcase_29 AC 390 ms
8,464 KB
testcase_30 AC 433 ms
8,864 KB
testcase_31 AC 16 ms
8,052 KB
testcase_32 AC 16 ms
8,228 KB
testcase_33 AC 16 ms
8,068 KB
testcase_34 AC 46 ms
8,044 KB
testcase_35 AC 127 ms
8,660 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

mod = 998244353

def matrix_tree_theory(mat, N):
    ret = 1
    for i in range(N-1):
        if mat[i][i] == 0:
            for j in range(i+1, N-1):
                if mat[j][i]:
                    for k in range(N-1):
                        mat[i][k], mat[j][k] = mat[j][k], mat[i][k]
                        ret = (-ret + mod) % mod
                    break

        if mat[i][i] == 0:
            return 0

        for j in range(i+1, N-1):
            if mat[j][i]:
                mul = (mat[j][i] * pow(mat[i][i], -1, mod)) % mod
                
                for k in range(i, N-1):
                    mat[j][k] = (mat[j][k] -  (mul * mat[i][k]) % mod + mod) % mod
    
    for i in range(N-1):
        ret = (ret * mat[i][i]) % mod

    return ret

N, K = map(int, input().split())
G = [[[0 for _ in range(N)] for _ in range(N)] for _ in range(K)]
for i in range(K):
    t = int(input())
    for j in range(t):
        a, b = map(int, input().split())
        a -= 1
        b -= 1
        G[i][a][b] = (G[i][a][b] -1 + mod) % mod
        G[i][b][a] = (G[i][b][a] -1 + mod) % mod

for i in range(K):
    for j in range(N):
        sum = 0
        for k in range(N):
            sum += G[i][j][k]
            sum %= mod
        sum = (-sum + mod) % mod
        G[i][j][j] = sum
ans = 0
for i in range(1<<K):
    mat = [[0 for _ in range(N)] for _ in range(N)]
    cnt = 0
    for j in range(K):
        if i & 1<<j:
            cnt += 1
            for a in range(N):
                for b in range(N):
                    mat[a][b] += G[j][a][b]
                    mat[a][b] %= mod

    if (K-cnt) % 2 == 0:
        ans += matrix_tree_theory(mat, N)
        ans %= mod
    else:
        ans -= matrix_tree_theory(mat, N)
        ans = (ans + mod) % mod

print(ans)
0