注意到线段树一个节点只有满了才会被用到,那时再建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. html5 +css3 点击后水波纹扩散效果 兼容移动端

    <!DOCTYPE html> <html> <head> <meta http-equiv="Content-Type" content ...

  2. Java中创建线程的三种方法以及区别

    Java使用Thread类代表线程,所有的线程对象都必须是Thread类或其子类的实例.Java可以用三种方式来创建线程,如下所示: 1)继承Thread类创建线程 2)实现Runnable接口创建线 ...

  3. 详细讲解安全升级MySQL的方法

    MySQL升级是非常必要的. 我们在Percona Support上列出了关于MySQL升级最佳实践的各种问题.这篇文章推荐了一些不同情况下升级MySQL的方法. 为什么MySQL升级是必须的? 原因 ...

  4. Java 图片转字节流 实现 图片->字节流(字符串)->图片

    //该方法实现图片转String 参数为图片的路径 可以是file.toString()得到public String testUpload(String path) { try { String s ...

  5. 基本控件文档-UISlider属性

    CHENYILONG Blog 基本控件文档-UISlider属性 Fullscreen     UISlide属性技术博客http://www.cnblogs.com/ChenYilong/ 新浪微 ...

  6. phpcms 仿站小结

    1.title<title>{if isset($SEO['title']) && !empty($SEO['title'])}{$SEO['title']}{/if}{$ ...

  7. JavaScript基本运算

    JavaScript基本运算 JavaScript:其实它的基本运算和我们数学的基本运算类似的. --------------------------------------------------- ...

  8. post请求远程url 报错“基础连接已经关闭...Authentication.AuthenticationException...远程证书无效”解决方案

    当我们有时用代码编写post请求url远程地址会报“基础连接已经关闭: 未能为 SSL/TLS 安全通道建立信任关系. ---> System.Security.Authentication.A ...

  9. JavaScript常用数组方法

    JavaScript数组方法有以下这些: forEach() 方法对数组的每一个元素执行一次提供的函数. map() 方法创建一个新数组,其结果是该数组都执行一次函数,原函数保持不变. filter( ...

  10. Linux信息搜集

    ## 1.取证工具 - LiME 内存获取工具 - volatility 内存分析工具 ## 2.机器信息收集 #sysinfo 16 # # 查看当前登录用户 who > who.txt # ...