BZOJ 3533 sdoi 2014 向量集
设(x,y)为Q的查询点,分类讨论如下:
1、y>0: 最大化a*x+b*y,维护一个上凸壳三分即可
2、y<0:最大化a*x+b*y 维护一个下凸壳三分即可
我们考虑对时间建出一棵线段树
对于每个区间,如果满了就做出两个凸壳
总时间复杂度是O(n*log^2n)
之后我们考虑查询,每个区间最多被分解为log(n)个区间
在每个区间的凸壳上三分最优解即可
至于优化,可以设定一个阈值,当区间长度小于阈值时不用做凸壳,查询时直接暴力就可以了
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
using namespace std; typedef long long LL;
const int maxn=400010;
const LL oo=1LL<<60;
char s,Q;
int n,cnt,a,b,T,TA,TB;
LL ans; struct Point{
int x,y;
Point(int x=0,int y=0):x(x),y(y){}
}p[maxn],tmp[maxn],now,A[4000010],B[4000010];
typedef Point Vector;
Vector operator -(const Point &A,const Point &B){return Vector(A.x-B.x,A.y-B.y);}
Vector operator +(const Point &A,const Point &B){return Vector(A.x+B.x,A.y+B.y);}
inline LL Dot(const Point &A,const Point &B){return 1LL*A.x*B.x+1LL*A.y*B.y;}
inline LL Cross(const Point &A,const Point &B){return 1LL*A.x*B.y-1LL*A.y*B.x;} struct Seg_Tree{
int L,R;
int AL,AR;
int BL,BR;
}t[maxn<<2]; bool cmp1(const Point &a,const Point &b){
if(a.x==b.x)return a.y>b.y;
return a.x<b.x;
}
bool cmp2(const Point &a,const Point &b){
if(a.x==b.x)return a.y<b.y;
return a.x<b.x;
}
inline void decode(int &x){x=(s=='E'?x:x=x^(ans&0x7fffffff));}
inline void Get(char &ch){
ch=getchar();
while(ch<'!')ch=getchar();
}
inline void read(int &num){
num=0;int f=1;char ch;Get(ch);
if(ch=='-')f=-1,ch=getchar();
while(ch>='0'&&ch<='9')num=num*10+ch-'0',ch=getchar();
num*=f;
}
void build(int o,int L,int R){
t[o].L=L;t[o].R=R;
if(L==R)return;
int mid=(L+R)>>1;
build(o<<1,L,mid);
build(o<<1|1,mid+1,R);
}
void add(int o){
int L=t[o].L,R=t[o].R;
int mid=(L+R)>>1;
if(R==cnt&&R-L>=70){
T=0;
for(int i=L;i<=R;++i)tmp[++T]=p[i];
sort(tmp+1,tmp+T+1,cmp1);
A[++TA]=tmp[1];t[o].AL=TA;
for(int i=2;i<=T;++i){
if(tmp[i].x!=tmp[i-1].x){
while(TA>t[o].AL&&Cross(tmp[i]-A[TA],A[TA]-A[TA-1])<=0)TA--;
A[++TA]=tmp[i];
}
}t[o].AR=TA;
sort(tmp+1,tmp+T+1,cmp2);
B[++TB]=tmp[1];t[o].BL=TB;
for(int i=2;i<=T;++i){
if(tmp[i].x!=tmp[i-1].x){
while(TB>t[o].BL&&Cross(tmp[i]-B[TB],B[TB]-B[TB-1])>=0)TB--;
B[++TB]=tmp[i];
}
}t[o].BR=TB;
}
if(L==R)return;
if(cnt<=mid)add(o<<1);
else add(o<<1|1);
}
LL Get_A(int L,int R){
while(R-L>=3){
int m1=(L+L+R)/3,m2=(R+R+L)/3;
if(Dot(now,A[m1])>Dot(now,A[m2]))R=m2;
else L=m1;
}
LL ans=-oo;
for(int i=L;i<=R;++i)ans=max(ans,Dot(now,A[i]));
return ans;
}
LL Get_B(int L,int R){
while(R-L>=3){
int m1=(L+L+R)/3,m2=(R+R+L)/3;
if(Dot(now,B[m1])>Dot(now,B[m2]))R=m2;
else L=m1;
}
LL ans=-oo;
for(int i=L;i<=R;++i)ans=max(ans,Dot(now,B[i]));
return ans;
}
LL Get_ask(int o){
LL ans=-oo;
int L=t[o].L,R=t[o].R;
int mid=(L+R)>>1;
if(a<=L&&R<=b){
if(R-L<70)for(int i=L;i<=R;++i)ans=max(ans,Dot(now,p[i]));
else if(now.y>0)ans=Get_A(t[o].AL,t[o].AR);
else ans=Get_B(t[o].BL,t[o].BR);
return ans;
}
if(b<=mid)return Get_ask(o<<1);
else if(a>mid)return Get_ask(o<<1|1);
else return max(Get_ask(o<<1),Get_ask(o<<1|1));
} int main(){
read(n);Get(s);build(1,1,n);
for(int i=1;i<=n;++i){
Get(Q);
if(Q=='A'){
++cnt;
read(p[cnt].x);read(p[cnt].y);
decode(p[cnt].x);decode(p[cnt].y);
add(1);
}else{
read(now.x);read(now.y);
decode(now.x);decode(now.y);
read(a);read(b);
decode(a);decode(b);
ans=Get_ask(1);
printf("%lld\n",ans);
}
}return 0;
}
BZOJ 3533 sdoi 2014 向量集的更多相关文章
- SDOI 2014 向量集
[SDOI2014]向量集 题目描述 维护一个向量集合,在线支持以下操作: - "A x y (|x|,|y| < =10^8)":加入向量(x,y); - " Q ...
- [BZOJ 3530][Sdoi 2014]数数
阿拉~好像最近总是做到 AC 自动机的题目呢喵~ 题目的算法似乎马上就能猜到的样子…… AC 自动机 + 数位 dp 先暴力转移出 f[i][j] :表示从 AC 自动机上第 j 号节点走 i 步且不 ...
- BZOJ 3533: [Sdoi2014]向量集( 线段树 + 三分 )
答案一定是在凸壳上的(y>0上凸壳, y<0下凸壳). 线段树维护, 至多N次询问, 每次询问影响O(logN)数量级的线段树结点, 每个结点O(logN)暴力建凸壳, 然后O(logN) ...
- [BZOJ 3123] [SDOI 2013]森林(可持久化线段树+并查集+启发式合并)
[BZOJ 3123] [SDOI 2013]森林(可持久化线段树+启发式合并) 题面 给出一个n个节点m条边的森林,每个节点都有一个权值.有两种操作: Q x y k查询点x到点y路径上所有的权值中 ...
- 【SDOI2014】向量集
[SDOI2014]向量集 题目描述 我们分析一波: 假设我们询问\((A,B)\),\(x_i>x_j\)若 \[ A\cdot x_i+B\cdot y_i>A\cdot x_j+B\ ...
- 「SDOI2014」向量集 解题报告
「SDOI2014」向量集 维护一个向量集合,在线支持以下操作: A x y :加入向量 \((x, y)\): Q x y l r:询问第 \(L\) 个到第 \(R\) 个加入的向量与向量 \(( ...
- [BZOJ 1879][SDOI 2009]Bill的挑战 题解(状压DP)
[BZOJ 1879][SDOI 2009]Bill的挑战 Description Solution 1.考虑状压的方式. 方案1:如果我们把每一个字符串压起来,用一个布尔数组表示与每一个字母的匹配关 ...
- [BZOJ 2299][HAOI 2011]向量 题解(裴蜀定理)
[BZOJ 2299][HAOI 2011]向量 Description 给你一对数a,b,你可以任意使用(a,b), (a,-b), (-a,b), (-a,-b), (b,a), (b,-a), ...
- BZOJ3533 [Sdoi2014]向量集 【线段树 + 凸包 + 三分】
题目链接 BZOJ3533 题解 我们设询问的向量为\((x_0,y_0)\),参与乘积的向量为\((x,y)\) 则有 \[ \begin{aligned} ans &= x_0x + y_ ...
随机推荐
- JS 截取字符串函数
一.函数:split() 功能:使用一个指定的分隔符把一个字符串分割存储到数组 例子: str=”jpg|bmp|gif|ico|png”; arr=theString.split(”|”); //a ...
- 抽象类[abstract]_C#
抽象类(abstract) abstract修饰符可以和类.方法.属性.索引器及事件一起使用.在类声明中使用abstract修饰符以指示某个类只能是其它类的基类.标记为抽象或包含在抽象类中的成员必须通 ...
- Emit
http://www.cnblogs.com/zhuweisky/archive/2008/09/20/1294666.html http://www.cnblogs.com/xiaoxiangfei ...
- 服务器 tfs不提供 TeamFoundation服务。基础连接已经关闭
服务器 tfs(服务器名或url)不提供 TeamFoundation服务.基础连接已经关闭,发送时发生错误.TFS突然间连接不上到,到服务器上配置团队项目的组成员资格提示这样的错误,客户端连接的时候 ...
- 关于linux下零散的东西 --慢慢补充
一.截图 ,使用Shift+Ctrl+PrtSc就可以截图. 二.tar命令参数 c:表示压缩 x:表示解压 z:表示gzip的方式解/压缩 j:表示bzip2的方式解/压缩 三.串口终端ke ...
- java.imageIo给图片添加水印
最近项目在做一个商城项目, 项目上的图片要添加水印①,添加图片水印;②:添加文字水印; 一下提供下个方法,希望大家可以用得着: package com.blogs.image; import java ...
- mysql table readonly
if the table does not have primary key or unique non-nullable defined, then MySql workbench could no ...
- php入门单引号与双引号区别
[1]单引号和双引号在处理变量的时候做法: 括在双引号内的变量会解释出值,但是括在单引号内则不做处理,直接输出: <?php $var = 'my name is huige'; echo &q ...
- 1080. Graduate Admission (30)
时间限制 200 ms 内存限制 65536 kB 代码长度限制 16000 B 判题程序 Standard 作者 CHEN, Yue It is said that in 2013, there w ...
- ^M的前世今生
^M 是ascii中的'\r', 回车符,是16进制的0x0D,8进制的015,十进制的13. 回想计算机发展,旧的打字机使用两个字符来开始一个新行.首先是一个字符命令使打印头移回开始位置(回车),然 ...