#include using namespace std; template struct HLPP { const T INF = numeric_limits::max(); struct edge { int to, rev; T f; }; int s = MAXN - 1, t = MAXN - 2; vector adj[MAXN]; vector lst[MAXN], gap[MAXN]; T excess[MAXN]; int highest, height[MAXN], cnt[MAXN], work; void addEdge(int from, int to, int f, bool isDirected = true) { adj[from].push_back({to, adj[to].size(), f}); adj[to].push_back({from, adj[from].size() - 1, isDirected ? 0 : f}); } void updHeight(int v, int nh) { work++; if (height[v] != MAXN) cnt[height[v]]--; height[v] = nh; if (nh == MAXN) return; cnt[nh]++, highest = nh; gap[nh].push_back(v); if (excess[v] > 0) lst[nh].push_back(v); } void globalRelabel() { work = 0; fill_n(height, MAXN, MAXN); fill_n(cnt, MAXN, 0); for (int i = 0; i < highest; i++) lst[i].clear(), gap[i].clear(); height[t] = 0; queue q({t}); while (!q.empty()) { int v = q.front(); q.pop(); for (auto &e : adj[v]) if (height[e.to] == MAXN && adj[e.to][e.rev].f > 0) q.push(e.to), updHeight(e.to, height[v] + 1); highest = height[v]; } } void push(int v, edge &e) { if (excess[e.to] == 0) lst[height[e.to]].push_back(e.to); T df = min(excess[v], e.f); e.f -= df, adj[e.to][e.rev].f += df; excess[v] -= df, excess[e.to] += df; } void discharge(int v) { int nh = MAXN; for (auto &e : adj[v]) { if (e.f > 0) { if (height[v] == height[e.to] + 1) { push(v, e); if (excess[v] <= 0) return; } else nh = min(nh, height[e.to] + 1); } } if (cnt[height[v]] > 1) updHeight(v, nh); else { for (int i = height[v]; i <= highest; i++) { for (auto j : gap[i]) updHeight(j, MAXN); gap[i].clear(); } } } T calc(int heur_n = MAXN) { fill_n(excess, MAXN, 0); excess[s] = INF, excess[t] = -INF; globalRelabel(); for (auto &e : adj[s]) push(s, e); for (; highest >= 0; highest--) { while (!lst[highest].empty()) { int v = lst[highest].back(); lst[highest].pop_back(); discharge(v); if (work > 4 * heur_n) globalRelabel(); } } return excess[t] + INF; } }; int n, m; int G[305][305]; int R[305]; int C[305]; HLPP<100000, long long> D; int main() { scanf("%d %d", &n, &m); for(int i = 0;i < n;i++){ for(int j = 0;j < m;j++){ scanf("%d", &G[i][j]); } } for(int i = 0;i < n;i++){ scanf("%d", &R[i]); } for(int i = 0;i < m;i++){ scanf("%d", &C[i]); } for(int i = 0;i < n;i++){ for(int j = 0;j < m;j++){ D.addEdge(i*m+j, D.t, G[i][j]); } } long long ans = 0; for(int i = 0;i < n;i++){ D.addEdge(D.s, n*m+i, R[i]); ans += R[i]; for(int j = 0;j < m;j++){ D.addEdge(n*m+i, i*m+j, R[i]); } } for(int i = 0;i < m;i++){ D.addEdge(D.s, n*m+n+i, C[i]); ans += C[i]; for(int j = 0;j < n;j++){ D.addEdge(n*m+n+i, j*m+i, C[i]); } } ans -= D.calc(); cout << ans << endl; return 0; }