BZOJ3262: 陌上花开

BZOJ3262

先按第一维排序，然后以第二维为序做CDQ分治，第三维用树状数组维护。
先考虑没有相同属性的花，即任意两朵花A，B之间不存在相互包含关系。
那么正常做就好了。
再考虑加入相同属性的花对答案有什么影响。
那么若A>B，对A的贡献应为numB,numi表示第i种花的个数。
那么对于同一个集合内的相同的花，每一个的答案应+=numi−1
结。

Warning!
排序里面不能偷懒啊。。少写是会粗事的！

【代码】

#include <cstdio>#include <iostream>#include <algorithm>#include <cstring>#include <queue>#include <cmath>#define N 100005#define M 200005#define INF 1<<29using namespace std;typedef long long ll;int read(){    int x=0,f=1;char ch=getchar();    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}    while(isdigit(ch)){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}    return x*f;}int n,K,tot,numa;int ans[N],Ans[N],num[N];class Flower{    public:    int A,B,C;}Flo[N];bool operator <(Flower a,Flower b){    return a.A<b.A||(a.A==b.A&&a.B<b.B)||(a.A==b.A&&a.B==b.B&&a.C<b.C);}bool operator ==(Flower a,Flower b){    return a.A==b.A&&a.B==b.B&&a.C==b.C;}class Query{    public:        int x,y,pos,w;    Query(){}    Query(int xx,int yy,int pp,int ww){        x=xx,y=yy,pos=pp,w=ww;    }}Q[N],tmp[N]; bool operator <(Query a,Query b){    return a.x<b.x||(a.x==b.x&&a.y<b.y)||(a.x==b.x&&a.y==b.y&&a.pos<b.pos);}namespace Szsz{    int szsz[M];    int lowbit(int x){return x&-x;}    void Sum_Up(int x,int y){        for(int i=x;i<=K;i+=lowbit(i)) szsz[i]+=y;    }    int Query(int x){        int rtn=0;        for(int i=x;i;i-=lowbit(i)) rtn+=szsz[i];        return rtn;    }    void Clear(int x){        for(int i=x;i<=K&&szsz[i];i+=lowbit(i)) szsz[i]=0;    }}void CDQ(int l,int r){    if(l==r) return;    int mid=l+r>>1;CDQ(l,mid);CDQ(mid+1,r);    int p=l,q=mid+1,o=0;    while(p<=mid&&q<=r) {        if(Q[p]<Q[q]) {            Szsz::Sum_Up(Q[p].y,Q[p].w);            tmp[o++]=Q[p++];        }        else {            ans[Q[q].pos]+=Szsz::Query(Q[q].y);            tmp[o++]=Q[q++];        }    }    while(p<=mid) {        tmp[o++]=Q[p++];    }    while(q<=r) {        ans[Q[q].pos]+=Szsz::Query(Q[q].y);        tmp[o++]=Q[q++];    }    for(int i=0;i<o;i++)        Szsz::Clear(tmp[i].y),Q[i+l]=tmp[i];}int main(){    n=read();K=read();    for(int i=1;i<=n;i++)        Flo[i].A=read(),Flo[i].B=read(),Flo[i].C=read();    sort(Flo+1,Flo+1+n);    int tot=0,ss=1;    for(int i=1;i<=n;i++)    {        if(Flo[i]==Flo[i+1]) {ss++;continue;}        Q[++tot]=Query(Flo[i].B,Flo[i].C,tot,ss);        num[tot]=ss;ss=1;    }    CDQ(1,tot);    for(int i=1;i<=tot;i++) ans[i]+=num[i]-1;    for(int i=1;i<=tot;i++) Ans[ans[i]]+=num[i];    for(int i=0;i<n;i++) printf("%d\n",Ans[i]);    return 0;}