P3309-[SDOI2014]向量集【线段树,凸壳】

正题

题目链接:https://www.luogu.com.cn/problem/P3309


题目大意

nnn个操作

  1. 在序列末尾加入一个向量(x,y)(x,y)(x,y)
  2. 询问加入的第l∼rl\sim rlr个向量中的一个向量和(x,y)(x,y)(x,y)的点积最大值

强制在线,点积的定义为x1x2+y1y2x_1x_2+y_1y_2x1x2+y1y2


解题思路

如果对于一个(x,y)(x,y)(x,y)对于两个(x1,y1)(x_1,y_1)(x1,y1)(x2,y2)(x_2,y_2)(x2,y2)如果后者更大那么有
x2x+y2y>x1x+y1y⇒yx≤x2−x1y2−y1x_2x+y_2y>x_1x+y_1y\Rightarrow \frac{y}{x}\leq \frac{x_2-x_1}{y_2-y_1}x2x+y2y>x1x+y1yxyy2y1x2x1
好像和斜率有关,可以维护凸壳来做,因为yyy可能是负数,如果是负数的时候就要求的是下凸壳了,所以两个凸壳都要维护。

因为强制在线所以上不了传统艺能CDQ\text{CDQ}CDQ

那怎么动态维护区间凸壳,平衡树支持动态插入但不支持区间问题。所以考虑线段树,因为一个位置修改了之后就不会再修改,而且是从左往右加的,可以利用这个性质。

每次我们修改一个位置后,如果一个区间[L,R][L,R][L,R]的节点内已经插入了R−L+1R-L+1RL+1个向量(也就是都插完了)的话就直接把它的两个儿子的凸壳合并起来。

然后询问的时候分成lognlog\ nlog n个区间询问的答案取最大值就好了。

合并凸壳的是用归并排序的话时间复杂度O(nlog⁡n)O(n\log n)O(nlogn)


