Codeforces 1221 F Choose a Square

题面

不知道大佬们怎么想的，反正我看到这种区间包含性质的并且score只和包含的区间与询问区间挂钩的题，马上就想到了扫描线23333

虽然革命方向无比正确，但却因为SB错误交了5次才 A。

WA第一发：考虑到我们扫描线的话是要坐标离散化的，然而这种做法是默认正方形内至少有一个点。可所有点权都是负值就尴尬了，这个时候需要找一个空点来把答案替代成0。虽然我考虑到了这种情况，但是第一次找空点的时候不仅写麻烦了(我是用 num[i]+1<num[i+1] 来判断两个坐标之间是否有空隙),而且漏了最大坐标后面的。。。。不过如果你观察题面描述仔细的话，你是会发现输出要求的 坐标是可以 >1e9的，这就是这题的坑比之处，答案<0就可以直接输出 0 \n 1000000001 1000000001 1000000001 1000000001 了。。

WA第二发：我开的vector的下标是离散化之后的坐标，显然是要开两倍空间的，忘开了。。。

WA第三发：又发现线段树的下标也是离散化之后的坐标，显然空间是要开n的八倍的，wori

WA第四发：操作线段树的时候有个 long long 写成int了，gan

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

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

/*
    注意数组空间啊！
	虽然只有n个点，不过离散化之后最多涉及 2n 个坐标值
	g[i]是用来存右端点在i的点的链表，所以空间开2n
	因为线段树的下标是离散后的坐标次序，所以空间要开 8n 而不是 4n
*/

vector<int> g[N*2];
int n,a[N],b[N],num[N*2],ky,c[N],L,R,pos[N*8],w,le=1,ri,P;
ll ans=-1e18,mx[N*8],ad[N*8],now;

inline void maintain(int o){
	pos[o]=mx[rc]>mx[lc]?pos[rc]:pos[lc];
	mx[o]=max(mx[lc],mx[rc]);
}

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

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

void build(int o,int l,int r){
	if(l==r){ pos[o]=mx[o]=num[l]; return;}
	build(lc,l,mid),build(rc,mid+1,r);
	maintain(o);
}

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);

	maintain(o);
}

void query(int o,int l,int r){
	if(l>=le&&r<=ri){ if(mx[o]>mx[P]) P=o; return;}
	// P在query后是最大的线段树上的节点，需要转换到具体的数轴上的点
	pushdown(o);

	if(le<=mid) query(lc,l,mid);
	if(ri>mid) query(rc,mid+1,r);
}

inline void solve(){
	build(1,1,ky),mx[0]=-1e18;
	for(int i=1;i<=ky;i++){
		for(int j:g[i]) ri=a[j],w=c[j],update(1,1,ky);
		P=0,ri=i,query(1,1,ky);
		if(mx[P]-num[i]>ans) ans=mx[P]-num[i],R=num[i],L=pos[P];
	}
}

int main(){
	n=read(),num[0]=-1;
	for(int i=1;i<=n;i++){
	    a[i]=read(),b[i]=read(),c[i]=read();
		num[++ky]=a[i],num[++ky]=b[i];
		if(a[i]>b[i]) swap(a[i],b[i]);
	}

	sort(num+1,num+ky+1),ky=unique(num+1,num+ky+1)-num-1;
	for(int i=1;i<=n;i++){
		a[i]=lower_bound(num+1,num+ky+1,a[i])-num;
		b[i]=lower_bound(num+1,num+ky+1,b[i])-num;
		g[b[i]].pb(i);
	}

	solve();
	if(ans<=0) ans=0,L=R=1e9+1;

	cout<<ans<<endl;
	printf("%d %d %d %d\n",L,L,R,R);

	return 0;
}

/*
    Σc[] + j - i
*/