#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; using ll = long long; template auto range(T s, T e) { return views::iota(s, max(s, e)); } template auto range(T n) { return range(0, n); } template void take(vector& vec, int n) { vec.resize(n); for (int i = 0; i < n; ++i) cin >> vec[i]; } template void sout(const Args &...args) { ((cout << args << ' '), ...); } template void soutn(const Args &...args) { ((cout << args << ' '), ...); cout << '\n'; } template struct In2 { T1 a; T2 b; friend std::istream& operator>>(std::istream& is, In2& obj) { T1 t1; T2 t2; is >> t1 >> t2; obj = {t1, t2}; return is; } }; template struct In3 { T1 a; T2 b; T3 c; friend std::istream& operator>>(std::istream& is, In3& obj) { T1 t1; T2 t2; T3 t3; is >> t1 >> t2 >> t3; obj = {t1, t2, t3}; return is; } }; template struct In4 { T1 a; T2 b; T3 c; T4 d; friend std::istream& operator>>(std::istream& is, In4& obj) { T1 t1; T2 t2; T3 t3; T4 t4; is >> t1 >> t2 >> t3 >> t4; obj = {t1, t2, t3, t4}; return is; } }; #ifdef LOCAL #include #else #define dump(...) ; #endif namespace { ll n, w, h; vector> va; void read() { cin >> n >> w >> h; take(va, n); } ll sol(ll t, vector vx) { if (t == 1) return 0; map slope; auto eval = [&](ll p) -> ll { ll tot = 0; for (ll x : vx) { ll d = abs(x - p); tot += min(d, t - d); } return tot; }; ll val0 = eval(0); ll val1 = eval(1); for (ll x : vx) { slope[x] += 2; ll half = t / 2; if (t % 2) { slope[x - half - 1]--; slope[x - half]--; slope[x + half]--; slope[x + half + 1]--; } else { slope[x - half] -= 2; slope[x + half] -= 2; } } dump(t, vx); dump(slope); ll val = val0, s = val1 - val0, x0 = 0; ll res = val; for (auto [x, sdiff] : slope) { if (x <= 0) continue; if (x >= t) break; val += (x - x0) * s; res = max(res, val); s += sdiff; x0 = x; dump(x, "=>", val, "slope=", s); } return res; } ll run() { vector vx, vy; for (auto [x, y] : va) { vx.push_back(x - 1); vy.push_back(y - 1); } ll res = sol(w, vx); res += sol(h, vy); return res; } } // namespace template void exec(F f) { if constexpr (std::is_same_v) f(); else cout << f() << endl; } int main(int argc, char** argv) { cerr << fixed << setprecision(12); cout << fixed << setprecision(12); int testcase = 1; if (argc > 1) testcase = atoi(argv[1]); while (testcase--) { read(); } exec(run); }