bzoj 3262: 陌上花开 -- CDQ分治

Posted 2020-10-04 lkhll--qaz

3262: 陌上花开

Time Limit: 20 Sec  Memory Limit: 256 MB

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

#include<map>
#include<cmath>
#include<queue>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define inf 1000000007
#define ll long long
#define N 200010
inline int rd()
{
    int x=0,f=1;char ch=getchar();
    while(ch<‘0‘||ch>‘9‘){if(ch==‘-‘)f=-1;ch=getchar();}
    while(ch>=‘0‘&&ch<=‘9‘){x=x*10+ch-‘0‘;ch=getchar();}
    return x*f;
}
struct qaz{int a,b,c,s,ans;}a[N],p[N];
bool cmp(qaz a,qaz b){return a.b==b.b&&a.a==b.a?a.c<b.c:(a.a==b.a?a.b<b.b:a.a<b.a);}
bool cmp2(qaz a,qaz b){return a.b==b.b?a.c<b.c:a.b<b.b;}
int n,k,nn,sum[N];
void add(int x,int v){for(int i=x;i<=k;i+=i&(-i)) sum[i]+=v;}
int fd(int x){int tt=0;for(int i=x;i;i-=i&(-i)) tt+=sum[i];return tt;}
void sol(int l,int r)
{
    if(l==r) return;
    int mid=l+r>>1,i=l,j=mid+1;
    sol(l,mid);sol(mid+1,r);
    sort(p+l,p+mid+1,cmp2);
    sort(p+mid+1,p+r+1,cmp2);
    while(j<=r)
    {
        while(i<=mid&&p[i].b<=p[j].b) add(p[i].c,p[i].s),i++;
        p[j].ans+=fd(p[j].c);j++;
    }
    for(j=l;j<i;j++) add(p[j].c,-p[j].s);
}
int ans[N];
int main()
{
    nn=rd();k=rd();
    for(int i=1;i<=nn;i++)
    {
        a[i].a=rd();a[i].b=rd();a[i].c=rd();
    }
    sort(a+1,a+nn+1,cmp);
    for(int i=1,ji=0;i<=nn;i++)
    {
        ji++;
        if(a[i].a!=a[i+1].a||a[i].b!=a[i+1].b||a[i].c!=a[i+1].c)
        {
            p[++n]=a[i];
            p[n].s=ji;ji=0;
        }
    }
    sol(1,n);
    for(int i=1;i<=n;i++) ans[p[i].ans+p[i].s-1]+=p[i].s;
    for(int i=0;i<nn;i++) printf("%d\n",ans[i]);
    return 0;
}