传送门

做法:每个卫星分别用连到左边圆与x轴交点的线的斜率和连到右边交点的线旋转90度的斜率可以表示成一个区间,问题转化成支持加/删区间和询问其中两个区间是否有交以及它们的交是否被其他区间包含。我一开始写了个O(nlog^2n)的cdq分治,判这个交的部分被包含次数是否超过2,后来dalao告诉我,直接线段树维护每个左端点对应的右端点最远在哪,每个左端点开一个堆支持删除就好了,询问的时候删掉两个区间后面再加回来,这样时间复杂度为O(nlogn)。还有为避免炸精度,直接用向量表示斜率。另外两个程序里都有许多卡常(cdq还是卡常过的),跑的飞快。

cdq:

#include<cstdio>
#include<algorithm>
using namespace std;
char B[<<],*S=B;
inline int read()
{
int x,f=;char c;
while((c=*S++)<''||c>'')if(c=='-')f=;
for(x=c-'';(c=*S++)>=''&&c<='';)x=x*+c-'';
return f?x:-x;
}
#define MN 500000
#define N 524288
struct vec
{
int x,y;
friend bool operator<(const vec&a,const vec&b){return 1LL*a.x*b.y<1LL*a.y*b.x;}
friend bool operator==(const vec&a,const vec&b){return 1LL*a.x*b.y==1LL*a.y*b.x;}
};
vec rot(const vec&a){return (vec){a.y,-a.x};}
struct work{int t,z;vec x,y;}w[MN+],q[MN+];
bool cmp(const work&a,const work&b){return a.x==b.x?a.t>b.t:a.x<b.x;}
vec l[MN+],r[MN+],c[MN+];
int cnt,ans[MN+],t[N*+],qn,sa[MN+],sb[MN+];
void add(int k,int x){for(k+=N;k;k>>=)t[k]+=x;}
int query(int l,int r)
{
int res=;
for(l+=N-,r+=N+;l^r^;l>>=,r>>=)
{
if(~l&)res+=t[l+];
if( r&)res+=t[r-];
}
return res;
}
void solve(int l,int r)
{
int mid=l+r>>,i;
if(l<r)solve(l,mid),solve(mid+,r);
if(sa[mid]==sa[l-]||sb[r]==sb[mid])return;
for(qn=,i=l;i<=mid;++i)if(w[i].t>)q[++qn]=w[i];
for(;i<=r;++i)if(!w[i].t)q[++qn]=w[i];
sort(q+,q+qn+,cmp);
for(i=;i<=qn;++i)
if(q[i].t)add(lower_bound(c+,c+cnt+,q[i].y)-c,q[i].z);
else ans[q[i].z]+=query(lower_bound(c+,c+cnt+,q[i].y)-c,cnt);
for(i=;i<=qn;++i)if(q[i].t)add(lower_bound(c+,c+cnt+,q[i].y)-c,-q[i].z);
}
int main()
{
B[fread(B,,<<,stdin)]=;
int R=read(),n=read(),i,t,x,y;
for(i=;i<=n;++i)
{
t=read();x=read();sa[i]=sa[i-];sb[i]=sb[i-];
if(t==)
{
l[++cnt]=(vec){x+R,y=read()};
c[cnt]=r[cnt]=(vec){x-R,y};
w[i]=(work){,,l[cnt],r[cnt]},++sa[i];
}
if(t==)w[i]=(work){,-,l[x],r[x]},++sa[i];
if(t==)
{
if(l[y=read()]<l[x])x^=y^=x^=y;
if(rot(r[x])<l[y])w[i]=(work){-,,,},ans[i]=;
else w[i]=(work){,i,l[y],r[x]<r[y]?r[x]:r[y]},++sb[i];
}
}
sort(c+,c+cnt+);cnt=unique(c+,c+cnt+)-c-;solve(,n);
for(i=;i<=n;++i)if(ans[i])puts(ans[i]>?"NO":"YES");
}

线段树:

