CCPC－Wannafly & Comet OJ 夏季欢乐赛（2019）D

题面

 

    一开始想到一个 O(N^2) 作法，先把x排序，而后顺次枚举x最大的点，看向前最多能够保留多少点 (也就是先无论正方形的上下长度限制，先考虑左右的限制)。而后再对这些点作一遍相似的。。(等等这么一说感受怎么成 O(N^2 logN)的了23333)。

    但反正这样确定是不行的，数据太大了，咱们要作一些优化。

    咱们保留限制左右边界的这个过程，看能不能优化一下限制y坐标找答案的过程。
    考虑一颗线段树，从左向右第i个叶子节点表明若是上边界在纵坐标i，此时的答案是多少 (注意，右边界是在枚举的过程当中不断右移的，因此这颗线段树是不断变化的)。

    显然，加入纵坐标为y的点以后，第 y 到 第 y+k 的叶子的答案都会 +1；删点反之。

    因而咱们把x坐标排序后单调扫一遍，边扫边维护线段树最大值&更新答案就能够了。

    不过如今还有一点小问题:线段树的叶子节点有多大 2*10^9多个。

    直接写的话由于涉及的边界叶子很少，因此动态开点线段树应该是能够过的；但咱们想要时空复杂度更低的算法，那么考虑一下是否能够离散化呢？

    像这种叶子过多的线段树想要离散化的话，必须能够证实出其正确性：好比这个题，咱们称一个没有被任何y[i]与y[i]+k覆盖到的纵坐标为空点，其他为非空点。若是离散化的话，咱们就默认最后的答案必定是在非空点处取到。事实也的确是这样的，考虑一个空点，若是它要被区间+1覆盖到，它两边最近的非空点确定也会被覆盖到；区间-1的话可能当作以前那个区间+1被撤销了，因此不用考虑。

 

    恩，这样就能够很舒服的AC，代码也很好写。

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int N=100005;
#define lc (o<<1)
#define mid (l+r>>1)
#define rc ((o<<1)|1)

inline int read(){
    int x=0; char ch=getchar();
    for(;!isdigit(ch);ch=getchar());
    for(;isdigit(ch);ch=getchar()) x=x*10+ch-'0';
    return x;
}

struct node{
	int x,y,z;
	bool operator <(const node &u)const{
		return x<u.x;
	}
}a[N];

int n,k,num[N*2],ky,ans,w;
int mx[N*8+1],tag[N*8+1],le,ri;

inline void Get(int o,int ad){ tag[o]+=ad,mx[o]+=ad;}

inline void pushdown(int o){ Get(lc,tag[o]),Get(rc,tag[o]),tag[o]=0;}

void update(int o,int l,int r){
	if(l>=le&&r<=ri){ Get(o,w); return;}
	pushdown(o);
	
	if(le<=mid) update(lc,l,mid);
	if(ri>mid) update(rc,mid+1,r);
	
	mx[o]=max(mx[lc],mx[rc]); 
}

inline void solve(){
	sort(a+1,a+n+1);
	
	for(int i=1;i<=n;i++) num[++ky]=a[i].y,num[++ky]=a[i].z;
	
	sort(num+1,num+ky+1),ky=unique(num+1,num+ky+1)-num-1;
	
	for(int i=1;i<=n;i++){
		a[i].y=lower_bound(num+1,num+ky+1,a[i].y)-num;
		a[i].z=lower_bound(num+1,num+ky+1,a[i].z)-num;		
	}
	
	for(int i=1,l=1;i<=n;ans=max(ans,mx[1]),i++){
		while(a[i].x-a[l].x>k) le=a[l].y,ri=a[l].z,w=-1,update(1,1,ky),l++;
		le=a[i].y,ri=a[i].z,w=1,update(1,1,ky);
	} 
}

int main(){
	n=read(),k=read();
	for(int i=1;i<=n;i++)
	    a[i].x=read(),a[i].y=read(),a[i].z=a[i].y+k;
	    
	solve(),printf("%d\n",ans);
	return 0;
}