ural1469
题解:
从左往右加入每一个点
判断一下和,pre,nxt是否相交
删除得时候也要判断
代码:
#pragma GCC optimize(2)
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
const int N=;
typedef double ld;
const ld eps=1e-,C=0.18479956785822313167427314019291;
struct nd{int id,ls,rs,rk,fa;ld key;}t[N];
struct se{ld x1,y1,x2,y2,k,b;}s[N];
struct Tp{ld sg,x;int id;}p[N];
int n,m,rt,date[N];
int cmp(const Tp &a,const Tp &b){return (a.x<b.x);}
void rotate(ld &x,ld &y)
{
ld p=x,q=y;
x=cos(C)*p-sin(C)*q;
y=sin(C)*p+cos(C)*q;
}
int zig(int d)
{
int tmp=t[d].ls;
t[d].ls=t[tmp].rs;t[tmp].rs=d;
t[tmp].fa=t[d].fa;
t[d].fa=tmp;
if (t[d].ls) t[t[d].ls].fa=d;
return tmp;
}
int zag(int d)
{
int tmp=t[d].rs;
t[d].rs=t[tmp].ls;t[tmp].ls=d;
t[tmp].fa=t[d].fa;
t[d].fa=tmp;
if (t[d].rs) t[t[d].rs].fa=d;
return tmp;
}
int ins(int d,ld y,ld x,int id)
{
if (!d)
{
date[id]=++m;
t[m].key=y;t[m].id=id;
t[m].fa=;t[m].rk=rand();
return m;
}
ld _y=s[t[d].id].k*x+s[t[d].id].b;
if (fabs(y-_y)<=eps)
{
printf("YES\n%d %d\n",id+,t[d].id+);
exit();
}
if (y<_y)
{
t[d].ls=ins(t[d].ls,y,x,id);
t[t[d].ls].fa=d;
if (t[t[d].ls].rk<t[d].rk)return zig(d);
else return d;
}
t[d].rs=ins(t[d].rs,y,x,id);
t[t[d].rs].fa=d;
if (t[t[d].rs].rk<t[d].rk)return zag(d);
else return d;
}
int del(int d,ld y,ld x)
{
ld _y=s[t[d].id].k*x+s[t[d].id].b;
if (fabs(y-_y)<=eps)
{
if (!(t[d].ls+t[d].rs)) return ;
if ((t[d].ls)&&((!t[d].rs)||(t[t[d].ls].rk<t[t[d].rs].rk)))d=zig(d);
else d=zag(d);
return del(d,y,x);
}
if (y<_y)t[d].ls=del(t[d].ls,y,x);
else t[d].rs=del(t[d].rs,y,x);
return d;
}
int pre(int d)
{
if (t[d].ls)
{
for (d=t[d].ls;t[d].rs;d=t[d].rs);
return d;
}
while (d)
{
if (t[t[d].fa].rs==d) return t[d].fa;
d=t[d].fa;
}
return -;
}
int nxt(int d)
{
if (t[d].rs)
{
for (d=t[d].rs;t[d].ls;d=t[d].ls);
return d;
}
while (d)
{
if (t[t[d].fa].ls==d) return t[d].fa;
d=t[d].fa;
}
return -;
}
int det(ld x1,ld y1,ld x2,ld y2,ld x3,ld y3)
{
ld res=(x2-x1)*(y3-y1)-(x3-x1)*(y2-y1);
if (fabs(res)<=eps) return ;
if (res>) return ;
else return -;
}
int cr(int a,int b)
{
if ((a==-)||(b==-)) return ;
a=t[a].id;b=t[b].id;
if ((s[a].x1+eps>=min(s[b].x1,s[b].x2))&&(s[a].x1-eps<=max(s[b].x1,s[b].x2))
&&(!det(s[b].x1,s[b].y1,s[b].x2,s[b].y2,s[a].x1,s[a].y1)))return ;
if ((s[a].x2+eps>=min(s[b].x1,s[b].x2))&&(s[a].x2-eps<=max(s[b].x1,s[b].x2))
&&(!det(s[b].x1,s[b].y1,s[b].x2,s[b].y2,s[a].x2,s[a].y2)))return ;
if ((s[b].x1+eps>=min(s[a].x1,s[a].x2))&&(s[b].x1-eps<=max(s[a].x1,s[a].x2))
&&(!det(s[a].x1,s[a].y1,s[a].x2,s[a].y2,s[b].x1,s[b].y1)))return ;
if ((s[b].x2+eps>=min(s[a].x1,s[a].x2))&&(s[b].x2-eps<=max(s[a].x1,s[a].x2))
&&(!det(s[a].x1,s[a].y1,s[a].x2,s[a].y2,s[b].x2,s[b].y2)))return ;
if ((det(s[a].x1,s[a].y1,s[a].x2,s[a].y2,s[b].x1,s[b].y1)*
det(s[a].x1,s[a].y1,s[a].x2,s[a].y2,s[b].x2,s[b].y2)<)&&
(det(s[b].x1,s[b].y1,s[b].x2,s[b].y2,s[a].x1,s[a].y1)*
det(s[b].x1,s[b].y1,s[b].x2,s[b].y2,s[a].x2,s[a].y2)<))return ;
else return ;
}
int main()
{
scanf("%d",&n);
for (int i=;i<n;i++)
{
scanf("%lf%lf%lf%lf",&s[i].x1,&s[i].y1,&s[i].x2,&s[i].y2);
rotate(s[i].x1,s[i].y1);rotate(s[i].x2,s[i].y2);
if (s[i].x1>s[i].x2) swap(s[i].x1,s[i].x2),swap(s[i].y1,s[i].y2);
s[i].k=(s[i].y2-s[i].y1)/(s[i].x2-s[i].x1);
s[i].b=s[i].y1-s[i].k*s[i].x1;
p[i].id=p[i+n].id=i;
p[i].sg=;p[i+n].sg=;
p[i].x=s[i].x1;p[i+n].x=s[i].x2;
}
sort(p,p+n*,cmp);
for (int i=;i<n*;i++)
if (p[i].sg)
{
rt=ins(rt,s[p[i].id].y1,p[i].x,p[i].id);
if (cr(date[p[i].id],pre(m)))
{
printf("YES\n%d %d",p[i].id+,t[pre(m)].id+);
return ;
}
if (cr(date[p[i].id],nxt(m)))
{
printf("YES\n%d %d",p[i].id+,t[nxt(m)].id+);
return ;
}
}
else
{
if (cr(pre(date[p[i].id]),nxt(date[p[i].id])))
{
printf("YES\n%d %d",t[pre(date[p[i].id])].id+,
t[nxt(date[p[i].id])].id+);
return ;
}
rt=del(rt,s[p[i].id].y2,p[i].x);
}
puts("NO");
return ;
}
ural1469的更多相关文章
随机推荐
- UVA 10462 Is There A Second Way Left?(次小生成树&Prim&Kruskal)题解
思路: Prim: 这道题目中有重边 Prim可以先加一个sec数组来保存重边的次小边,这样不会影响到最小生成树,在算次小生成树时要同时判断次小边(不需判断是否在MST中) Kruskal: Krus ...
- Sql Server 创建表添加说明
http://bbs.csdn.net/topics/340184487 在此感谢 提供参考 CREATE TABLE ToPayFee ( Id INT IDENTITY(1,1) PRIMA ...
- ovs-ofctl: s1 is not a bridge or a socket 解决方法
参考: ovs-vsctl: Error detected while setting up bridge ovs-ofctl: s1 is not a bridge or a socket 解决方法 ...
- Jenkins+Ant/maven+Svn实现自动化部署,编译,运行,测试结果自动邮件通知
Jenkins+Ant+Svn实现自动化部署,编译,运行,测试结果自动邮件通知
- jquery 封装插件
如今做web开发,jquery 几乎是必不可少的,就连vs神器在2010版本开始将Jquery 及ui 内置web项目里了.至于使用jquery好处这里就不再赘述了,用过的都知道.今天我们来讨论下jq ...
- 微信小程序之可滚动视图 scroll-view 的使用注意
微信小程序之可滚动视图 scroll-view 使用注意: 1.scroll-view 中的需要滑动的元素不可以用 float 浮动: 2.scroll-view 中的包裹需要滑动的元素的大盒子用 d ...
- c++ 指定长度容器元素的拷贝(copy_n)
#include <iostream> // cout #include <algorithm> // copy #include <vector> ...
- testNG 学习笔记 Day 3 常用的断言
TestNG中最常用的一个断言类是Assert.java,里面有多个静态方法,这个类我们习惯叫硬断言.对应的还有一个软断言的类,叫SoftAssert.java,这个类是需要创建实例对象,才能调用相关 ...
- zookeeper篇-如何修改源码
提一个问题先 zxid有64位,分成两部分: 高32位是Leader的epoch:选举时钟,每次选出新的Leader,epoch累加1 低32位是在这轮epoch内的事务id:对于用户的每一次更新操作 ...
- Cglib方法实现动态代理
除了使用JDK方式产生动态代理外,Java还给我们提供了另外一种产生动态代理的方法,那就是使用cglib. cglib是这样实现动态代理的: · ①.针对类来实现代理 · ②对指定目标类产生一个子类 ...