Hdu 5618 CDQ分治

最新推荐文章于 2020-06-13 16:10:58 发布

MrBird_to_fly

最新推荐文章于 2020-06-13 16:10:58 发布

阅读量461

点赞数

CC 4.0 BY-SA版权

分类专栏：离线乱搞

本文链接：https://blog.youkuaiyun.com/MrBird_to_fly/article/details/70148838

离线乱搞专栏收录该内容

3 篇文章

订阅专栏

本文通过CDQ分治策略解决了一个三维偏序问题，并详细展示了使用线段树进行区间更新和查询的过程。作者分享了从调试到最终实现的曲折经历，强调了面对复杂错误时的调试技巧。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

上次写填坑帖顺手学了一下markdown然后顺手用了一下发现无序列表一直跪然后就很不爽直接一个月没写blog。。。

然后这个题就是个傻逼三维偏序啊。。。但是我怎么写OLE了啊。。。菜到随便改改都是bug啊。。。
以后这种奇怪的错误。。。直接认为WA。。。可能会好一点啊。。。

CDQ就TM是分治一下。。。然后考虑两部分之间的影响。。。但是我TM现在才会啊。。。我好菜菜啊。。。

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;

int n;

int ans[105000];

struct node{
    int x, y, z;
    int id;
}a[105000], b[105000];

bool cmp(node a, node b){
    if(a.x!=b.x)return a.x<b.x;
    if(a.y!=b.y)return a.y<b.y;
    return a.z<b.z;
}

bool cmp1(node a, node b){
    if(a.y!=b.y)return a.y<b.y;
    return a.z<b.z;
}

int c[105000];

int lowbit(int x){
    return x&(-x);
}

void update(int p){
    while(p<=100000){
        c[p]++;
        p+=lowbit(p);
    }
}

int query(int p){
    int sum=0;
    while(p){
        sum+=c[p];
        p-=lowbit(p);
    }
    return sum;
}

void reset(int p){
    while(p<=100000){
        c[p]=0;
        p+=lowbit(p);
    }
}

void cdq(int l, int r){
    if(l==r)return;
    int m=(l+r)/2;
    cdq(l, m);
    cdq(m+1, r);
    for(int i=l;i<=r;i++){
        b[i]=a[i];
    }
    sort(b+l, b+m+1, cmp1);
    sort(b+m+1, b+r+1, cmp1);
    int j=l;
    for(int i=m+1;i<=r;i++){
        while(j<=m&&b[j].y<=b[i].y){
            update(b[j++].z);
        }
        ans[b[i].id]+=query(b[i].z);
    }
    for(int i=l;i<=m;i++)reset(b[i].z);
    return;
}

int main()
{
    int T;
    scanf("%d", &T);
    while(T--){
        scanf("%d", &n);
        for(int i=1;i<=n;i++){
            scanf("%d%d%d", &a[i].x, &a[i].y, &a[i].z);
            a[i].id=i;
        }
        sort(a+1, a+1+n, cmp);
        memset(c, 0, sizeof(c));
        memset(ans, 0, sizeof(ans));
        cdq(1, n);
        for(int i=n-1;i>=1;i--){
            if(a[i].x==a[i+1].x&&a[i].y==a[i+1].y&&a[i].z==a[i+1].z)ans[a[i].id]=ans[a[i+1].id];
        }
        for(int i=1;i<=n;i++){
            printf("%d\n", ans[i]);
        }
    }
    return 0;
}