code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#define ll long long
using namespace std;
struct point{ll x,y;point(ll xx=0,ll yy=0){x=xx;y=yy;return;}
}z;
point operator+(point x,point y)
{return point(x.x+y.x,x.y+y.y);}
point operator-(point x,point y)
{return point(x.x-y.x,x.y-y.y);}
ll operator^(point x,point y)
{return x.x*y.y-x.y*y.x;}
ll operator*(point x,point y)
{return x.x*y.x+x.y*y.y;}
bool operator<(point x,point y)
{return (x.x==y.x)?x.y<y.y:x.x<y.x;}const ll N=4e5+10;
char pe[3];
ll n,num,siz[N<<2];
vector<point> v[N<<2][2],tmp;void Make(ll x){ll ls=x*2,rs=x*2+1;for(ll k=0;k<2;k++){ll i=0,j=0,l1=v[ls][k].size()-1,l2=v[rs][k].size()-1;tmp.clear();while(i<=l1||j<=l2){if(i>l1||(j<=l2&&v[rs][k][j]<v[ls][k][i]))tmp.push_back(v[rs][k][j]),j++;else tmp.push_back(v[ls][k][i]),i++;}ll cnt=0;for(ll i=0;i<tmp.size();i++){while(cnt>1&&((v[x][k][cnt-1]-v[x][k][cnt-2])^(tmp[i]-v[x][k][cnt-1]))>=0)v[x][k].pop_back(),cnt--;v[x][k].push_back(tmp[i]);cnt++;}}return;
}
ll Calc(ll x,point p){ll f=0;if(p.y<0)p=z-p,f^=1;ll l=0,r=v[x][f].size()-2;while(l<=r){ll mid=(l+r)>>1;point tmp=v[x][f][mid+1]-v[x][f][mid];tmp.x*=-1;if(p.x*tmp.x>=p.y*tmp.y)r=mid-1;else l=mid+1;}return p*v[x][f][l];
}
void Change(ll x,ll L,ll R,ll pos,point p){if(L==R){v[x][0].push_back(p);v[x][1].push_back(z-p);return;}ll mid=(L+R)>>1;siz[x]++;if(pos<=mid)Change(x*2,L,mid,pos,p);else Change(x*2+1,mid+1,R,pos,p);if(siz[x]==R-L+1)Make(x);
}
ll Ask(ll x,ll L,ll R,ll l,ll r,point p){if(L==l&&R==r)return Calc(x,p);ll mid=(L+R)>>1;if(r<=mid)return Ask(x*2,L,mid,l,r,p);if(l>mid)return Ask(x*2+1,mid+1,R,l,r,p);return max(Ask(x*2,L,mid,l,mid,p),Ask(x*2+1,mid+1,R,mid+1,r,p));
}
void dc(ll &x,ll lastans) {if(pe[0]=='E')return;x=x^(lastans&0x7fffffff);return;
}
signed main()
{scanf("%lld%s",&n,pe);ll last=0;for(ll i=1;i<=n;i++){char op[3];ll x,y,l,r;scanf("%s%lld%lld",op,&x,&y);dc(x,last);dc(y,last);if(op[0]=='A'){++num;Change(1,1,n,num,point(x,y));}else{scanf("%lld%lld",&l,&r);dc(l,last);dc(r,last);printf("%lld\n",last=Ask(1,1,n,l,r,point(x,y)));}}return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/319723.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

CanalSharp-mysql数据库binlog的增量订阅消费组件Canal的.NET客户端

一.前言CanalSharp是阿里巴巴开源项目mysql数据库binlog的增量订阅&消费组件 Canal 的.NET客户端&#xff0c;关于什么是 Canal&#xff1f;又能做什么&#xff1f;我会在后文为大家一一介绍。CanalSharp 这个项目&#xff0c;是由我和 WithLin(主要贡献) 完成&#xff0c;…

Codeforces Round #657 (Div. 2)

A. Acacius and String 爆零&#xff01;太菜了&#xff0c;下来终于把A题代码调AC了 #define IO ios::sync_with_stdio(false);cin.tie();cout.tie(0) #include<iostream> #include<algorithm> #include<cstring> #include<string> using namespace…

【树链剖分】软件管理(luogu 2146/金牌导航 树链剖分-2)

正题 luogu 2146 金牌导航 树链剖分-2 题目大意 有若干软件&#xff0c;除了软件0&#xff0c;所有软件都依赖且只依赖于另外一个软件 当要删除一个软件时&#xff0c;所有依赖于该软件的软件都要删掉 当安装一个软件时&#xff0c;该软件依赖的软件都要安装 问你每次操作…

【每日一题】7月14日题目精讲—压缩

来源&#xff1a;牛客网&#xff1a; 时间限制&#xff1a;C/C 1秒&#xff0c;其他语言2秒 空间限制&#xff1a;C/C 262144K&#xff0c;其他语言524288K 64bit IO Format: %lld题目描述 给一个由小写字母组成的字符串&#xff0c;我们可以用一种简单的方法来压缩其中的重复…

P3306-[SDOI2013]随机数生成器【BSGS】

正题 题目链接:https://www.luogu.com.cn/problem/P3306 题目大意 给出一个p,a,b,x1,tp,a,b,x_1,tp,a,b,x1​,t&#xff0c;有xiaxi−1bx_iax_{i-1}bxi​axi−1​b 求一个最小的nnn使得xntx_ntxn​t 解题思路 下标缩一下先变成x0x_0x0​会更好算一点&#xff0c;只考虑x0x_0…

【DevOps+LIVE】直播复盘 – DevOps能力成长模型2018首发

昨天晚上小编和DevOpsDays中国核心组织者刘征老师和张乐老师一起嗨了一把&#xff0c;本来规划进行1小时的直播结果延长到了2个小时才结束&#xff0c;观众人数在前半段一直处于增长状态&#xff0c;直到开播一小时到达顶峰。小伙伴也非常给力&#xff0c;一共提出了40多个问题…

Codeforces Round #658 (Div. 2)

A - Common Subsequence 最短相同子序列长度肯定为1&#xff0c;如果一个元素都不相等之间不存在相同子序列 #define IO ios::sync_with_stdio(false);cin.tie();cout.tie(0) #define debug(x) cout<<#x<<": "<<x<<" " #include…

【树链剖分】染色(luogu 2486/金牌导航 树链剖分-3)

正题 luogu 2486 金牌导航 树链剖分-3 题目大意 给你一棵树&#xff0c;让你进行以下操作&#xff1a; 1.把一条路径染上一个颜色 2.查询一条路径上有多少个颜色段 解题思路 用树链剖分把问题转化为链上问题 然后维护一下左右端点颜色和颜色总数就好了 代码 #include<c…

小小粉刷匠

来源&#xff1a;牛客网&#xff1a; 时间限制&#xff1a;C/C 1秒&#xff0c;其他语言2秒 空间限制&#xff1a;C/C 32768K&#xff0c;其他语言65536K 64bit IO Format: %lld题目描述 “lalala,我是一个快乐的粉刷匠”,小名一边快活地唱着歌,一边开心地刷着墙",兴致突…

P5363-[SDOI2019]移动金币【阶梯博弈,dp,组合数学】

正题 题目链接:https://www.luogu.com.cn/problem/P5363 题目大意 1n1\times n1n的网格上有mmm个硬币&#xff0c;两个人轮流向前移动一个硬币但是不能超过前一个硬币&#xff0c;无法移动者输。 求有多少种情况先手必胜。 解题思路 竟然有我会的题&#xff0c;我感动 位置做…

M-SOLUTIONS Programming Contest 2020总结

A - Kyu in AtCoder 直接模拟 #define IO ios::sync_with_stdio(false);cin.tie();cout.tie(0) #define debug(x) cout<<#x<<": "<<x<<" " #include<iostream> #include<algorithm> using namespace std; int main(…

Emit动态代理.NetCore迁移之旅

前言】前面我们介绍了Aop 从静态代理到动态代理&#xff1a;https://www.cnblogs.com/7tiny/p/9657451.html我们在.NetFramework平台下使用微软提供的Emit技术实现了动态代理类的生成。但是.NetCore作为微软.Net平台的春天&#xff0c;如果类库光支持.NetFramework&#xff0c;…

【LCT】弹飞绵羊(luogu 3203/金牌导航 LCT-2)

正题 luogu 3203 金牌导航 LCT-2 题目大意 给你n个格子&#xff0c;当你在第i个格子时&#xff0c;可以往后跳aia_iai​格&#xff0c;让你进行几下操作&#xff1a; 1.修改第i个数 2.查询在第i个格子跳多少下会跳出界 解题思路 往后跳相当于连接格子&#xff0c;由此建立一…

【每日一题】7月15日题目精讲—生日快乐

【每日一题】7月15日题目精讲—生日快乐 [SCOI2009]生日快乐 时间限制&#xff1a;C/C 1秒&#xff0c;其他语言2秒 空间限制&#xff1a;C/C 262144K&#xff0c;其他语言524288K 64bit IO Format: %lld题目描述 windy的生日到了&#xff0c;为了庆祝生日&#xff0c;他的朋友…

P1712-[NOI2016]区间【线段树,尺取法】

正题 题目链接:https://www.luogu.com.cn/problem/P1712 题目大意 nnn个区间&#xff0c;求出其中mmm个区间使得它们有覆盖同一个点且最长区间长度减去最短长度最小。 解题思路 因为是最接近的mmm个&#xff0c;考虑一种叫尺取法的做法。 先把区间按照长度排序&#xff0c;每…

Codeforces Round #659 (Div. 2)

日常爆零掉分wa A. Common Prefixes 根据前一个字符串构造后一个字符串&#xff0c;从哪不同就从哪换 #define IO ios::sync_with_stdio(false);cin.tie();cout.tie(0) #define debug(x) cout<<#x<<": "<<x<<" " #include<i…

asp.net core结合Gitlab-CI实现自动化部署

一、前言在之前的文章中写过k8sJenkinsGitLab-自动化部署asp.net core项目 的topic&#xff0c;这次讲解一下gitlab-ci的CI/CD&#xff0c;说实话&#xff0c;自动化部署是在是非常的舒服&#xff0c;而且gitlab-ci很适合中小项目的部署&#xff0c;大项目推荐使用k8s。阅读目录…

【LCT】城市旅行(luogu 4842/金牌导航 LCT-3)

正题 luogu 4842 金牌导航 LCT-3 题目大意 给你一棵树&#xff0c;让你进行一些操作&#xff1a; 1.删除一条边 2.连接一条边 3.给一条路径上的点加上x 4.给出一条路径&#xff0c;在该路径选取两个点&#xff0c;求这两个点之间路径的权值和的期望值 解题思路 该树可以用LC…

P6329-[模板]点分树 | 震波

正题 题目链接:https://www.luogu.com.cn/problem/P6329 解题思路 给出nnn个点的一棵树&#xff0c;每个点有权值&#xff0c;有mmm次操作 修改一个点xxx的权值为yyy询问距离点xxx不超过kkk的所有点点权和 解题思路 点分树的模板题&#xff0c;先点分治构造出点分树&#x…

Codeforces Round #660 (Div. 2)

A - Captain Flint and Crew Recruitment 刚开始还想筛法求质数&#xff0c;最后发现是个脑筋急转弯 #define IO ios::sync_with_stdio(false);cin.tie();cout.tie(0) #define debug(x) cout<<#x<<": "<<x<<" " #include<ios…