博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
2019 Multi-University Training Contest 4 - K-th Closest Distance
阅读量:4951 次
发布时间:2019-06-11

本文共 2861 字,大约阅读时间需要 9 分钟。

主席树 + 二分答案

对于这种区间内的值域问题一般用主席树进行求解。

因为数据范围只有1e6,所以不用离散化,直接建树即可。

题目要求找到区间内离p第k近的数,可以想到,这个问题具有单调性(某个区间长度有大于k个值,那么比这个区间长度更长的比如也有大于k个值)

于是我们可以二分答案这个距离,枚举离p距离为mid的范围(max(1, p - mid), min(p + mid, 1e6)

这样我们每次主席树上询问这个值域范围内的个数是否大于等于k判断枚举的mid是否可行,最后mid停下来的最小值就是答案。

#include 
#define INF 0x3f3f3f3f#define full(a, b) memset(a, b, sizeof a)#define __fastIn ios::sync_with_stdio(false), cin.tie(0)#define pb push_backusing namespace std;typedef long long LL;inline int lowbit(int x){ return x & (-x); }inline int read(){ int ret = 0, w = 0; char ch = 0; while(!isdigit(ch)){ w |= ch == '-', ch = getchar(); } while(isdigit(ch)){ ret = (ret << 3) + (ret << 1) + (ch ^ 48); ch = getchar(); } return w ? -ret : ret;}inline int lcm(int a, int b){ return a / __gcd(a, b) * b; }template
inline A fpow(A x, B p, C lyd){ A ans = 1; for(; p; p >>= 1, x = 1LL * x * x % lyd)if(p & 1)ans = 1LL * x * ans % lyd; return ans;}const int N = 200005;const int M = 1000005;int _, n, m, tot, a[N], lc[M*20], rc[M*20], tree[M*20], root[N];int buildTree(int l, int r){ int cur = ++ tot; if(l == r) return cur; int mid = (l + r) >> 1; buildTree(l, mid); buildTree(mid + 1, r); return cur;}int insert(int rt, int l, int r, int val){ int cur = ++ tot; tree[cur] = tree[rt] + 1, lc[cur] = lc[rt], rc[cur] = rc[rt]; if(l == r) return cur; int mid = (l + r) >> 1; if(val <= mid) lc[cur] = insert(lc[rt], l, mid, val); else rc[cur] = insert(rc[rt], mid + 1, r, val); return cur;}int query(int a, int b, int l, int r, int ql, int qr){ if(l == ql && r == qr){ return tree[b] - tree[a]; } int mid = (l + r) >> 1; if(qr <= mid) return query(lc[a], lc[b], l, mid, ql, qr); else if(ql > mid) return query(rc[a], rc[b], mid + 1, r, ql, qr); return query(lc[a], lc[b], l, mid, ql, mid) + query(rc[a], rc[b], mid + 1, r, mid + 1, qr);}bool calc(int l, int r, int p, int k, int mid){ int ret = query(root[l - 1], root[r], 1, M, max(1, p - mid), min(p + mid, M)); return ret >= k;}int main(){ //freopen("data.txt", "r", stdin); for(_ = read(); _; _ --){ tot = 0; n = read(), m = read(); for(int i = 1; i <= n; i ++) a[i] = read(); root[0] = buildTree(1, M); for(int i = 1; i <= n; i ++){ root[i] = insert(root[i - 1], 1, M, a[i]); } int ans = 0; while(m --){ int L = read(), R = read(), p = read(), k = read(); L ^= ans, R ^= ans, p ^= ans, k ^= ans; int l = 0, r = 1e6; while(l < r){ int mid = (l + r) >> 1; if(calc(L, R, p, k, mid)) r = mid; else l = mid + 1; } printf("%d\n", l); ans = l; } } return 0;}

转载于:https://www.cnblogs.com/onionQAQ/p/11280424.html

你可能感兴趣的文章
OS笔记047代理传值和block传值
查看>>
Android应用程序与SurfaceFlinger服务的连接过程分析
查看>>
coco2dx服务器简单例子
查看>>
Java回顾之多线程
查看>>
sqlite
查看>>
机电行业如何进行信息化建设
查看>>
Windows Azure Platform Introduction (4) Windows Azure架构
查看>>
【转】chrome developer tool 调试技巧
查看>>
mahout运行测试与kmeans算法解析
查看>>
互相给一巴掌器
查看>>
Android SDK环境变量配置
查看>>
VM10虚拟机安装图解
查看>>
9、总线
查看>>
Git 笔记 - section 1
查看>>
JZOJ 4.1 B组 俄罗斯方块
查看>>
HDU6409 没有兄弟的舞会
查看>>
2018 Multi-University Training Contest 10 - TeaTree
查看>>
HDU6205 card card card
查看>>
2018 Multi-University Training Contest 10 - Count
查看>>
HDU6198 number number number
查看>>