BZOJ_P3262 陌上花开(CDQ分治+树状数组)

BZOJ传送门

Time Limit: 20 Sec Memory Limit: 256 MB
Submit: 1045 Solved: 475
[Submit][Status][Discuss]
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的每级花的数量。

Sample Input

10 3
3 3 3
2 3 3
2 3 1
3 1 1
3 1 2
1 3 1
1 1 2
1 2 2
1 3 2
1 2 1

Sample Output

3
1
3
0
1
0
1
0
0
1

HINT

1 <= N <= 100,000, 1 <= K <= 200,000

Source

树套树 CDQ分治

Sol:
CDQ分治裸题，三维偏序，CDQ分治s，排序c，树状数组搞m
这道题我调了半天QuQ，处理相同的情况比较麻烦，先对s排序，统计一下三个全相同的，对s重新编号，然后按c排序，分治s，用树状数组维护一下，就跟上一道题很像了，s就是id！

#include<cstdio>
#include<algorithm>
using namespace std;
#define N 100005
#define K 200005
inline int in(int x=0,char ch=getchar()){while(ch>'9'||ch<'0') ch=getchar();
    while(ch>='0'&&ch<='9') x=(x<<3)+(x<<1)+ch-'0',ch=getchar();return x;}
struct Node{int s,c,m,w,v;}d[N],tmp[N];
bool operator == (const Node &a,const Node &b){return a.s==b.s&&a.c==b.c&&a.m==b.m;}
inline int cmps(const Node &a,const Node &b){
    if(a.s==b.s&&a.c==b.c) return a.m<b.m;if(a.s==b.s) return a.c<b.c;return a.s<b.s;}
inline int cmpc(const Node &a,const Node &b){
    if(a.m==b.m&&a.c==b.c) return a.s<b.s;if(a.c==b.c) return a.m<b.m;return a.c<b.c;}
int n,k,deep,tot;int ans[N],c[K],g[K];
inline void Add(int x,int v){for(;x<=k;x+=x&-x){if(g[x]!=deep) g[x]=deep,c[x]=0;c[x]+=v;}}
inline int Sum(int x,int res=0){for(;x;x-=x&-x)if(g[x]==deep) res+=c[x];return res;}
void Solve(int l,int r){
    if(l==r){d[l].v+=d[l].w-1;return;};int mid=(l+r)>>1,h=l,t=mid+1;
    for(int i=l;i<=r;i++) if(d[i].s<=mid) tmp[h++]=d[i];else tmp[t++]=d[i];
    for(int i=l;i<=r;i++) d[i]=tmp[i];
    Solve(l,mid);deep++;
    for(int i=mid+1,j=l;i<=r;i++){
        for(;j<=mid&&d[j].c<=d[i].c;j++) Add(d[j].m,d[j].w);
        d[i].v+=Sum(d[i].m);
    }
    Solve(mid+1,r);h=l,t=mid+1;
    for(int i=l;i<=r;i++) if(h<=mid&&(t>r||cmpc(d[h],d[t]))) tmp[i]=d[h++];else tmp[i]=d[t++];
    for(int i=l;i<=r;i++) d[i]=tmp[i];
}
int main(){
    n=in(),k=in();
    for(int i=1;i<=n;i++) d[i].s=in(),d[i].c=in(),d[i].m=in(),d[i].w=1;
    sort(d+1,d+n+1,cmps);
    for(int i=1;i<=n;i++) if(i!=1&&d[i]==d[i-1]) d[tot].w++;else d[++tot]=d[i];
    for(int i=1;i<=tot;i++) d[i].s=i;
    sort(d+1,d+tot+1,cmpc);Solve(1,tot);
//  for(int i=1;i<=tot;i++) printf("%d%c",d[i].v," \n"[i==tot]);
//  for(int i=1;i<=tot;i++) printf("%d%c",d[i].w," \n"[i==tot]);
//  putchar('\n');
    for(int i=1;i<=tot;i++) ans[d[i].v]+=d[i].w;
    for(int i=0;i<n;i++) printf("%d\n",ans[i]);
    return 0;
}