算法复习——cdq分治

题目：

Description

有n朵花，每朵花有三个属性：花形(s)、颜色(c)、气味(m)，又三个整数表示。现要对每朵花评级，一朵花的级别是它拥有的美丽能超过的花的数量。定义一朵花A比另一朵花B要美丽，当且仅当Sa>=Sb,Ca>=Cb,Ma>=Mb。显然，两朵花可能有同样的属性。需要统计出评出每个等级的花的数量。

Input

第一行为N,K (1 <= N <= 100,000, 1 <= K <= 200,000 ), 分别表示花的数量和最大属性值。

以下N行，每行三个整数si, ci, mi (1 <= si, ci, mi <= K)，表示第i朵花的属性

Output

包含N行，分别表示评级为0...N-1的每级花的数量。

Sample Input

10 3
3 3 3
2 3 3
2 3 1
3 1 1
3 1 2
1 3 1
1 1 2
1 2 2
1 3 2
1 2 1

Sample Output

3
1
3
0
1
0
1
0
0
1

HINT

1 <= N <= 100,000, 1 <= K <= 200,000

Source

树套树 CDQ分治

题解：

CDQ分治解三维偏序的模板题，以下是关于cdq分治的大概介绍

CDQ的核心就是分为左右两组后，算左边对右边的贡献····

但CDQ分治运用最广泛的是解决偏序的问题···比如上面的陌上花开就是解决关于(a,b,c)的三维偏序问题····

解决这类问题的核心是降维···一个一个忽略每一维影响···比如这道题，先以a优先排序，从而只用管b和c···然后在分成左右两组时每组以b优先排序··进而最后利用树状数组解决c

详细见：http://www.cnblogs.com/mlystdcall/p/6219421.html

不得不说真神奇·····

代码：

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<ctime>
#include<cctype>
#include<cstring>
#include<string>
#include<algorithm>
using namespace std;
const int N=1e5+;
const int M=2e5+;
inline int R()
{
  char c;int f=,i=;
  for(c=getchar();(c<''||c>'')&&c!='-';c=getchar());
  if(c=='-')  i=-,c=getchar();
  for(;c<=''&&c>='';c=getchar())
    f=(f<<)+(f<<)+c-'';
  return f*i;
}
struct node
{
  int x,y,z,cnt,ans;
}f[N],a[N];
bool cmp(node a,node b)
{
  if(a.x==b.x&&a.y==b.y)  return a.z<b.z;
  else if(a.x==b.x)  return a.y<b.y;
  else return a.x<b.x;
}
int m,n,k,tot,tree[M],ans[N];
inline bool operator < (node a,node b)
{
  if(a.y==b.y)  return a.z<b.z;
  else return a.y<b.y;
}
inline void insert(int u,int w)
{
  for(int i=u;i<=k;i+=(i&(-i)))
    tree[i]+=w;
}
inline int query(int u)
{
  int temp=;
  for(int i=u;i>;i-=(i&(-i)))
    temp+=tree[i];
  return temp;
}
inline void solve(int l,int r)
{
  if(l==r)  return;
  int mid=(l+r)/,i,j;
  solve(l,mid),solve(mid+,r);
  sort(a+l,a+mid+);sort(a+mid+,a+r+);
  i=l,j=mid+;
  while(j<=r)
  {
    while(i<=mid&&a[i].y<=a[j].y)
    {
      insert(a[i].z,a[i].cnt);i++;
    }
    a[j].ans+=query(a[j].z);j++;
  }
  for(j=l;j<i;j++)  insert(a[j].z,-a[j].cnt);
}
int main()
{
  //freopen("a.in","r",stdin);
  m=R(),k=R();
  for(int i=;i<=m;i++)
  {
    f[i].x=R();f[i].y=R();f[i].z=R();
  }
  sort(f+,f+m+,cmp);
  a[++n]=f[];tot=;
  a[].cnt=;
  for(int i=;i<=m;i++)
  {
    if(f[i].x!=f[i-].x||f[i].y!=f[i-].y||f[i].z!=f[i-].z)  ++n,tot=;
    else tot++;
    a[n]=f[i],a[n].cnt=tot;
  }
  solve(,n);
  for(int i=;i<=n;i++)
    ans[a[i].ans+a[i].cnt-]+=a[i].cnt;
  for(int i=;i<m;i++)
    printf("%d\n",ans[i]);
  return ;
}