#include #include #define rep(i,n) for(int i=0;i vi; typedef vector vl; typedef vector> vvi; typedef vector> vvl; typedef long double ld; typedef pair P; ostream& operator<<(ostream& os, const modint& a) {os << a.val(); return os;} template ostream& operator<<(ostream& os, const static_modint& a) {os << a.val(); return os;} template istream& operator>>(istream& is, vector& v){int n = v.size(); assert(n > 0); rep(i, n) is >> v[i]; return is;} template ostream& operator<<(ostream& os, const pair& p){os << p.first << ' ' << p.second; return os;} template ostream& operator<<(ostream& os, const vector& v){int n = v.size(); rep(i, n) os << v[i] << (i == n - 1 ? "\n" : " "); return os;} template ostream& operator<<(ostream& os, const vector>& v){int n = v.size(); rep(i, n) os << v[i] << (i == n - 1 ? "\n" : ""); return os;} template struct Edge_Dijkstra{ int from, to; T cost; Edge_Dijkstra(int from, int to, T cost) : from(from), to(to), cost(cost) {}; }; template struct Dijkstra{ const T INF = 1001001001001001; int n, m; vector initialized; vector> E; vector> G; map> dist; map> idx; Dijkstra(int _n) : n(_n), m(0), initialized(n, false), G(n){} void add_edge(int from, int to, T cost){ Edge_Dijkstra e(from, to, cost); E.push_back(e); G[from].emplace_back(m); m++; } void calc(int s){ initialized[s] = true; dist[s] = vector(n, INF); idx[s] = vector(n, -1); priority_queue, vector>, greater>> pq; pq.emplace(0, s, -1); while(pq.size()){ auto [cost, from, index] = pq.top(); pq.pop(); if(dist[s][from] <= cost) continue; dist[s][from] = cost; idx[s][from] = index; for(int index : G[from]){ int to = E[index].to; T cost_plus = E[index].cost; if(dist[s][to] <= cost + cost_plus) continue; pq.emplace(cost + cost_plus, to, index); } } } int farthest(int s){ if(!initialized[s]) calc(s); int idx = 0; rep(i, n) if(dist[s][i] > dist[s][idx]) idx = i; return idx; } T get_dist(int s, int t){ if(!initialized[s]) calc(s); return dist[s][t]; } vi restore(int s, int t){ if(!initialized[s]) calc(s); if(dist[s][t] == INF) return vi(0); vi res; while(idx[s][t] != -1){ auto e = E[idx[s][t]]; res.push_back(idx[s][t]); t = e.from; } reverse(res.begin(), res.end()); return res; } }; int main(){ int h, w; cin >> h >> w; vvi a(h, vi(h + w)); Dijkstra graph(h * (h + w)); rep(x, h) rep(y, h + w){ int from = x * (h + w) + y; if(y < h + w - 1){ int to = x * (h + w) + y + 1; graph.add_edge(from, to, 1); graph.add_edge(to, from, 1); } if(x < h - 1){ int to = (x + 1) * (h + w) + y; graph.add_edge(from, to, 1); graph.add_edge(to, from, 1); } if(y < h + w - 1 and x < h - 1){ int to = (x + 1) * (h + w) + y + 1; graph.add_edge(from, to, 1); graph.add_edge(to, from, 1); } } rep(i, w){ int x, y; cin >> x >> y; x--; y--; a[x][y]++; } int cnt = 0; long long ans = 0; rep(y, h + w) for(int x = h - 1; x >= 0; x--){ if(a[x][y] > 0){ ans += graph.get_dist(x * (h + w) + y, cnt); cnt++; } } cout << ans << "\n"; return 0; }