【BZOJ4445】[SCOI2015]小凸想跑步(半平面交)

题面

BZOJ

洛谷

题解

首先把点给设出来,\(A(x_a,y_a),B(x_b,y_b),C(x_c,y_c),D(x_d,y_d),P(x,y)\)

然后我们考虑\(S_\Delta ABP<S_\Delta CDP\)什么情况下满足。

根据点积来求面积,得到:

\[(x_a-x,y_a-y)\times(x_b-x,y_b-y)<(x_c-x,y_c-y)\times(x_d-x,y_d-y)
\]

这个东西左边拆开之后得到:

\[\begin{aligned}
&\ \ \ \ (x_a-x)(y_b-y)-(x_b-x)(y_a-y)\\
&=x(y_a-y_b)-y(x_a-x_b)+x_ay_b-x_by_a
\end{aligned}\]

右侧类似。然后就可以移项,得到:

\[x(y_a-y_b-y_c+y_d)-y(x_a-x_b-x_c+x_d)+x_ay_b-x_by_a-x_cy_d+x_dy_c<0
\]

那么对于相邻的两个点和\(0,1\)两个点进行一次比较,每次都可以得到一个半平面,最后求解这个半平面交就是结果了。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
#define MAX 100100
const double eps=1e-7;
inline int read()
{
int x=0;bool t=false;char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')t=true,ch=getchar();
while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
return t?-x:x;
}
struct Node{double x,y;}p[MAX],Qp[MAX<<1];
Node operator+(Node a,Node b){return (Node){a.x+b.x,a.y+b.y};}
Node operator-(Node a,Node b){return (Node){a.x-b.x,a.y-b.y};}
Node operator*(Node a,double b){return (Node){a.x*b,a.y*b};}
double operator*(Node a,Node b){return a.x*b.x+a.y*b.y;}
double cross(Node a,Node b){return a.x*b.y-a.y*b.x;}
struct Line{Node p,v;double alpha;}S[MAX<<1],Q[MAX<<1];int tot;
bool cmp(Line a,Line b){return a.alpha<b.alpha;}
Node Intersection(Line a,Line b)
{
Node c=b.p-a.p;
double t=cross(b.v,c)/cross(b.v,a.v);
return a.p+a.v*t;
}
int zero(double x)
{
if(fabs(x)<eps)return 0;
return x>0?1:-1;
}
bool Left(Node a,Line b){return zero(cross(b.v,a-b.p))>0;}
int n;double area;
int l,r;
bool HalfPlaneIntersection()
{
for(int i=1;i<=tot;++i)S[i].alpha=atan2(S[i].v.y,S[i].v.x);
sort(&S[1],&S[tot+1],cmp);
Q[l=r=1]=S[1];
for(int i=2;i<=tot;++i)
{
while(l<r&&!Left(Qp[r-1],S[i]))--r;
while(l<r&&!Left(Qp[l],S[i]))++l;
if(zero(cross(Q[r].v,S[i].v))==0)
Q[r]=Left(Q[r].p,S[i])?Q[r]:S[i];
else Q[++r]=S[i];
if(l<r)Qp[r-1]=Intersection(Q[r],Q[r-1]);
}
while(l<r&&!Left(Qp[r-1],Q[l]))--r;
return (r-l)>1;
}
int main()
{
n=read();
for(int i=0;i<n;++i)p[i].x=read(),p[i].y=read();p[n]=p[0];
for(int i=0;i<n;++i)S[++tot]=(Line){p[i],p[i+1]-p[i]};
for(int i=1;i<n;++i)area+=fabs(cross(p[i]-p[0],p[i+1]-p[0]));
for(int i=1;i<n;++i)
{
double a=p[0].y-p[1].y-p[i].y+p[i+1].y;
double b=-p[0].x+p[1].x+p[i].x-p[i+1].x;
double c=cross(p[0],p[1])-cross(p[i],p[i+1]);
Line d;
if(fabs(a)>eps)d=(Line){(Node){-c/a,0},(Node){-b,a}};
else d=(Line){(Node){0,-c/b},(Node){-b,a}};
S[++tot]=d;
}
HalfPlaneIntersection();
double ans=0;Qp[r]=Intersection(Q[l],Q[r]);
for(int i=l+1;i<r;++i)ans+=fabs(cross(Qp[i]-Qp[l],Qp[i+1]-Qp[l]));
printf("%.4lf\n",ans/area);
return 0;
}

