BZOJ4527: K-D-Sequence 线段树

别人家的题解。

#include<bits/stdc++.h>
#define N 200005
#define M (l+r>>1)
#define P (k<<1)
#define S (k<<1|1)
#define K l,r,k
#define L l,M,P
#define R M+1,r,S
#define Z \
int l=0,int r=n-1,int k=1
using namespace std;
int n,z,d;
typedef int seg[N<<2];
seg s[2],t[2],e;
void apply(int i,int d,Z){
	s[i][k]=d;
	t[i][k]=d+l;
	e[k]=d+t[i^1][k];
}
void devolve(Z){
	for(int i=0;i!=2;++i)
		if(s[i][k]!=INT_MIN){
			apply(i,s[i][k],L);
			apply(i,s[i][k],R);
			s[i][k]=INT_MIN;
		}
}
void update(int k){
	for(int i=0;i!=2;++i)
		t[i][k]=min(
		t[i][P],t[i][S]);
	e[k]=min(e[P],e[S]);
}
void change(
int i,int d,int s,int t,Z){
	if(s==l&&t==r)
		apply(i,d,K);
	else{
		devolve(K);
		if(t<=M)
			change(i,d,s,t,L);
		else if(s>M)
			change(i,d,s,t,R);
		else{
			change(i,d,s,M,L);
			change(i,d,M+1,t,R);
		}
		update(k);
	}
}
int query(int s,int t,Z){
	static int d;
	return s==l&&t==r
	&&e[k]>z?-1:l==r?l
	:(devolve(K),
	t<=M?query(s,t,L)
	:s>M?query(s,t,R)
	:~(d=query(s,M,L))
	?d:query(M+1,t,R));
}
typedef int ds[N];
ds f,a,q,u,v,p;
typedef pair<int,int> vec;
vec b;
void foo(){
	int s=*min_element(a,a+n);
	for(int i=0;i!=n;++i)
		a[i]-=s;
	set<int> t;
	for(int i=0,j=0;i!=n;++i){
		while(t.count(a[i])
		||a[i]%d!=a[j]%d)
			t.erase(a[j++]);
		f[i]=j;
		t.insert(a[i]);
	}
	for(int i=0,j=0,k=0;
	k!=n;++k,++z){
		a[u[i+1]=k]/=d;
		while(i&&a[p[i]]>a[k])
			--i;
		p[++i]=v[j+1]=k;
		change(0,-a[k],u[i],k);
		while(j&&a[q[j]]<a[k])
			--j;
		q[++j]=k;
		change(1,a[k],v[j],k);
		int s=query(f[k],k);
		b=max(b,vec(k-s+1,-s));
	}
}
void bar(){
	for(int i=0;i!=n;++i){
		int j=i;
		while(i+1!=n&&a[i]==a[i+1])
			++i;
		b=max(b,vec(i-j+1,-j));
	}
}
int main(){
	scanf("%d%d%d",&n,&z,&d);
	for(int i=0;i!=n;++i)
		scanf("%d",a+i);
	(d?foo:bar)();
	printf("%d %d\n",1-b.second,
	b.first-b.second);
}