【bzoj4311】向量 线段树对时间分治+STL-vector维护凸包
题目描述
输入
输出
样例输入
5
1 3 3
1 1 4
3 3 3
2 1
3 3 3
样例输出
18
15
题解
线段树对时间分治+STL-vector维护凸包
首先需要知道向量 $(x_1,y_1)$ 和 $(x_2,y_2)$ 的点积(数量积)为 $x_1x_2+y_1y_2$ 。
问题转化为给出一堆 $x$ 、$y$ 和 $a$ 、$b$ ,求 $ax+by$ 的最大值。
设 $ax+by=c$ ,则整理可得 $y=-\frac abx+\frac cb$ ,要让 $c$ 最大即让截距最大。因此答案一定在上凸壳上取到。
由于有删除操作,因此使用线段树对时间分治,把一个向量出现的时间段分为线段树上的log个。
然后考虑怎么统计答案:由于平衡树维护凸包的时间复杂度时均摊的,因此不能在线段树上进行插入与恢复的操作。考虑到答案所在的每一段互不影响,因此可以对于线段树的每个节点维护凸包,查询时在每一段上进行二分,取最大值即为答案。
因此使用STL-vector,对于线段树的每个节点,维护上凸壳;查询时在其到线段树根节点的每个节点的凸包上二分,每个节点求出最优解后再取最大值即为答案。
这样做的时间复杂度时 $O(n\log ^2n)$ 。
存在一种更优的解法:对每个询问按照 $-\frac ab$ 从大到小排序,这样决策就是按 $x$ 单调不降的了。每个节点维护当前决策位置,统计时不断判断下一个是否比当前的优即可。
时间复杂度 $O(n\log n)$ ,然而跑不过 $O(n\log^2n)$ 什么鬼。。。
#include <cstdio>
#include <vector>
#include <algorithm>
#define N 200010
#define lson l , mid , x << 1
#define rson mid + 1 , r , x << 1 | 1
using namespace std;
typedef long long ll;
struct point
{
ll x , y;
point() {}
point(ll a , ll b) {x = a , y = b;}
point operator-(const point &a)const {return point(x - a.x , y - a.y);}
ll operator*(const point &a)const {return x * a.y - y * a.x;}
};
struct data
{
point p;
int l , r;
}a[N] , b[N] , c[N];
vector<point> v[N << 2];
int pos[N << 2] , val[N] , ta , tb , tc;
ll ans[N];
bool cmp1(const data &a , const data &b)
{
return a.p.x == b.p.x ? a.p.y < b.p.y : a.p.x < b.p.x;
}
bool cmp2(const data &a , const data &b)
{
return a.p * b.p < 0;
}
inline ll calc(point a , point b)
{
return a.x * b.x + a.y * b.y;
}
void insert(int b , int e , point p , int l , int r , int x)
{
if(b <= l && r <= e)
{
while(v[x].size() > 1 && (p - v[x][v[x].size() - 1]) * (p - v[x][v[x].size() - 2]) <= 0) v[x].pop_back();
v[x].push_back(p);
return;
}
int mid = (l + r) >> 1;
if(b <= mid) insert(b , e , p , lson);
if(e > mid) insert(b , e , p , rson);
}
ll query(int k , point p , int l , int r , int x)
{
ll ans = 0;
if(v[x].size())
{
while(pos[x] < (int)v[x].size() - 1 && calc(p , v[x][pos[x] + 1]) >= calc(p , v[x][pos[x]])) pos[x] ++ ;
ans = calc(p , v[x][pos[x]]);
}
if(l == r) return ans;
int mid = (l + r) >> 1;
if(k <= mid) return max(ans , query(k , p , lson));
else return max(ans , query(k , p , rson));
}
int main()
{
int n , i , opt , k;
scanf("%d" , &n);
for(i = 1 ; i <= n ; i ++ )
{
scanf("%d" , &opt);
if(opt == 1) scanf("%lld%lld" , &a[i].p.x , &a[i].p.y) , a[i].l = i , a[i].r = n , val[++ta] = i;
else if(opt == 2) scanf("%d" , &k) , a[val[k]].r = i - 1 , a[i].r = -1;
else tc ++ , scanf("%lld%lld" , &c[tc].p.x , &c[tc].p.y) , c[tc].l = i;
}
for(i = 1 ; i <= n ; i ++ )
if(a[i].l)
b[++tb] = a[i];
sort(b + 1 , b + tb + 1 , cmp1);
for(i = 1 ; i <= tb ; i ++ ) insert(b[i].l , b[i].r , b[i].p , 1 , n , 1);
sort(c + 1 , c + tc + 1 , cmp2);
for(i = 1 ; i <= tc ; i ++ ) ans[c[i].l] = query(c[i].l , c[i].p , 1 , n , 1);
for(i = 1 ; i <= n ; i ++ )
if(!a[i].r)
printf("%lld\n" , ans[i]);
return 0;
}
【bzoj4311】向量 线段树对时间分治+STL-vector维护凸包的更多相关文章
- BZOJ_4311_向量_线段树按时间分治
BZOJ_4311_向量_CDQ分治+线段树按时间分治 Description 你要维护一个向量集合,支持以下操作: 1.插入一个向量(x,y) 2.删除插入的第i个向量 3.查询当前集合与(x,y) ...
- 【CF576E】Painting Edges 线段树按时间分治+并查集
[CF576E]Painting Edges 题意:给你一张n个点,m条边的无向图,每条边是k种颜色中的一种,满足所有颜色相同的边内部形成一个二分图.有q个询问,每次询问给出a,b代表将编号为a的边染 ...
- 【BZOJ-4184 】 Shallot 线段树按时间分治 + 线性基
4184: shallot Time Limit: 30 Sec Memory Limit: 128 MBSubmit: 356 Solved: 180[Submit][Status][Discu ...
- BZOJ_4025_二分图_线段树按时间分治+并查集
BZOJ_4025_二分图_线段树按时间分治+并查集 Description 神犇有一个n个节点的图.因为神犇是神犇,所以在T时间内一些边会出现后消失.神犇要求出每一时间段内这个图是否是二分图.这么简 ...
- BZOJ_4184_shallot_线段树按时间分治维护线性基
BZOJ_4184_shallot_线段树按时间分治维护线性基 Description 小苗去市场上买了一捆小葱苗,她突然一时兴起,于是她在每颗小葱苗上写上一个数字,然后把小葱叫过来玩游戏. 每个时刻 ...
- Bipartite Checking CodeForces - 813F (线段树按时间分治)
大意: 动态添边, 询问是否是二分图. 算是个线段树按时间分治入门题, 并查集维护每个点到根的奇偶性即可. #include <iostream> #include <sstream ...
- BZOJ 1492 货币兑换 cdq分治或平衡树维护凸包
题意:链接 方法:cdq分治或平衡树维护凸包 解析: 这道题我拒绝写平衡树的题解,我仅仅想说splay不要写挂,insert边界条件不要忘.del点的时候不要脑抽d错.有想写平衡树的去看140142或 ...
- codevs 1080 线段树练习 CDQ分治
codevs 1080 线段树练习 http://codevs.cn/problem/1080/ 时间限制: 1 s 空间限制: 128000 KB 题目描述 Description 一行N个 ...
- BZOJ4012[HNOI2015]开店——树链剖分+可持久化线段树/动态点分治+vector
题目描述 风见幽香有一个好朋友叫八云紫,她们经常一起看星星看月亮从诗词歌赋谈到 人生哲学.最近她们灵机一动,打算在幻想乡开一家小店来做生意赚点钱.这样的 想法当然非常好啦,但是她们也发现她们面临着一个 ...
随机推荐
- 初识Linux的感受与对它的印象——20155328张钰清
之前从未接触过虚拟机的我,由于这次寒假预备作业,稍稍地认识了一下Linux操作系统. 在自己笔记本上安装Linux操作系统 根据老师提供的<基于VirtualBox虚拟机安装Ubuntu图文教程 ...
- CH03 课下作业
CH03 课下作业 缓冲区溢出漏洞实验 缓冲区溢出攻击:通过往程序的缓冲区写超出其长度的内容,造成缓冲区的溢出,从而破坏程序的堆栈,造成程序崩溃或使程序转而执行其它指令,以达到攻击的目的. 缓冲区溢出 ...
- switchsharp
https://www.switchysharp.com/file/switchysharp-v1.10.4.zip
- thinkphp5 开启多语言
一.配置点击打开链接1.开启语言包功能'lang_switch_on' => true,2.支持的语言列表'lang_list' => ['zh-cn','en-us'],二.语言定义(默 ...
- 如何写一个简单的HTTP服务器(重做版)
最近几天用C++重新写了之前的HTTP服务器,对以前的代码进行改进.新的HTTP服务器采用Reactor模式,有多个线程并且每个线程有一个EventLoop,主程序将任务分发到每个线程,其中采用的是轮 ...
- fastdfs+nginx+image_filter安装与生成缩略图
fastdfs简介 类似google FS的一个轻量级分布式文件系统,纯C实现,支持linux.FreeBSD等UNIX系统: 只能通过API访问,不支持POXIS: 文件不分块存储,上传的文件和OS ...
- 修改Config文件
/// <summary> /// Config文件操作 /// </summary> public class Config { /// <summary> // ...
- cmake-index-3.11.4机翻
index next | CMake » git-stage git-master latest release 3.13 3.12 3.11.4 3.10 3.9 3.8 3.7 3.6 3.5 3 ...
- Numpy入门笔记第三天
__TITLE__ = "利用Numpy进行历史股价分析" __DATASOURCE__ = "ATAGURU" # CSV文件读取 import numpy ...
- JSON.stringify处理对象时的问题
1. JSON.stringify({entry_key: 'test', entry_detail: undefined}) 结果 为 "{"entry_key": & ...