【BZOJ4445】[SCOI2015]小凸想跑步(半平面交)的更多相关文章

  1. 【BZOJ4445】[Scoi2015]小凸想跑步 半平面交

    [BZOJ4445][Scoi2015]小凸想跑步 Description 小凸晚上喜欢到操场跑步,今天他跑完两圈之后,他玩起了这样一个游戏. 操场是个凸n边形,N个顶点按照逆时针从0-n-l编号.现 ...

  2. 4445: [Scoi2015]小凸想跑步 半平面交

    题目大意: http://www.lydsy.com/JudgeOnline/problem.php?id=4445 题解: 设点坐标,利用叉积可以解出当p坐标为\((x_p,y_p)\)时,与边i- ...

  3. bzoj 4445 小凸想跑步 - 半平面交

    题目传送门 vjudge的快速通道 bzoj的快速通道 题目大意 问在一个凸多边形内找一个点,连接这个点和所有顶点,使得与0号顶点,1号顶点构成的三角形是最小的概率. 假设点的位置是$(x, y)$, ...

  4. BZOJ4445: [Scoi2015]小凸想跑步

    裸半平面交. 记得把P0P1表示的半平面加进去,否则点可能在多边形外. #include<bits/stdc++.h> #define N 100009 using namespace s ...

  5. [bzoj4445] [SCOI2015]小凸想跑步 (半平面交)

    题意:凸包上一个点\(p\),使得\(p\)和点\(0,1\)组成的三角形面积最小 用叉积来求: \(p,i,i+1\)组成的三角形面积为: (\(\times\)为叉积) \((p_p-i)\tim ...

  6. BZOJ4445 SCOI2015小凸想跑步(半平面交)

    考虑怎样的点满足条件.设其为(xp,yp),则要满足(x0-xp,y0-yp)×(x1-xp,y1-yp)<=(xi-xp,yi-yp)×(xi+1-xp,yi+1-yp)对任意i成立.拆开式子 ...

  7. 2018.10.15 bzoj4445: [Scoi2015]小凸想跑步(半平面交)

    传送门 话说去年的省选计算几何难度跟前几年比起来根本不能做啊(虽然去年考的时候并没有学过计算几何) 这题就是推个式子然后上半平面交就做完了. 什么? 怎么推式子? 先把题目的概率转换成求出可行区域. ...

  8. BZOJ 4445 [Scoi2015]小凸想跑步:半平面交

    传送门 题意 小凸晚上喜欢到操场跑步,今天他跑完两圈之后,他玩起了这样一个游戏. 操场是个凸 $ n $ 边形,$ n $ 个顶点 $ P_i $ 按照逆时针从 $ 0 $ 至 $ n-1 $ 编号. ...

  9. [SCOI2015]小凸想跑步

    题目描述 小凸晚上喜欢到操场跑步,今天他跑完两圈之后,他玩起了这样一个游戏. 操场是个凸 n 边形, nn 个顶点按照逆时针从 0 ∼n−1 编号.现在小凸随机站在操场中的某个位置,标记为p点.将 p ...

随机推荐

  1. linux内核与分析 心得与体会

    作业目录: (1)计算机是如何工作的:http://www.cnblogs.com/20135335hs/p/5213394.html (2)操作系统是如何工作的:http://www.cnblogs ...

  2. Redis之字符串

    Redis 字符串 String 是redis最基本的类型,value 不仅可以是 String,也可以是数字. 使用 Strings 类型,可以完全实现目前 Memcached 的功能,并且效率更高 ...

  3. Gym 101147A The game of Osho(SG找规律)

    https://vjudge.net/problem/Gym-101147A 题意:给出G组数,每组数包括两个数B,N,两玩家轮流取数,使得N-num,num<=N并且num是N的整次幂.判断谁 ...

  4. UVa 11248 网络扩容(最大流(需要优化))

    https://vjudge.net/problem/UVA-11248 题意: 给定一个有向网络,每条边均有一个容量.问是否存在一个从点1到点N,流量为C的流.如果不存在,是否可以恰好修改一条弧的容 ...

  5. Linux中重定向--转载

    转:http://blog.csdn.net/songyang516/article/details/6758256 1重定向 1.1      重定向符号 >               输出 ...

  6. 爬虫模拟登陆之formdata表单数据

    首先HTTP协议是个无连接的协议,浏览器和服务器之间是以循环往复的请求回复来交互的,交互的形式是以文件形式来进行的.比如在chrome开发者工具network中看到了 每一行是一个文件,又文件大小啊, ...

  7. hiho一下 第二周 trie树

    Trie树 时间限制:10000ms 单点时限:1000ms 内存限制:256MB 描述 小Hi和小Ho是一对好朋友,出生在信息化社会的他们对编程产生了莫大的兴趣,他们约定好互相帮助,在编程的学习道路 ...

  8. [转]Cryengine渲染引擎剖析

    转篇Napoleon314 大牛的分析,排版好乱,见谅,耐心读,这是个好东西,注意看他自己的实现,是个技术狂人啊,Ogre焕发次时代的光芒啊~~~努力 ------------------------ ...

  9. maven3常用命令、java项目搭建、web项目搭建

    ------------------------------maven3常用命令--------------------------- 1.常用命令 1)创建一个Project mvn archety ...

  10. curl使用记录

    $header = array("Connection: Keep-Alive", "Accept: text/html,application/xhtml+xml,ap ...