BZOJ3262陌上花开 CDQ分治_BIT+Treap
Posted
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了BZOJ3262陌上花开 CDQ分治_BIT+Treap相关的知识,希望对你有一定的参考价值。
三个属性, 第一个属性用cdq分治处理, 以第一个属性为关键字从小到大排序, 那么考虑一朵花的等级, 只需考虑排在其前面的花的其他属性(特殊情况是有相同的花,根据题意,对一段相同的花,以排在最后的一朵花的答案为准), 第二三维可以用树状数组加Treap解决, 以每朵花第二属性数值作为位置(因为最大属性k < 2e5, 可以不用离散化, 直接用属性的数值对应树状数组中的下标), 树状数组的每个节点建一颗Treap, 这颗Treap里存的是相应区间里的花的第三个属性, 询问时类似于树状数组求前缀和, 依次将询问的花的位置的前面这一段分成不超过logk棵Treap对应的区间, 在这些区间里可以找出, 由于是以第二属性为位置插入到树状数组里的, 保证前面一段区间的第二属性都是不超过我们询问的花的第二属性, 只要在这不超过logk棵Treap里找出第三属性不超过询问的花的数量, 就是三个属性均不超过询问的花的三个属性的花的数量
#include<algorithm> #include<iostream> #include<cstring> #include<cstdio> using namespace std; const int maxn = 2e5 + 7, maxk = 2e5 + 7, maxnd = 1e6 + 7; void readin(int &ret) { ret = 0; int f = 1; char ch = getchar(); while(ch < ‘0‘ || ch > ‘9‘) {if(ch == ‘-‘) f = -1; ch = getchar();} while(ch >= ‘0‘ && ch <= ‘9‘) ret *= 10, ret += ch - ‘0‘, ch = getchar(); ret *= f; } struct flower { int s, c, m, id, ans; void read() { readin(s); readin(c); readin(m); } } f[maxn]; bool cmp(const flower& a, const flower& b) { return a.s < b.s || (a.s == b.s && a.c < b.c) || (a.s == b.s && a.c == b.c && a.m < b.m); } bool ssame(const flower& a, const flower& b) { return a.s == b.s && a.c == b.c && a.m == b.m; } int n, k, siz, ans[maxn], root[maxn], sz[maxnd], val[maxnd], rnd[maxnd], cnt[maxnd], c[maxnd][2]; void update(int k) { sz[k] = sz[c[k][0]] + sz[c[k][1]] + cnt[k]; } void rotate(int &k, int ch) { int t = c[k][ch]; c[k][ch] = c[t][ch ^ 1]; c[t][ch ^ 1] = k; sz[t] = sz[k]; update(k); k = t; } void insert(int &k, int x) { if(!k) { val[k = ++siz] = x; rnd[k] = rand(); sz[k] = cnt[k] = 1; c[k][0] = c[k][1] = 0; } else if(val[k] == x) sz[k]++, cnt[k]++; else { sz[k]++; int d = val[k] < x; insert(c[k][d], x); if(rnd[c[k][d]] < rnd[k]) rotate(k, d); } } int find(int k, int x) { if(!k) return 0; if(val[k] == x) return sz[c[k][0]] + cnt[k]; if(val[k] > x) return find(c[k][0], x); return sz[c[k][0]] + cnt[k] + find(c[k][1], x); } //find(k, x)找出节点k为根的子树中, 权值不大于x的节点数量 #define lowbit(x) (x&-x) void ins(int pos, int x) { while(pos <= k) { insert(root[pos], x); pos += lowbit(pos); } } int query(int pos, int x) { int ret = 0; while(pos > 0) { ret += find(root[pos], x); pos -= lowbit(pos); } return ret; } int top, stk[maxn], num[maxn]; int main() { readin(n); readin(k); for(int i = 1; i <= n; i++) f[i].read(); sort(f + 1, f + n + 1, cmp); top = 0; for(int i = 1; i <= n; i++) { if(ssame(f[i], f[i + 1])) stk[++top] = i;//如果有相同的花, 先存到数组里, 直到处理了最后一朵这种花再更新前面的答案 else { f[i].ans = query(f[i].c, f[i].m); while(top) f[stk[top--]].ans = f[i].ans; } ins(f[i].c, f[i].m);//询问了这朵花后再将这朵花插入树状数组中 } for(int i = 1; i <= n; i++) num[f[i].ans]++; for(int i = 0; i <= n - 1; i++) printf("%d\n", num[i]); return 0; }
以上是关于BZOJ3262陌上花开 CDQ分治_BIT+Treap的主要内容,如果未能解决你的问题,请参考以下文章