注意到线段树一个节点只有满了才会被用到,那时再建ConvexHull就行了。。。

#include <bits/stdc++.h>
using namespace std;
namespace my_useful_tools {
#define rep(_i, _k, _j) for(int _i = _k; _i <= _j; ++_i)
#define reu(_i, _k, _j) for(int _i = _k; _i < _j; ++_i)
#define red(_i, _k, _j) for(int _i = _k; _j <= _i; --_i)
#define foreach(_i, _s) for(typeof(_s.begin()) _i = _s.begin(); _i != _s.end(); ++_i)
#define pb push_back
#define mp make_pair
#define ipir pair<int, int>
#define ivec vector<int>
#define clr(t) memset(t, 0, sizeof t)
#define pse(t, v) memset(t, v, sizeof t)
#define brl puts("")
#define file(x) freopen(#x".in", "r", stdin), freopen(#x".out", "w", stdout)
#define file_hza freopen("input.txt", "r", stdin), freopen("output.txt", "w", stdout);
#define file_gen(x) freopen(#x".in", "w", stdout);
const int INF = 0x3f3f3f3f;
typedef long long LL;
typedef double DB;
inline void pc(char c) { putchar(c); }
template<class T> inline T gcd(T a, T b) { return b == ? a : gcd(b, a % b); }
template<class T> inline void W(T p) { if(p < ) pc('-'), p = -p; if(p / != ) W(p / ); pc('' + p % ); }
template<class T> inline void Wn(T p) { W(p), brl; } template<class T> inline void W(T a, T b) { W(a), pc(' '), W(b); }
template<class T> inline void Wn(T a, T b) { W(a), pc(' '), Wn(b); }
template<class T> inline void W(T a, T b, T c) { W(a), pc(' '), W(b), pc(' '), W(c); }
inline char gchar() { char ret = getchar(); for(; ret == '\n' || ret == '\r' || ret == ' '; ret = getchar()); return ret; }
template<class T> inline void fr(T&ret) { char c = ' '; int flag = ; for(c = getchar(); c != '-' && !('' <= c && c <= ''); c = getchar());
if(c == '-') flag = -, ret = ; else ret = c - ''; for(c = getchar(); '' <= c && c <= ''; c = getchar()) ret = ret * + c - '';
ret = ret * flag;
}
inline int fr() { int x; fr(x); return x; }
template<class T> inline void fr(T&a, T&b) { fr(a), fr(b); } template<class T> inline void fr(T&a, T&b, T&c) { fr(a), fr(b), fr(c); }
template<class T> inline T fast_pow(T base, T index, T mod = , T ret = ) {
for(; index; index >>= , base = base * base % mod) if(index & ) ret = ret * base % mod;
return ret;
}
};
using namespace my_useful_tools; const int maxNode = 1e7;
const int maxQueryNode = 4e5 + ;
const int maxSegNode = (maxQueryNode<<);
struct Point {
int x, y;
Point() {}
Point(int x, int y):x(x), y(y) {}
bool operator < (const Point&rhs) const {
return x < rhs.x || (x==rhs.x && y < rhs.y);
}
Point operator - (const Point&rhs) {
return Point(x-rhs.x, y-rhs.y);
}
} pool[maxNode], *loc = pool, P[maxQueryNode];
LL Dot(Point a, Point b) {
return (LL)a.x*b.x+(LL)a.y*b.y;
}
LL Cross(Point a, Point b) {
return (LL)a.x*b.y-(LL)a.y*b.x;
} struct ConvexHull {
Point *up, *dw;
int un, dn;
void buildConvexHull(Point*p, int sz) {
sort(p, p+sz);
up = loc;
int cnt = -;
reu(i, , sz) {
while (cnt > && Cross(p[i]-up[cnt-], up[cnt]-up[cnt-]) <= )
--cnt;
++cnt, up[cnt] = p[i];
}
loc += cnt+;
un = cnt+;
dw = loc;
cnt = -;
reu(i, , sz) {
while (cnt > && Cross(p[i]-dw[cnt-], dw[cnt]-dw[cnt-]) >= )
--cnt;
++cnt, dw[cnt] = p[i];
}
loc += cnt+;
dn = cnt+;
}
LL qMaxDot(Point p) {
LL ret = LLONG_MIN;
int l = , r = un-;
while (l <= r) {
if (r - l + <= ) {
for (register int i = l; i <= r; ++i) {
ret = max(ret, Dot(up[i], p));
}
break;
}
int mid = (l+r)>>;
if (Dot(up[mid], p) < Dot(up[mid+], p)) {
l = mid+;
} else {
r = mid;
}
}
l = , r = dn-;
while (l <= r) {
if (r - l + <= ) {
for (register int i = l; i <= r; ++i) {
ret = max(ret, Dot(dw[i], p));
}
break;
}
int mid = (l+r)>>;
if (Dot(dw[mid], p) < Dot(dw[mid+], p)) {
l = mid+;
} else {
r = mid;
}
}
// return ret;
// LL tmp = LLONG_MIN;
// for (int i = 0; i < un; ++i)
// tmp = max(tmp, Dot(up[i], p));
// for (int i = 0; i < dn; ++i)
// tmp = max(tmp, Dot(dw[i], p));
return ret;
}
} ch_pool[maxSegNode], *ch_loc = ch_pool; struct SegNode {
SegNode*lc, *rc;
ConvexHull*ch;
int l, r, sz;
SegNode() {
lc = rc = NULL;
ch = ch_loc++;
sz = ;
}
} sn_pool[maxSegNode], *sn_loc = sn_pool, *root; void buildSegTree(int l, int r, SegNode*&rt) {
if (rt == NULL) {
rt = sn_loc++;
rt->l = l, rt->r = r;
}
if (l == r) return ;
int mid = (l+r)>>;
buildSegTree(l, mid, rt->lc), buildSegTree(mid+, r, rt->rc);
} void addPoint(int pos, SegNode*rt) {
if (rt->l == rt->r) {
++rt->sz;
if (rt->sz == rt->r-rt->l+)
rt->ch->buildConvexHull(P+(rt->l), rt->sz);
return ;
}
int mid = (rt->l+rt->r)>>;
if (pos <= mid) addPoint(pos, rt->lc);
else addPoint(pos, rt->rc);
rt->sz = rt->lc->sz+rt->rc->sz;
if (rt->sz == rt->r-rt->l+)
rt->ch->buildConvexHull(P+(rt->l), rt->sz);
} Point qp;
LL Query(int ql, int qr, SegNode*rt) {
if (ql <= rt->l && rt->r <= qr) {
return rt->ch->qMaxDot(qp);
}
int mid = (rt->l+rt->r)>>;
LL ret = LLONG_MIN;
if (ql <= mid) ret = max(ret, Query(ql, qr, rt->lc));
if (mid < qr) ret = max(ret, Query(ql, qr, rt->rc));
return ret;
} int n;
bool e;
LL last_ans = 0LL; inline int decode(int x) {
if (e)
return x ^ (last_ans & 0x7fffffff);
else return x;
} void print_tree(SegNode*rt) {
if (rt == NULL) return ;
printf("%d %d %d\n", rt->l, rt->r, rt->sz);
print_tree(rt->lc);
print_tree(rt->rc);
} int main() {
fr(n), e=(gchar()!='E');
int x, y, ql, qr;
int cnt = ;
buildSegTree(, n, root);
while (n--) {
if (gchar() == 'A') {
x = decode(fr()), y = decode(fr());
P[++cnt] = Point(x, y);
addPoint(cnt, root);
} else {
x = decode(fr()), y = decode(fr());
ql = decode(fr()), qr = decode(fr());
// printf("%d %d %d %d\n", x, y, ql, qr);
qp = Point(x, y);
last_ans = Query(ql, qr, root);
printf("%lld\n", last_ans);
}
// print_tree(root);
} return ;
}

