洛谷 P3810 【模板】三维偏序（陌上花开）（CDQ分治）

传送门

---

真是一道毒瘤题（弄了接近一下午+一晚上）

解题思路

先对所有的点按照a-b-c的优先值从小到大排序，然后去重（注意要记录每一个有几个重复的，因为条件是小于等于，所以重复的对答案也有影响），然后再按照关键字b归并排序，排序过程中，用树状数组记录c值，每一次合并时如果是右半部分的b值较小，答案就加上先前的所有c值小于这个c值的个数。

由于用memset时间复杂度过大，所以进行完操作后用for循环对所有更改过的点进行修改。

然后因为答案问的不是总数，所以用一系列数组来维护各种值，及其恶心。

然后注意k值实际上就是树状数组的边界，千万不能用n，否则会wa三个点。

然后才能顺利AC。

AC代码

#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstdio>
#include<cstring>
using namespace std;
const int maxn=100005;
int n,k,ans[maxn],anss[maxn],numm[maxn];
struct node{
    int id,a,b,c,num;
}p[maxn],q[maxn];
bool cmp(const node &x,const node &y){
    if(x.a!=y.a){
        return x.a<y.a;
    }
    if(x.b!=y.b){
        return x.b<y.b;
    }
    return x.c<y.c;
}
inline int lowbit(int x){
    return x&-x;
}
int c[2*maxn];
void update(int x,int v){
    for(int i=x;i<=k;i+=lowbit(i)){
        c[i]+=v;
    }
}
int query(int x){
    int res=0;
    for(int i=x;i>0;i-=lowbit(i)){
        res+=c[i];
    }
    return res;
}
void divide(int l,int r){
    if(l==r) return;
    int mid=(l+r)>>1;
    int p1=l,p2=mid+1,tot=0;
    divide(l,mid);
    divide(mid+1,r);
    while(p1<=mid||p2<=r){
        if(p1<=mid&&(p2>r||p[p1].b<=p[p2].b)){
            q[++tot]=p[p1++];
            update(q[tot].c,q[tot].num);
        }else{
            q[++tot]=p[p2++];
            ans[q[tot].id]+=query(q[tot].c);
        }
    }
    for(int i=l;i<=mid;i++){
        update(p[i].c,-p[i].num);
    }
    for(int i=1;i<=tot;i++){
        p[l+i-1]=q[i];
    }
}
int main()
{
    cin>>n>>k;
    for(int i=1;i<=n;i++){
        scanf("%d%d%d",&p[i].a,&p[i].b,&p[i].c);
    }
    sort(p+1,p+n+1,cmp);
    int tot=0;
    for(int i=1;i<=n;){
        int j=i,cnt=0;
        while(j<=n&&p[i].a==p[j].a&&p[i].b==p[j].b&&p[i].c==p[j].c){
            p[i].num++;
            cnt++;
            j++;
        }
        q[++tot]=p[i];
        numm[tot]=cnt;
        i=j;
    }
    int nn=n;
    n=tot;
    memset(p,0,sizeof(p));
    for(int i=1;i<=tot;i++) p[i]=q[i],p[i].id=i;
    divide(1,n);
    for(int i=1;i<=n;i++){
        ans[i]+=numm[i]-1;
    }
    for(int i=1;i<=n;i++){
        anss[ans[i]]+=numm[i];
    }
    for(int i=0;i<nn;i++) printf("%d
",anss[i]);
    return 0;
}

 寒假作业已经放弃了qwq