#include<cstdio>
#include<algorithm>
#include<queue>
using namespace std;
inline int read()
{
int x,f=;char c;
while((c=getchar())<''||c>'')if(c=='-')f=;
for(x=c-'';(c=getchar())>=''&&c<='';)x=x*+c-'';
return f?x:-x;
}
#define MN 500000
#define N 524288
#define mp(x,y) make_pair(x,y)
struct vec
{
int x,y,z;
friend bool operator<(const vec&a,const vec&b)
{return a.z||(!b.z&&1LL*a.x*b.y<1LL*a.y*b.x);}
friend bool operator==(const vec&a,const vec&b){return 1LL*a.x*b.y==1LL*a.y*b.x;}
}l[MN+],r[MN+],c[MN+],T[N*+];
priority_queue< pair<vec,int> > p[MN+];
int cnt,u[MN+],t[MN+],x[MN+],y[MN+],lp[MN+];
void change(int k,const vec&x){for(T[k+=N]=x;k>>=;)T[k]=max(T[k<<],T[k<<|]);}
vec query(int l,int r)
{
vec res=(vec){,,};
for(l+=N-,r+=N+;l^r^;l>>=,r>>=)
{
if(~l&)res=max(res,T[l+]);
if( r&)res=max(res,T[r-]);
}
return res;
}
void ins(int x)
{
int k=lp[x]?lp[x]:lp[x]=lower_bound(c+,c+cnt+,l[x])-c;
u[x]=;p[k].push(mp(r[x],x));
if(p[k].top()==mp(r[x],x))change(k,r[x]);
}
void del(int x)
{
int k=lp[x]?lp[x]:lp[x]=lower_bound(c+,c+cnt+,l[x])-c,s=;
for(u[x]=;u[p[k].top().second];++s)p[k].pop();
if(s)change(k,p[k].top().first);
}
int main()
{
int R=read(),n=read(),i,k;
for(i=;i<=n;++i)
{
t[i]=read();x[i]=read();if(t[i]!=)y[i]=read();
if(t[i]==)++cnt,c[cnt]=l[cnt]=(vec){x[i]+R,y[i],},r[cnt]=(vec){y[i],R-x[i],},x[i]=cnt;
}
sort(c+,c+cnt+);cnt=unique(c+,c+cnt+)-c-;
for(i=;i<*N;++i)T[i]=(vec){,,};
for(i=;i<=cnt;++i)p[i].push(mp(T[],));
for(i=;i<=n;++i)
{
if(t[i]==)ins(x[i]);
if(t[i]==)del(x[i]);
if(t[i]==)
{
if(l[y[i]]<l[x[i]])swap(x[i],y[i]);
if(r[x[i]]<l[y[i]]){puts("NO");continue;}
k=lp[y[i]]?lp[y[i]]:lp[y[i]]=lower_bound(c+,c+cnt+,l[y[i]])-c;
del(x[i]);del(y[i]);
puts(query(,k)<(r[x[i]]<r[y[i]]?r[x[i]]:r[y[i]])?"YES":"NO");
ins(x[i]);ins(y[i]);
}
}
}

