#include #include #include #include #include #include #include #include #include #include int KthBit(const int x, const int k) { return x >> k & 1; } // // できる限り枝刈りを行っている。 template int Solve(const std::vector& rc, const std::vector>& a) { const int n = a.size(), m = a.front().size(); assert(std::ssize(rc) == n + m); // ビットごと 2 グループに分ける std::vector group(n + m, std::bitset()); for (int bit = 0; bit < D; ++bit) { const int criteria = KthBit(a.front().front(), bit); for (int i = 1; i < n; ++i) { const int is_diff = (KthBit(a[i].front(), bit) != criteria ? 1 : 0); for (int j = 1; j < m; ++j) { if (KthBit(a.front()[j], bit) ^ is_diff != KthBit(a[i][j], bit)) { return -1; } } if (is_diff == 1) group[i].set(bit); } if (criteria == 1) group[n].set(bit); for (int j = 1; j < m; ++j) { const int is_diff = (KthBit(a.front()[j], bit) != criteria ? 1 : 0); for (int i = 1; i < n; ++i) { if (KthBit(a[i].front(), bit) ^ is_diff != KthBit(a[i][j], bit)) { return -1; } } if (KthBit(a.front()[j], bit) == 1) group[j + n].set(bit); } } const auto Operation = [&rc, &group]( const int bit, const int choice, const std::vector& or_012, const std::vector& or_01, const std::vector& or_12) -> std::tuple, std::vector, std::vector> { int cost = 0; std::vector next_or_012, next_or_01, next_or_12; for (const int i : or_012) { if (group[i][bit] == choice) { // 1 if (KthBit(rc[i], bit) == 0) { return {false, 0, std::vector{}, std::vector{}, std::vector{}}; } ++cost; next_or_12.emplace_back(i); } else { // 0 (KthBit(rc[i], bit) == 1 ? next_or_01 : next_or_012).emplace_back(i); } } for (const int i : or_01) { if (group[i][bit] == choice) { ++cost; } else { next_or_01.emplace_back(i); } } for (const int i : or_12) { const int bit_of_i = (group[i][bit] == choice ? 1 : 0); const int bit_of_rc = KthBit(rc[i], bit); if (bit_of_i > bit_of_rc) { ++cost; } else if (bit_of_i == bit_of_rc) { next_or_12.emplace_back(i); } } return {true, cost, next_or_012, next_or_01, next_or_12}; }; const int argmin = std::distance(rc.begin(), std::ranges::min_element(rc)); const int need_to_choose = // これより上位のビットは自由に選べない std::bit_width(static_cast(rc[argmin])); int current_cost = 0; std::vector or_012(n + m), or_01, or_12; std::iota(or_012.begin(), or_012.end(), 0); for (int bit = D - 1; bit >= need_to_choose; --bit) { const int choice = group[argmin][bit] ^ 1; bool is_valid; int added_cost; std::tie(is_valid, added_cost, or_012, or_01, or_12) = Operation(bit, choice, or_012, or_01, or_12); if (!is_valid) return -1; current_cost += added_cost; } int min_cost = (n + m) * 2 + 1; // 上位からグループを決定する const auto dfs = [Operation, &min_cost]( auto dfs, const int bit, const int current_cost, const std::vector& or_012, const std::vector& or_01, const std::vector& or_12) -> void { if ((or_012.empty() && or_01.empty() && or_12.empty()) || bit < 0) { min_cost = std::min(min_cost, current_cost); return; } std::array next_cost{}; std::array, 2> next_or_012{}, next_or_01{}, next_or_12{}; for (const int choice : std::array{0, 1}) { std::tie(std::ignore, next_cost[choice], next_or_012[choice], next_or_01[choice], next_or_12[choice]) = Operation(bit, choice, or_012, or_01, or_12); next_cost[choice] += current_cost; } const int earlier = (next_cost[0] + next_or_012[0].size() + next_or_01[0].size() < next_cost[1] + next_or_012[1].size() + next_or_01[1].size() ? 0 : 1); for (const int parity : std::array{0, 1}) { const int choice = parity ^ earlier; if (next_cost[choice] >= min_cost) continue; // 枝刈り dfs(dfs, bit - 1, next_cost[choice], next_or_012[choice], next_or_01[choice], next_or_12[choice]); } }; dfs(dfs, need_to_choose - 1, current_cost, or_012, or_01, or_12); return min_cost; } int main() { constexpr int kMaxT = 50000, kMaxNM = 50000, kDigits = 30; int t; std::cin >> t; assert(1 <= t && t <= kMaxT); while (t--) { int n, m; std::cin >> n >> m; assert(1 <= n && 1 <= m && n * m <= kMaxNM); std::vector rc(n + m); for (int i = 0; i < n + m; ++i) { std::cin >> rc[i]; assert(0 <= rc[i] && rc[i] < (1 << kDigits)); } std::vector a(n, std::vector(m, 0)); for (int i = 0; i < n; ++i) { for (int j = 0; j < m; ++j) { std::cin >> a[i][j]; assert(0 <= a[i][j] && a[i][j] < (1 << kDigits)); } } std::cout << Solve(rc, a) << '\n'; } return 0; }