結果

問題 No.2907 Business Revealing Dora Tiles
ユーザー 👑 NachiaNachia
提出日時 2024-08-11 17:20:20
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 4,799 bytes
コンパイル時間 5,169 ms
コンパイル使用メモリ 82,188 KB
実行使用メモリ 97,312 KB
最終ジャッジ日時 2024-09-23 12:56:02
合計ジャッジ時間 42,653 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 49 ms
72,136 KB
testcase_01 AC 52 ms
64,956 KB
testcase_02 AC 51 ms
64,828 KB
testcase_03 AC 84 ms
77,068 KB
testcase_04 AC 92 ms
76,832 KB
testcase_05 AC 88 ms
77,020 KB
testcase_06 AC 144 ms
78,116 KB
testcase_07 AC 69 ms
65,676 KB
testcase_08 AC 56 ms
68,748 KB
testcase_09 AC 98 ms
76,736 KB
testcase_10 AC 69 ms
74,816 KB
testcase_11 AC 695 ms
84,676 KB
testcase_12 AC 80 ms
76,860 KB
testcase_13 AC 77 ms
76,668 KB
testcase_14 AC 99 ms
77,068 KB
testcase_15 AC 104 ms
76,888 KB
testcase_16 AC 57 ms
69,088 KB
testcase_17 AC 80 ms
76,820 KB
testcase_18 AC 289 ms
79,332 KB
testcase_19 AC 89 ms
76,960 KB
testcase_20 AC 71 ms
75,112 KB
testcase_21 AC 757 ms
82,760 KB
testcase_22 AC 184 ms
78,392 KB
testcase_23 AC 420 ms
82,756 KB
testcase_24 AC 905 ms
86,204 KB
testcase_25 AC 904 ms
85,816 KB
testcase_26 AC 949 ms
86,712 KB
testcase_27 AC 1,101 ms
86,444 KB
testcase_28 AC 895 ms
86,076 KB
testcase_29 AC 956 ms
86,200 KB
testcase_30 AC 1,038 ms
85,968 KB
testcase_31 AC 955 ms
85,940 KB
testcase_32 AC 943 ms
86,716 KB
testcase_33 AC 952 ms
85,928 KB
testcase_34 AC 911 ms
86,448 KB
testcase_35 AC 963 ms
86,028 KB
testcase_36 AC 929 ms
85,920 KB
testcase_37 AC 947 ms
86,244 KB
testcase_38 AC 913 ms
86,200 KB
testcase_39 AC 929 ms
85,944 KB
testcase_40 AC 937 ms
86,252 KB
testcase_41 AC 972 ms
86,076 KB
testcase_42 AC 953 ms
85,928 KB
testcase_43 AC 988 ms
85,804 KB
testcase_44 AC 2,691 ms
97,312 KB
testcase_45 TLE -
testcase_46 -- -
testcase_47 -- -
testcase_48 -- -
testcase_49 -- -
testcase_50 -- -
testcase_51 -- -
testcase_52 -- -
testcase_53 -- -
testcase_54 -- -
testcase_55 -- -
testcase_56 -- -
testcase_57 -- -
testcase_58 -- -
testcase_59 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

'''
translated https://yukicoder.me/submissions/998882
'''

nim_precalc1 = []
nim_inv_precalc1 = []

def nim_fill_precalc1() :
    global nim_precalc1
    nim_precalc1 = [0] * (2 ** 16)
    nim_precalc1[(1 << 8) ^ 1] = 1
    dd = 1
    while dd < 8 :
        d = 1 << dd
        c = d >> 1
        for a0 in range(d) :
            for a1 in range(d) :
                if (a0 | a1) != 0 :
                    for b0 in range(d) :
                        for b1 in range(d) :
                            if (b0 | b1) != 0 :
                                buf = 0
                                buf ^= nim_precalc1[(a1 << 8) ^ b1]
                                buf ^= nim_precalc1[(a1 << 8) ^ b0]
                                buf ^= nim_precalc1[(a0 << 8) ^ b1]
                                buf <<= dd
                                buf ^= nim_precalc1[(c << 8) ^ nim_precalc1[(a1 << 8) ^ b1]]
                                buf ^= nim_precalc1[(a0 << 8) ^ b0]
                                nim_precalc1[(((a1 << dd) ^ a0) << 8) ^ ((b1 << dd) ^ b0)] = buf
        dd *= 2

