結果

問題 No.2274 三角彩色
ユーザー ecotteaecottea
提出日時 2023-03-02 17:10:21
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
AC  
実行時間 1,288 ms / 2,000 ms
コード長 1,308 bytes
コンパイル時間 725 ms
コンパイル使用メモリ 11,988 KB
実行使用メモリ 13,568 KB
最終ジャッジ日時 2023-10-17 13:37:45
合計ジャッジ時間 7,522 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 31 ms
10,268 KB
testcase_01 AC 28 ms
10,268 KB
testcase_02 AC 29 ms
10,268 KB
testcase_03 AC 28 ms
10,268 KB
testcase_04 AC 29 ms
10,268 KB
testcase_05 AC 29 ms
10,268 KB
testcase_06 AC 28 ms
10,268 KB
testcase_07 AC 28 ms
10,268 KB
testcase_08 AC 29 ms
10,268 KB
testcase_09 AC 29 ms
10,272 KB
testcase_10 AC 29 ms
10,272 KB
testcase_11 AC 31 ms
10,288 KB
testcase_12 AC 31 ms
10,296 KB
testcase_13 AC 41 ms
10,388 KB
testcase_14 AC 48 ms
10,496 KB
testcase_15 AC 568 ms
13,372 KB
testcase_16 AC 557 ms
13,372 KB
testcase_17 AC 550 ms
13,568 KB
testcase_18 AC 1,221 ms
13,488 KB
testcase_19 AC 1,288 ms
13,492 KB
testcase_20 AC 1,169 ms
13,372 KB
testcase_21 AC 53 ms
11,056 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import bisect

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

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 = [[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 = [[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