很明显是一道半平面交的题。

先说一下半平面交的步骤:

1.用点向法(点+向量)表示直线

2.极角排序,若极角相同,按相对位置排序。

3.去重,极角相同的保留更优的

4.枚举边维护双端队列

5.求答案

1就不说了,2中的极角可以用atan2(y,x)来求,因为atan2精度要高



双端队列的原因是新加的一条边对头和尾都有影响,如图:

如何去判断:只要判断线head和线head+1,的交点p与新的一条线的位置关系就可以

至于交点的求法:先见图:

求\(p_1v_1,p_2v_2\)的交点\(p_0\)

设\(p_0=p_2+kv_2\ \ u=p_2-p_1\)

\(S_1=u\times v_1,S_2=v_1\times v_2,k=S_1/S2\)

所以\(p_0=p_2+kv_2\)

\(S_1\)为\(u\)与\(v_1\)的面积,\(S_2\)为\(v_1\)与\(v_2\)的面积,按比例求得\(k\)再乘一下就求出\(p_0\)

最后统计答案

细节见代码:

#include<iostream>
#include<cstdio>
#include<string>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<queue>
#include<stack>
#include<set>
#include<bitset>
#include<vector>
#include<cstdlib>
#define QAQ int
#define TAT long long
#define OwO bool
#define ORZ double
#define F(i,j,n) for(QAQ i=j;i<=n;++i)
#define E(i,j,n) for(QAQ i=j;i>=n;--i)
#define MES(i,j) memset(i,j,sizeof(i))
#define MEC(i,j) memcpy(i,j,sizeof(j)) using namespace std;
const int N=505;
const double eps=1e-8; int n;
struct Point{
double x,y; friend Point operator - (Point a,Point b){
Point t;
t.x=a.x-b.x;t.y=a.y-b.y;
return t;
}
friend Point operator + (Point a,Point b){
Point t;
t.x=a.x+b.x;t.y=a.y+b.y;
return t;
}
friend double operator * (Point a,Point b){
return a.x*b.x+a.y*b.y;
}
friend double operator ^ (Point a,Point b){
return a.x*b.y-b.x*a.y;
}
friend Point operator * (double k,Point b){
Point t;
t.x=k*b.x;t.y=k*b.y;
return t;
}
}b[N];
int sign(double x){
return fabs(x)<=eps ? 0 : (x>0 ? 1 : -1);
}
struct Line{
Point p,v;
double poa;
friend OwO operator < (Line x,Line y){
return sign(x.poa-y.poa)==0 ? sign((x.v-x.p) ^ (y.v-x.p)) >0 : sign(x.poa-y.poa)<0;
//因为我是向量左侧求交,所以极角相同时靠左的更优,把优的放在后面,方便之后的操作 ,可以画图体会一下
}
}a[N],q[N];
int js,cnt,head,tail;
double ans; Point inter(Line a,Line b){//求交点
Point p1=a.p,p2=b.p,v1=a.v,v2=b.v;
v1=v1-p1;v2=v2-p2;
Point u=p2-p1;
Point p=p2+((u^v1)/(v1^v2))*v2;
return p;
} OwO pd(Line i,Line j,Line k){
Point p=inter(i,j);
return sign((k.v-k.p) ^ (p-k.p))<0;
} void Half_Plane(){
sort(a+1,a+js+1);//排序
F(i,1,js) {
if(sign(a[i].poa-a[i-1].poa)!=0) cnt++;
a[cnt]=a[i];//因为排过序,即使极角相同,后面的也比前面的优
}
head=1;tail=0;
q[++tail]=a[1];q[++tail]=a[2];
F(i,3,cnt){
while(head<tail&&pd(q[tail-1],q[tail],a[i])) tail--;//维护双端队列
while(head<tail&&pd(q[head+1],q[head],a[i])) head++;
q[++tail]=a[i];
}
while(head<tail&&pd(q[tail-1],q[tail],q[head])) tail--;
while(head<tail&&pd(q[head+1],q[head],q[tail])) head++;
q[tail+1]=q[head];
js=0;
F(i,head,tail) b[++js]=inter(q[i],q[i+1]);
} int main(){
scanf("%d",&n);
F(i,1,n){
int k;
scanf("%d",&k);
F(j,1,k) scanf("%lf%lf",&b[j].x,&b[j].y);
b[k+1]=b[1];
F(j,1,k) a[++js].p=b[j],a[js].v=b[j+1];
}
F(i,1,js) a[i].poa=atan2(a[i].v.y-a[i].p.y,a[i].v.x-a[i].p.x);
Half_Plane();
b[js+1]=b[1];
if(js>2) F(i,1,js) ans+=(b[i]^b[i+1]);
ans=fabs(ans)/2.0;
printf("%.3lf\n",ans);
return 0;
}