sdoi2014-向量集-线段树-二分斜率的更多相关文章

  1. BZOJ 3533: [Sdoi2014]向量集( 线段树 + 三分 )

    答案一定是在凸壳上的(y>0上凸壳, y<0下凸壳). 线段树维护, 至多N次询问, 每次询问影响O(logN)数量级的线段树结点, 每个结点O(logN)暴力建凸壳, 然后O(logN) ...

  2. 【bzoj3533】[Sdoi2014]向量集 线段树+STL-vector维护凸包

    题目描述 维护一个向量集合,在线支持以下操作:"A x y (|x|,|y| < =10^8)":加入向量(x,y);"Q x y l r (|x|,|y| < ...

  3. BZOJ3533:[SDOI2014]向量集(线段树,三分,凸包)

    Description 维护一个向量集合,在线支持以下操作: "A x y (|x|,|y| < =10^8)":加入向量(x,y); " Q x y l r (| ...

  4. bzoj 3533 [Sdoi2014]向量集 线段树+凸包+三分(+动态开数组) 好题

    题目大意 维护一个向量集合,在线支持以下操作: "A x y (|x|,|y| < =10^8)":加入向量(x,y); "Q x y l r (|x|,|y| & ...

  5. bzoj 3533: [Sdoi2014]向量集 线段树维护凸包

    题目大意: http://www.lydsy.com/JudgeOnline/problem.php?id=3533 题解: 首先我们把这些向量都平移到原点.这样我们就发现: 对于每次询问所得到的an ...

  6. [SDOI2014][BZOJ3533] 向量集 [线段树+凸包]

    题面 BZOJ传送门 思路 首先当然是推式子 对于一个询问点$(x_0,y_0$和给定向量$(x_1,y_1)$来说,点积这么表达: $A=x_0x_1+y_0y_1$ 首先肯定是考虑大小关系:$x_ ...

  7. bzoj4311向量(线段树分治+斜率优化)

    第二道线段树分治. 首先设当前向量是(x,y),剩余有两个不同的向量(u1,v1)(u2,v2),假设u1>u2,则移项可得,若(u1,v1)优于(u2,v2),则-x/y>(v1-v2) ...

  8. 【BZOJ4311】向量(线段树分治,斜率优化)

    [BZOJ4311]向量(线段树分治,斜率优化) 题面 BZOJ 题解 先考虑对于给定的向量集,如何求解和当前向量的最大内积. 设当前向量\((x,y)\),有两个不同的向量\((u1,v1),(u2 ...

  9. bzoj4399 魔法少女LJJ 线段树合并+线段树二分+并查集

    题目传送门 https://lydsy.com/JudgeOnline/problem.php?id=4399 题解 毒瘤题 \(9\) 种操作还有支持动态图的连通性 仔细读题 $ c<=7$. ...

随机推荐

  1. samba和nginx服务

    samba和nginx服务 1.s配置amba samba的功能: samba是一个网络服务器,用于Linux和Windows之间共享文件. 1.1配置环境 关闭防火墙和selinux systemc ...

  2. 解析Job,bpmn文件的小项目总结

    1.在使用String类中split(String regex)切割字符串abcd.job遇得到job字符串时,直接使用split("."),导致数组超出界限错误 原因:得到的数组 ...

  3. 科学计算三维可视化---Mayavi入门(Mayavi管线)

    一:Mayavi管线 mlab.show_pipeline() #显示管线层级,来打开管线对话框 (一)管线中的对象scene Mayavi Scene:处于树的最顶层的对象,他表示场景,配置界面中可 ...

  4. PHP常见数组方法和函数

    current();当前游标指向的数组单元值 next();下一个数组单元值 end()最后一个 reset()复位 prev()把数组指针往前一位 写法:$arr=array('a','b','c' ...

  5. 手把手教你使用koa2

    简介 koa是由express的原班人马打造的web框架.但是相对于express,koa的性能要更高,因为koa通过使用aysnc函数,帮你丢弃回调函数,并有力的增强了错误处理.而且koa没有绑定任 ...

  6. 2015/11/6用Python写游戏,pygame入门(6):控制大量的对象

    昨天我们已经实现了这个游戏的三个基本类. 但是现在它还是没办法做成一个适合玩的游戏,毕竟只有一架敌机的游戏是很乏味的.所以,我们需要好多子弹,也需要好多敌机. 所以,我们要创建list,这个list存 ...

  7. Android的音频解码原来是直接调用的本地C方法直接通过硬件解码

    Android就是披着JAVA外衣的C啊~音频解码原来是直接调用的本地C方法直接通过硬件解码的,JAVA和C的字节数组存放模式不同(java是大端,C根据不同平台不同),不同格式需要转化以后才能用. ...

  8. Asp.Net使用加密cookie代替session验证用户登录状态 源码分享

    首先 session 和 cache 拥有各自的优势而存在.  他们的优劣就不在这里讨论了. 本实例仅存储用户id于用户名,对于多级权限的架构,可以自行修改增加权限字段   本实例采用vs2010编写 ...

  9. Linux系统接入小区宽带

    jollywing(jollywing@foxmail.com) 安装 rp-pppoe 今天去联通营业厅开通了家庭宽带,回到家就搜索怎么用Linux接入小区宽带,发现大多数人都选择用PPPOE拨号上 ...

  10. Django 基础命令