BNU校赛总决赛J 小白兔小灰兔 相交计算几何模板
J 小白兔小灰兔
空间限制:C/C++ 32768K,其他语言65536K
Special Judge, 64bit IO Format: %lld
题目描述
老山羊伯伯在地里收白菜,小白兔和小灰兔看见了就一起来帮忙。
他们干了半天,终于干完了。
羊伯伯:小灰兔,这车白菜送给你!
小灰兔:谢谢羊伯伯!
羊伯伯:小白兔,我也送你一车白菜!
小白兔:我不要白菜!给我一包白菜种子吧!
羊伯伯:好!给你~
小白兔:谢谢羊伯伯~
小灰兔把白菜拉到家里之后,就跟大家梦想中的生活一样,躺着啥都不干,吃吃吃,吃了玩~(好想一辈子都这样啊~小灰兔心想。)
小白兔把种子拿回去,打算开始勤劳地种白菜。然而他发现不是所有土地都能用来种白菜,只有被阳光照到的地方可以种白菜。
小白兔生活的星球可以看作二维平面中的一个简单多边形,太阳可以看作一个点。小白兔想知道,这个星球上一共有长度多少的土地可以用来种白菜。
输入描述:
第一行是一个正整数T(≤ 20),表示测试数据组数,
对于每组测试数据,
第一行是一个整数n(3≤ n≤ 10),表示简单多边形的顶点数,
接下来n行,每行是两个整数xi,yi(-10≤ xi,yi≤10),按照逆时针绕向给出简单多边形的n个顶点的坐标,
最后一行是两个整数x,y(-10≤ x,y≤10),表示太阳的坐标,保证太阳在多边形外且不在多边形任意一条边所在直线上。
输出描述:
对于每组测试数据,输出一行,包含一个实数,表示可以用来种白菜的土地的总长度,要求相对误差不超过,
也就是说,令输出结果为a,标准答案为b,若满足,则输出结果会被认为是正确答案。
输入例子:
2
4
0 0
2 0
2 2
0 2
4 1
4
0 0
2 0
2 2
0 2
4 4
输出例子:
2.000000000000
4.000000000000
-->
输入
2
4
0 0
2 0
2 2
0 2
4 1
4
0 0
2 0
2 2
0 2
4 4
输出
2.000000000000
4.000000000000

