#include #include #include #include #include #include #include #include #include #include #include using namespace std; typedef long long ll; const int INF = INT_MAX; struct Edge { int to; int cap; int rev; Edge(int to = -1, int cap = -1, int rev = -1) { this->to = to; this->cap = cap; this->rev = rev; } }; const int MAX_N = 100; const int MAX_V = MAX_N * MAX_N + 10; class MaxFlow { public: int level[MAX_V]; int iter[MAX_V]; vector G[MAX_V]; void add_edge(int from, int to, int cap) { G[from].push_back(Edge(to, cap, G[to].size())); G[to].push_back(Edge(from, 0, G[from].size() - 1)); } void bfs(int s) { memset(level, -1, sizeof(level)); queue que; level[s] = 0; que.push(s); while (!que.empty()) { int v = que.front(); que.pop(); for (int i = 0; i < (int) G[v].size(); ++i) { Edge *e = &G[v][i]; if (e->cap > 0 && level[e->to] < 0) { level[e->to] = level[v] + 1; que.push(e->to); } } } } int dfs(int v, int t, int f) { if (v == t) return f; for (int &i = iter[v]; i < (int) G[v].size(); ++i) { Edge *e = &G[v][i]; if (e->cap > 0 && level[v] < level[e->to]) { int d = dfs(e->to, t, min(f, e->cap)); if (d > 0) { e->cap -= d; G[e->to][e->rev].cap += d; return d; } } } return 0; } int max_flow(int s, int t) { int flow = 0; for (;;) { bfs(s); if (level[t] < 0) return flow; memset(iter, 0, sizeof(iter)); int f; while ((f = dfs(s, t, INF)) > 0) { flow += f; } } } }; int main() { int H, W; cin >> H >> W; map> S; for (int y = 0; y < H; ++y) { for (int x = 0; x < W; ++x) { int a; cin >> a; int z = y * W + x; S[a].push_back(z); } } int s = H + W + 1; int t = H + W + 2; int ans = 0; for (auto &[a, positions] : S) { MaxFlow mf; bool checked[H + W]; memset(checked, false, sizeof(checked)); if (a == 0) continue; for (int z : positions) { int y = z / W; int x = z % W; if (not checked[y]) { checked[y] = true; mf.add_edge(s, y, 1); } if (not checked[H + x]) { checked[H + x] = true; mf.add_edge(H + x, t, 1); } mf.add_edge(y, H + x, 1); } ans += mf.max_flow(s, t); } cout << ans << endl; return 0; }