题面

题解

李超线段树

为了与机房大佬 HYJ 同步伐

学习笔记请移步 yyb的博客

Code

#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdio>
#include <cmath>
const int N = 100005;
const int mod = 39989;
#define val(id, x) (k[id] * x + b[id])
using namespace std; int Q, cnt, ans;
struct node
{
int id;
double k, b;
} t[N << 2];
double k[N], b[N]; template < typename T >
inline T read()
{
T x = 0, w = 1; char c = getchar();
while(c < '0' || c > '9') { if(c == '-') w = -1; c = getchar(); }
while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
return x * w;
} void pushtag(int p, int l, int r, int id)
{
if(!t[p].id) { t[p] = (node) { id, k[id], b[id] }; return; }
int now = t[p].id;
double l1 = val(now, l), r1 = val(now, r), l2 = val(id, l), r2 = val(id, r);
if(l2 <= l1 && r2 <= r1) return;
if(l1 < l2 && r1 < r2) { t[p] = (node) { id, k[id], b[id] }; return; }
int mid = (l + r) >> 1; double x = (b[id] - b[now]) / (k[now] - k[id]);
if(x <= mid)
{
if(l1 <= l2) pushtag(p << 1, l, mid, id);
else pushtag(p << 1, l, mid, now), t[p] = (node) { id, k[id], b[id] };
}
else
{
if(l1 <= l2) pushtag(p << 1 | 1, mid + 1, r, now), t[p] = (node) { id, k[id], b[id] };
else pushtag(p << 1 | 1, mid + 1, r, id);
}
} void modify(int p, int l, int r, int ql, int qr, int id)
{
if(ql <= l && r <= qr) { pushtag(p, l, r, id); return; }
int mid = (l + r) >> 1;
if(ql <= mid) modify(p << 1, l, mid, ql, qr, id);
if(mid < qr) modify(p << 1 | 1, mid + 1, r, ql, qr, id);
} void check(int &a, int c, int x)
{
double y0 = val(a, x), y1 = val(c, x);
if(y1 > y0 || (fabs(y1 - y0) < 1e-7 && a > c)) a = c;
} int query(int p, int l, int r, int x)
{
if(l == r) return t[p].id;
int mid = (l + r) >> 1, res = t[p].id;
if(x <= mid) check(res, query(p << 1, l, mid, x), x);
else check(res, query(p << 1 | 1, mid + 1, r, x), x);
return res;
} int main()
{
/* freopen("cpp.in", "r", stdin);
freopen("cpp.out", "w", stdout);
*/ Q = read <int> ();
int opt, x0, x1, y0, y1;
while(Q--)
{
opt = read <int> ();
if(!opt)
{
// printf("%d\n", x0);
x0 = read <int> (), x0 = ((x0 + ans - 1) % mod + 1);
printf("%d\n", ans = query(1, 1, mod, x0));
}
else
{
x0 = read <int> (), y0 = read <int> (), x1 = read <int> (), y1 = read <int> ();
x0 = ((x0 + ans - 1) % mod + 1), y0 = ((y0 + ans - 1) % 1000000000 + 1);
x1 = ((x1 + ans - 1) % mod + 1), y1 = ((y1 + ans - 1) % 1000000000 + 1);
if(x0 > x1) swap(x0, x1), swap(y0, y1);
// printf("%d %d\n", x0, x1);
k[++cnt] = 1.0 * (y1 - y0) / (x1 - x0), b[cnt] = 1.0 * y0 - k[cnt] * x0;
modify(1, 1, mod, x0, x1, cnt);
}
}
return 0;
}
/*
3
1 3 5 7 7
1 3 9 7 3
0 5
*/

