結果
問題 | No.803 Very Limited Xor Subset |
ユーザー | tpyneriver |
提出日時 | 2019-05-30 13:34:05 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 95 ms / 2,000 ms |
コード長 | 3,678 bytes |
コンパイル時間 | 170 ms |
コンパイル使用メモリ | 82,284 KB |
実行使用メモリ | 76,684 KB |
最終ジャッジ日時 | 2024-09-17 16:06:24 |
合計ジャッジ時間 | 4,921 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 40 ms
52,816 KB |
testcase_01 | AC | 39 ms
53,636 KB |
testcase_02 | AC | 39 ms
53,820 KB |
testcase_03 | AC | 40 ms
53,584 KB |
testcase_04 | AC | 40 ms
54,316 KB |
testcase_05 | AC | 39 ms
53,248 KB |
testcase_06 | AC | 40 ms
53,972 KB |
testcase_07 | AC | 40 ms
53,596 KB |
testcase_08 | AC | 40 ms
54,044 KB |
testcase_09 | AC | 40 ms
54,108 KB |
testcase_10 | AC | 40 ms
54,032 KB |
testcase_11 | AC | 42 ms
54,340 KB |
testcase_12 | AC | 41 ms
54,536 KB |
testcase_13 | AC | 40 ms
53,664 KB |
testcase_14 | AC | 88 ms
76,316 KB |
testcase_15 | AC | 87 ms
76,352 KB |
testcase_16 | AC | 90 ms
76,560 KB |
testcase_17 | AC | 88 ms
76,344 KB |
testcase_18 | AC | 91 ms
76,584 KB |
testcase_19 | AC | 90 ms
76,348 KB |
testcase_20 | AC | 89 ms
76,228 KB |
testcase_21 | AC | 88 ms
76,192 KB |
testcase_22 | AC | 92 ms
76,192 KB |
testcase_23 | AC | 92 ms
76,300 KB |
testcase_24 | AC | 91 ms
76,312 KB |
testcase_25 | AC | 93 ms
76,684 KB |
testcase_26 | AC | 92 ms
76,272 KB |
testcase_27 | AC | 90 ms
76,272 KB |
testcase_28 | AC | 91 ms
76,608 KB |
testcase_29 | AC | 92 ms
76,272 KB |
testcase_30 | AC | 92 ms
76,560 KB |
testcase_31 | AC | 93 ms
76,212 KB |
testcase_32 | AC | 92 ms
76,204 KB |
testcase_33 | AC | 94 ms
76,440 KB |
testcase_34 | AC | 59 ms
67,828 KB |
testcase_35 | AC | 91 ms
76,248 KB |
testcase_36 | AC | 80 ms
75,920 KB |
testcase_37 | AC | 82 ms
75,688 KB |
testcase_38 | AC | 46 ms
60,984 KB |
testcase_39 | AC | 47 ms
61,776 KB |
testcase_40 | AC | 93 ms
76,260 KB |
testcase_41 | AC | 95 ms
76,416 KB |
testcase_42 | AC | 93 ms
76,420 KB |
testcase_43 | AC | 90 ms
76,196 KB |
testcase_44 | AC | 40 ms
53,364 KB |
testcase_45 | AC | 39 ms
53,864 KB |
testcase_46 | AC | 52 ms
65,244 KB |
ソースコード
class GE(): #mod2 なら Aは1重リストで渡して, max(bit_length) をMに渡してあげた方が速い def __init__(self, A, mod, M = 0): self.N = len(A) if M: self.M = M else: self.M = len(A[0]) self.mod = mod self.A = A[:] self.uptri = None self.A2 = None self.pivot = [] self.R = None def ut(self): if self.mod != 2: if self.uptri is not None: return self.uptri self.uptri = self.A[:] c = 0 for i in range(self.N): if i + c >= self.M: break while self.uptri[i][i+c] == 0: for j in range(i+1, self.N): if self.uptri[j][i+c]: self.uptri[i], self.uptri[j] = self.uptri[j], self.uptri[i] break else: c += 1 if i + c == self.M: break else: self.pivot.append((i, i+c)) t = pow(self.uptri[i][i+c], self.mod - 2, self.mod) for j in range(i+1, self.N): tj = t * self.uptri[j][i+c] self.uptri[j][i+c:] = [(aj - tj*ai) % self.mod for ai, aj in zip(self.uptri[i][i+c:], self.uptri[j][i+c:])] continue break for pi, pj in self.pivot[::-1]: t = pow(self.uptri[pi][pj], self.mod-2, self.mod) self.uptri[pi][pj:] = [(ai * t) % self.mod for ai in self.uptri[pi][pj:]] for i in range(pi-1, -1, -1): ti = self.uptri[i][pj] self.uptri[i][pj:] = [(ai - api*ti) % self.mod for ai, api in zip(self.uptri[i][pj:], self.uptri[pi][pj:])] return self.uptri else: if self.A2 is not None: return self.A2 self.A2 = self.A[:] c = 0 for i in range(self.N): if i + c >= self.M: break while not self.A2[i] & 2**(i+c): for j in range(i+1, self.N): if self.A2[j] & 2**(i+c): self.A2[i], self.A2[j] = self.A2[j], self.A2[i] break else: c += 1 if i + c == self.M: break else: self.pivot.append((i, i+c)) for j in range(i+1, self.N): if self.A2[j] & 2**(i+c): self.A2[j] ^= self.A2[i] continue break for pi, pj in self.pivot[::-1]: for i in range(pi-1, -1, -1): if 2**pj & self.A2[i]: self.A2[i] ^= self.A2[pi] return self.A2 def rank(self): if self.R is not None: return self.R self.ut() self.R = len(self.pivot) return self.R import sys N, M, X = map(int, input().split()) mod = 10**9+7 A = list(map(int,input().split())) + [X] AM = [sum(2**j if 2**i & a else 0 for j, a in enumerate(A)) for i in range(30)] for _ in range(M): t, l, r = map(int, sys.stdin.readline().split()) AM.append(2**r - 2**(l-1) + 2**N*t) T = GE(AM, 2, N+1) ans = 0 for u in T.ut(): if u == 2**N: break else: ans = pow(2, N-T.rank(), mod) print(ans)