BZOJ3165: [Heoi2013]Segment(李超线段树)
题意
Sol
李超线段树板子题。具体原理就不讲了。
一开始自己yy着写差点写自闭都快把叉积搬出来了。。。
后来看了下litble的写法才发现原来可以写的这么清晰简洁Orz
#include<bits/stdc++.h>
#define pdd pair<double, double>
#define MP make_pair
#define fi first
#define se second
using namespace std;
const int MAXN = 1e6 + 10, Lim = 1e9;
template <typename A, typename B> inline bool chmin(A &a, B b){if(a > b) {a = b; return 1;} return 0;}
template <typename A, typename B> inline bool chmax(A &a, B b){if(a < b) {a = b; return 1;} return 0;}
inline int read() {
char c = getchar(); int x = 0, f = 1;
while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();}
while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
return x * f;
}
int N = 39989, M;
int ls[MAXN], rs[MAXN], root, cnt, tot;
pdd mx[MAXN];
struct Line {
double k, b;
int id;
}s[MAXN];
pdd get(int x0, int y0, int x1, int y1) {
double k = (double) (y1 - y0) / (x1 - x0),
b = (double) y0 - k * x0;
return {k, b};
}
double calc(Line line, int x) {
return line.k * x + line.b;
}
double GetPoint(Line a, Line b) {
return (b.b - a.b) / (a.k - b.k);
}
pdd ret;
void Query(int k, int l, int r, int p) {//fi: val se: id
if(chmax(ret.fi, calc(s[k], p))) ret.se = s[k].id;
if(l == r) return ;
int mid = l + r >> 1;
if(p <= mid) Query(ls[k], l, mid, p);
else Query(rs[k], mid + 1, r, p);
}
void Modify(int &k, int l, int r, int ql, int qr, Line seg) {
if(!k) k = ++tot;
int mid = l + r >> 1;
if(ql <= l && r <= qr) {
if(!s[k].id) {s[k] = seg; return ;}
int p = GetPoint(s[k], seg);
int pl = calc(s[k], l), pr = calc(s[k], r), nl = calc(seg, l), nr = calc(seg, r);
if(pl > nl && pr > nr) return ;
if(pl < nl && pr < nr) {s[k] = seg; return ;}
if(pl < nl) {
if(p > mid) Modify(rs[k], mid + 1, r, mid + 1, r, s[k]), s[k] = seg;
else Modify(ls[k], l, mid, l, mid, seg);
} else {
if(p > mid) Modify(rs[k], mid + 1, r, mid + 1, r, seg);
else Modify(ls[k], l, mid, l, mid, s[k]), s[k] = seg;
}
return ;
}
if(l == r) return ;
if(ql <= mid) Modify(ls[k], l, mid, ql, qr, seg);
if(qr > mid) Modify(rs[k], mid + 1, r, ql, qr, seg);
}
signed main() {
M = read();
for(int i = 1, lastans = 0; i <= M; i++) {
int opt = read();
if(!opt) {
int k = read(), x = (k + lastans - 1) % 39989 + 1;
ret.fi = 0; ret.se = 0;
Query(root, 1, N, x);
printf("%d\n", lastans = (mx[x].fi > ret.fi ? mx[x].se : ret.se));
} else {
int x0 = (read() + lastans - 1) % 39989 + 1, y0 = (read() + lastans - 1) % Lim + 1,
x1 = (read() + lastans - 1) % 39989 + 1, y1 = (read() + lastans - 1) % Lim + 1;
if(x0 > x1) swap(x0, x1), swap(y0, y1);
if(x0 == x1 && chmax(mx[x0].fi, max(y0, y1))) mx[x0].se = i;
pdd li = get(x0, y0, x1, y1);
Modify(root, 1, N, x0, x1, {li.fi, li.se, ++cnt});
}
}
return 0;
}
BZOJ3165: [Heoi2013]Segment(李超线段树)的更多相关文章
- BZOJ3165[Heoi2013]Segment——李超线段树
题目描述 要求在平面直角坐标系下维护两个操作: 1.在平面上加入一条线段.记第i条被插入的线段的标号为i. 2.给定一个数k,询问与直线 x = k相交的线段中,交点最靠上的线段的编号. 输入 第一行 ...
- 【BZOJ-3165】Segment 李超线段树(标记永久化)
3165: [Heoi2013]Segment Time Limit: 40 Sec Memory Limit: 256 MBSubmit: 368 Solved: 148[Submit][Sta ...
- Luogu P4097 [HEOI2013]Segment 李超线段树
题目链接 \(Click\) \(Here\) 李超线段树的模板.但是因为我实在太\(Naive\)了,想象不到实现方法. 看代码就能懂的东西,放在这里用于复习. #include <bits/ ...
- 2019.02.11 bzoj3165: [Heoi2013]Segment(线段树)
传送门 题意简述:要求支持两种操作: 插入一条线段. 询问与直线x=kx=kx=k相交的线段中,交点最靠上的线段的编号. 思路: 直接上李超线段树即可. 代码: #include<bits/st ...
- 【BZOJ 3165】 [Heoi2013]Segment 李超线段树
所谓李超线段树就是解决此题一类的问题(线段覆盖查询点最大(小)),把原本计算几何的题目变成了简单的线段树,巧妙地结合了线段树的标记永久化与标记下传,在不考虑精度误差的影响下,打法应该是这样的. #in ...
- P4097 [HEOI2013]Segment 李超线段树
$ \color{#0066ff}{ 题目描述 }$ 要求在平面直角坐标系下维护两个操作: 在平面上加入一条线段.记第 i 条被插入的线段的标号为 i 给定一个数 k,询问与直线 x = k 相交的线 ...
- BZOJ.3165.[HEOI2013]Segment(李超线段树)
BZOJ 洛谷 对于线段,依旧是存斜率即可. 表示精度误差一点都不需要管啊/托腮 就我一个人看成了mod(10^9+1)吗.. //4248kb 892ms #include <cstdio&g ...
- Segment 李超线段树
题目大意: 要求在平面直角坐标系下维护两个操作: 1.在平面上加入一条线段.记第 i 条被插入的线段的标号为 i 2.给定一个数 k,询问与直线 x = k 相交的线段中,交点最靠上的线段的编号. 若 ...
- 【洛谷P4097】Segment 李超线段树
题目大意:维护一个二维平面,给定若干条线段,支持询问任意整数横坐标处对应的纵坐标最靠上的线段的 id,相同高度取 id 值较小的,强制在线. 题解:初步学习了李超线段树.李超线段树的核心思想在于通过标 ...
随机推荐
- Swift5 语言参考(三) 类型
在Swift中,有两种类型:命名类型和复合类型.一个名为类型是当它的定义可以给出一个特定名称的类型.命名类型包括类,结构,枚举和协议.例如,名为的用户定义类的实例MyClass具有该类型MyClass ...
- [Leetcode]下一个更大元素II
题目 给定一个循环数组(最后一个元素的下一个元素是数组的第一个元素),输出每个元素的下一个更大元素.数字 x 的下一个更大的元素是按数组遍历顺序,这个数字之后的第一个比它更大的数,这意味着你应该循环地 ...
- POJ 2707
#include<iostream> #include<stdio.h> #include<algorithm> using namespace std; int ...
- 课程一(Neural Networks and Deep Learning),第三周(Shallow neural networks)—— 1、两层神经网络的单样本向量化表示与多样本向量化表示
如上图所示的两层神经网络, 单样本向量化: ...
- crontab命令使用文档.txt
基本格式 : * * * * * command 分 时 日 月 周 命令 第1列表示分钟1-59 每分钟用*或者 */1表示 第2列表示小时1-23(0表示0点) 第3列表示日期1-31 第4列表示 ...
- Java模式—静态代理模式
静态代理模式(Proxy):为其他对象提供一种代理以控制对这个对象的访问,提供“真实对象”的代表,在访问对象时引入一定程度的间接性,这种间接性可以附加多种用途. 代理模式的主要作用是为其他对象 ...
- 在商城系统中使用设计模式----策略模式之在spring中使用策略模式
1.前言: 这是策略模式在spring中的使用,对策略模式不了解对同学可以移步在商城中简单对使用策略模式. 2.问题: 在策略模式中,我们创建表示各种策略的对象和一个行为,随着策略对象改变而改变的 c ...
- 安装python-devel开发包
1.概述 有时在安装某些软件的时候,会报错: Error: must have python development packages -devel, python2.-devel, python2. ...
- JavaScript核心基础语法
1 什么是JavaScript? 是一种嵌入在网页中的程序段. 是一种解释型语言,被浏览器解释执行. 由Netscape发明,ECMA(欧洲计算机制造商协会)将其标准化. JavaScript借用了J ...
- JS pop push unshift shift的作用与区别
白话JS中数组方法pop push unshift shift的作用与区别,通过本文,你大概能知道这四种数组方法的基本使用与大致区别. 首先,这四种方法会直接修改数组,请先记住这一点. 我们先把pop ...