#line 1 "library-cpp/other/template.hpp" // clang-format off #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; using uint = unsigned int; using ll = long long; using ull = unsigned long long; using i128 = __int128_t; using ld = long double; using pii = pair; using pll = pair; template using maxheap = priority_queue; template using minheap = priority_queue, greater>; template using vec = vector; template using vvec = vector>; #define OVERLOAD_REP(_1, _2, _3, name, ...) name #define REP0(n) for (auto minato = decay_t{}; minato < (n); ++minato) #define REP1(i, n) for (auto i = decay_t{}; (i) < (n); (i)++) #define REP2(i, l, r) for (auto i = (l); (i) < (r); (i)++) #define rep(...) OVERLOAD_REP(__VA_ARGS__, REP2, REP1, REP0)(__VA_ARGS__) #define OVERLOAD_RREP(_1, _2, _3, name, ...) name #define RREP1(i, n) for (auto i = (n) - 1; (i) >= decay_t{}; (i)--) #define RREP2(i, l, r) for (auto i = (r) - 1; (i) >= (l); (i)--) #define rrep(...) OVERLOAD_RREP(__VA_ARGS__, RREP2, RREP1)(__VA_ARGS__) #define ALL(x) begin(x), end(x) template int SZ(const Container& v) { return int(v.size()); } template void UNIQUE(vector& v) { v.erase(unique(v.begin(), v.end()), v.end()); } template T MAX(const vector& v) { return *max_element(v.begin(), v.end()); } template T MIN(const vector& v) { return *min_element(v.begin(), v.end()); } template T SUM(const vector& v) { return accumulate(v.begin(), v.end(), T(0)); } template T ABS(T x) { return max(x, -x); } template constexpr T POW(T x, ull n) { T ret = 1; while (n > 0) { if (n & 1) ret *= x; x *= x; n >>= 1; } return ret; } template constexpr T POW(T x, ull n, T mod) { T ret = 1; while (n > 0) { if (n & 1) ret = ret * x % mod; x = x * x % mod; n >>= 1; } return ret; } constexpr long long floor_div(long long n, long long d) { return n / d - ((n ^ d) < 0 && n % d); } template bool chmax(T1& a, T2 b) { if (a < b) { a = b; return true; } return false; } template bool chmin(T1& a, T2 b) { if (a > b) { a = b; return true; } return false; } constexpr int topbit(ull x) { return 63 - std::countl_zero(x); } constexpr int botbit(ull x) { return std::countr_zero(x); } constexpr int popcnt(ull x) { return std::popcount(x); } constexpr int kthbit(ull x, int k) { return (x >> k) & 1; } constexpr int popparity(ull x) { return __builtin_parityll(x); } constexpr int parity_sign(ull x) { return (x & 1) ? -1 : 1; } constexpr long long TEN(int x) { return x == 0 ? 1 : TEN(x - 1) * 10; } template void rearrange(const vector& id) { (void)id; } template void rearrange_exec(const vector& id, vector& v) { vector w(v.size()); for (size_t i = 0; i < id.size(); i++) { w[i] = v[id[i]]; } v.swap(w); } template void rearrange(const vector& id, Head& a, Tail& ...tail) { rearrange_exec(id, a); rearrange(id, tail...); } istream& operator>>(istream& is, __int128_t& x) { x = 0; string s; is >> s; int n = int(s.size()), it = 0; if (s[0] == '-') it++; for (; it < n; it++) x = (x * 10 + s[it] - '0'); if (s[0] == '-') x = -x; return is; } ostream& operator<<(ostream& os, __int128_t x) { if (x == 0) return os << 0; if (x < 0) os << '-', x = -x; deque deq; while (x) deq.emplace_front(x % 10), x /= 10; for (int e : deq) os << e; return os; } template vector &operator++(vector& v) { for (auto& e : v) { e++; } return v;} template vector operator++(vector& v, int) { auto res = v; for (auto& e : v) { e++; } return res; } template vector &operator--(vector& v) { for (auto& e : v) { e--; } return v; } template vector operator--(vector& v, int) { auto res = v; for (auto& e : v) { e--; } return res; } template pair operator-(const pair& x) { return pair(-x.first, -x.second); } template pair operator-(const pair& x, const pair& y) { return pair(x.first - y.first, x.second - y.second); } template pair operator+(const pair& x, const pair& y) { return pair(x.first + y.first, x.second + y.second); } template pair operator+=(pair& l, const pair& r) { return l = l + r; } template pair operator-=(pair& l, const pair& r) { return l = l - r; } constexpr char ln = '\n'; const string YESNO[2] = {"NO", "YES"}; const string YesNo[2] = {"No", "Yes"}; void YES(bool t = true) { cout << YESNO[t] << "\n"; } void NO(bool t = 1) { YES(!t); } void Yes(bool t = true) { cout << YesNo[t] << "\n"; } void No(bool t = 1) { Yes(!t); } template void drop(T x) { cout << x << "\n"; exit(0); } #define INT(...) \ int __VA_ARGS__; \ IN(__VA_ARGS__) #define LL(...) \ ll __VA_ARGS__; \ IN(__VA_ARGS__) #define STR(...) \ string __VA_ARGS__; \ IN(__VA_ARGS__) #define CHR(...) \ char __VA_ARGS__; \ IN(__VA_ARGS__) #define LDB(...) \ long double __VA_ARGS__; \ IN(__VA_ARGS__) #define VEC(type, name, size) \ vector name(size); \ IN(name) #define VEC2(type, name1, name2, size) \ vector name1(size), name2(size); \ for (int i = 0; i < size; i++) IN(name1[i], name2[i]) #define VEC3(type, name1, name2, name3, size) \ vector name1(size), name2(size), name3(size); \ for (int i = 0; i < size; i++) IN(name1[i], name2[i], name3[i]) #define VEC4(type, name1, name2, name3, name4, size) \ vector name1(size), name2(size), name3(size), name4(size); \ for (int i = 0; i < size; i++) IN(name1[i], name2[i], name3[i], name4[i]); #define VV(type, name, N, M) \ vector> name(N, vector(M)); \ IN(name) template void scan(T& a) { cin >> a; } template void scan(vector& a) { for (auto& i : a) scan(i); } void IN() {} template void IN(Head& head, Tail&... tail) { scan(head); IN(tail...); } std::vector> read_graph(int N, int M) { std::vector> g(N); for (int i = 0; i < M; i++) { int a, b; cin >> a >> b; a--, b--; g[a].emplace_back(b); g[b].emplace_back(a); } return g; } std::vector> read_tree(int N) { return read_graph(N, N - 1); } void print() { cout << "\n"; } template void print(const vector& v) { for (auto it = v.begin(); it != v.end(); ++it) { if (it != v.begin()) { cout << " "; } cout << *it; } print(); } template void print(const T& x, const Args& ... args) { cout << x; if (sizeof...(Args)) cout << " "; print(args...); } #ifdef MINATO_LOCAL template ostream& operator<<(ostream& os, pair p); template ostream& operator<<(ostream& os, tuple t); template ostream& operator<<(ostream& os, vector v); template ostream& operator<<(ostream& os, array a); template enable_if_t>, ostream>& operator<<(ostream& os, T(&a)[N]); template ostream& operator<<(ostream& os, set s); template ostream& operator<<(ostream& os, map mp); template ostream& operator<<(ostream& os, multiset s); template ostream& operator<<(ostream& os, queue q); template ostream& operator<<(ostream& os, deque q); template ostream& operator<<(ostream& os, priority_queue q); template ostream& operator<<(ostream& os, priority_queue, greater> q); template ostream& operator<<(ostream& os, pair p) { return os << "(" << p.first << ", " << p.second << ")"; } template void debug_tuple(ostream& os, TUPLE _) { (void)os; (void)_; } template void debug_tuple(ostream &os, TUPLE t) { os << (N == 0 ? "" : ", ") << get(t); debug_tuple(os, t); } template ostream& operator<<(ostream& os, tuple t) { os << "("; debug_tuple<0, tuple, Args...>(os, t); return os << ")"; } string debug_delim(int& i) { return i++ == 0 ? "" : ", "; } #define debug_embrace(x) { int i = 0; os << "{"; { x } return os << "}"; } template ostream& operator<<(ostream& os, vector v) { debug_embrace( for (T e : v) { os << debug_delim(i) << e; } ) } template ostream& operator<<(ostream& os, array a) { debug_embrace( for (T e : a) { os << debug_delim(i) << e; } ) } template enable_if_t>, ostream>& operator<<(ostream& os, T(&a)[N]) { debug_embrace( for (T e : a) { os << debug_delim(i) << e; } ) } template ostream& operator<<(ostream& os, set s) { debug_embrace( for (Key e : s) { os << debug_delim(i) << e; }) } template ostream& operator<<(ostream& os, map mp) { debug_embrace( for (auto e : mp) { os << debug_delim(i) << e; }) } template ostream& operator<<(ostream& os, multiset s) { debug_embrace( for (Key e : s) { os << debug_delim(i) << e; }) } template ostream& operator<<(ostream& os, queue q) { debug_embrace( for (; !q.empty(); q.pop()) { os << debug_delim(i) << q.front(); } ) } template ostream& operator<<(ostream& os, deque q) { debug_embrace( for (T e : q) { os << debug_delim(i) << e; } ) } template ostream& operator<<(ostream& os, priority_queue q) { debug_embrace( for (; !q.empty(); q.pop()) { os << debug_delim(i) << q.top(); } ) } template ostream& operator<<(ostream& os, priority_queue, greater> q) { debug_embrace( for (; !q.empty(); q.pop()) { os << debug_delim(i) << q.top(); } ) } void debug_out() { cerr << endl; } template void debug_out(const T& x, const Args& ... args) { cerr << " " << x; debug_out(args...); } #define debug(...) cerr << __LINE__ << " : [" << #__VA_ARGS__ << "] =", debug_out(__VA_ARGS__) #else #define debug(...) (void(0)) #endif /////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// // clang-format on #line 2 "library-cpp/flow/ProjectSelectionProblem.hpp" #line 6 "library-cpp/flow/ProjectSelectionProblem.hpp" #include template struct ProjectSelectionProblem { struct Edge { int from, to; Cap cost; Edge(int from, int to, Cap cost) : from(from), to(to), cost(cost) { } }; int n, mf_n, s, t; Cap base; std::vector edges; atcoder::mf_graph g; ProjectSelectionProblem() { } ProjectSelectionProblem(int n_) : n(n_), mf_n(n_ + 2), s(n_), t(n_ + 1), base(0) { } /** * 頂点 u が赤色かつ頂点 v が青色だとペナルティを受ける */ void if_u_red_and_v_blue_then_penalty(int u, int v, Cap penalty) { assert(0 <= u && u < n && 0 <= v && v < n && u != v && penalty >= 0); edges.emplace_back(u, v, penalty); } /** * 頂点 v が赤色になるとペナルティを受ける */ void red_penalty(int v, Cap penalty) { assert(0 <= v && v < n && penalty >= 0); edges.emplace_back(v, t, penalty); } /** * 頂点 v が青色になるとペナルティを受ける */ void blue_penalty(int v, Cap penalty) { assert(0 <= v && v < n && penalty >= 0); edges.emplace_back(s, v, penalty); } /** * 頂点 u と v の色が異なるとペナルティを受ける */ void difference_penalty(int u, int v, Cap penalty) { assert(0 <= u && u < n && 0 <= v && v < n && u != v && penalty >= 0); edges.emplace_back(u, v, penalty); edges.emplace_back(v, u, penalty); } /** * 頂点 u が赤色になると頂点 v も赤色になる */ void if_u_red_then_v_red(int u, int v) { assert(0 <= u && u < n && 0 <= v && v < n && u != v); edges.emplace_back(u, v, std::numeric_limits::max()); } /** * 頂点 u が青色になると頂点 v も青色になる */ void if_u_blue_then_v_blue(int u, int v) { assert(0 <= u && u < n && 0 <= v && v < n && u != v); edges.emplace_back(v, u, std::numeric_limits::max()); } /** * 頂点 v が赤色になる */ void force_red(int v) { assert(0 <= v && v < n); edges.emplace_back(s, v, std::numeric_limits::max()); } /** * 頂点 v が青色になる */ void force_blue(int v) { assert(0 <= v && v < n); edges.emplace_back(v, t, std::numeric_limits::max()); } /** * 頂点 u と v が同じ色になる */ void force_same(int u, int v) { assert(0 <= u && u < n && 0 <= v && v < n && u != v); edges.emplace_back(u, v, std::numeric_limits::max()); edges.emplace_back(v, u, std::numeric_limits::max()); } /** * 頂点 v が赤色になると報酬を得る */ void red_reward(int v, Cap reward) { assert(0 <= v && v < n && reward >= 0); base += reward; edges.emplace_back(s, v, reward); } /** * 頂点 v が青色になると報酬を得る */ void blue_reward(int v, Cap reward) { assert(0 <= v && v < n && reward >= 0); base += reward; edges.emplace_back(v, t, reward); } /** * 頂点 u と v が同じ色になると報酬を得る */ void same_reward(int u, int v, Cap reward) { assert(0 <= u && u < n && 0 <= v && v < n && u != v && reward >= 0); base += reward; edges.emplace_back(u, v, reward); edges.emplace_back(v, u, reward); } /** * 赤色の頂点が存在するとペナルティを受ける */ void any_red_peanlty(const std::vector& vs, Cap penalty) { int u = mf_n++; for (int v : vs) { assert(0 <= v && v < n); edges.emplace_back(v, u, std::numeric_limits::max()); } edges.emplace_back(u, t, penalty); } /** * 青色の頂点が存在するとペナルティを受ける */ void any_blue_penalty(const std::vector& vs, Cap penalty) { int u = mf_n++; for (int v : vs) { assert(0 <= v && v < n); edges.emplace_back(u, v, std::numeric_limits::max()); } edges.emplace_back(s, u, penalty); } /** * すべての頂点が赤になると報酬を得る */ void all_red_reward(const std::vector& vs, Cap reward) { base += reward; int u = mf_n++; for (int v : vs) { assert(0 <= v && v < n); edges.emplace_back(u, v, std::numeric_limits::max()); } edges.emplace_back(s, u, reward); } /** * すべての頂点が青になると報酬を得る */ void all_blue_reward(const std::vector& vs, Cap reward) { base += reward; int u = mf_n++; for (int v : vs) { assert(0 <= v && v < n); edges.emplace_back(v, u, std::numeric_limits::max()); } edges.emplace_back(u, t, reward); } /** * 報酬の最大値を求める */ Cap solve() { g = atcoder::mf_graph(mf_n); for (const auto& [u, v, cost] : edges) { g.add_edge(u, v, cost); } return base - g.flow(s, t); } std::vector red_selection() { auto selection = g.min_cut(s); selection.resize(n); return selection; } }; #line 3 "A.cpp" void solve() { INT(H, W); VV(ll, G, H, W); VEC(ll, R, H); VEC(ll, C, W); ProjectSelectionProblem psp(H + W); rep(i, H) { ll sum = 0; rep(j, W) { sum += G[i][j]; } ll val = R[i] - sum; if (val > 0) { psp.red_reward(i, val); } else if (val < 0) { psp.red_penalty(i, -val); } } rep(j, W) { ll sum = 0; rep(i, H) { sum += G[i][j]; } ll val = C[j] - sum; if (val > 0) { psp.red_reward(H + j, val); } else if (val < 0) { psp.red_penalty(H + j, -val); } } rep(i, H) rep(j, W) { psp.all_red_reward({i, H + j}, G[i][j]); } ll ans = psp.solve(); print(ans); } int main() { cin.tie(nullptr); ios::sync_with_stdio(false); cout << fixed << setprecision(20); cerr << fixed << setprecision(7); int T = 1; // cin >> T; for (int test_case = 1; test_case <= T; test_case++) { // debug(test_case); solve(); } }