Description

给出一个N个顶点的简单多边形,对于每个顶点,假如它和原点连成的线段只在这个顶点处和多边形相交,就称为满足要求的顶点。你的任务是输出所有满足要求的顶点编号。

Input

第一行一个正整数N。下面N行每行两个不超过106的正整数,依次表示每个顶点的坐标。顶点按照输入顺序用正整数1..N编号,并且顶点保证按照顺时针或逆时针顺序给出。

Output

第一行一个正整数M,表示满足要求的顶点个数。第二行M个正整数,按照升序给出满足要求的顶点编号。
 
对所有点极角排序并离散化,特判极角相同的情况,建线段树,将所有边插到线段树上,最后查询每个点方向的射线与多边形最近的交点。
按极角做扫描线也可以在相同复杂度内实现,但细节可能更多。
#include<bits/stdc++.h>
typedef long long i64;
int min(int a,int b){return a<b?a:b;}
int max(int a,int b){return a>b?a:b;}
const int N=;
char buf[],*ptr=buf;
int _(){
int x=;
while(*ptr<)++ptr;
while(*ptr>)x=x*+*ptr++-;
return x;
}
struct pos{
double x,y;
double val()const{return x+y;}
}xs[N];
double operator*(pos a,pos b){return a.x*b.y-a.y*b.x;}
pos operator-(pos a,pos b){return (pos){a.x-b.x,a.y-b.y};}
int n;
struct vertex{
pos p;
int id,x;
bool operator<(vertex w)const{
double d=p*w.p;
return d?d<:p.val()<w.p.val();
}
}vs[N];
int ps[N],xp=;
bool un[N];
struct itv{
pos a,b;
int l,r;
double operator()(pos&w)const{
return w.val()*(b*a)/(w*(a-b));
}
}is[];
int mx;
void ins(itv&a,itv&b){
if(!a.r)a=b;
else{
int l=max(a.l,b.l),r=min(a.r,b.r);
if(a(xs[l])+a(xs[r])>b(xs[l])+b(xs[r]))a=b;
}
}
void ins(vertex a,vertex b){
itv w=(itv){a.p,b.p,a.x,b.x};
for(int l=a.x+mx-,r=b.x+mx+;r-l!=;l>>=,r>>=){
if(~l&)ins(is[l^],w);
if(r&)ins(is[r^],w);
}
}
int as[N],ap=;
int main(){
fread(buf,,sizeof(buf),stdin);
n=_();
for(mx=;mx<=n+;mx<<=);
for(int i=;i<=n;++i){
vs[i].p.x=_();
vs[i].p.y=_();
vs[i].id=i;
}
std::sort(vs+,vs+n+);
for(int i=,j=;i<=n;i=j){
for(++j;j<=n&&vs[i].p*vs[j].p==;un[vs[j].id]=,++j);
xs[++xp]=vs[i].p;
for(int k=i;k<j;++k)vs[k].x=xp;
}
for(int i=;i<=n;++i)ps[vs[i].id]=i;
ps[]=ps[n];
for(int i=;i<=n;++i){
int a=ps[i-],b=ps[i];
if(a>b)std::swap(a,b);
if(vs[a].x<vs[b].x)ins(vs[a],vs[b]);
}
for(int i=;i<=n;++i)if(!un[vs[i].id]){
for(int w=vs[i].x+mx;w;w>>=)if(is[w].r&&is[w](vs[i].p)<vs[i].p.val()-1e-){
un[vs[i].id]=;
break;
}
}
for(int i=;i<=n;++i)if(!un[i])as[ap++]=i;
printf("%d\n",ap);
for(int i=;i<ap;++i)printf("%d ",as[i]);
return ;
}

