bzoj 3262

题意：给你一些三维上的点，对于每个点，统计三个坐标都小于等于该点的点数。

如果点的范围在300以内，可以用三维树状数组搞，但这题坐标范围太大。

考虑将所有点按照x坐标排序，从左到右，相当于在一个二维平面上插入点，并询问某个点左下方的点数，而后者可以按时间分治，在O(nloglog)复杂度内搞定。（其实可以把x轴看成时间轴）

注意相同点的处理。

（感觉对时间分治就是：将左区间和右区间的修改对其后面的询问的贡献处理掉，合并时就只有左区间的修改和右区间的询问，这就达到了以一个log的复杂度将"修改 询问 修改 修改 询问...”的问题变成了“修改 修改 修改 询问 询问”）

 #include <cstdio>
 #include <algorithm>
 #define N 100010
 #define K 200010
 using namespace std;

 struct Trid {
     int x, y, z;
     int cnt;
     int xx;
     Trid(){}
     Trid( int x, int y, int z ):x(x),y(y),z(z),cnt(),xx(x){}
     bool operator==( const Trid & o ) const {
         return x==o.x && y==o.y && z==o.z;
     }
     bool operator<( const Trid &b ) const {
         if( x!=b.x ) return x<b.x;
         if( y!=b.y ) return y<b.y;
         return z<b.z;
     }
 };

 int n, k;
 Trid trid[N];
 int bit[K];
 int ans[N];

 void modify( int pos, int v ) {
     for( int i=pos; i<=k; i+=i&-i )
         bit[i] += v;
 }
 int query( int pos ) {
     int rt=;
     for( int i=pos; i; i-=i&-i )
         rt += bit[i];
     return rt;
 }
 bool cmp_yzx( const Trid &a, const Trid &b ) {
     if( a.y!=b.y ) return a.y<b.y;
     if( a.z!=b.z ) return a.z<b.z;
     return a.x<b.x;
 }
 void cdq( int lf, int rg ) {
     if( lf==rg ) return;
     int mid=(lf+rg)>>;
     cdq(lf,mid);
     cdq(mid+,rg);
     sort( trid+lf, trid+rg+, cmp_yzx );
     for( int i=lf; i<=rg; i++ ) {
         if( trid[i].x<=mid )
             modify( trid[i].z, + );
         if( trid[i].x>mid )
             trid[i].cnt += query( trid[i].z );
     }
     for( int i=lf; i<=rg; i++ )
         if( trid[i].x<=mid )
             modify( trid[i].z, - );
 }
 int main() {
     scanf( "%d%d", &n, &k );
     for( int i=,x,y,z; i<=n; i++ ) {
         scanf( "%d%d%d", &x, &y, &z );
         trid[i] = Trid( x, y, z );
     }
     sort( trid+, trid++n );
     for( int i=,j; i<=n; i++ ) {
         for( j=i; j<=n && trid[j]==trid[i]; j++ );
         for( int k=i; k<j; k++ )
             trid[k].cnt += j-k-;
         i = j-;
     }
     for( int i=; i<=n; i++ )
         trid[i].x = i;
     cdq( , n );
     for( int i=; i<=n; i++ ) {
         ans[trid[i].cnt]++;
     }
     for( int i=; i<n; i++ )
         printf( "%d\n", ans[i] );
 }