bzoj3262 陌上花开
做数据结构体搜到了这道题
然而轩神告诉我这题应该用CDQ分治（smg 不会）
于是我用了树状数组套treap
首先考虑如果只有两维：第一维可以排序，然后第二维用树状数组$log_2(n)$插入、查询；
那么有三维的话就在树状数组的每个点上挂一颗treap，插入和查询就是$log^2_2(n)$的；
（还有就是要注意三维都相同的情况要取后面的最大值）
1 #include<algorithm> 2 #include<iostream> 3 #include<cstdlib> 4 #include<cstring> 5 #include<cstdio> 6 #include<string> 7 #include<cmath> 8 #include<ctime> 9 #include<queue> 10 #include<stack> 11 #include<map> 12 #include<set> 13 #define rre(i,r,l) for(int i=(r);i>=(l);i--) 14 #define re(i,l,r) for(int i=(l);i<=(r);i++) 15 #define Clear(a,b) memset(a,b,sizeof(a)) 16 #define inout(x) printf("%d",(x)) 17 #define douin(x) scanf("%lf",&x) 18 #define strin(x) scanf("%s",(x)) 19 #define LLin(x) scanf("%lld",&x) 20 #define op operator 21 #define CSC main 22 typedef unsigned long long ULL; 23 typedef const int cint; 24 typedef long long LL; 25 using namespace std; 26 void inin(int &ret) 27 { 28 ret=0;int f=0;char ch=getchar(); 29 while(ch<'0'||ch>'9'){if(ch=='-')f=1;ch=getchar();} 30 while(ch>='0'&&ch<='9')ret*=10,ret+=ch-'0',ch=getchar(); 31 ret=f?-ret:ret; 32 } 33 struct data 34 { 35 int a,b,c; 36 bool op < (const data &rhs)const {return a==rhs.a?b==rhs.b?c<rhs.c:b<rhs.b:a<rhs.a;} 37 }da[100010]; 38 int ch[5000050][2],w[5000050],c[5000050],s[5000050],r[5000050],ed; 39 struct treap 40 { 41 int root; 42 void maintain(int k){if(k)s[k]=c[k]+s[ch[k][0]]+s[ch[k][1]];} 43 void rotate(int &k,int d){int p=ch[k][d^1];ch[k][d^1]=ch[p][d];ch[p][d]=k;maintain(k);maintain(p);k=p;} 44 void add(int &k,int x) 45 { 46 if(!k) 47 { 48 k=++ed,r[k]=rand(); 49 w[k]=x,s[k]=c[k]=1; 50 ch[k][0]=ch[k][1]=0; 51 return ; 52 } 53 s[k]++; 54 if(w[k]==x){c[k]++;return ;} 55 int d=x>w[k];add(ch[k][d],x); 56 } 57 int findrank(int k,int x) 58 { 59 if(!k)return 0; 60 int pp=(ch[k][0]?s[ch[k][0]]:0); 61 if(x==w[k])return pp+c[k]; 62 if(x<w[k])return findrank(ch[k][0],x); 63 else return pp+c[k]+findrank(ch[k][1],x); 64 } 65 }t[200020]; 66 int ans[100010]; 67 int n,m; 68 inline int lowbit(int x){return x&-x;} 69 void add(int k,int x) 70 { 71 while(k<=m) 72 { 73 t[k].add(t[k].root,x); 74 k+=lowbit(k); 75 } 76 } 77 int query(int k,int x) 78 { 79 int ret=0; 80 while(k) 81 { 82 ret+=t[k].findrank(t[k].root,x); 83 k-=lowbit(k); 84 } 85 return ret; 86 } 87 int sum[100010]; 88 int CSC() 89 { 90 freopen("in.in","r",stdin); 91 inin(n);inin(m); 92 re(i,1,n)inin(da[i].a),inin(da[i].b),inin(da[i].c); 93 sort(da+1,da+n+1); 94 re(i,1,n) 95 { 96 if(da[i].a==da[i+1].a&&da[i].b==da[i+1].b&&da[i].c==da[i+1].c) 97 { 98 add(da[i].b,da[i].c); 99 continue; 100 } 101 sum[i]=query(da[i].b,da[i].c); 102 add(da[i].b,da[i].c); 103 } 104 rre(i,n,1) 105 { 106 if(da[i].a==da[i+1].a&&da[i].b==da[i+1].b&&da[i].c==da[i+1].c) 107 sum[i]=sum[i+1]; 108 ans[sum[i]]++; 109 } 110 re(i,0,n-1)printf("%d\n",ans[i]); 111 return 0; 112 }
