[BZOJ2961] 共点圆 [cdq分治+凸包]
题面
思路
首先考虑一个点$(x_0,y_0)$什么时候在一个圆$(x_1,y_1,\sqrt{x_12+y_12})$内
显然有:$x_12+y_12\geq (x_0-x_1)2+(y_0-y_1)2$
化简:$2x_0x_1+2y_0y_1\geq x_02+y_02$
所有含$x_1,y_1$的项挪到同一边,除掉一个$2y_0$(假设它是正的),得到:
$y_1\geq -\frac{x_0}{y_0}x_1+\frac{x_02+y_02}{2y_0}$
如果是负的:
$y_1\leq -\frac{x_0}{y_0}x_1+\frac{x_02+y_02}{2y_0}$
DUANG!半平面来了
那么现在的询问变成了:给定一个半平面,问是不是所有的点都在这个半平面的上方(或者下方)
显然,我们如果维护了所有输入节点的上下凸包,这个问题就迎刃而解了
众所周知,维护动态上下凸壳可以用$set$或者平衡树做到$O(n\log n)$
然而博主并不想写这种东西
所以他写了非常沙雕的cdq分治23333
分治
分治开始之前先按照所有询问点的斜率排个序(非询问点不管)
我们对时间顺序分治
进入分治后,首先分治左区间,返回按照横坐标排好序的左区间所有点
求出左区间中所有非询问的点的上下凸壳
然后对右边的所有询问,因为一开始排好序了,所以直接在凸壳上顺次双指针过去
注意:$y_0 > 0$的时候所有点在直线上方,用的是下凸包,反之亦然
最后分治右区间,按照x坐标归并排序,合并左右区间
就是个很模板的cdq斜率分治
Code
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cassert>
#include<cmath>
#define sqr(x) ((x)*(x))
#define eps 1e-10
#define ll long long
using namespace std;
inline int read(){
int re=0,flag=1;char ch=getchar();
while(!isdigit(ch)){
if(ch=='-') flag=-1;
ch=getchar();
}
while(isdigit(ch)) re=(re<<1)+(re<<3)+ch-'0',ch=getchar();
return re*flag;
}
struct node{
double x,y,k;
int id,qid,op;
inline friend double slope(const node &a,const node &b){return ((fabs(a.x-b.x)<eps)?(1e30):((a.y-b.y)/(a.x-b.x)));}
inline friend double dis(const node &a,const node &b){return sqrt(sqr(a.x-b.x)+sqr(a.y-b.y));}
inline friend double rad(const node &a){return sqrt(sqr(a.x)+sqr(a.y));}
}a[500010],tmp[500010],q1[500010],q2[500010];//q1 up q2 down
int n,top1,top2;bool ans[500010];
inline bool cmp(const node &a,const node &b){return a.k<b.k;}
void solve(int l,int r){
if(l==r) return;
int i,t1,t2,mid=(l+r)>>1;
t1=l;t2=mid+1;
for(i=l;i<=r;i++){
if(a[i].id<=mid) tmp[t1++]=a[i];
else tmp[t2++]=a[i];
}
memcpy(a+l,tmp+l,sizeof(node)*(r-l+1));
solve(l,mid);
top1=top2=0;
for(i=l;i<=mid;i++){
if(a[i].op) continue;
while(top1>1&&slope(q1[top1-1],q1[top1])<slope(q1[top1],a[i])+eps) top1--;
q1[++top1]=a[i];
while(top2>1&&slope(q2[top2-1],q2[top2])+eps>slope(q2[top2],a[i])) top2--;
q2[++top2]=a[i];
}
t1=t2=1;
for(i=mid+1;i<=r;i++){
if(!a[i].op) continue;
if(a[i].y>0){
while(t2<top2&&slope(q2[t2],q2[t2+1])<a[i].k) t2++;
if(t2<=top2) ans[a[i].qid]&=(dis(a[i],q2[t2])<rad(q2[t2]));
}
else{
while(t1<top1&&slope(q1[top1-1],q1[top1])<a[i].k) top1--;
if(t1<=top1) ans[a[i].qid]&=(dis(a[i],q1[top1])<rad(q1[top1]));
}
}
solve(mid+1,r);
t1=l;t2=mid+1;
for(i=l;i<=r;i++){
if(t2==r+1||(t1<=mid&&a[t1].x<a[t2].x)) tmp[i]=a[t1++];
else tmp[i]=a[t2++];
}
memcpy(a+l,tmp+l,sizeof(node)*(r-l+1));
}
int main(){
n=read();int i,flag=0,cntq=0;
for(i=1;i<=n;i++){
scanf("%d%lf%lf",&a[i].op,&a[i].x,&a[i].y);
a[i].id=i;
if(a[i].op){
a[i].qid=++cntq;
if(flag) ans[cntq]=1;
if(a[i].y) a[i].k=-a[i].x/a[i].y;
else a[i].k=1e30;
}
else flag=1;
}
sort(a+1,a+n+1,cmp);
solve(1,n);
for(i=1;i<=cntq;i++) puts(ans[i]?"Yes":"No");
}
[BZOJ2961] 共点圆 [cdq分治+凸包]的更多相关文章
- BZOJ2961 共点圆[CDQ分治]
题面 bzoj 其实就是推一下圆的式子 长成这个样子 假设要查询的点是(x, y) 某个圆心是(p, q) \((x - p)^2 + (y - q)^2 \leq p^2 + q^2\) 变成 \( ...
- BZOJ2961: 共点圆(CDQ分治+凸包)
题面 传送门 题解 这题解法真是多啊--据说可以圆反演转化为动态插入半平面并判断给定点是否在半平面交中,或者化一下改成给定点判断是否所有点都在某一个半平面内-- 鉴于圆反演我也不会,这里讲一下直接推的 ...
- bzoj2961 共点圆 (CDQ分治, 凸包)
/* 可以发现可行的圆心相对于我们要查询的点是在一个半平面上, 然后我们要做的就是动态维护凸壳然后用这个半平面去切它 看看是否是在合法的那一面 然后cdq分治就可以了 代码基本是抄的, */ #inc ...
- bzoj2961 共点圆 bzoj 4140
题解: 比较水的一道题 首先我们化简一下式子发现是维护xxo+yyo的最值 显然是用凸包来做 我们可以直接用支持插入删除的凸包 也是nlogn的 因为没有强制在线,我们也可以cdq,考虑前面一半对答案 ...
- Bzoj2149拆迁队:cdq分治 凸包
国际惯例的题面:我们考虑大力DP.首先重新定义代价为:1e13*选择数量-(总高度+总补偿).这样我们只需要一个long long就能维护.然后重新定义高度为heighti - i,这样我们能选择高度 ...
- BZOJ2961: 共点圆
好久没发了 CDQ分治,具体做法见XHR的论文… /************************************************************** Problem: 29 ...
- [BZOJ2961]共点圆-[凸包+cdq分治]
Description 传送门 Solution 考虑对于每一个点: 设圆的坐标为(x,y),点的坐标为(x0,y0).依题意得,当一个点在圆里,需要满足(x-x0)2+(y-y0)2<=x2+ ...
- bzoj 2961 共点圆 cdq+凸包+三分
题目大意 两种操作 1)插入一个过原点的圆 2)询问一个点是否在所有的圆中 分析 在圆中则在半径范围内 设圆心 \(x,y\) 查询点\(x_0,y_0\) 则\(\sqrt{(x-x_0)^2+(y ...
- 【bzoj2961】共点圆 k-d树
更新:此题我的代码设置eps=1e-8会WA,现在改为1e-9貌似T了 此题网上的大部分做法是cdq分治+凸包,然而我觉得太烦了,于是自己口胡了一个k-d树做法: 加入一个圆$(x,y)$,直接在k- ...
随机推荐
- Hive支持行级update、delete时遇到的问题
Hive从0.14版本开始支持事务和行级更新,但缺省是不支持的,需要一些附加的配置.要想支持行级insert.update.delete,需要配置Hive支持事务.(行级的insert好像不配置也能运 ...
- unity3d 计时功能舒爽解决方案
上次也写了一篇计时功能的博客 今天这篇文章和上次的文章实现思路不一样,结果一样 上篇文章地址:http://www.cnblogs.com/shenggege/p/4251123.html 思路决定一 ...
- 「题目代码」P1054~P1059(Java)
P1054 猴子吃桃 import java.util.*; import java.io.*; import java.math.BigInteger; import java.lang.Chara ...
- Ubuntu目录与权限
Ubuntu目录 / /bin /sbin /boot /etc /mnt /home d :directory - :file b :block 磁盘以块进行 l :link Ubuntu权限 U ...
- 图像质量评价指标之 PSNR 和 SSIM
1. PSNR (Peak Signal-to-Noise Ratio) 峰值信噪比 给定一个大小为 \(m×n\) 的干净图像 \(I\) 和噪声图像 \(K\),均方误差 \((MSE)\) 定义 ...
- Fluent Python: Classmethod vs Staticmethod
Fluent Python一书9.4节比较了 Classmethod 和 Staticmethod 两个装饰器的区别: 给出的结论是一个非常有用(Classmethod), 一个不太有用(Static ...
- Python中的list
list的创建 1 字面量 >>>L = [1, 2, 3] [1, 2, 3] 2 通过iterable可迭代对象,比如str对象,range对象,map对象 >>&g ...
- protected、public、private
一.protected成员 1. 受保护的成员的可访问性 对于一个类的protected成员,①该类的用户(如类对象)不能访问它,②该类的派生类的成员(及其友元)可以访问它. 派生类的成员及其友元不能 ...
- Thunder团队——文案+美工
团队名称:Thunder 组长:王航 成员:李传康.代秋彤.邹双黛.苗威.宋雨.胡佑蓉.杨梓瑞 项目名称:爱阅app 需求概述: 现在市面上有很多手机阅读器,但是基本上不是收费就是广告满天飞.基于这种 ...
- Java 变量和输入输出
一些重要知识 一个源文件里只能有一个public类,其它类数量不限.文件名与public类名相同 JAVA程序严格区分大小写 JAVA应用程序的执行入口是main方法固定写法:public stati ...