[CQOI2006]凸多边形的更多相关文章

  1. 【BZOJ 2618】 2618: [Cqoi2006]凸多边形 (半平面交)

    2618: [Cqoi2006]凸多边形 Description 逆时针给出n个凸多边形的顶点坐标,求它们交的面积.例如n=2时,两个凸多边形如下图: 则相交部分的面积为5.233. Input 第一 ...

  2. bzoj 2618 2618: [Cqoi2006]凸多边形(半平面交)

    2618: [Cqoi2006]凸多边形 Time Limit: 5 Sec  Memory Limit: 128 MBSubmit: 656  Solved: 340[Submit][Status] ...

  3. bzoj 2618: [Cqoi2006]凸多边形 [半平面交]

    2618: [Cqoi2006]凸多边形 半平面交 注意一开始多边形边界不要太大... #include <iostream> #include <cstdio> #inclu ...

  4. 【BZOJ2618】[CQOI2006]凸多边形(半平面交)

    [BZOJ2618][CQOI2006]凸多边形(半平面交) 题面 BZOJ 洛谷 题解 这个东西就是要求凸多边形的边所形成的半平面交. 那么就是一个半平面交模板题了. 这里写的是平方的做法. #in ...

  5. 2018.07.04 BZOJ 2618 Cqoi2006凸多边形(半平面交)

    2618: [Cqoi2006]凸多边形 Time Limit: 5 Sec Memory Limit: 128 MB Description 逆时针给出n个凸多边形的顶点坐标,求它们交的面积.例如n ...

  6. 洛谷 P4196 [CQOI2006]凸多边形 (半平面交)

    题目链接:P4196 [CQOI2006]凸多边形 题意 给定 \(n\) 个凸多边形,求它们相交的面积. 思路 半平面交 半平面交的模板题. 代码 #include <bits/stdc++. ...

  7. bzoj2618: [Cqoi2006]凸多边形

    Description 逆时针给出n个凸多边形的顶点坐标,求它们交的面积.例如n=2时,两个凸多边形如下图: 则相交部分的面积为5.233. Input 第一行有一个整数n,表示凸多边形的个数,以下依 ...

  8. BZOJ2618[Cqoi2006]凸多边形——半平面交

    题目描述 逆时针给出n个凸多边形的顶点坐标,求它们交的面积.例如n=2时,两个凸多边形如下图: 则相交部分的面积为5.233. 输入 第一行有一个整数n,表示凸多边形的个数,以下依次描述各个多边形.第 ...

  9. LG4196 [CQOI2006]凸多边形

    题意 题目描述 逆时针给出n个凸多边形的顶点坐标,求它们交的面积.例如n=2时,两个凸多边形如下图: 则相交部分的面积为5.233. 输入输出格式 输入格式: 第一行有一个整数n,表示凸多边形的个数, ...

  10. P4196 [CQOI2006]凸多边形 半平面交

    \(\color{#0066ff}{题目描述}\) 逆时针给出n个凸多边形的顶点坐标,求它们交的面积.例如n=2时,两个凸多边形如下图: 则相交部分的面积为5.233. \(\color{#0066f ...

随机推荐

  1. Spark算子--mapPartitions和mapPartitionsWithIndex

    mapPartitions--Transformation类算子 代码示例 result   mapPartitionsWithIndex--Transformation类算子 代码示例 result ...

  2. mysql 远程连接数据库的二种方法

    一.连接远程数据库: 1.显示密码 如:MySQL 连接远程数据库(192.168.5.116),端口"3306",用户名为"root",密码"123 ...

  3. vue2.0集成百度UE编辑器,上传图片报错!!!

    我这边配置进去之后,界面加载,文本输入都没有问题,就是上传图片会有问题 这张图, 左边红色框框 就是目录结构咯, 右边红色框框 就是各种网上教程给出的第一个路径配置对吧, 下面的就是绿色 服务器接口配 ...

  4. 本地如何使用phpstudy环境搭建多站点

    http://jingyan.baidu.com/article/e52e36154227ef40c70c5147.html 平时在开发项目的时候, 多个项目同时开发的时候会遇到都得放到根目录才能正常 ...

  5. 利用光场进行深度图估计(Depth Estimation)算法之二——匹配算法

    光场相机由于能够捕获相机内部光线的强度和方向而得到整个光场,可以实现重聚焦(refocus)和视角变换等功能.进而可以进行深度估计获取深度图,前面说过利用重聚焦的图像进行深度估计,今天说一下利用不同视 ...

  6. @RequestMapping 相关 spring

    *      * @param request HttpServletRequest     * @param delList 削除Idエスト     * @return 削除結果     * @th ...

  7. HTTP常用方法

    GET : 获取资源 get方法用来请求访问已被URI识别的资源. 请求 GET /index.html HTTP/1.1 HOST:www.baidu.com 响应 返回index.html的页面资 ...

  8. (1-2)SpringCloud:服务的消费者rest+ribbon

    服务发现的任务是由Eureka客户端完成,而服务的消费任务由Ribbon完成.Ribbon是一个基于HTTP和TCP的客户端负载据衡器,它可以通过客户端中配置ribbonServerList服务端列表 ...

  9. mybatis-枚举类型的typeHandler&自定义枚举类型typeHandler

    MyBatis内部提供了两个转化枚举类型的typeHandler给我们使用. org.apache.ibatis.type.EnumTypeHandler 是使用枚举字符串名称作为参数传递的 org. ...

  10. 【javascript】onblur调用函数失效问题记录

    在处理页面上一个失焦效果时遇到了  xxxx is not a function 这样的错误,但函数确实有定义成功 ,经过测试发现是 id 值与函数名一致导致的问题,但是产生的原因没有找到详细的解释 ...