BZOJ 3262: 陌上花開 [CDQ分治 三維偏序]

Description

有n朵花，每朵花有三個屬性：花形(s)、顏色(c)、氣味(m)，又三個整數表示。現要對每朵花評級，一朵花的級別是它擁有的美麗能超過的花的數量。定義一朵花A比另一朵花B要美麗，當且僅當Sa>=Sb,Ca>=Cb,Ma>=Mb。顯然，兩朵花可能有同樣的屬性。需要統計出評出每個等級的花的數量。

Input

第一行為N,K (1 <= N <= 100,000, 1 <= K <= 200,000 ), 分別表示花的數量和最大屬性值。

以下N行，每行三個整數si, ci, mi (1 <= si, ci, mi <= K)，表示第i朵花的屬性

Output

包含N行，分別表示評級為0...N-1的每級花的數量。

---

太坑了！！！！！！！！！

各種修改比較操作各種修改細節最后發現是因為我去重之后把$n$改變了然而人家輸出需要原來的$n$!!!!!!!!!!!!!!!!!

 

裸的三維偏序

$a$排序，$b$ $CDQ$分治，$c$樹狀數組

一開始$a$排序需要考慮$a,b,c$整體排序，然后分治里$b$排序要用$\le$，或者$b,c$整體$\le$，注意這里就不關$a$什么事了，分治已經保證用左面更新右面了

注意需要去重，一個元素代表多個

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
using namespace std;
typedef long long ll;
const int N=2e5+5;
inline int read(){
    char c=getchar();int x=0,f=1;
    while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
    while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
    return x*f;
}
int n,maxVal;
struct Operation{
    int a,b,c,w;
    int f;//ans
    bool operator <(const Operation &r)const{
        //return a<r.a || (a==r.a&&b<r.b) || (a==r.a&&b==r.b&&c<r.c);
        return (a==r.a&&b==r.b) ? c<r.c : (a==r.a?b<r.b:a<r.a);
    }
}a[N],t[N];
int c[N];
inline int lowbit(int x){return x&-x;}
inline void add(int p,int v){for(;p<=maxVal;p+=lowbit(p)) c[p]+=v;}
inline int sum(int p){
    int re=0;
    for(;p;p-=lowbit(p)) re+=c[p];
    return re;
}
int ans[N];
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,j=mid+1,p=l;
    while(i<=mid||j<=r){
        if(j>r||(i<=mid&&a[i].b<=a[j].b)) add(a[i].c,a[i].w),t[p++]=a[i++];
        else a[j].f+=sum(a[j].c),t[p++]=a[j++];
    }
    for(int i=l;i<=mid;i++) add(a[i].c,-a[i].w);
    for(int i=l;i<=r;i++) a[i]=t[i];
}
int main(){
    freopen("in","r",stdin);
    n=read();maxVal=read();
    for(int i=1;i<=n;i++) 
        a[i].a=read(),a[i].b=read(),a[i].c=read(),a[i].w=1;
    sort(a+1,a+1+n);
    
    int p=1;
    for(int i=2;i<=n;i++){
        if(a[i].a==a[p].a&&a[i].b==a[p].b&&a[i].c==a[p].c) a[p].w++;
        else a[++p]=a[i];
    }
    int nnnnn=n;
    n=p;
    
    CDQ(1,n);
    for(int i=1;i<=n;i++) ans[a[i].f+a[i].w-1]+=a[i].w;
    for(int i=0;i<=nnnnn-1;i++) printf("%d\n",ans[i]);
}