博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
洛谷P4169 天使玩偶 (算竞进阶习题)
阅读量:6758 次
发布时间:2019-06-26

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

CDQ分治

分成四个方向讨论最小值,把所有坐标全部离线处理。

把左边按x轴排序,保证x的顺序,然后树状数组维护每个方向需要的最值。。

然后CDQ分治。。必须手动撤销树状数组的修改,保证分治的时间复杂度。

#include 
#define INF 0x3f3f3f3f#define full(a, b) memset(a, b, sizeof a)#define FAST_IO ios::sync_with_stdio(false), cin.tie(0), cout.tie(0)using namespace std;typedef long long ll;inline int lowbit(int x){ return x & (-x); }inline int read(){ int X = 0, w = 0; char ch = 0; while(!isdigit(ch)) { w |= ch == '-'; ch = getchar(); } while(isdigit(ch)) X = (X << 3) + (X << 1) + (ch ^ 48), ch = getchar(); return w ? -X : X;}inline int gcd(int a, int b){ return b ? gcd(b, a % b) : a; }inline int lcm(int a, int b){ return a / gcd(a, b) * b; }template
inline T max(T x, T y, T z){ return max(max(x, y), z); }template
inline T min(T x, T y, T z){ return min(min(x, y), z); }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 = 1000005;int n, m, tot, t[N], k, ans[N];struct rec { int x, y, z; bool operator < (const rec &rhs) const { return x < rhs.x; }}a[N], b[N];inline void add(int index, int val){ for(; index < tot; index += lowbit(index)) t[index] = max(t[index], val);}inline int query(int index){ int ret = -(1 << 30); for(; index; index -= lowbit(index)) ret = max(ret, t[index]); return ret;}inline void calc1(){ for(int i = 1; i <= k; i ++){ int y = b[i].y; int temp = b[i].x + b[i].y; if(a[b[i].z].z == 1) add(y, temp); else ans[b[i].z] = min(ans[b[i].z], abs(temp - query(y))); } for(int i = 1; i <= k; i ++){ int y = b[i].y; if(a[b[i].z].z == 1){ for(int j = y; j < tot; j += lowbit(j)) t[j] = -(1 << 30); } }}inline void calc2(){ for(int i = 1; i <= k; i ++){ int y = b[i].y; int temp = b[i].x - b[i].y; if(a[b[i].z].z == 1) add(tot - y, temp); else ans[b[i].z] = min(ans[b[i].z], abs(temp - query(tot - y))); } for(int i = 1; i <= k; i ++){ int y = b[i].y; if(a[b[i].z].z == 1){ for(int j = tot - y; j < tot; j += lowbit(j)) t[j] = -(1 << 30); } }}inline void calc3(){ for(int i = k; i >= 1; i --){ int y = b[i].y; int temp = -b[i].x - b[i].y; if(a[b[i].z].z == 1) add(tot - y, temp); else ans[b[i].z] = min(ans[b[i].z], abs(temp - query(tot - y))); } for(int i = k; i >= 1; i --){ int y = b[i].y; if(a[b[i].z].z == 1){ for(int j = tot - y; j < tot; j += lowbit(j)) t[j] = -(1 << 30); } }}inline void calc4(){ for(int i = k; i >= 1; i --){ int y = b[i].y; int temp = -b[i].x + b[i].y; if(a[b[i].z].z == 1) add(y, temp); else ans[b[i].z] = min(ans[b[i].z], abs(temp - query(y))); } for(int i = k; i >= 1; i --){ int y = b[i].y; if(a[b[i].z].z == 1){ for(int j = y; j < tot; j += lowbit(j)) t[j] = -(1 << 30); } }}void cdqDiv(int l, int r){ if(l == r) return; int mid = (l + r) >> 1; cdqDiv(l, mid), cdqDiv(mid + 1, r); k = 0; for(int i = l; i <= r; i ++){ if((i <= mid && a[i].z == 1) || (i > mid && a[i].z == 2)) b[++k] = a[i], b[k].z = i; } sort(b + 1, b + k + 1); calc1(), calc2(); calc3(), calc4();}int main(){ full(t, 0xcf), full(ans, INF); n = read(), m = read(); m += n; for(int i = 1; i <= n; i ++){ a[i].x = read(), a[i].y = read() + 1, a[i].z = 1; } for(int i = n + 1; i <= m; i ++){ a[i].z = read(), a[i].x = read(), a[i].y = read() + 1; } for(int i = 1; i <= m; i ++) tot = max(tot, a[i].y); tot ++; cdqDiv(1, m); for(int i = 1; i <= m; i ++){ if(a[i].z == 2) printf("%d\n", ans[i]); } return 0;}

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

你可能感兴趣的文章
存储虚拟化技术之解读
查看>>
X5Pro惊艳双子塔,vivo国际化渐入佳境
查看>>
(运维)VMware-vShield 5.0 安装与部署
查看>>
报表,是件容易的事吗?
查看>>
技术,技术人员,谁是风,谁是草
查看>>
Android应用程序内部启动Activity过程(startActivity)的源代码分析
查看>>
iOS网络编程-ASIHTTPRequest框架同步请求
查看>>
Gartner:2011年北美MSS市场分析
查看>>
一个网上找到的,在Grid中嵌套Grid的示例:Nested Grids Example
查看>>
IT群侠传第四回大鱼小虾
查看>>
《Python从小白到大牛》第9章 数据结构
查看>>
Xcode中四种build for 的区别
查看>>
酷客多小程序百城宣讲会-嵊州站完美落幕
查看>>
搞机年代,ivvi用“爱情”细分市场
查看>>
学员问答之3-View桌面问题
查看>>
思科路由器开机进入 miniIOS 的原因分析
查看>>
卢松松:性格决定网站风格
查看>>
Redis 数据结构与内存管理策略(上)
查看>>
Management Console 工具管理类软件通用开发框架(开放源码)
查看>>
Gnome 3.2 发布计划及新功能
查看>>