[题解] [HEOI2013] Segment的更多相关文章

  1. Bzoj 3165 [Heoi2013]Segment题解

    3165: [Heoi2013]Segment Time Limit: 40 Sec  Memory Limit: 256 MBSubmit: 668  Solved: 276[Submit][Sta ...

  2. 【BZOJ3165】[HEOI2013]Segment(李超线段树)

    [BZOJ3165][HEOI2013]Segment(李超线段树) 题面 BZOJ 洛谷 题解 似乎还是模板题QwQ #include<iostream> #include<cst ...

  3. bzoj 3165: [Heoi2013]Segment 动态凸壳

    3165: [Heoi2013]Segment Time Limit: 40 Sec  Memory Limit: 256 MBSubmit: 202  Solved: 89[Submit][Stat ...

  4. 洛谷 P4097 [HEOI2013]Segment 解题报告

    P4097 [HEOI2013]Segment 题目描述 要求在平面直角坐标系下维护两个操作: 在平面上加入一条线段.记第 \(i\) 条被插入的线段的标号为 \(i\) 给定一个数 \(k\),询问 ...

  5. BZOJ 3165: [Heoi2013]Segment

    3165: [Heoi2013]Segment Time Limit: 40 Sec  Memory Limit: 256 MBSubmit: 465  Solved: 187[Submit][Sta ...

  6. 【题解】Luogu P4097 [HEOI2013]Segment

    原题传送门 这珂以说是李超线段树的模板题 按着题意写就行了,时间复杂度为\(O(n\log^2n)\) #include <bits/stdc++.h> #define N 40005 # ...

  7. BZOJ3165 & 洛谷4097:[HEOI2013]Segment——题解

    https://www.lydsy.com/JudgeOnline/problem.php?id=3165 https://www.luogu.org/problemnew/show/P4097 要求 ...

  8. HEOI2013 Segment

    传说中的“李超树”. 大意:给你若干线段,试求横坐标x上的最上方一条线段的编号.无则输出零. 解:用线段树维护. 插入的时候保存自己这个区间上可能成为最大值的线段,被抛弃的则看情况下放. 查询时从最底 ...

  9. 洛谷P4097 [HEOI2013]Segment(李超线段树)

    题面 传送门 题解 调得咱自闭了-- 不难发现这就是个李超线段树,不过因为这里加入的是线段而不是直线,所以得把线段在线段树上对应区间内拆开之后再执行李超线段树的操作,那么复杂度就是\(O(n\log^ ...

随机推荐

  1. C#特性 详解

    一:Conditional:条件特性,预定义了一个条件方法. 使用方法: [Conditional("DEBUG")] public void test() { MessageBo ...

  2. JS中this和call

    首先来了解一下JS中this的原理: 要访问自己的属性就必须使用 this.属性名 1.this总是指向它的直接调用者: var a={ user:'Artimis', fn:function(){ ...

  3. XML文件解析之SAX解析

    使用DOM解析的时候是需要把文档的所有内容读入内存然后建立一个DOM树结构,然后通过DOM提供的接口来实现XML文件的解析,如果文件比较小的时候肯定是很方便的.但是如果是XML文件很大的话,那么这种方 ...

  4. python数据写入Excel表格

    from openpyxl import Workbook def main(): sheet_name = "表名1" row_count = 6 # 行数 info_resul ...

  5. SmartBinding实现DataSet与ListView的绑定及同步显示

    kbmMW 5.10.10发布了,这个版本解决了我提出的问题,当对DataSet增删记录时,ListView能够同步显示.下面看看具体的实现代码. 为了解决上面的问题,作者为IkbmMWBinding ...

  6. Xen 虚拟化技术

    Xen 是一种开源的.属于类型1(裸金属虚拟化,Baremetal Hypervisor)的虚拟化技术,它使多个同样操作系统或不同操作系统的虚拟机运行在同一个物理主机节点上成为可能并实现. Xen 是 ...

  7. Centos7下安装Elasticsearch 5.6.6

    环境 因为elasticsearch是用java编写的,所以需要先安装JDK ES 5,安装需要 JDK 8 以上ES 6.5,安装需要 JDK 11 以上ES 7.2.1,内置了 JDK 12 安装 ...

  8. Vue下的package.json

    { "name": "demo2", "version": "1.0.0", "description&quo ...

  9. No result defined for action com.java.test.Action.HelloAction and result index

    Struts中配置action访问出错: Struts Problem Report Struts has detected an unhandled exception: Messages: No ...

  10. python3 random

    一.random 1.生成伪随机数 2.伪随机数是可预测的,严格意义上不具有随机性质,通常用数学公式的方法(比如统计分布,平方取中等)获得 3.正如数列需要有首项,产生伪随机数需要一个初值用来计算整个 ...