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 }

 

posted @ 2016-01-22 12:53  HugeGun  阅读(738)  评论(0编辑  收藏  举报