結果

問題 No.2435 Order All Company
コンテスト
ユーザー yuyu_5510
提出日時 2023-08-15 00:30:04
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
AC  
実行時間 931 ms / 2,000 ms
コード長 1,782 bytes
コンパイル時間 137 ms
コンパイル使用メモリ 12,672 KB
実行使用メモリ 11,520 KB
最終ジャッジ日時 2024-11-23 22:52:29
合計ジャッジ時間 14,033 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 36
権限があれば一括ダウンロードができます

ソースコード

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