bzoj 1537 bus|cdq分治

一看就知道是离散化+j排序降维+线段树（树状数组）求前缀最大值

炜神看到了这道题，吼道一定有不用数据结构的方法，我又想了一天，确实有。

然后cdq分治被我大材小用了。想到cdq应该就能yy出来了，可是快排依赖症的我归并调了好久。

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>

#define md
#define ll long long
#define inf (int) 1e9
#define eps 1e-8
#define N 100010
using namespace std;
struct P{ int x,y,w,f;} p[N],q1[N],q2[N],rank[N];
bool operator <= (P a,P b) { return a.x==b.x?a.y<=b.y:a.x<=b.x;}
bool operator < (P a,P b) { return a.y==b.y?a.x<b.x:a.y<b.y;}
bool cmpx(P a,P b) { return a.x==b.x?a.y<b.y:a.x<b.x;}
bool cmpy(P a,P b) { return a.y==b.y?a.x<b.x:a.y<b.y;}
void solve(int l,int r)
{
if (l==r) { p[l].f=max(p[l].f,p[l].w); return;}
P mid=rank[(l+r)>>1];
//------------------------------------------------------
int w1=0,w2=0;
for (int i=l;i<=r;i++)
{
if (p[i]<=mid) q1[++w1]=p[i]; else q2[++w2]=p[i];
}
int w=l-1;
for (int i=1;i<=w1;i++) p[++w]=q1[i];
for (int i=1;i<=w2;i++) p[++w]=q2[i];
solve(l,l+w1-1);
//------------------------------------------------------
w=l-1;
for (int i=1;i<=w1;i++) q1[i]=p[++w];
for (int i=1;i<=w2;i++) q2[i]=p[++w];
int a=1,b=1,mx=0;
for (int i=l;i<=r;i++)
{
if (b>w2||(a<=w1&&q1[a]<q2[b])) { p[i]=q1[a]; a++;} else { p[i]=q2[b]; b++;}
}
for (int i=l;i<=r;i++)
{
if (p[i]<=mid) mx=max(mx,p[i].f); else p[i].f=max(p[i].f,mx+p[i].w);
}
//-----------------------------------------------------
w1=0,w2=0;
for (int i=l;i<=r;i++)
{
if (p[i]<=mid) q1[++w1]=p[i]; else q2[++w2]=p[i];
}
w=l-1;
for (int i=1;i<=w1;i++) p[++w]=q1[i];
for (int i=1;i<=w2;i++) p[++w]=q2[i];
solve(l+w1,r);
//-------------------------------------------------------
w=l-1;
for (int i=1;i<=w1;i++) q1[i]=p[++w];
for (int i=1;i<=w2;i++) q2[i]=p[++w];
a=1,b=1;
for (int i=l;i<=r;i++)
{
if (b>w2||(a<=w1&&q1[a]<q2[b])) { p[i]=q1[a]; a++;} else { p[i]=q2[b]; b++;}
}
}

int main()
{
int n;
scanf("%d%d%d",&n,&n,&n);
for (int i=1;i<=n;i++)
{
scanf("%d%d%d",&p[i].x,&p[i].y,&p[i].w);
}
sort(p+1,p+n+1,cmpx);
for (int i=1;i<=n;i++) rank[i]=p[i];
sort(p+1,p+n+1,cmpy);
solve(1,n);
int ans=0;
for (int i=1;i<=n;i++) ans=max(ans,p[i].f);
printf("%d\n",ans);
return 0;
}