#include using namespace atcoder; #include using namespace internal; #include using namespace std; #include #pragma GCC target("avx2") #pragma GCC optimize("O3") #pragma GCC optimize("unroll-loops") #define rep(i, n) for (int i = 0; i < (int)(n); i++) #define rep2(i,a,b) for (int i = (int)(a); i < (int)(b); i++) #define all(v) v.begin(),v.end() #define inc(x,l,r) ((l)<=(x)&&(x)<(r)) #define Unique(x) sort(all(x)), x.erase(unique(all(x)), x.end()) #define pcnt __builtin_popcountll typedef long long ll; #define int ll using ld = long double; using vi = vector; using vs = vector; using P = pair; using vp = vector

; using Bint = boost::multiprecision::cpp_int; template bool chmax(T1 &a, const T2 b) {if (a < b) {a = b; return true;} else return false; } template bool chmin(T1 &a, const T2 b) {if (a > b) {a = b; return true;} else return false; } template using priority_queue_greater = priority_queue, greater>; template ostream &operator<< (ostream &os, const pair &p){os << p.first <<" "<> (istream &is, modint1000000007 &m){ll in;is>>in;m=in;return is;} ostream &operator<< (ostream &os, const modint998244353 &m){os << m.val();return os;} istream &operator>> (istream &is, modint998244353 &m){ll in;is>>in;m=in;return is;} template istream &operator>>(istream& is,vector &v){for(T &in:v)is>>in;return is;} template void input(T&... a){(cin>> ... >> a);} #ifdef LOCAL template ostream &operator<<(ostream &os,const vector &v){os<<"\x1b[32m";rep(i,v.size())os< int print(T& a){cout << "\x1b[32m"<< a<< '\n' << "\x1b[0m";return 0;} template int print(const T&a, const Ts&... b){cout << "\x1b[32m" << a;(cout<<...<<(cout<<' ',b));cout<<'\n' << "\x1b[0m";return 0;} #else template ostream &operator<<(ostream &os,const vector &v){rep(i,v.size())os< int print(T& a){cout << a<< '\n';return 0;} template int print(const T&a, const Ts&... b){cout << a;(cout<<...<<(cout<<' ',b));cout<<'\n';return 0;} #endif #define VI(v,n) vi v(n); input(v) #define INT(...) int __VA_ARGS__; input(__VA_ARGS__) #define STR(...) string __VA_ARGS__; input(__VA_ARGS__) #define CHAR(...) char __VA_ARGS__; input(__VA_ARGS__) int sign(ll x){return x>0?1:x<0?-1:0;} ll ceil(ll x,ll y){assert(y!=0);if(sign(x)==sign(y))return (x+y-1)/y;return -((-x/y));} ll floor(ll x,ll y){assert(y!=0);if(sign(x)==sign(y))return x/y;if(y<0)x*=-1,y*=-1;return x/y-(x%y<0);} ll abs(ll x,ll y){return abs(x-y);} ll bit(int n){return 1ll< bool ins(string s,T t){return s.find(t)!=string::npos;} P operator+ (const P &p, const P &q){ return P{p.first+q.first,p.second+q.second};} P operator- (const P &p, const P &q){ return P{p.first-q.first,p.second-q.second};} int yesno(bool ok,string y="Yes",string n="No"){ cout<<(ok?y:n)< bit, sum; BitVector(){} BitVector(int len) : len(len), blocks((len+w-1)/w){ bit.assign(blocks, 0); sum.assign(blocks, 0); } void set(const int k, const int b = 1){ if(b) bit[k / w] |= 1U << rem(k); else bit[k / w] &= ~(1U << rem(k)); } void build(){ sum[0] = 0; for(int i = 1; i < blocks; i++){ sum[i] = sum[i - 1] + __builtin_popcount(bit[i - 1]); } } int operator[](const int k) const{ return bit[k / w] >> rem(k) & 1; } inline int rank(const int k) const{ return sum[k / w] + __builtin_popcount(bit[k / w] & ((1U << rem(k)) - 1)); } inline int rank(const int val, const int k) const{ return val ? rank(k) : k - rank(k); } #undef rem }; template struct WaveletMatrixComp { int len; BitVector matrix[maxlog]; int mid[maxlog]; vector d; WaveletMatrixComp(){} WaveletMatrixComp(const vector &data) : len(data.size()), d(data){ sort(d.begin(), d.end()); d.erase(unique(d.begin(), d.end()), d.end()); vector v(len); for(int i = 0; i < len; i++) v[i] = get_pos(data[i]); vector l(len), r(len); fill(matrix, matrix + maxlog, BitVector(len + 1)); for(int level = maxlog - 1; level >= 0; level--){ int left = 0, right = 0; for(int i = 0; i < len; i++){ if(v[i] >> level & 1){ matrix[level].set(i); r[right++] = v[i]; }else{ l[left++] = v[i]; } } mid[level] = left; matrix[level].build(); v.swap(l); copy(r.begin(), r.begin() + right, v.begin() + left); } } inline int get_pos(const T x) const{ return lower_bound(d.begin(), d.end(), x) - d.begin(); } inline int succ(const int f, const int l, const int level) const{ return matrix[level].rank(f, l) + mid[level] * f; } //find v[k] T access(int k) const{ assert(0 <= k && k < len); T res = 0; for(int level = maxlog - 1; level >= 0; level--){ const bool f = matrix[level][k]; if(f) res |= T(1) << level; k = matrix[level].rank(f, k) + mid[level]*f; } return d[res]; } T operator[](const int k) const{ return access(k); } // count i (0 <= i < r) && v[i] == x int rank(const T x, int r) const{ assert(0 <= r && r <= len); const int pos = get_pos(x); if(pos >= (int)d.size() || d[pos] != x) return 0; int l = 0; for(int level = maxlog - 1; level >= 0; level--){ l = succ(pos >> level & 1, l, level); r = succ(pos >> level & 1, r, level); } return r - l; } //k-th value v[l,r) T kth_smallest(int l, int r, int k) const{ assert(0 <= l && l <= r && r <= len); assert(0 <= k && k < r - l); T res = 0; for(int level = maxlog - 1; level >= 0; level--){ int cnt = matrix[level].rank(0, r) - matrix[level].rank(0, l); bool f = cnt <= k; if(f){ res |= T(1) << level; k -= cnt; } l = succ(f, l, level); r = succ(f, r, level); } return d[res]; } T kth_largest(int l, int r, int k) const{ return kth_smallest(l, r, r-l-k-1); } // count i (l <= i < r), (v[i] < upper) int range_freq(int l, int r, const T upper) const{ const int pos = get_pos(upper); int res = 0; for(int level = maxlog - 1; level >= 0; level--){ bool f = pos >> level & 1; if(f) res += matrix[level].rank(0, r) - matrix[level].rank(0, l); l = succ(f, l, level); r = succ(f, r, level); } return res; } // count i (l <= i < r), (lower <= v[i] < upper) int range_freq(int l, int r, const T lower, const T upper) const{ assert(0 <= l && l <= r && r <= len); return range_freq(l, r, upper) - range_freq(l, r, lower); } }; template struct static_count2d{ vector> ps; WaveletMatrixComp mat; bool built = false; static_count2d(){}; void add_point(T x,T y){ ps.emplace_back(x,y); built = false; } void build(){ sort(ps.begin(),ps.end()); vector a(ps.size()); rep(i,ps.size())a[i]=ps[i].second; mat = WaveletMatrixComp(a); built = true; } int count(int xl,int xr,int yl,int yr){ assert(built); auto lt = lower_bound(ps.begin(),ps.end(),make_pair(xl,T(-INF))) - ps.begin(); auto rt = lower_bound(ps.begin(),ps.end(),make_pair(xr,T(-INF))) - ps.begin(); return mat.range_freq(lt,rt,yl,yr); } }; signed main() { cin.tie(0); ios_base::sync_with_stdio(false); cout << fixed << setprecision(20); INT(n); vi u(n),v(n); static_count2d st; rep(i,n){ INT(x,y); u[i] = x+y; v[i] = x-y; st.add_point(u[i],v[i]); } st.build(); int res = INF; rep(i,n){ int mx,mn; { int ok = 4e5+10; int ng = 0; auto f = [&](int x){ return st.count(u[i]-x,u[i]+x+1,v[i]-x,v[i]+x+1)>=n; }; while(abs(ng-ok)>1){ int mid = ng+(ok-ng)/2; if(f(mid))ok=mid; else ng = mid; } mx = ok; } { int ok = 4e5+10; int ng = 0; auto f = [&](int x){ return st.count(u[i]-x,u[i]+x+1,v[i]-x,v[i]+x+1)>=2; }; while(abs(ng-ok)>1){ int mid = ng+(ok-ng)/2; if(f(mid))ok=mid; else ng = mid; } mn = ok; } chmin(res,mx-mn); } print(res); return 0; }