[Codeforces]856E - Satellites的更多相关文章

  1. python爬虫学习(5) —— 扒一下codeforces题面

    上一次我们拿学校的URP做了个小小的demo.... 其实我们还可以把每个学生的证件照爬下来做成一个证件照校花校草评比 另外也可以写一个物理实验自动选课... 但是出于多种原因,,还是绕开这些敏感话题 ...

  2. 【Codeforces 738D】Sea Battle(贪心)

    http://codeforces.com/contest/738/problem/D Galya is playing one-dimensional Sea Battle on a 1 × n g ...

  3. 【Codeforces 738C】Road to Cinema

    http://codeforces.com/contest/738/problem/C Vasya is currently at a car rental service, and he wants ...

  4. 【Codeforces 738A】Interview with Oleg

    http://codeforces.com/contest/738/problem/A Polycarp has interviewed Oleg and has written the interv ...

  5. CodeForces - 662A Gambling Nim

    http://codeforces.com/problemset/problem/662/A 题目大意: 给定n(n <= 500000)张卡片,每张卡片的两个面都写有数字,每个面都有0.5的概 ...

  6. CodeForces - 274B Zero Tree

    http://codeforces.com/problemset/problem/274/B 题目大意: 给定你一颗树,每个点上有权值. 现在你每次取出这颗树的一颗子树(即点集和边集均是原图的子集的连 ...

  7. CodeForces - 261B Maxim and Restaurant

    http://codeforces.com/problemset/problem/261/B 题目大意:给定n个数a1-an(n<=50,ai<=50),随机打乱后,记Si=a1+a2+a ...

  8. CodeForces - 696B Puzzles

    http://codeforces.com/problemset/problem/696/B 题目大意: 这是一颗有n个点的树,你从根开始游走,每当你第一次到达一个点时,把这个点的权记为(你已经到过不 ...

  9. CodeForces - 148D Bag of mice

    http://codeforces.com/problemset/problem/148/D 题目大意: 原来袋子里有w只白鼠和b只黑鼠 龙和王妃轮流从袋子里抓老鼠.谁先抓到白色老鼠谁就赢. 王妃每次 ...

随机推荐

  1. beta冲刺7-咸鱼

    前言:最后一篇惹.明天就是正式交差了.有点慌-- 昨天的未完成: 用户试用+测评 输入部分的正则式判定 今天的工作: 登陆界面修改 我的社团显示效果优化 部分信息注册后锁定无法修改 其他部分功能优化 ...

  2. Beta 第三天

    今天遇到的困难: 组员对github极度的不适应 github的版本控制和协同化编程确实操作起来需要一定的熟练度,我们缺乏这种熟练度 Android Studio版本不一致项目难以打开的问题仍然无法解 ...

  3. 20155306 2017-2018-1《信息安全系统设计》第二周课堂测试以及myod的实现

    20155306 2017-2018-1<信息安全系统设计>第二周课堂测试以及myod的实现 第二周课堂测验: (注:前两项在课堂已提交,在此不做详解) 第一项: 每个.c一个文件,每个. ...

  4. 冲刺NO.8

    Alpha冲刺第八天 站立式会议 项目进展 项目稳步进行,项目的基础部分如基本信息管理,信用信息管理等部分已相对比较完善. 问题困难 技术困难在短期内很难发生质的变化,而本项目由于选择了队员不太熟悉的 ...

  5. css变化代码2

    <!DOCTYPE html><html>    <head>        <meta charset="utf-8" />   ...

  6. 织梦cms/dedecms清理冗余废弃未引用图片方法

    原理描述: 在原有织梦后台菜单中增加"清理冗余图片按钮",实现清理冗余图片的功能. 操作步骤: 1. 打开后台dede\sys_sql_query.php代码 在该文件中搜索如下代 ...

  7. python 常用算法学习(2)

    一,算法定义 算法(Algorithm)是指解题方案的准确而完整的描述,是一系列解决问题的清晰指令,算法代表着用系统的方法描述解决问题的策略机制.也就是说,能够对一定规范的输入,在有限时间内获得所要求 ...

  8. JAVA_SE基础——4.path的临时配置&Classpath的配置

    这次,我来写下关于path的临时配置的心的 我来说个有可能的实例:如果你去到别人的电脑 又想写代码 又不想改乱别人的path配置的话  再说别人愿意你在别人的电脑上瞎配吗? 那该怎么办呢? 那没问题 ...

  9. 自定义SpringBoot启动banner

    序: springboot启动的时候会有一个启动logo似的东西,如图,这个logo似的东西叫做banner,本文小计修改此banner显示与关闭banner.没什么用,有兴趣可以玩玩-- 正文: 自 ...

  10. 阿里云API网关(17)签名算法

    网关指南: https://help.aliyun.com/document_detail/29487.html?spm=5176.doc48835.6.550.23Oqbl 网关控制台: https ...