BZOJ3262: 陌上花开
题意概述
有N朵花,对于每一朵花,有三个属性:s,c,m. 当且仅当si>sj,ci>cj,mi>mj,有花i比花j美丽. 一朵花的评级为比其他花更美丽的数量(不包括自己),输出评级为0 N−1的花的数量.
题目分析:
初学CDQ分治,写了一下午.orz……
这是一个三维偏序的问题,可以用CDQ来搞一搞. 当然此题当中,一朵花既是修改操作,又是询问操作. 先以一维s从小到大排序,开始CDQ分治. 对于每一次完成两个子问题后,合并时,先按照c从小到大排序,第三维m用树状数组维护,对于当前花,若是左区间的花,树状数组add操作,否则,查询m比其小的花的个数.
(千万要注意三个属性均相同的花的处理,一种方法是先去重,再计算)
代码:
//Continue#include<cstdio>#include<iostream>#include<algorithm>using namespace std;const int maxn=100000+10;int read() { char ch=getchar();int ret=0; while(ch<'0'||ch>'9') ch=getchar(); while(ch>='0'&&ch<='9') ret=ret*10+ch-'0',ch=getchar(); return ret;}struct Flower { int s,c,m,k,cnt,ans;//k:0 add操作,1 sum操作 //cnt记录这种花有多少,ans指的是这种花的评级 Flower(int cnt=0,int ans=0){} bool Operator < (const Flower& rhs) const { return c<rhs.c||(c==rhs.c&&(m<rhs.m||(m==rhs.m&&s<rhs.s))); } bool operator == (const Flower& rhs) const { return s==rhs.s&&c==rhs.c&&m==rhs.m; } void input() { s=read();c=read();m=read(); }}t[maxn],f[maxn];#define lowbit(x) (x&-x)int bit[maxn],N,K;void add(int d,int v) { for(;d<=K;d+=lowbit(d)) bit[d]+=v;}int sum(int d,int ret=0) { for(;d;d-=lowbit(d)) ret+=bit[d]; return ret;}#define mid ((l+r)>>1)int cmp(const Flower& a,const Flower& b) { return a.s<b.s||(a.s==b.s&&(a.c<b.c||(a.c==b.c&&a.m<b.m)));}void solve(int l,int r) { if(l==r) return ; //先递归解决左右子区间内部更新 solve(l,mid); solve(mid+1,r); //此时解决左右区间间的更新 for(int i=l;i<=mid;i++) f[i].k=0;//左区间修改 for(int i=mid+1;i<=r;i++) f[i].k=1;//右区间查询 sort(f+l,f+r+1);//按照c排序 for(int i=l;i<=r;i++) if(f[i].k) f[i].ans+=sum(f[i].m);//右区间查询 else add(f[i].m,f[i].cnt);//左区间修改 //注意在用分治求解的时候不能引入序列长度的时间复杂度,否则会造成时间复杂度退化 for(int i=l;i<=r;i++) if(!f[i].k) add(f[i].m,-f[i].cnt);//采用倒着修改回去的方法 }int ans[maxn],tot;int main() { N=read();K=read(); for(int i=1;i<=N;i++) t[i].input(); sort(t+1,t+N+1,cmp); for(int i=1;i<=N;i++)//先去重操作 if(i==1||!(t[i]==t[i-1])) f[++tot]=t[i],f[tot].cnt=1; else ++f[tot].cnt; for(int i=1;i<=tot;i++) f[i].ans=f[i].cnt-1;//一开始评级赋值为个数-1 sort(f+1,f+tot+1,cmp);//按s排序 solve(1,tot); for(int i=1;i<=tot;i++) ans[f[i].ans]+=f[i].cnt; for(int i=0;i<N;i++) PRintf("%d/n",ans[i]); return 0;}
