CDQ分治

最新推荐文章于 2024-02-26 17:43:38 发布

Vicente.

最新推荐文章于 2024-02-26 17:43:38 发布

阅读量211

点赞数

本文链接：https://blog.csdn.net/qq_42447015/article/details/99992042

版权

CDQ分治

多维偏序问题

二维偏序问题

HDU 1541 Stars

给定 $n$ 个元素，第 $i$ 个元素有 $a_i 、 b_i$ 两个属性，设 $f (i)$ 表示满足 $a_j\leq a_i$ 且 $b_j\leq b_i$ 的 j 的数量。

对于 $d\in [0,n]$ ，求满足 $f (i) = d$ 的数量。

首先，我们可以把两个元素抽象成一个点 $(a, b)$ ，那么我们就是求一个矩形中有多少个点。 (1<=N<=15000, 0<=X,Y<=32000)

比如我们要求这个矩形内有多少个点：

方法一：按照x坐标排序，y坐标离散一下，从左向右跑一下树状数组插入，在第i个位置找一下在1~i-1位置上纵坐标小于i的纵坐标的点的个数。
用二维偏序CDQ解决树状数组问题

P3374 【模板】树状数组 1
在这篇博客里有很详细的讲解，说的很清晰。

通过这个题第一次理解了CDQ的查询的左右端点分开存储使用的神奇操作。

#include<iostream>
#include<cstdio>
#include<cstring>
#define MogeKo qwq
using namespace std;
const int maxn = 500005*3;

int n,m,cnt,cqry,opt,x,y,ans[maxn];
struct node{
    int type,id,val;
    bool operator < (const node & x) const{
        if(id != x.id)return id < x.id;
        else return type < x.type;
    }
}q[maxn],tem[maxn];

void cdq(int L,int R){
    if(L == R) return;
    int mid = L+R>>1;
    cdq(L,mid),cdq(mid+1,R);
    int t1 = L,t2 = mid+1;
    int sum = 0;
    for(int i = L;i <= R;i++){
        if( (t1 <= mid && q[t1]<q[t2]) || t2 > R){
            if(q[t1].type == 1) sum += q[t1].val;
            tem[i] = q[t1++];
        }
        else{
            if(q[t2].type == 2) ans[q[t2].val] -= sum;
            if(q[t2].type == 3) ans[q[t2].val] += sum;
            tem[i] = q[t2++];
        }
    }
    for(int i = L;i <= R;i++) q[i] = tem[i];
}

int main(){
    scanf("%d%d",&n,&m);
    for(int i = 1;i <= n;i++){
        cnt++;
        scanf("%d",&y);
        q[cnt].type = 1;
        q[cnt].id = i;
        q[cnt].val = y;
    }
    for(int i = 1;i <= m;i++){
        scanf("%d%d%d",&opt,&x,&y);
        if(opt == 1){
            q[++cnt].type = 1;
            q[cnt].id = x;
            q[cnt].val = y;
        }
        if(opt == 2){
            cqry++;
            q[++cnt].type = 2;
            q[cnt].id = x-1;
            q[cnt].val = cqry;
            q[++cnt].type = 3;
            q[cnt].id = y;
            q[cnt].val = cqry;
        }
    }
    cdq(1,cnt);
    for(int i = 1;i <= cqry;i++)
        printf("%d\n",ans[i]);
    return 0;
}