不知道大佬们怎么想的,反正我看到这种区间包含性质的并且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 */