bzoj 1537 bus|cdq分治

CDQ分治优化算法

最新推荐文章于 2022-10-05 14:28:42 发布

原创最新推荐文章于 2022-10-05 14:28:42 发布 · 463 阅读

0 ·

CC 4.0 BY-SA版权

乱搞专栏收录该内容

5 篇文章

订阅专栏

本文介绍了一种使用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;
}