【BZOJ3262】陌上花开(树套树)

时间:2024-08-13 13:03:32

【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;
}