题面

题解

加入一条线段,可以把它转化为在[L,R]区间内加一条线 y=ax+b (如果原线段与y轴平行,就相当于在{x1}处加一条线 y=max(y1,y2))

我们可以把它加到线段树上,线段树上每个点存一个区间内贯穿整个区间的一条对答案有贡献的线段(因为是贯穿整个区间,所以存一个 {a,b} 表示 y=ax+b 就行了)。

一条线段是这么加上去的:

查询要遍历路径上的所有点:

接下来,最大的问题就是,它会产生冲突!

当我们加一条线段到线段树上,而原来节点上就已经有了一条线段怎么办?

首先,如果其中一条线段凌驾于另一条之上,肯定直接判了嘛,

但是,出问题的就是两条线段交叉,

而线段树上的节点存不了“V”字形,也维护不了(除了写平衡树的大佬)

这里就要用到一个“启发式下传”的思想

一般的暴力想法就是先把原来的线段先下传了,再把新线段放上去:

稍微有点想法的人会先把“V”字形处理出来,并把两条线段拆开下传(虽然并没有什么*用):

这时候,较聪明的人就会想,为什么原来的地方要空着呢?把其中随便一条线段留在原处,对答案也没有影响,但是却快得多:

最后,深谙暴力优化精髓的人会想,数据肯定会卡每次下传的线段长度只比原先段短一点的情况,便把“V”字形两条线段中短的一条下传,

于是便优化出了正解:

我们会发现,每次一旦冲突,便会开始下传,每次下传的线段长度至少会减半,所以,下传复杂度为O(log),加线段复杂度就为O(log^2)

查询复杂度没变,总复杂度就为O(mlog^2)

CODE

#include<map>
#include<queue>
#include<cmath>
#include<vector>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define MAXN 50005
#define MAXM 100005
#define LL long long
#define DB double
#define ENDL putchar('\n')
#define lowbit(x) ((-x)&(x))
//#define int LL
//#pragma GCC optimize(2)
LL read() {
LL f = 1,x = 0;char s = getchar();
while(s < '0' || s > '9') {if(s=='-')f=-f;s = getchar();}
while(s >= '0' && s <= '9') {x=x*10+(s-'0');s=getchar();}
return f*x;
}
const int MOD = 998244353;
int n,m,i,j,s,o,k;
struct it{
DB x,y;
it(){x = y = 0;}
it(DB X,DB Y){x=X;y=Y;}
};
it cg(it a,it b) {
DB k = (a.y - b.y) / (a.x - b.x);
DB B = a.y - a.x * k;
return it(k,B);
}
it jd(it a,it b) {
DB xx = (b.y - a.y) / (a.x - b.x);
return it(xx,a.x*xx+a.y);
}
bool abov(it a,it st) {
DB y2 = a.x * st.x + st.y;
return a.y >= y2;
}
bool operator < (it a,it b) {return a.y < b.y;}
bool operator > (it a,it b) {return b < a;}
bool operator <= (it a,it b) {return a.y <= b.y;}
bool operator >= (it a,it b) {return b <= a;}
struct tr{
int l,r;
it st;
tr(){l=r=0;st=it(0,-1e9);}
}tre[MAXM<<4];
void maketree(int a,int l,int r) {
tre[a].l = l;tre[a].r = r;
if(l < r) {
int mid = (l+r)>>1;
maketree(a<<1,l,mid);
maketree(a<<1|1,mid+1,r);
}
return ;
}
void addtree(int a,int l,int r,it st) {
if(tre[a].l > r || tre[a].r < l) return ;
if(tre[a].l >= l && tre[a].r <= r) {
DB ll = (DB)tre[a].l,rr = (DB)tre[a].r;
it l1 = it(ll,ll*tre[a].st.x+tre[a].st.y);
it r1 = it(rr,rr*tre[a].st.x+tre[a].st.y);
it l2 = it(ll,ll*st.x+st.y),r2 = it(rr,rr*st.x+st.y);
if(l1 >= l2 && r1 >= r2) return ;
if(l2 >= l1 && r2 >= r1) {
tre[a].st = st; return ;
}
it md = jd(st,tre[a].st); int mid = (tre[a].l+tre[a].r)>>1;
it ls = (l1 >= l2 ? tre[a].st:st),rs = (r1 >= r2 ? tre[a].st:st);
if(md.x <= (DB)mid) tre[a].st = rs,addtree(a<<1,l,r,ls);
else tre[a].st = ls,addtree(a<<1|1,l,r,rs);
return ;
}
addtree(a<<1,l,r,st);addtree(a<<1|1,l,r,st);
return ;
}
DB query(int a,int ad) {
if(tre[a].l > ad || tre[a].r < ad) return -1e9;
if(tre[a].l == tre[a].r) {
return ad*tre[a].st.x + tre[a].st.y;
}
return max(ad*tre[a].st.x+tre[a].st.y,max(query(a<<1,ad),query(a<<1|1,ad)));
}
int main() {
n = read();m = read();
maketree(1,1,100000);
for(int i = 1;i <= n;i ++) {
it a,b;
a.x = (DB)read();
a.y = (DB)read();
b.x = (DB)read();
b.y = (DB)read();
if(a.x > b.x) swap(a,b);
it st;
if(a.x == b.x) st = it(0,max(a.y,b.y));
else st = cg(a,b);
int ll = (int)a.x,rr = (int)b.x;
addtree(1,ll,rr,st);
}
for(int i = 1;i <= m;i ++) {
k = read();
if(!k) {
it a,b;
a.x = (DB)read();
a.y = (DB)read();
b.x = (DB)read();
b.y = (DB)read();
if(a.x > b.x) swap(a,b);
it st;
if(a.x == b.x) st = it(0,max(a.y,b.y));
else st = cg(a,b);
int ll = (int)a.x,rr = (int)b.x;
addtree(1,ll,rr,st);
}
else {
int x0 = read();
DB ans = query(1,x0);
printf("%.4f\n",ans < -1e8 ? 0.0:ans);
}
}
return 0;
}

