【三维偏序】【分块】bzoj3262 陌上花开

裸的三维偏序。 对x坐标排序，y、z坐标分块。复杂度O(n*sqrt(n*log(n)))。代码很短。

 #include<cstdio>
 #include<cmath>
 #include<algorithm>
 #include<vector>
 using namespace std;
 struct Point{int x,y,z,num;void Read(){scanf("%d%d%d",&x,&y,&z);}}p[];
 bool operator < (const Point &a,const Point &b){return a.x<b.x;}
 bool cmp (const Point &a,const Point &b){return a.y<b.y;}
 vector<int>b[];
 vector<Point>a[];
 int n,rank[],m,head,maxv[],sum;
 void makeblock()
 {
     int sz=(int)sqrt((double)n*(log((double)n)/log(2.0))); if(!sz) sz=;
     for(sum=;sum*sz<n;sum++)
       {
         int R=sum*sz;
         for(int i=(sum-)*sz+;i<=R;i++) p[i].num=sum;
         maxv[sum]=p[R].y;
       }
     for(int i=(sum-)*sz+;i<=n;i++) p[i].num=sum;
     maxv[sum]=p[n].y;
 }
 void insert(const Point &U)
 {
     b[U.num].insert(upper_bound(b[U.num].begin(),b[U.num].end(),U.z),U.z);
     a[U.num].push_back(U);
 }
 int query(const Point &U)
 {
     int cnt=,i;
     for(i=;i<=sum && maxv[i]<=U.y;++i)
       cnt+=upper_bound(b[i].begin(),b[i].end(),U.z)-b[i].begin();
     for(vector<Point>::iterator it=a[i].begin();it!=a[i].end();++it)
       if((*it).z<=U.z&&(*it).y<=U.y) ++cnt;
     return cnt;
 }
 int main()
 {
     scanf("%d%d",&n,&m);
     for(int i=;i<=n;i++) p[i].Read();
     sort(p+,p+n+,cmp); makeblock();
     sort(p+,p+n+);
     for(int i=;i<=n;i++)
       {
           if(p[i].x!=p[i-].x) head=i;
         if(p[i].x!=p[i+].x)
           {
               for(int j=head;j<=i;j++) insert(p[j]);
               for(int j=head;j<=i;j++) ++rank[query(p[j])-];
           }
       }
     for(int i=;i<n;i++) printf("%d\n",rank[i]);
     return ;
 }