bzoj3262 陌上花开【解法一】

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的每级花的数量。

解法二树套树见【这里】。
CDQ分治，先按照$s$排序，显然只会有前面的对后面的产生影响。对一个区间$[l,r]$分治处理之后得到分别按照$c$排序的$[l,mid]$和$[mid+1,r]$。这里只需要考虑前一半对后一半的影响。按照$c$归并的同时按照$m$维护树状数组前缀和，更新答案。
注意处理有相同元素的情况。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int rd()
{
    int x=0;
    char c=getchar();
    while (c<'0'||c>'9') c=getchar();
    while (c>='0'&&c<='9')
    {
        x=x*10+c-'0';
        c=getchar();
    }
    return x;
}
struct str
{
    int x,y,z,ans,cnt;
    void read()
    {
        x=rd();
        y=rd();
        z=rd();
        cnt=1;
    }
    bool operator == (const str &ss) const
    {
        return x==ss.x&&y==ss.y&&z==ss.z;
    }
    bool operator < (const str &ss) const
    {
        if (x!=ss.x) return x<ss.x;
        if (y!=ss.y) return y<ss.y;
        return z<ss.z;
    }
}a[100010],b[100010];
int s[200010],ans[100010],n,k;
void add(int p,int x)
{
    for (;p<=k;p+=p&-p) s[p]+=x;
}
int qry(int p)
{
    int ret=0;
    for (;p;p-=p&-p) ret+=s[p];
    return ret;
}
void solve(int l,int r)
{
    if (l==r) return;
    int mid=l+r>>1;
    solve(l,mid);
    solve(mid+1,r);
    for (int i=l,j=mid+1,clo=l;i<=mid||j<=r;)
    {
        if (i<=mid&&(j>r||a[i].y<=a[j].y))
        {
            b[clo++]=a[i];
            add(a[i].z,a[i].cnt);
            i++;
        }
        else
        {
            a[j].ans+=qry(a[j].z);
            b[clo++]=a[j];
            j++;
        }
    }
    for (int i=l;i<=mid;i++) add(a[i].z,-a[i].cnt);
    for (int i=l;i<=r;i++) a[i]=b[i];
}
int main()
{
    int nn;
    nn=rd();
    k=rd();
    for (int i=1;i<=nn;i++) b[i].read();
    sort(b+1,b+nn+1);
    for (int i=1;i<=nn;i++)
        if (!(b[i]==b[i-1])) a[++n]=b[i];
        else a[n].cnt++;
    solve(1,n);
    for (int i=1;i<=n;i++)
        ans[a[i].ans+a[i].cnt-1]+=a[i].cnt;
    for (int i=0;i<nn;i++) printf("%d\n",ans[i]);
}