#include<bits/stdc++.h>
#define clr(x) memset(x,0,sizeof(x))
#define clr_1(x) memset(x,-1,sizeof(x))
#define mod 1000000007
#define INF 0x3f3f3f3f
#define LL long long
#define pb push_back
#define pbk pop_back
#define ls(i) (i<<1)
#define rs(i) (i<<1|1)
#define mp make_pair
using namespace std;
typedef double db;
const int N=5e5+;
const db eps=1e-;
int equzero(db t)//带精度大小判断
{
if(t>eps) return ;
if(t<-eps) return -;
return ;
}
struct Point
{
db x,y;
Point(){}
Point(db _x,db _y):x(_x),y(_y) {}
Point operator +(const Point &b) const //+
{
return Point(x+b.x,y+b.y);
}
Point operator -(const Point &b) const//-
{
return Point(x-b.x,y-b.y);
}
Point operator *(const db &b) const// 放大
{
return Point(x*b,y*b);
}
Point operator /(const db &b) const//缩小
{
return Point (x/b,y/b);
}
db dot(const Point &b) const//点积
{
return x*b.x+y*b.y;
}
db det(const Point &b) const//叉积
{
return x*b.y-b.x*y;
}
bool operator ==(const Point &b) const
{
return equzero(x-b.x)== && equzero(y-b.y)==;
}
bool operator !=(const Point &b) const
{
return !(Point(x,y)==b);
}
db mo() //模
{
return sqrt(x*x+y*y);
}
}pt[N],u,v,mid;
struct Line
{
Point a,b;
Line (){}
Line (Point _a,Point _b):a(_a),b(_b) {};
bool quickrej(const Line &t) const//快速排斥判断
{
return equzero(min(a.x,b.x)-max(t.a.x,t.b.x))<= && equzero(min(t.a.x,t.b.x)-max(a.x,b.x))<= && equzero(min(a.y,b.y)-max(t.a.y,t.b.y))<= && equzero(min(t.a.y,t.b.y)-max(a.y,b.y))<=;
}
db len() const
{
return (b-a).mo();
}
bool stra(const Line &t) const //跨立判断
{
return equzero((a-t.a).det(t.b-t.a)*(b-t.a).det(t.b-t.a))< && equzero((t.a-a).det(b-a)*(t.b-a).det(b-a))< ;
}
bool gfxj(const Line &t) const //规范相交
{
return quickrej(t) && stra(t);
}
bool on(const Point &t) const //点在线段上
{
Line p=Line(a,b);
if(p.a.x>p.b.x)
swap(p.a,p.b);
if(equzero(t.x-p.a.x)< || equzero(t.x-p.b.x)>)
return false;
if(equzero((p.a-t).det(p.b-t))==)
return true;
else
return false;
}
bool kddxj(const Line &t) const //可多点相交
{
return quickrej(t) && (stra(t) || on(t.a) || on(t.b) || t.on(a) || t.on(b));
}
bool para(const Line &t) const // 平行或在一条直线上
{
return equzero((a.y-b.y)*(t.a.x-t.b.x)-(t.a.y-t.b.y)*(a.x-b.x))== ;
}
db operator & (const Line &t) const //返回切点占直线比例,距a比例
{
if(equzero((a-b).det(t.a-t.b))==)return -INF;
db tk=((a-t.a).det(t.a-t.b))/((a-b).det(t.a-t.b));
return equzero(tk)>= && equzero(tk-)<= ? tk : -INF;
}
}le[N],lu,lv;
int n,m,T;
vector<db> cp;
db ans,ins,rate;
bool flag;
int main()
{
scanf("%d",&T);
while(T--)
{
scanf("%d",&n);
for(int i=;i<=n;i++)
scanf("%lf%lf",&pt[i].x,&pt[i].y);
for(int i=;i<n;i++)
le[i]=Line(pt[i],pt[(i+)%n]);
ans=;
for(int i=;i<n;i++)
{
cp.clear();
for(int j=;j<n;j++)
{
ins=le[i]&Line(pt[n],pt[j]);
if(ins>=-) cp.pb(ins);
}
sort(cp.begin(),cp.end());
rate=;
for(int j=;j<cp.size()-;j++)
{
lu=Line(pt[n],le[i].a+(le[i].b-le[i].a)*((cp[j]+cp[j+])/));
flag=;
for(int k=;k<n;k++)
{
if(lu.gfxj(le[k]))
{
flag=;
break;
}
}
if(!flag)
{
rate+=cp[j+]-cp[j];
}
}
ans+=rate*le[i].len();
}
printf("%.12f\n",ans);
}
return ;
}
BNU校赛总决赛J 小白兔小灰兔 相交计算几何模板的更多相关文章
- 2018BNU校赛总决赛
题解是qls的题解我就懒得写了23333 A塞特斯玛斯塔 时间限制:C/C++ 1秒,其他语言2秒 空间限制:C/C++ 32768K,其他语言65536K 64bit IO Format: %lld ...
- 牛客网暑期ACM多校训练营(第三场) J Distance to Work 计算几何求圆与多边形相交面积模板
链接:https://www.nowcoder.com/acm/contest/141/J来源:牛客网 Eddy has graduated from college. Currently, he i ...
- 2014哈商大ICPC/ACM校赛解题报告
被debug邀请去參加校赛,哎,被虐..我对不起工大.. 由于本人不搞ACM,算法处于HelloWorld水准.. 虽然题目除了鸟不拉屎星人之外都非常水,但我能做到这个程度,全然是超水平发挥了.. 数 ...
- 2016 华南师大ACM校赛 SCNUCPC 非官方题解
我要举报本次校赛出题人的消极出题!!! 官方题解请戳:http://3.scnuacm2015.sinaapp.com/?p=89(其实就是一堆代码没有题解) A. 树链剖分数据结构板题 题目大意:我 ...
- 2014上半年acm总结(1)(入门+校赛)
大一下学期才开始了acm,不得不说有一点迟,但是acm确实使我的生活充实了很多,,不至于像以前一样经常没事干= = 上学期的颓废使我的c语言学的渣的一笔..靠考前突击才基本掌握了语法 寒假突然醒悟, ...
- HZNU第十二届校赛赛后补题
愉快的校赛翻皮水! 题解 A 温暖的签到,注意用gets #include <map> #include <set> #include <ctime> #inclu ...
- 牛客网多校赛第9场 E-Music Game【概率期望】【逆元】
链接:https://www.nowcoder.com/acm/contest/147/E 来源:牛客网 时间限制:C/C++ 1秒,其他语言2秒 空间限制:C/C++ 262144K,其他语言524 ...
- ZJU 17th 校赛
第一次参加校赛,和小伙伴们拿了7个气球,还是挺开心的. 简单记个流水账吧. A:判断出INF的情况后 暴力模拟即可. INF的情况有x=1 || y=1 || (x==2 && y= ...
- 长春理工大学第十四届程序设计竞赛(重现赛)J.Printout
链接:https://ac.nowcoder.com/acm/contest/912/J 题意: 小r为了打校赛,他打算去打字社打印一份包含世界上所有算法的模板. 到了打字社,小r一看价格:总打印页数 ...
随机推荐
- flex实例(阮一峰)
Flex 布局教程:实例篇 作者: 阮一峰 日期: 2015年7月14日 上一篇文章介绍了Flex布局的语法,今天介绍常见布局的Flex写法. 你会看到,不管是什么布局,Flex往往都可以几行命令 ...
- NYOJ 1012 RMQ with Shifts (线段树)
题目链接 In the traditional RMQ (Range Minimum Query) problem, we have a static array A. Then for each q ...
- C# 动态调取 soap 接口
调用示例 string url = "http://localhost:8080/server/PatientService.asmx"; Hashtable ht = new H ...
- Windows降权
使用invoke-tokenmanipulation进行降权 枚举所有令牌 PS C:\Users\SMC> Get-ExecutionPolicy Restricted PS C:\Users ...
- 巅峰极客第二场CTF部分writeup
word-MISC 微信回答问题+word字体里. sqli-WEB 注册个admin空格即可,长字符截断. 晚上把后续的写出来.现在睡觉
- python并发编程之threading线程(一)
进程是系统进行资源分配最小单元,线程是进程的一个实体,是CPU调度和分派的基本单位,它是比进程更小的能独立运行的基本单位.进程在执行过程中拥有独立的内存单元,而多个线程共享内存等资源. 系列文章 py ...
- virsh 命令最新整理。 每个“;”之后是正解
1,migrate --domain --destURL --dname --live(热迁移) migrate lf 192.168.16.3 dname 2,managedsave domain ...
- 修改 firefox accesskey 的快捷键
Chrome中,如果设置了 accesskey 的话,可以通过 Alt + 快捷键 来之直接跳转的.但在Firefox 中,可能是为了防止于菜单的快捷键冲突,所以设置了 Shift + Alt + 快 ...
- 教你如何修改FireFox打开新标签页(NewTab Page)的行列数
FireFox的打开新建标签页(即NewTab Page)默认只能显示3x3个网站缩略图,这9个自定义的网站,非常方便快捷,什么hao123的弱爆了,本人从未用过此类导航网站,曾经用过的也只是abou ...
- mac 安装mongodb与常用操作
1.安装 brew update brew install mongodb 2.启动mongo mongod --config /usr/local/etc/mongod.conf 3.启动 mong ...