【BZOJ】【3262】陌上花开

树套树

---

　　orz zyf

　　这题的思路……算是让我了解到了树套树的一种用途吧

三维。。。第一维排序，第二维树状数组，第三维treap

具体实现就是每个树状数组的节点保存一颗treap，然后就可以查询了。

　　好神啊……

　　树套树可以方便的进行特殊的区间求和，大多数满足区间加法的运算都可以用树套树来搞，比如这题的对所有第二维坐标小于当前点（相当于一段前缀）求第三维坐标小于当前点的和。= =啊……就是一个特殊的前缀和嘛……

 

　　还有就是这题其实有个K的……表示坐标的范围……所以树状数组要一直到k，而不是n

 /**************************************************************
     Problem: 3262
     User: Tunix
     Language: C++
     Result: Accepted
     Time:6728 ms
     Memory:99128 kb
 ****************************************************************/
 
 //BZOJ 3262
 #include<cmath>
 #include<vector>
 #include<cstdio>
 #include<cstring>
 #include<cstdlib>
 #include<iostream>
 #include<algorithm>
 #define rep(i,n) for(int i=0;i<n;++i)
 #define F(i,j,n) for(int i=j;i<=n;++i)
 #define D(i,j,n) for(int i=j;i>=n;--i)
 #define pb push_back
 #define CC(a,b) memset(a,b,sizeof(a))
 using namespace std;
 int getint(){
     int v=,sign=; char ch=getchar();
     while(!isdigit(ch)) {if(ch=='-') sign=-; ch=getchar();}
     while(isdigit(ch))  {v=v*+ch-''; ch=getchar();}
     return v*sign;
 }
 const int N=,M=,INF=~0u>>;
 const double eps=1e-;
 /*******************template********************/
 int n,k,tot,l[M],r[M],s[M],rnd[M],w[M],v[M];
 #define L l[x]
 #define R r[x]
 inline void Push_up(int x){ s[x]=s[L]+s[R]+w[x]; }
 inline void zig(int &x){int t=L; L=r[t]; r[t]=x; s[t]=s[x]; Push_up(x); x=t;}
 inline void zag(int &x){int t=R; R=l[t]; l[t]=x; s[t]=s[x]; Push_up(x); x=t;}
 void ins(int &x,int num,int z){
     if (!x){
         x=++tot; v[x]=num; s[x]=w[x]=z; L=R=; rnd[x]=rand(); return;
     }
     s[x]+=z;
     if (v[x]==num) w[x]+=z;
     else if(num<v[x]){
         ins(L,num,z); if(rnd[L]<rnd[x]) zig(x);
     }else{
         ins(R,num,z); if(rnd[R]<rnd[x]) zag(x);
     }
 }
 int rank(int x,int num){
     if (!x) return ;
     if (v[x]==num) return s[L]+w[x];
     else if(num<v[x]) return rank(L,num);
     else return s[L]+w[x]+rank(R,num);
 }
 #undef L
 #undef R
 /*********************Treap*********************/
 int rt[N];
 void update(int x,int y,int z){
     for(int i=x;i<=k;i+=i&-i) ins(rt[i],y,z);
 }
 int query(int x,int y){
     int ans=;
     for(int i=x;i;i-=i&-i) ans+=rank(rt[i],y);
     return ans;
 }
 /*******************Fenwick*********************/
 struct data{
     int s,c,m;
 }a[N];
 bool cmp(data a,data b){
     if (a.s==b.s){
         if (a.c==b.c) return a.m<b.m;
         return a.c<b.c;
     }
     return a.s<b.s;
 }
 int ans[N];
 int main(){
 #ifndef ONLINE_JUDGE
     freopen("3262.in","r",stdin);
     freopen("3262.out","w",stdout);
 #endif
     n=getint(); k=getint();
     F(i,,n){
         a[i].s=getint();
         a[i].c=getint();
         a[i].m=getint();
     }
     sort(a+,a+n+,cmp);
     int num=;
     F(i,,n){
         if (a[i].s==a[i+].s &&
             a[i].c==a[i+].c &&
             a[i].m==a[i+].m) num++;
         else{
             ans[query(a[i].c,a[i].m)+num-]+=num;
             update(a[i].c,a[i].m,num);
             num=;
         }
     }
     rep(i,n) printf("%d\n",ans[i]);
     return ;
 }