【头疼的模板】三维偏序（陌上花开）

CDQ分治 + 树状数组 的模板题

这道题是zyd在济南三期qbxt讲bitset的时候提到的，不过我好像当时秒了正解的样子。。。

\(===========================分割线================================\)

什么是CDQ分治：

其实cdq分治的思想与应用都能被很简单地描述——它是用来解决各种区间段转移问题\(x -> y (x < y)\)的一种算法。
我们用f[x]表示位置x转移之后的结果，用solve(l,r)来传递完全限制在[l,r]范围内的状态转移，并且转移\(a->b\)一定有\(a<b\)
那么对于solve(l,r)

1.可以把这些问题排成一个序列，用一个区间[L,R]表示。
2.分。递归处理左边区间[L,M]和右边区间[M+1,R]的问题。
3.治。合并两个子问题，同时考虑到[L,M]内的修改对[M+1,R]内的查询产生的影响。即，用左边的子问题帮助解决右边的子问题

这样CDQ分治就做完啦

那CDQ有什么好处呢？

CDQ分治是我们处理各类问题的重要武器
它的优势在于可以顶替复杂的高级数据结构，而且常数比较小；
缺点在于必须离线操作。

然后我们考虑下这道题的解法：

   我们考虑一下，这道题的要求是求三维偏序，但是好像CDQ分治只能求二维偏序的样子qwq。。。
   所以我们先在CDQ分治之前，把其中1维变成有序的，然后再分治下去，这时在区间\([l,mid]\)关于区间\([mid + 1,r]\)就不存在某一维的逆序了，所以就只剩下2维，于是就将三维偏序成功的转化为可以用CDQ分治来解的二维偏序问题。
   这个时候来一波树状数组求逆序对的操作搞一下二维偏序，就可以把跨过中线的，左边更新右边的情况计算出来,然后左右两边递归处理就好了。。。

注意：只计算左边的操作对右边的询问的贡献！！！

时间复杂度 : 分治 + 树状数组 = \(O(nlog^2 n)\)

CODE:

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

using namespace std;
const int N = 200100;

struct QAQ {
    int x,y;
    int z,p;
    inline bool operator == (const QAQ &b) const {
        return x == b.x && y == b.y && z == b.z;
    }
    inline bool operator < (const QAQ &b) const {
        return y < b.y || (y == b.y && z < b.z);
    }
}a[N],c[N];
int m,n,ans[N],opt[N];

inline int lowbit(int x){
    return (x & (-x));
}
inline bool cmp(QAQ a,QAQ b) {
    return a.x < b.x || (a.x == b.x && (a.y < b.y || (a.y == b.y && a.z < b.z)));
}
void print(int x,int y) {
    for( ; x <= m ; x += lowbit(x)) {
        opt[x] += y;
    //  x += lowbit(x);
    }
}
int gets(int x) {
    int sum = 0;
    for(; x ; x-= lowbit(x)) {
        sum += opt[x];
        //x -= lowbit(x);
    }
    return sum;
}
void clean(int x) {
    for(;x <= m;x += lowbit(x)) {
        opt[x] = 0;
    }
}
void CDQ(int l,int r) {
    if(l >= r) return ;
    int mid = (l + r) >> 1;
    CDQ(l,mid);
    CDQ(mid+1,r);
    int i = l; 
    for(int j = mid+1 ; j <= r ; j++) {
        while(i <= mid && a[i].y <= a[j].y) {
            print(a[i].z,1);
            i++;
        }
        a[j].p += gets(a[j].z);
    }
    for(int i = l ; i <= mid ; i++) clean(a[i].z);
    for(int i = l ; i <= r; i++) c[i] = a[i];
    i = l; int j = mid+1;
    for(int k = l ; k <= r ; k++) {
        if(j > r || ((i <= mid && c[i] < c[j]))) {
            a[k] = c[i++];
        }else a[k] = c[j++];
    }
}

int main() {
    scanf("%d%d",&n,&m);
    for(int i = 1 ;i <= n ; i++) {
        scanf("%d%d%d",&a[i].x,&a[i].y,&a[i].z);
    }
    sort(a+1,a+n+1,cmp);
    CDQ(1,n);
    sort(a+1,a+n+1,cmp);
    int i = 1;
    while(i <= n) {
        int j = i ,Max = 0;
        while(j <= n && a[i] == a[j])
            Max = max(Max,a[j++].p);
            ans[Max] += j-i;
            i = j;
    }
    for(int i = 0 ; i < n ; i++)
        printf("%d\n",ans[i]);
    return 0;
}

转载于:https://www.cnblogs.com/Repulser/p/9630278.html