poj2761&&poj2104 主席树(静态区间第K大)
2015-08-19 15:47
387 查看
搞了两天参考大牛的博客终于写除了自己的主席树但是一直挖,在离散化时关于是否有相同的数,如果不考虑的话就莫名其妙的ac了,考虑的话就一直wa。
主席树:
好高大上的名字,说白了就是以序列的所有前缀或者后缀建立起来的线段树,同时在建树的过程中还用了上一个线段树的树根。
线段树里记录了其前缀序列里出现的值的次数。
离散后标记为1~n,要建立n个线段树从root[1]到root
对于1~n中的第K小
从根节点开始,如果左儿子保存的数的数量大于K递归处理左儿子,寻找左儿子中第k小的数,否则在右儿子中寻找第k-(左儿子中元素个数)
对于[L, R]区间中的第k小
假设我们知道了[1,R]区间中P结点的值域中所含的元素个数,也知道[1,L-1]区间中P结点的值域中所包含的元素个数,显然用第一个个数减去第二个个数,就可以得到[L,R]区间中的元素个数。
在这里贴下老板的博客:
http://blog.csdn.net/neko01/article/details/39326365
代码如下:
主席树:
好高大上的名字,说白了就是以序列的所有前缀或者后缀建立起来的线段树,同时在建树的过程中还用了上一个线段树的树根。
线段树里记录了其前缀序列里出现的值的次数。
离散后标记为1~n,要建立n个线段树从root[1]到root
对于1~n中的第K小
从根节点开始,如果左儿子保存的数的数量大于K递归处理左儿子,寻找左儿子中第k小的数,否则在右儿子中寻找第k-(左儿子中元素个数)
对于[L, R]区间中的第k小
假设我们知道了[1,R]区间中P结点的值域中所含的元素个数,也知道[1,L-1]区间中P结点的值域中所包含的元素个数,显然用第一个个数减去第二个个数,就可以得到[L,R]区间中的元素个数。
在这里贴下老板的博客:
http://blog.csdn.net/neko01/article/details/39326365
代码如下:
[code]#include<cstdio> #include<cstring> #include<iostream> #include<algorithm> #include<cmath> #include<cstdlib> #define ls(i) T[(i)].ls #define rs(i) T[(i)].rs #define w(i) T[(i)].w using namespace std; int n,m; const int maxn = 100000 + 5; int sta[maxn], p[maxn], root[maxn], b[maxn]; struct Tree{ int ls, rs, w; Tree(){ls = rs = w = 0;} }T[maxn*20]; bool cmp(int i, int j) { return sta[i] < sta[j]; } int siz; void Insert(int &i, int l, int r, int x) { T[++siz] = T[i]; i = siz; w(i)++; if (l == r) return; int mid = (l + r) >> 1; if (x <= mid) Insert(ls(i), l, mid, x); else Insert(rs(i), mid + 1, r, x); } int query(int i, int j, int l, int r, int k) { if (l == r) return l; int mid = (l + r) >> 1; int t = w(ls(j)) - w(ls(i)); if (k <= t) return query(ls(i), ls(j), l, mid, k); else return query(rs(i), rs(j), mid+1, r, k-t); } int main() { //freopen("in", "r", stdin); while (cin >> n >> m) { for (int i = 1; i <= n; i++) { scanf("%d", &sta[i]); p[i] = i; } sort(p+1, p+1+n, cmp); for (int i = 1; i <= n; i++) b[p[i]] = i;//离散化 root[0] = 0; siz = 0; for (int i = 1; i <= n; i++) { root[i] = root[i-1]; Insert(root[i], 1, n, b[i]); } for (int i = 0; i < m; i++) { int x, y, k; scanf("%d%d%d", &x, &y, &k); int t = query(root[x-1], root[y], 1, n, k); printf("%d\n", sta[p[t]]); } } return 0; }
相关文章推荐
- Servlet 容器如何处理请求资源路径
- c++ primer plus 第七章《编程题7.13.10》
- 桌面支持--adobe reader设置记录上次阅读位置
- SQL大量数据查询的优化
- Quartz的任务的临时启动和暂停和恢复
- poj-2531 Network Saboteur DFS
- 进程/线程同步的方式和机制,进程间通信
- SPFA算法
- mysql general log 查看mysql 执行历史
- 【转】Python调用C函数
- 收起键盘
- 软件测试方法-边界值分析法
- ios 对齐属性
- 桌面支持--512-Rear chassis fan not detected
- 桌面支持--不懂不要乱动-尤其是别人的东西
- java向文件中追加内容与读写文件内容源码
- hdu 5396 Expression 2015多校联合训练赛#9 动态规划 dp
- .Net程序员安卓学习之路2:访问网络API
- MVC Razor模板引擎 @RenderBody、@RenderPage、@RenderSection及Html.RenderPartial、Html.RenderAction
- 齐次坐标