def nim_inv_precalc() :
    global nim_inv_precalc1
    nim_inv_precalc1 = [0] * 256
    for i in range(256) :
        for j in range(256) :
            if nim_precalc1[(i << 8) ^ j] == 1 :
                nim_inv_precalc1[i] = j
                break


def nim_product_full(a, b, d = 6) :
    if a == 0 or b == 0 :
        return 0
    if d == 3 :
        return nim_precalc1[(a << 8) ^ b]
    d -= 1
    lm = (1 << (1 << d)) - 1
    us = (1 << d)
    buf = 0
    a1b1 = nim_product_full(a >> us, b >> us, d)
    a2b2 = nim_product_full(a & lm, b & lm, d)
    aabb = nim_product_full((a & lm) ^ (a >> us), (b & lm) ^ (b >> us), d)
    buf ^= (aabb ^ a2b2)
    buf <<= us
    buf ^= a2b2
    buf ^= nim_product_full(1 << (us - 1), a1b1, d)
    return buf

def nim_inv_full(a, d = 6) :
    if a < 256 :
        return nim_inv_precalc1[a]
    p = 1 << (d - 1)
    a_h = a >> p
    a_l = a - (a_h << p)
    half_inv = nim_inv_full(nim_product_full(a_h ^ a_l, a_l, d) ^ nim_product_full(nim_product_full(a_h, a_h, d - 1), 1 << (p - 1)), d - 1)
    return (nim_product_full(half_inv, a_h, d) << p) ^ nim_product_full(half_inv, a_h ^ a_l, d)

def nim_product(a, b) :
    return nim_product_full(a, b)
def nim_inv(a) :
    return nim_inv_full(a)

nim_fill_precalc1()
nim_inv_precalc()

def popcount(a) :
    ans = 0
    for i in range(64) :
        ans += (a >> i) & 1
    return ans

def testcase() :
    n, t = map(int,input().split())
    a = [list(map(int,input().split())) for _ in range(t)]
    for i in range(t) :
        for j in range(n) :
            a[i][j] -= 1
    y = 0
    for x in reversed(range(0, n)) :
        if y < t :
            for tt in range(y, t) :
                if a[tt][x] != 0 :
                    a[y], a[tt] = a[tt], a[y]
                    break
            if a[y][x] == 0 :
                continue
            inv_b = nim_inv(a[y][x])
            for tt in range(n) :
                a[y][tt] = nim_product(a[y][tt], inv_b)
            for yy in range(t) :
                if y != yy :
                    times = a[yy][x]
                    for j in range(n) :
                        a[yy][j] ^= nim_product(a[y][j], times)
            y += 1
    ranks = [0] * (1 << n)
    def most_significant(x) :
        f = len(x) - 1
        while f >= 0 and x[f] == 0 :
            f -= 1
        return f
    def dfs(q, offseti, offsetd, z) :
        if z == 1 :
            ranks[offseti] = offsetd
            if q[0][0] :
                offsetd += 64
            ranks[offseti + 1] = offsetd
            return
        q1 = [[i for i in j] for j in q] # deep copy
        d1 = offsetd
        if q1[-1][-1] != 0 :
            d1 += 64
        q1.pop()
        for qq in q1 :
            qq.pop()
        dfs(q1, offseti + (1 << (z-1)), d1, z-1)
        
        q1 = q
        for qq in q1 :
            qq.pop()

        q1x = q1[-1]
        q1.pop()
        b = most_significant(q1x)

        if b >= 0 :
            inv_b = nim_inv(q1x[b])
            for t in range(b+1) :
                q1x[t] = nim_product(q1x[t], inv_b)
            for t in range(b+1, len(q1)) :
                times = q1[t][b]
                for i in range(len(q1x)) :
                    q1[t][i] ^= nim_product(q1x[i], times)
            q1[b], q1x = q1x, q1[b]
    
        dfs(q1, offseti, offsetd, z-1)
    
    qinit = [[0] * n for _ in range(n)]
    for i in range(y) :
        qinit[most_significant(a[i])] = a[i]
    dfs(qinit, 0, 0, n)
    w = [0] * (1 << n)
    for i in range(1 << n) :
        w[i] = pow(2, 64 * popcount(i) - ranks[i], 998244353)

    ans = 0
    for i in range(1 << n) :
        ans += w[i] * (-1 if (n - popcount(i)) % 2 == 1 else 1)
    print(ans % 998244353)

testcase()
0