题目描述

你要维护一个向量集合,支持以下操作:
1.插入一个向量(x,y)
2.删除插入的第i个向量
3.查询当前集合与(x,y)点积的最大值是多少。如果当前是空集输出0

输入

第一行输入一个整数n,表示操作个数
接下来n行,每行先是一个整数t表示类型,如果t=1,输入向量
(x,y);如果t=2,输入id表示删除第id个向量;否则输入(x,y),查询
与向量(x,y)点积最大值是多少。
保证一个向量只会被删除一次,不会删没有插入过的向量

输出

对于每条t=3的询问,输出一个答案

样例输入

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维护凸包的更多相关文章

  1. BZOJ_4311_向量_线段树按时间分治

    BZOJ_4311_向量_CDQ分治+线段树按时间分治 Description 你要维护一个向量集合,支持以下操作: 1.插入一个向量(x,y) 2.删除插入的第i个向量 3.查询当前集合与(x,y) ...

  2. 【CF576E】Painting Edges 线段树按时间分治+并查集

    [CF576E]Painting Edges 题意:给你一张n个点,m条边的无向图,每条边是k种颜色中的一种,满足所有颜色相同的边内部形成一个二分图.有q个询问,每次询问给出a,b代表将编号为a的边染 ...

  3. 【BZOJ-4184 】 Shallot 线段树按时间分治 + 线性基

    4184: shallot Time Limit: 30 Sec  Memory Limit: 128 MBSubmit: 356  Solved: 180[Submit][Status][Discu ...

  4. BZOJ_4025_二分图_线段树按时间分治+并查集

    BZOJ_4025_二分图_线段树按时间分治+并查集 Description 神犇有一个n个节点的图.因为神犇是神犇,所以在T时间内一些边会出现后消失.神犇要求出每一时间段内这个图是否是二分图.这么简 ...

  5. BZOJ_4184_shallot_线段树按时间分治维护线性基

    BZOJ_4184_shallot_线段树按时间分治维护线性基 Description 小苗去市场上买了一捆小葱苗,她突然一时兴起,于是她在每颗小葱苗上写上一个数字,然后把小葱叫过来玩游戏. 每个时刻 ...

  6. Bipartite Checking CodeForces - 813F (线段树按时间分治)

    大意: 动态添边, 询问是否是二分图. 算是个线段树按时间分治入门题, 并查集维护每个点到根的奇偶性即可. #include <iostream> #include <sstream ...

  7. BZOJ 1492 货币兑换 cdq分治或平衡树维护凸包

    题意:链接 方法:cdq分治或平衡树维护凸包 解析: 这道题我拒绝写平衡树的题解,我仅仅想说splay不要写挂,insert边界条件不要忘.del点的时候不要脑抽d错.有想写平衡树的去看140142或 ...

  8. codevs 1080 线段树练习 CDQ分治

    codevs 1080 线段树练习 http://codevs.cn/problem/1080/  时间限制: 1 s  空间限制: 128000 KB   题目描述 Description 一行N个 ...

  9. BZOJ4012[HNOI2015]开店——树链剖分+可持久化线段树/动态点分治+vector

    题目描述 风见幽香有一个好朋友叫八云紫,她们经常一起看星星看月亮从诗词歌赋谈到 人生哲学.最近她们灵机一动,打算在幻想乡开一家小店来做生意赚点钱.这样的 想法当然非常好啦,但是她们也发现她们面临着一个 ...

随机推荐

  1. 初识Linux的感受与对它的印象——20155328张钰清

    之前从未接触过虚拟机的我,由于这次寒假预备作业,稍稍地认识了一下Linux操作系统. 在自己笔记本上安装Linux操作系统 根据老师提供的<基于VirtualBox虚拟机安装Ubuntu图文教程 ...

  2. CH03 课下作业

    CH03 课下作业 缓冲区溢出漏洞实验 缓冲区溢出攻击:通过往程序的缓冲区写超出其长度的内容,造成缓冲区的溢出,从而破坏程序的堆栈,造成程序崩溃或使程序转而执行其它指令,以达到攻击的目的. 缓冲区溢出 ...

  3. switchsharp

    https://www.switchysharp.com/file/switchysharp-v1.10.4.zip

  4. thinkphp5 开启多语言

    一.配置点击打开链接1.开启语言包功能'lang_switch_on' => true,2.支持的语言列表'lang_list' => ['zh-cn','en-us'],二.语言定义(默 ...

  5. 如何写一个简单的HTTP服务器(重做版)

    最近几天用C++重新写了之前的HTTP服务器,对以前的代码进行改进.新的HTTP服务器采用Reactor模式,有多个线程并且每个线程有一个EventLoop,主程序将任务分发到每个线程,其中采用的是轮 ...

  6. fastdfs+nginx+image_filter安装与生成缩略图

    fastdfs简介 类似google FS的一个轻量级分布式文件系统,纯C实现,支持linux.FreeBSD等UNIX系统: 只能通过API访问,不支持POXIS: 文件不分块存储,上传的文件和OS ...

  7. 修改Config文件

    /// <summary> /// Config文件操作 /// </summary> public class Config { /// <summary> // ...

  8. 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 ...

  9. Numpy入门笔记第三天

    __TITLE__ = "利用Numpy进行历史股价分析" __DATASOURCE__ = "ATAGURU" # CSV文件读取 import numpy ...

  10. JSON.stringify处理对象时的问题

    1. JSON.stringify({entry_key: 'test', entry_detail: undefined}) 结果 为 "{"entry_key": & ...