結果

問題 No.2274 三角彩色
コンテスト
ユーザー ecottea
提出日時 2023-03-02 17:18:53
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 172 ms / 2,000 ms
コード長 1,609 bytes
コンパイル時間 252 ms
コンパイル使用メモリ 82,048 KB
実行使用メモリ 80,512 KB
最終ジャッジ日時 2024-09-17 11:39:48
合計ジャッジ時間 2,847 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 22
権限があれば一括ダウンロードができます

ソースコード

diff # 

import bisect

# リスト a を座標圧縮したものを返す．
def coordinate_compression(a):
    n = len(a)
    x = sorted(list(set(a)))
    b = [0] * n
    for i in range(n):
        b[i] = bisect.bisect_left(x, a[i])
    return b

# 0-1 行列 a の階数を返す．
def bit_matrix_rank(a):
    n, m = len(a), len(a[0])

    i, j = 0, 0

    while i < n and j < m:
        i2 = i
        while i2 < n and not a[i2][j]:
            i2 += 1
        
        if i2 == n:
            j += 1
            continue

        a[i], a[i2] = a[i2], a[i]

        for i2 in range(n):
            if i2 == i:
                continue

            if a[i2][j]:
                for j2 in range(m):
                    a[i2][j2] ^= a[i][j2]
        
        i += 1
        j += 1
    
    return i


N, M, B, Q = map(int, input().split())

ijs = []
for _ in range(M):  
    i, j = map(int, input().split())
    ijs += [i, j]

# グラフを座標圧縮する．
ijs_cp = coordinate_compression(ijs)
h = max(ijs_cp) + 1

# mat : G の接続行列（行：頂点，列：辺）
mat = [[0 for j in range(M)] for i in range(h)]
for m in range(M):
    for t in range(2):
        mat[ijs_cp[2 * m + t]][m] = 1

# mat2 : クエリを並べた行列（行：クエリ，列：辺）
mat2 = [[0 for j in range(M)] for i in range(Q)]
for q in range(Q):
    m = list(map(int, input().split()))
    for t in range(3):
        mat2[q][m[t]] = 1

# それぞれの行列の階数を求める．
r = bit_matrix_rank(mat)
r2 = bit_matrix_rank(mat2)

resT = pow(2, r2, B)
resS = (pow(2, M - r, B) - resT + B) % B

print(resS, resT)
0