首先我们枚举横着一刀切在哪,然后考虑此时竖着一刀切在哪里最好。我们令横线上方的点记为0,下方的点记为1,竖线左边的点记为l,右边的点记为r,则随着竖线的右移,显然有
max(l0,l1)
m
a
x
(
l
0
,
l
1
)
逐渐变大,
max(r0,r1)
m
a
x
(
r
0
,
r
1
)
逐渐变小。
记
ans=max(max(l0,l1),max(r0,r1))
a
n
s
=
m
a
x
(
m
a
x
(
l
0
,
l
1
)
,
m
a
x
(
r
0
,
r
1
)
)
,则我们一定有ans先减小后增大,是个单峰函数,我们可以在线段树上二分一下找到这个峰值点,更新答案。
复杂度 O(nlogn) O ( n l o g n )
upd:这个单峰函数的最值可能在两端,需要取min。但是大概因为枚举了横线,就算此次求出的竖线不是当前横线最优的,那最优的也一定会在枚举某个横线时出现。
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 100010
inline char gc(){
static char buf[1<<16],*S,*T;
if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
return *S++;
}
inline int read(){
int x=0,f=1;char ch=gc();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
return x*f;
}
int n,num[N*10],ans=inf,X0,X1;
struct node{
int x0,x1;
}tr[N*10*4];
struct P{
int x,y;
friend bool operator<(P a,P b){return a.y<b.y;}
}a[N];
inline void pushup(int p){
tr[p].x0=tr[p<<1].x0+tr[p<<1|1].x0;
tr[p].x1=tr[p<<1].x1+tr[p<<1|1].x1;
}
inline void build(int p,int l,int r){
if(l==r){tr[p].x0=num[l];tr[p].x1=0;return;}
int mid=l+r>>1;build(p<<1,l,mid);build(p<<1|1,mid+1,r);pushup(p);
}
inline void change(int p,int l,int r,int x){
if(l==r){tr[p].x0--;tr[p].x1++;return;}
int mid=l+r>>1;
if(x<=mid) change(p<<1,l,mid,x);
else change(p<<1|1,mid+1,r,x);pushup(p);
}
inline int ask(int p,int l,int r,int x0,int x1){//找峰值
if(l==r) return l;
int mid=l+r>>1;int sum0=x0+tr[p<<1].x0,sum1=x1+tr[p<<1].x1;
if(max(sum0,sum1)>=max(tr[1].x0-sum0,tr[1].x1-sum1))
return ask(p<<1,l,mid,x0,x1);
return ask(p<<1|1,mid+1,r,sum0,sum1);
}
inline void qsum(int p,int l,int r,int x,int y){
if(x>y) return;
if(x<=l&&r<=y){X0+=tr[p].x0;X1+=tr[p].x1;return;}
int mid=l+r>>1;
if(x<=mid) qsum(p<<1,l,mid,x,y);
if(y>mid) qsum(p<<1|1,mid+1,r,x,y);
}
int main(){
// freopen("a.in","r",stdin);
n=read();
for(int i=1;i<=n;++i) a[i].x=read(),a[i].y=read(),num[a[i].x]++;
sort(a+1,a+n+1);build(1,1,1e6);
for(int i=1,j;i<=n;i=j){//枚举横向切在a[i].y上面
j=i;
while(j<=n&&a[j].y==a[i].y) change(1,1,1e6,a[j].x),++j;
int pos=ask(1,1,1e6,0,0);
X0=0;X1=0;qsum(1,1,1e6,1,pos);
ans=min(ans,max(X0,X1));
X0=0;X1=0;qsum(1,1,1e6,1,pos-1);
ans=min(ans,max(tr[1].x0-X0,tr[1].x1-X1));
}printf("%d\n",ans);
return 0;
}