結果
問題 | No.2505 matriX cOnstRuction |
ユーザー | rniya |
提出日時 | 2023-09-10 17:48:40 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 304 ms / 2,500 ms |
コード長 | 3,854 bytes |
コンパイル時間 | 3,916 ms |
コンパイル使用メモリ | 216,720 KB |
実行使用メモリ | 79,912 KB |
最終ジャッジ日時 | 2023-10-13 18:19:15 |
合計ジャッジ時間 | 10,187 ms |
ジャッジサーバーID (参考情報) |
judge12 / judge14 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
4,356 KB |
testcase_01 | AC | 25 ms
4,348 KB |
testcase_02 | AC | 20 ms
4,348 KB |
testcase_03 | AC | 20 ms
4,352 KB |
testcase_04 | AC | 21 ms
4,476 KB |
testcase_05 | AC | 21 ms
4,348 KB |
testcase_06 | AC | 27 ms
4,348 KB |
testcase_07 | AC | 13 ms
4,352 KB |
testcase_08 | AC | 13 ms
4,348 KB |
testcase_09 | AC | 13 ms
4,352 KB |
testcase_10 | AC | 13 ms
4,352 KB |
testcase_11 | AC | 13 ms
4,348 KB |
testcase_12 | AC | 13 ms
4,352 KB |
testcase_13 | AC | 13 ms
4,352 KB |
testcase_14 | AC | 13 ms
4,352 KB |
testcase_15 | AC | 13 ms
4,352 KB |
testcase_16 | AC | 13 ms
4,348 KB |
testcase_17 | AC | 13 ms
4,352 KB |
testcase_18 | AC | 13 ms
4,352 KB |
testcase_19 | AC | 13 ms
4,352 KB |
testcase_20 | AC | 14 ms
4,348 KB |
testcase_21 | AC | 14 ms
4,352 KB |
testcase_22 | AC | 13 ms
4,352 KB |
testcase_23 | AC | 13 ms
4,352 KB |
testcase_24 | AC | 13 ms
4,352 KB |
testcase_25 | AC | 13 ms
4,352 KB |
testcase_26 | AC | 12 ms
4,352 KB |
testcase_27 | AC | 13 ms
4,352 KB |
testcase_28 | AC | 13 ms
4,348 KB |
testcase_29 | AC | 13 ms
4,348 KB |
testcase_30 | AC | 10 ms
4,348 KB |
testcase_31 | AC | 9 ms
4,348 KB |
testcase_32 | AC | 9 ms
4,348 KB |
testcase_33 | AC | 9 ms
4,352 KB |
testcase_34 | AC | 15 ms
5,332 KB |
testcase_35 | AC | 9 ms
4,428 KB |
testcase_36 | AC | 15 ms
6,776 KB |
testcase_37 | AC | 13 ms
5,464 KB |
testcase_38 | AC | 19 ms
8,300 KB |
testcase_39 | AC | 41 ms
20,608 KB |
testcase_40 | AC | 17 ms
8,164 KB |
testcase_41 | AC | 249 ms
73,628 KB |
testcase_42 | AC | 19 ms
8,852 KB |
testcase_43 | AC | 226 ms
79,864 KB |
testcase_44 | AC | 216 ms
74,760 KB |
testcase_45 | AC | 227 ms
79,772 KB |
testcase_46 | AC | 210 ms
73,384 KB |
testcase_47 | AC | 222 ms
79,912 KB |
testcase_48 | AC | 96 ms
40,728 KB |
testcase_49 | AC | 222 ms
74,012 KB |
testcase_50 | AC | 223 ms
74,296 KB |
testcase_51 | AC | 304 ms
4,352 KB |
testcase_52 | AC | 183 ms
79,644 KB |
testcase_53 | AC | 172 ms
73,964 KB |
testcase_54 | AC | 175 ms
78,888 KB |
testcase_55 | AC | 16 ms
4,348 KB |
testcase_56 | AC | 15 ms
4,356 KB |
testcase_57 | AC | 16 ms
4,352 KB |
testcase_58 | AC | 19 ms
5,432 KB |
testcase_59 | AC | 22 ms
6,736 KB |
testcase_60 | AC | 22 ms
7,268 KB |
testcase_61 | AC | 15 ms
4,348 KB |
testcase_62 | AC | 15 ms
4,352 KB |
testcase_63 | AC | 10 ms
4,368 KB |
ソースコード
#include <bits/stdc++.h> #ifdef LOCAL #include <debug.hpp> #else #define debug(...) void(0) #endif using namespace std; typedef long long ll; #define all(x) begin(x), end(x) constexpr int INF = (1 << 30) - 1; constexpr long long IINF = (1LL << 60) - 1; constexpr int dx[4] = {1, 0, -1, 0}, dy[4] = {0, 1, 0, -1}; template <class T> istream& operator>>(istream& is, vector<T>& v) { for (auto& x : v) is >> x; return is; } template <class T> ostream& operator<<(ostream& os, const vector<T>& v) { auto sep = ""; for (const auto& x : v) os << exchange(sep, " ") << x; return os; } template <class T, class U = T> bool chmin(T& x, U&& y) { return y < x and (x = forward<U>(y), true); } template <class T, class U = T> bool chmax(T& x, U&& y) { return x < y and (x = forward<U>(y), true); } template <class T> void mkuni(vector<T>& v) { sort(begin(v), end(v)); v.erase(unique(begin(v), end(v)), end(v)); } template <class T> int lwb(const vector<T>& v, const T& x) { return lower_bound(begin(v), end(v), x) - begin(v); } int topbit(signed t) { return t == 0 ? -1 : 31 - __builtin_clz(t); } const int MAX_LOG = 30; /* 上から桁ごとに考えて良い r, c のいずれか 1 つの値を決めると全て決まる r[0] を決めるとする r[0] ^ (a[0][0] ^ a[i][0]) = r[i] r[0] ^ a[0][i] = c[i] r[0] = r[i] ^ (a[0][0] ^ a[i][0]) r[0] = c[i] ^ a[0][i] a[i][j] = r[i] ^ c[j] = a[0][0] ^ a[i][0] ^ a[0][j] r[i] = r[0] ^ x が 0, R[i], 2^k より大きくなったらそれぞれ 1, 1, inf のコスト */ struct node { array<int, 2> nxt; ll w; node() : nxt{-1, -1}, w(0) {} }; void solve() { int n, m; cin >> n >> m; vector<int> R(n), C(m); cin >> R >> C; vector a(n, vector<int>(m)); cin >> a; vector b(n, vector<int>(m)); vector<int> r(n), c(m); for (int i = 0; i < n; i++) { for (int j = 0; j < m; j++) { if (a[0][0] ^ a[i][0] ^ a[0][j] ^ a[i][j]) { cout << -1 << '\n'; return; } } } vector<tuple<int, int, int>> cs; for (int i = 0; i < n; i++) { int x = a[0][0] ^ a[i][0]; if (R[i] == 0) cs.emplace_back(x, 0, INF); else { cs.emplace_back(x, 0, 1); cs.emplace_back(x, R[i], 1); cs.emplace_back(x, (1 << (topbit(R[i]) + 1)) - 1, INF); } } for (int j = 0; j < m; j++) { int x = a[0][j]; if (C[j] == 0) cs.emplace_back(x, 0, INF); else { cs.emplace_back(x, 0, 1); cs.emplace_back(x, C[j], 1); cs.emplace_back(x, (1 << (topbit(C[j]) + 1)) - 1, INF); } } vector<node> nodes(1); auto add = [&](int x, int y, int z) { for (int i = MAX_LOG - 1, cur = 0; i >= 0; i--) { int f = x >> i & 1, g = y >> i & 1; for (int j = 0; j < 2; j++) { if (nodes[cur].nxt[j] == -1) { nodes[cur].nxt[j] = nodes.size(); nodes.emplace_back(); } } if (g == 0) { int ch = nodes[cur].nxt[f ^ g ^ 1]; nodes[ch].w += z; } cur = nodes[cur].nxt[f ^ g]; } }; for (auto [x, y, z] : cs) add(x, y, z); ll ans = INF; auto dfs = [&](auto self, int v, ll d) -> void { d += nodes[v].w; bool leaf = true; for (int i = 0; i < 2; i++) { if (nodes[v].nxt[i] == -1) continue; self(self, nodes[v].nxt[i], d); leaf = false; } if (leaf) ans = min(ans, d); }; dfs(dfs, 0, 0); cout << (ans == INF ? -1 : ans) << '\n'; } int main() { ios::sync_with_stdio(false); cin.tie(nullptr); int T; cin >> T; for (; T--;) solve(); return 0; }