BZOJ3262: 陌上花开

最新推荐文章于 2020-01-13 20:23:43 发布

CR1SceNT

最新推荐文章于 2020-01-13 20:23:43 发布

阅读量449

点赞数

分类专栏： CDQ分治

本文链接：https://blog.csdn.net/Ep1C_HeReT1c/article/details/72655164

版权

CDQ分治专栏收录该内容

6 篇文章 0 订阅

订阅专栏

BZOJ3262

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

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

【代码】

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <queue>
#include <cmath>
#define N 100005
#define M 200005
#define INF 1<<29
using 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;
}