D. Ball-(CDQ分治)

D. Ball

题目大意:

n n n ( a , b , c ) (a,b,c) (a,b,c),类型的对。问有多少个对存在严格大于关系。

思路:

多维偏序先排序降维。
C D Q CDQ CDQ时习惯用左区间更新右区间,故排序按降序排序。
第二维采用归并排序的方法。
题目的特殊地方在于:对单个对只用判断是否存在不需要记录个数
第三维的树状数组可以省去。

思考分治时左区间如何更新右区间。
分治时对二维排序的同时可以对第三维取最值,当取到右区间的值时直接和最值比较即可。

但会有一个很蛋疼的问题,一维不能保证严格大于,只是大于等于。
上述分治策略在两区间存在一维相同情况下,最值会偏大。

在分治时合并区间后一维降序不能保证,但区间与区间之间保证降序。
故右区间一维只可能和左区间最小值相同。
记录两个最值一个包含左区间最小,一个不包含即可。

具体看代码。

Code:

#include <iostream>
#include <map>
#include <set>
#include <queue>
#include <stack>
#include <algorithm>
#include <vector>
#include <string>
#include <iomanip>
#include <cmath>
#include <ctime>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <climits>
//#include <unordered_map>
#define guo312 std::ios::sync_with_stdio(false), cin.tie(0), cout.tie(0)
#define ll long long
#define Inf LONG_LONG_MAX
#define inf INT_MAX
#define endl "\n"
#define PI 3.1415926535898
using namespace std;
const int N=5e5+10;

int n;
struct PW{
	int ai,bi,ci,ans;	
}a[N];

bool cmp_ai(PW &s1,PW &s2){
	return s1.ai>s2.ai;
}

bool cmp_bi(PW &s1,PW &s2){
	return s1.bi>s2.bi;
}

int lg[N];
void init(){
	lg[1] = 0, lg[2] = 1;  
    for(int i = 3 ; i < N ; i ++)
    lg[i] = lg[i >> 1] + 1; 
}

int minn[N][20];
void work(){  // ST表
	init();
	for(int j = 0 ; j <= 19 ; j ++){
        for(int i = 1 ; i + (1 << j) - 1 <= n; i ++){
            if(!j) minn[i][j] = a[i].ai;
            else minn[i][j] = min(minn[i][j - 1],minn[i + (1 << j - 1)][j - 1]);
        }
    } 
}

int query(int l ,int r){
    int k = lg[r - l + 1];
    return min(minn[l][k], minn[r + 1 - (1 << k)][k]);
}

void CDQ(int l,int r){
	if(l==r) return;
	int mid=l+r>>1;
	CDQ(l,mid); CDQ(mid+1,r);
	int minn=query(l,mid),maxn1=0,maxn2=0;
	for(int i=mid+1,j=l;i<=r;i++){
		while(j<=mid&&a[j].bi>a[i].bi){  // 保证二维严格大于
			maxn1=max(maxn1,a[j].ci);
			if(a[j].ai!=minn){
				maxn2=max(maxn2,a[j].ci);
			}
			j++;
		}
		if(a[i].ai==minn){
			if(maxn2>a[i].ci){
				a[i].ans=1;
			}
		}
		else if(maxn1>a[i].ci){
			a[i].ans=1;
		}
	}
	sort(a+l,a+r+1,cmp_bi);
}

int main(){
guo312;
	cin>>n;
	for(int i=1;i<=n;i++){
		cin>>a[i].ai;
	}
	for(int i=1;i<=n;i++){
		cin>>a[i].bi;
	}
	for(int i=1;i<=n;i++){
		cin>>a[i].ci;
	}
	sort(a+1,a+1+n,cmp_ai);
	work(); CDQ(1,n); int ans=0;
	for(int i=1;i<=n;i++){
		if(a[i].ans) ans++;
	}
	cout<<ans;
	return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

要用bug来打败bug

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值