【BZOJ3262】陌上花开（树套树）

题面

对于权限题，我这种苦逼肯定是从别的OJ上搞的对不对？？？
CJOJ
洛谷

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

题解

然而我太菜了，不会CDQ分治

智商不够数据结构来补 ——– 菊开
所以，当然就是数据结构来补智商

很显然的，傻逼都想得到的第一维当然是sort呀
第二维，显然就是要找前面的比他小的个数
那要维护区间，要搞值域，当然是树状数组或者线段树了
第三维要找这段区间内比他小的点的个数，当然是平衡数了

于是，平衡树只会Splay的我就打了一个树状数组套Splay,慢的飞起
我还是太菜啦。。。。
orz球爷，成为了继就老爷后，第二个替罪羊树dalao
这题我改了好久。。。
原因是Splay在Insert之后忘记Spaly
然后就没有pushup，所以size挂了。。。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
#define MAX 100010
inline int read()
{
    register int x=0,t=1;
    register char ch=getchar();
    while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
    if(ch=='-'){t=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-48;ch=getchar();}
    return x*t;
}
int ans[MAX];
struct seq
{
    int a,b,c;
}s[MAX];
int K;
inline int lowbit(int x){return x&(-x);}
struct Node
{
    int ch[2];
    int v,f,size,cnt;
}t[MAX<<5];
int root[MAX<<1];
int n,tot;
inline void pushup(int x)
{
    t[x].size=t[t[x].ch[0]].size+t[t[x].ch[1]].size+t[x].cnt;
}
void rotate(int x)
{
    int y=t[x].f,z=t[y].f;
    int k=t[y].ch[1]==x;
    t[z].ch[t[z].ch[1]==y]=x;t[x].f=z;
    t[y].ch[k]=t[x].ch[k^1];t[t[x].ch[k^1]].f=y;
    t[x].ch[k^1]=y;t[y].f=x;
    pushup(y);pushup(x);
}
void Splay(int x,int goal)
{
    while(t[x].f!=goal)
    {
        int y=t[x].f,z=t[y].f;
        if(z!=goal)
            (t[y].ch[0]==x)^(t[z].ch[0]==y)?rotate(x):rotate(y);
        rotate(x);
    }
    if(goal<=K)root[goal]=x;
}
void insert(int x,int k,int dd)
{
    int u=root[k],f=k;
    while(u&&t[u].v!=x)
        f=u,u=t[u].ch[t[u].v<x];
    if(u>K)
    {
        t[u].cnt+=dd;t[u].size+=dd;
    }
    else
    {
        u=++tot;
        if(f>K)
            t[f].ch[t[f].v<x]=u;
        t[u].f=f;t[u].v=x;t[u].size=dd;t[u].cnt=dd;
    }
    Splay(u,k);
}
int get(int x,int k)
{
    int u=root[k],re=0;
    while(u)
    {
        if(x>=t[u].v)re+=t[t[u].ch[0]].size+t[u].cnt;
        if(x==t[u].v)break;
        u=t[u].ch[t[u].v<x];
    }
    return re;
}
void Add(int x,int w,int dd)
{
    while(x<=K)
    {
        insert(w,x,dd);
        x+=lowbit(x);
    }
}
int getsum(int x,int w)
{
    int re=0;
    while(x)
    {
        re+=get(w,x);
        x-=lowbit(x);
    }
    return re;
}
inline bool operator<(seq a,seq b)
{
    if(a.a!=b.a)
        return a.a<b.a;
    if(a.b!=b.b)
        return a.b<b.b;
    return a.c<b.c;
}
inline bool operator==(seq a,seq b)
{
    return a.a==b.a&&a.b==b.b&&a.c==b.c;
}
int main()
{
    n=read();tot=K=read();
    for(int i=1;i<=n;++i)s[i].a=read(),s[i].b=read(),s[i].c=read();
    sort(&s[1],&s[n+1]);
    int ccc=0;
    for(int i=1;i<=n;++i)
    {
        ccc++;
        if(s[i]==s[i+1])continue;
        ans[getsum(s[i].b,s[i].c)+ccc-1]+=ccc;
        Add(s[i].b,s[i].c,ccc);
        ccc=0;
    }
    for(int i=0;i<n;++i)
        printf("%d\n",ans[i]);
    return 0;
}