#pragma GCC optimize("O3") //#pragma GCC target("avx2") //#pragma GCC optimize("unroll-loops") #include using namespace std; typedef long long ll; typedef unsigned long long ull; typedef long double ld; typedef vector vi; typedef vector vll; typedef vector vld; typedef pair Pll; typedef pair Pin; ll INF = 1e16; int inf = 1e9; #define ALL(x) (x).begin(), (x).end() #define FOR(i, m, n) for (ll i = (m); i < (n); ++i) #define REVFOR(i, m, n) for (ll i = (n - 1); i >= (m); --i) #define REP(i, n) FOR(i, 0, n) #define REVREP(i, n) REVFOR(i, 0, n) #define fi first #define se second #define pb push_back #define mp make_pair #define eb emplace_back #define bcnt __builtin_popcountll #ifdef LOCAL #include #define debug(...) cerr << "[" << #__VA_ARGS__ << "]: ", d_err(__VA_ARGS__); #else #define debug(...) 83; #endif void d_err() { cerr << endl; } template void d_err(H h, T... t) { cerr << h << " "; d_err(t...); } template void print(T x) { cout << x << "\n"; } template void print(vector& x) { int N = x.size(); REP(i, N) { if (i > 0) cout << " "; cout << x[i]; } cout << "\n"; } // Dinic template struct Dinic { struct Graph { int to, rev; T cap; Graph(int t, int r, const T &c) : to(t), rev(r), cap(c) {}; }; vector> g; vector level; vector iter; int sz; Dinic(int N) : sz(N) { g.resize(N); iter.resize(N); level.resize(N); } void add_edge(int from, int to, T cap) { g[from].eb(to, (int) g[to].size(), cap); g[to].eb(from, (int) g[from].size() - 1, 0); } T operator()(int start, int end) { T flow = 0; while(true) { bfs(start); if (level[end] < 0) break; iter.assign(sz, 0); T f; while(f = dfs(start, end, inf) > 0) { flow += f; } } return flow; } private: void bfs(int s) { level.assign(sz, -1); level[s] = 0; queue q; q.push(s); while(!q.empty()) { int v = q.front(); q.pop(); for(auto &e: g[v]) { int t = e.to; if (e.cap > 0 && level[t] < 0) { level[t] = level[v] + 1; q.push(t); } } } } T dfs(int v, int t, T flow) { if (v == t || flow == 0) return flow; for (int &i = iter[v]; i < g[v].size(); ++i) { auto &e = g[v][i]; if(e.cap <= 0 || level[v] >= level[e.to]) continue; if (T f = dfs(e.to, t, min(flow, e.cap)) > 0) { e.cap -= f; g[e.to][e.rev].cap += f; return f; } } return 0; } }; int main(){ cin.tie(0); ios_base::sync_with_stdio(false); cout << fixed << setprecision(20); int h, w; cin >> h >> w; vector a(h, vll(w)); REP(i, h) REP(j, w) cin >> a[i][j]; vector> v(1e6); REP(i, h) REP(j, w) { if (a[i][j] == 0) continue; v[a[i][j]].eb(i, j); } ll ans = 0; for (auto &e: v) { if (e.empty()) continue; map mi, mj; for(auto &e2: e) { if (mi.find(e2.fi) == mi.end()) mi[e2.fi] = mi.size(); if (mj.find(e2.se) == mj.end()) mj[e2.se] = mj.size(); } Dinic g(mi.size() + mj.size() + 2); for (auto &e2: e) { g.add_edge(mi[e2.fi], mi.size() + mj[e2.se], 1); } for (auto &p: mi) { g.add_edge(mi.size() + mj.size(), p.se, 1); } for (auto &p: mj) { g.add_edge(mi.size() + p.se, mi.size() + mj.size() + 1, 1); } ans += g(mi.size() + mj.size(), mi.size() + mj.size() + 1); } print(ans); }