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

【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
",ans[i]);
	return 0;
}