//#define _GLIBCXX_DEBUG #include using namespace std; #define rep(i, n) for(int i=0; i; using vs = vector; using vi = vector; using vvi = vector; template using PQ = priority_queue; template using PQG = priority_queue, greater>; const int INF = 0xccccccc; const ll LINF = 0xcccccccccccccccLL; template inline bool chmax(T1 &a, T2 b) {return a < b && (a = b, true);} template inline bool chmin(T1 &a, T2 b) {return a > b && (a = b, true);} template istream &operator>>(istream &is, pair &p) { return is >> p.first >> p.second;} template ostream &operator<<(ostream &os, const pair &p) { return os << p.first << ' ' << p.second;} struct edge{ int to, rev; int64_t cap; edge(int t=0, int64_t c=0, int r=0):to(t), cap(c), rev(r) {} }; struct max_flow { vector > G; vector level, iter; int n; max_flow(int n=0):n(n), G(n), level(n), iter(n) {} inline int add(int from, int to, int64_t cap) { G[from].emplace_back(to, cap, G[to].size()); G[to].emplace_back(from, 0, G[from].size()-1); return G[from].size()-1; } void bfs(int s) { for(int i = 0; i < n; i++) level[i] = -1; queue que; level[s] = 0; que.emplace(s); while(!que.empty()) { int now = que.front(); que.pop(); for(edge &e:G[now]) { if(e.cap > 0 and level[e.to] < 0) { level[e.to] = level[now] + 1; que.emplace(e.to); } } } } int64_t dfs(int v, int t, int64_t f) { if(v == t) return f; for(int &i = iter[v]; i < G[v].size(); i++) { edge &e = G[v][i]; if(e.cap > 0 and level[v] < level[e.to]) { int64_t 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; } int64_t Dinic(int s, int t) { int64_t flow = 0; while(true) { bfs(s); if(level[t] < 0) return flow; for(int i = 0; i < n; i++) iter[i] = 0; int64_t f; while((f = dfs(s, t, INT64_MAX)) > 0) flow += f; } } }; //head #define N 500010 vector

_[N]; int main() { ios::sync_with_stdio(false); cin.tie(0); int h, w; cin >> h >> w; rep(i, h) rep(j, w) { int a; cin >> a; _[a].emplace_back(i, j); } int ans = 0; for(int i = N-1; i > 0; i--) if(not _[i].empty()) { vi lst; for(auto [a, b]:_[i]) { lst.emplace_back(a); lst.emplace_back(b+N); } sort(all(lst)); lst.erase(unique(all(lst)), lst.end()); const int s = lst.size(), t = s+1; max_flow mfl(t+1); for(auto [a, b]:_[i]) { int ia = lower_bound(all(lst), a)-lst.begin(); int ib = lower_bound(all(lst), b+N)-lst.begin(); mfl.add(ia, ib, 1); } rep(i, s) { if(lst[i] >= N) mfl.add(i, t, 1); else mfl.add(s, i, 1); } ans += mfl.Dinic(s, t); } cout << ans << endl; }