【BZOJ】4056: [Ctsc2015]shallot
题意
在线、可持久化地维护一条二维平面上的折线,支持查询与任意一条直线的交点个数。
点的个数和操作个数小于\(10^5\)
分析
一条折线可以用一个序列表示,可持久化序列考虑用可持久化treap。
如何判断交点?如果有交点,那么一定与包含这个折线的矩阵有交点。
题解
所以我们可持久化treap一下即可,虽然这个复杂度很不靠谱,纯rp算法。
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int Lim=20000000, M=1e5+10, oo=~0u>>1;
struct node *null;
struct node {
int x, y, s, mx[2], mn[2];
node *c[2], *ch[2];
void up() {
s=c[0]->s+c[1]->s+1;
mx[0]=max(x, max(c[0]->mx[0], c[1]->mx[0]));
mn[0]=min(x, min(c[0]->mn[0], c[1]->mn[0]));
mx[1]=max(y, max(c[0]->mx[1], c[1]->mx[1]));
mn[1]=min(y, min(c[0]->mn[1], c[1]->mn[1]));
if(c[0]!=null) ch[0]=c[0]->ch[0]; else ch[0]=this;
if(c[1]!=null) ch[1]=c[1]->ch[1]; else ch[1]=this;
}
void init(int _x, int _y) {
x=mx[0]=mn[0]=_x;
y=mx[1]=mn[1]=_y;
c[0]=c[1]=null;
ch[0]=ch[1]=this;
s=1;
}
}Po[Lim], *iT=Po, *t[M];
node *newnode(int x, int y) {
iT->init(x, y);
return iT++;
}
void init() {
null=iT++;
null->init(0, 0);
null->mx[0]=null->mx[1]=-oo;
null->mn[0]=null->mn[1]=oo;
null->s=0;
}
node *build(int l, int r) {
if(l>r) {
return null;
}
node *ret;
int x, y, mid=(l+r)>>1;
node *le=build(l, mid-1);
scanf("%d%d", &x, &y);
node *ri=build(mid+1, r);
ret=newnode(x, y);
ret->c[0]=le;
ret->c[1]=ri;
ret->up();
return ret;
}
node *update(int p, int X, int Y, node *x) {
node *y;
if(x->s==p) {
y=newnode(X, Y);
y->c[0]=x;
y->up();
return y;
}
y=iT++;
*y=*x;
if(y->c[0]->s>=p) {
y->c[0]=update(p, X, Y, y->c[0]);
}
else {
y->c[1]=update(p-y->c[0]->s-1, X, Y, y->c[1]);
}
y->up();
return y;
}
ll cross(int x0, int y0, int x, int y) {
return (ll)x0*y-(ll)x*y0;
}
bool jiao(int x, int y, int xx, int yy, int x0, int y0, int X, int Y) {
ll a=cross(x-x0, y-y0, X, Y),
b=cross(xx-x0, yy-y0, X, Y);
a=a>=0?a>0:-1;
b=b>=0?b>0:-1;
return a*b<=0;
}
bool jiao(node *x, int x0, int y0, int X, int Y) {
return jiao(x->mn[0], x->mx[1], x->mn[0], x->mn[1], x0, y0, X, Y) ||
jiao(x->mn[0], x->mn[1], x->mx[0], x->mn[1], x0, y0, X, Y) ||
jiao(x->mx[0], x->mn[1], x->mx[0], x->mx[1], x0, y0, X, Y) ||
jiao(x->mx[0], x->mx[1], x->mn[0], x->mx[1], x0, y0, X, Y);
}
int query(int x0, int y0, int X, int Y, node *x) {
if(x->s<=1 || !jiao(x, x0, y0, X, Y)) {
return 0;
}
int ret=0;
if(x->c[0]!=null && jiao(x->x, x->y, x->c[0]->ch[1]->x, x->c[0]->ch[1]->y, x0, y0, X, Y)) {
++ret;
}
if(x->c[1]!=null && jiao(x->x, x->y, x->c[1]->ch[0]->x, x->c[1]->ch[0]->y, x0, y0, X, Y)) {
++ret;
}
return query(x0, y0, X, Y, x->c[0])+query(x0, y0, X, Y, x->c[1])+ret;
}
int main() {
init();
int n, m, cn, last=0;
scanf("%d%d%d", &n, &m, &cn);
t[0]=build(1, n);
for(int kk=1; kk<=m; ++kk) {
static char s[5];
int T, x0, y0, x, y;
scanf("%s", s);
if(s[0]=='H') {
scanf("%d%d%d%d%d", &T, &x0, &y0, &x, &y);
if(cn) {
x0^=last;
y0^=last;
x^=last;
y^=last;
}
t[kk]=t[T];
printf("%d\n", last=query(x0, y0, x, y, t[kk]));
}
else {
scanf("%d%d%d%d", &T, &x0, &x, &y);
if(cn) {
x^=last;
y^=last;
}
t[kk]=update(x0, x, y, t[T]);
}
}
return 0;
}
【BZOJ】4056: [Ctsc2015]shallot的更多相关文章
- 【BZOJ】3052: [wc2013]糖果公园
http://www.lydsy.com/JudgeOnline/problem.php?id=3052 题意:n个带颜色的点(m种),q次询问,每次询问x到y的路径上sum{w[次数]*v[颜色]} ...
- 【BZOJ】3319: 黑白树
http://www.lydsy.com/JudgeOnline/problem.php?id=3319 题意:给一棵n节点的树(n<=1e6),m个操作(m<=1e6),每次操作有两种: ...
- 【BZOJ】3319: 黑白树(并查集+特殊的技巧/-树链剖分+线段树)
http://www.lydsy.com/JudgeOnline/problem.php?id=3319 以为是模板题就复习了下hld............................. 然后n ...
- 【BZOJ】1013: [JSOI2008]球形空间产生器sphere
[BZOJ]1013: [JSOI2008]球形空间产生器sphere 题意:给n+1个n维的点的坐标,要你求出一个到这n+1个点距离相等的点的坐标: 思路:高斯消元即第i个点和第i+1个点处理出一个 ...
- 【BZOJ】1002:轮状病毒(基尔霍夫矩阵【附公式推导】或打表)
Description 轮状病毒有很多变种,所有轮状病毒的变种都是从一个轮状基产生的.一个N轮状基由圆环上N个不同的基原子和圆心处一个核原子构成的,2个原子之间的边表示这2个原子之间的信息通道.如下图 ...
- 【BZOJ】【3083】遥远的国度
树链剖分/dfs序 其实过了[BZOJ][4034][HAOI2015]T2以后就好搞了…… 链修改+子树查询+换根 其实静态树的换根直接树链剖分就可以搞了…… 因为其实只有一样变了:子树 如果roo ...
- 【BZOJ】【2434】【NOI2011】阿狸的打字机
AC自动机+DFS序+BIT 好题啊……orz PoPoQQQ 大爷 一道相似的题目:[BZOJ][3172][TJOI2013]单词 那道题也是在fail树上数有多少个点,只不过这题是在x的fail ...
- 【BZOJ】【2738】&【Tsinsen】【A1333】矩阵乘法
整体二分+树状数组 过了[BZOJ][2527][POI2011]Meteors以后这题就没那么难啦~ 关键是[从小到大]依次插入数字,然后整体二分每个查询的第k大是在第几次插入中被插入的……嗯大概就 ...
- 【BZOJ】【3170】【TJOI2103】松鼠聚会
切比雪夫距离+曼哈顿距离 题解:http://www.cnblogs.com/zyfzyf/p/4105456.html 其实应该先做这题再做[BZOJ][3210]花神的浇花集会的吧…… 我们发现d ...
随机推荐
- required - HTML5里的input标签的required属性提示文字修改
input 里面增加这样的语句: <input type="text" placeholder="您的姓名" required oninvalid=&qu ...
- Java transient 关键字
1)一旦变量被transient修饰,变量将不再是对象持久化的一部分,该变量内容在序列化后无法获得访问. 2)transient关键字只能修饰变量,而不能修饰方法和类.注意,本地变量是不能被trans ...
- mysql数据库史上最详细起步教程(1)
本文主要讲解mysql的操作,尽量保证步骤的详细与清晰,希望能帮到大家. 1.登录后进行数据库的创建:create database lf(数据库名); (一定要记住分号,mysql在语句的结束符就 ...
- 部署React+webpack工程的步骤
# 部署React+webpack工程的步骤ps:以Mac os系统做开发环境.因为npm现在使用灰常的慢,所以我使用淘宝境像cnpm. 1,准备工作: 先确保存已经安装了node.js: 2,文件部 ...
- DLL 生成与使用的全过程(2010-01-18 14:50:17)
转载自 水滴的博客http://blog.sina.com.cn/spiritofwater 个人学习用 转载▼ 分类: 技术 由dll导出的lib文件: 包含了每一个dll导出函数的符号名和 ...
- 【HEOI2012】采花 BZOJ2743
Description 萧芸斓是Z国的公主,平时的一大爱好是采花. 今天天气晴朗,阳光明媚,公主清晨便去了皇宫中新建的花园采花.花园足够大,容纳了n朵花,花有c种颜色(用整数1-c表示),且花是排成一 ...
- docker 初探之简单安装 ----Windows10
报错一 $ docker run hello-world Unable to find image 'hello-world:latest' locally Pulling repository do ...
- linux之netstat命令
netstat用于显示各种网络相关的信息,如网络连接,路由表接口状态(interface statistics), masquerade连接,多播成员(Multicast Memberships)等等 ...
- .net测试学习--理解.net测试选项
1.创建基于测试简单应用程序 (1)启动visual studio(有安装c#的) (2) 选择File|New project (3)创建一个C# project,名字和保存路径自己设定,假设取名 ...
- How to Fix GNOME License Not Accepted Issue on CentOS 7
This post assume that you have just finished the Gnome GUI installation on CentOS 7 by using “yum gr ...