「雅礼集训 2017 Day2」线段游戏(线段树懒标记“启发式下传”,李超树)的更多相关文章

  1. #6034. 「雅礼集训 2017 Day2」线段游戏 李超树

    #6034. 「雅礼集训 2017 Day2」线段游戏 内存限制:256 MiB时间限制:1000 ms标准输入输出 题目类型:传统评测方式:Special Judge 上传者: 匿名 提交提交记录统 ...

  2. 【loj6034】「雅礼集训 2017 Day2」线段游戏

    #6034. 「雅礼集训 2017 Day2」线段游戏 内存限制:256 MiB 时间限制:1000 ms 标准输入输出 题目类型:传统 评测方式:Special Judge 上传者: 匿名 题目描述 ...

  3. 「雅礼集训 2017 Day2」解题报告

    「雅礼集训 2017 Day2」水箱 我怎么知道这种题目都能构造树形结构. 根据高度构造一棵树,在树上倍增找到最大的小于约束条件高度的隔板,开一个 \(vector\) 记录一下,然后对于每个 \(v ...

  4. 【LOJ6029】「雅礼集训 2017 Day1」市场(线段树裸题)

    点此看题面 大致题意: 维护序列,支持区间加法,区间除法(向下取整),区间求\(min\)和区间求和. 线段树维护区间除法 区间加法.区间求\(min\)和区间求和都是线段树基本操作,因此略过不提. ...

  5. loj #6032. 「雅礼集训 2017 Day2」水箱 线段树优化DP转移

    $ \color{#0066ff}{ 题目描述 }$ 给出一个长度为 \(n\) 宽度为 \(1\) ,高度无限的水箱,有 \(n-1\) 个挡板将其分为 \(n\) 个 \(1 - 1\) 的小格, ...

  6. loj#6034 「雅礼集训 2017 Day2」线段游戏

    分析 区间李超树板子题 代码 #include<bits/stdc++.h> using namespace std; #define db double const int inf = ...

  7. LOJ#6032. 「雅礼集训 2017 Day2」水箱

    传送门 首先可以有一个平方复杂度的 \(DP\) 设 \(f_{i,j}\) 表示前面 \(i\) 个小格,高度为 \(j\) 的最大答案 令 \(h_i\) 表示隔板 \(i\) 的高度 当 \(j ...

  8. LOJ#6048. 「雅礼集训 2017 Day10」数列(线段树)

    题面 传送门 题解 我的做法似乎非常复杂啊-- 首先最长上升子序列长度就等于把它反过来再接到前面求一遍,比方说把\(2134\)变成\(43122134\),实际上变化之后的求一个最长上升子序列和方案 ...

  9. 「雅礼集训 2017 Day2」水箱

    题目链接 题意分析 我们用\(f[i][j]\)表示当前到达第\(i\)个位置水位高度为\(j\)的答案 如果那么\(h[i]\)为\(i\)和\(i+1\)之间的支柱高度 那么如果\(j≤h[i]\ ...

随机推荐

  1. CSP 2021 总结

    CSP 2021 总结 PJ 开题顺序:1342 应该先做 T2 ,导致我 T2 直接看错 T1.T3 T1 :直接推规律即可,考场的想法应该正确 T3 :好家伙直接 map 走起 T2 最崩溃的来了 ...

  2. ExtJS 布局-Absolute布局(Absolute layout)

    更新记录: 2022年5月31日 发布本篇 1.说明 使用xy配置项设置子组件在父容器中绝对位置,本质是将子组件的CSS的position设置为absolute,然后使用x和y配置项映射到CSS的to ...

  3. 『忘了再学』Shell流程控制 — 36、for循环介绍

    目录 1.for循环介绍 2.示例 语法一举例: 语法二举例: 3.for循环总结 4.练习:批量解压缩脚本 方式一:批量解压缩 方式二:批量解压缩 1.for循环介绍 for循环是固定循环,也就是在 ...

  4. SAP Grid control( ALV Grid 列表 自定义 按钮)

    ALV 列表和按钮 效果 源代码 PROGRAM bcalvc_tb_menu_with_def_but. *&&&&&&&&& ...

  5. 一个bug肝一周...忍不住提了issue

    导航 Socket.IO是什么 Socket.IO的应用场景 为什么选socket.io-client-java 实战案例 参考 本文首发于智客工坊-<socket.io客户端向webserve ...

  6. 聊聊 C++ 大一统的初始化运算符 {}

    一:背景 最近发现 C++ 中的类型初始化操作,没有 {} 运算符搞不定的,蛮有意思,今天我们就来逐一列一下各自的用法以及汇编展现,本来想分为 值类型 和 引用类型 两大块,但发现在 C++ 中没这种 ...

  7. linux 文件名乱码的文件无法删除

    1.通过ls -i命令获得文件的节点号 2.通过节点号删除 find -inum 节点号 -delete 这样就可以删除文件名乱码的文件

  8. 前 K 个高频元素问题

    前 K 个高频元素问题 作者:Grey 原文地址: 前 K 个高频元素问题 题目描述 LeetCode 347. Top K Frequent Elements 思路 第一步,针对数组元素封装一个数据 ...

  9. 虚拟机win7系统安装

    win7 x64虚拟机安装步骤 1.点击创建新虚拟机,选择典型 2.选择客户机操作系统:windows 7 x64 3.命名虚拟机 4.指定磁盘容量 5.点击完成 6.编辑此虚拟机设置 7.内存设置 ...

  10. vivado没用上的寄存器变量

    vivado中定义了但没用上的寄存器变量,在综合时会被移除,即没有综合出来.(如下cnt,虽然在y的过程块中用了cnt作为判断条件,但实际上cnt用了跟没用效果一样,所以综合时cnt_reg就被放弃了 ...