bzoj2856: [ceoi2012]Printed Circuit Board的更多相关文章

  1. Printed Circuit Board (board)

    Printed Circuit Board (board) 题目描述 给出一个N个顶点的简单多边形,对于每个顶点,假如它和原点连成的线段只在这个顶点处和多边形相交,就称为满足要求的顶点.你的任务是输出 ...

  2. ZOJ1648 Circuit Board 2017-04-18 20:31 34人阅读 评论(0) 收藏

    Circuit Board Time Limit: 2 Seconds      Memory Limit: 65536 KB On the circuit board, there are lots ...

  3. ZOJ 1648 Circuit Board(计算几何)

    Circuit Board Time Limit: 2 Seconds Memory Limit: 65536 KB On the circuit board, there are lots of c ...

  4. ZOJ1648 Circuit Board(线段相交)

    裸的判断线段相交

  5. kickstart2019 round_C B. Circuit Board

    思路: 这题应该不止一种解法,其中的一种可以看作是leetcode85https://www.cnblogs.com/wangyiming/p/11059176.html的加强版: 首先对于每一行,分 ...

  6. zoj 1648 Circuit Board

    题目:意思就是推断给定的几条线段是否有相交的. 方法:模版吧,有空在来细细学习. 代码: #include <iostream> #include <cstdio> using ...

  7. bzoj AC倒序

    Search GO 说明:输入题号直接进入相应题目,如需搜索含数字的题目,请在关键词前加单引号 Problem ID Title Source AC Submit Y 1000 A+B Problem ...

  8. COB(Chip On Board) 工艺技术

    COX(Chip On X) •X 基板:  PCB (Printed circuit board)  FPC (Flexible Printed Circuit)  Glass •导线焊接 球形焊接 ...

  9. 特邀美国EMC实战专家Mark来华授课

    “轻松搞定EMC-PCB和系统设计”课程介绍 本次课程特邀美国EMC领域权威专家Mark Montrose主讲,将涵盖满足产品电磁兼容性和信号完整性的基本原理.课程涉及多个领域,不仅仅针对PCB设计, ...

随机推荐

  1. php 中使用正则

    1.匹配一个由一个小写字母和一位数字组成的字符串,比如”z2″   用^[a-z][0-9]$ 2.当在一组方括号里使用^是,它表示“非”或“排除”的意思   比如要求第一个字符不能是数字:^[^0- ...

  2. 关于 transparent rgba display:none; opacity visiblity 关于em

    关于 transparent  rgba   display:none; opacity   visiblity   display 之后不会占位. 其余都会占位 opacity 还会继承,子元素也会 ...

  3. awk 相关的复习

    1. awk 引用外部变量: aa=666  echo "." | awk -v GET_A=$aa '{print GET_A}' . sort -n fuxi.awk |awk ...

  4. PTA——最大公约数和最小公倍数

    PTA 7-26 最大公约数和最小公倍数 #include<stdio.h> int main(){ int num1,num2,temp1,temp2,r; scanf("%d ...

  5. 【liunx】Linux下的压缩和解压缩命令——jar

    原文链接:http://blog.chinaunix.net/uid-692788-id-2681136.html JAR包是Java中所特有一种压缩文档,其实大家就可以把它理解为.zip包.当然也是 ...

  6. WebService 学习记录

    -------------------------------------------PS:这个WebService 服务必须一直开着,关闭就没法访问了 Web Service 教程 一.webser ...

  7. 手动制作openstack CentOS 镜像

    https://docs.openstack.org/image-guide/centos-image.html This example shows you how to install a Cen ...

  8. 通过阅读python subprocess源码尝试实现非阻塞读取stdout以及非阻塞wait

    http://blog.chinaunix.net/uid-23504396-id-4661783.html 执行subprocess的时候,执行不是问题最麻烦的是获取进程执行后的回显来确认是否正确执 ...

  9. makeinfo is missing on your system(转)

    ubunut14.04 make install 提示 makeinfo is missing on your system: 输入makeinfo后,提示没有安装该命令,然后提示可以安装texinf ...

  10. java 对象的创建

    jvm在创建对象时,首先判断该对象所对应的类是否已经被加载.链接和初始化,如果没有则先执行类的加载过程.类的加载过程见 java类加载. 类加载检查通过后,虚拟机在堆中为新生对象分配内存. 内存分配完 ...