[SHOI2012]信用卡凸包(计算几何)
/*
考验观察法??
可以发现最终答案等于所有作为圆心的点求出凸包的周长加上一个圆的周长
向量旋转
(x1, y1) 相较于 (x2, y2) 旋转角c 答案是
(dtx * cosc - dty * sinc + x2, dty * cosc + dtx * sinc + y2)
貌似我的凸包也不鲁棒耶
*/
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<iostream>
#include<cmath>
#include<queue>
#define ll long long
#define M 400010
#define mmp make_pair
using namespace std;
int read()
{
int nm = 0, f = 1;
char c = getchar();
for(; !isdigit(c); c = getchar()) if(c == '-') f = -1;
for(; isdigit(c); c = getchar()) nm = nm * 10 + c - '0';
return nm * f;
}
const double pi = acos(-1), eps = 1e-12;
struct Vec
{
double x, y;
}pt[M];
double operator * (const Vec &a, const Vec &b)
{
return a.x * b.y - a.y * b.x;
}
Vec operator + (Vec a, Vec b)
{
return Vec{a.x + b.x, a.y + b.y};
}
Vec operator - (Vec a, Vec b)
{
return Vec{a.x - b.x, a.y - b.y};
}
double cross(Vec a, Vec b, Vec c)
{
return (a - c) * (b - c);
}
bool cmp1(Vec a, Vec b)
{
return a.y == b.y ? a.x < b.x : a.y < b.y;
}
bool cmp(Vec a, Vec b)
{
return cross(a, b, pt[1]) >= 0;
}
Vec rotate(Vec a, double p)
{
return Vec{a.x * cos(p) - a.y * sin(p), a.y * cos(p) + a.x * sin(p)};
}
double diss(Vec a, Vec b)
{
double x = a.x - b.x, y = a.y - b.y;
return sqrt(x * x + y * y);
}
double r, h, l, a, b, ap;
int n, tp, sta[M];
int main()
{
n = read();
scanf("%lf%lf%lf", &l, &h, &r);
h = (h - r * 2) / 2.0, l = (l - r * 2) / 2.0;
for(int i = 1; i <= n; i++)
{
scanf("%lf%lf%lf", &a, &b, &ap);
pt[i] = rotate(Vec{h, l}, ap) + Vec{a, b};
pt[i + 2 * n] = rotate(Vec{-h, l}, ap) + Vec{a, b};
pt[i + 3 * n] = rotate(Vec{h, -l}, ap) + Vec{a, b};
pt[i + n] = rotate(Vec{-h, -l}, ap) + Vec{a, b};
}
n *= 4;
sort(pt + 1, pt + n + 1, cmp1);
sta[++tp] = 1;
sort(pt + 2, pt + n + 1, cmp);
sta[++tp] = 2;
sta[++tp] = 3;
for(int i = 4; i <= n; i++)
{
while(tp > 1 && cross(pt[i], pt[sta[tp]], pt[sta[tp - 1]]) >= 0) tp--;
sta[++tp] = i;
}
double ans = 0;
for(int i = 1; i < tp; i++) ans += diss(pt[sta[i]], pt[sta[i + 1]]);
ans += diss(pt[sta[1]], pt[sta[tp]]);
ans += r * pi * 2;
printf("%.2lf\n", ans);
return 0;
}
[SHOI2012]信用卡凸包(计算几何)的更多相关文章
- 【BZOJ2829】[SHOI2012]信用卡凸包(凸包)
[BZOJ2829][SHOI2012]信用卡凸包(凸包) 题面 BZOJ 洛谷 题解 既然圆角的半径都是一样的,而凸包的内角和恰好为\(360°\),所以只需要把圆角的圆心弄下来跑一个凸包,再额外加 ...
- [SHOI2012]信用卡凸包(凸包+直觉)
这个题还是比较有趣. 小心发现,大胆猜想,不用证明! 我们发现所谓的信用卡凸包上弧的长度总和就是圆的周长! 然后再加上每个长宽都减去圆的直径之后的长方形的凸包周长即可! #include<ios ...
- luogu P3829 [SHOI2012]信用卡凸包 凸包 点的旋转
LINK:信用卡凸包 当 R==0的时候显然是一个点的旋转 之后再求凸包即可. 这里先说点如何旋转 如果是根据原点旋转的话 经过一个繁杂的推导可以得到一个矩阵. [cosw,-sinw] [sinw, ...
- P3829 [SHOI2012]信用卡凸包
思路 注意到结果就是每个信用卡边上的四个圆心的凸包周长+一个圆的周长 然后就好做了 注意平行时把距离小的排在前面,栈中至少要有1个元素(top>1),凸包中如果存在叉积为0的点也要pop,否则可 ...
- [洛谷P3829][SHOI2012]信用卡凸包
题目大意:有$n$张一模一样的信用卡,每个角进行了圆滑处理,问这些卡组成的“凸包”的周长 题解:发现是圆滑处理的圆心围成的凸包加上一个圆周即可 卡点:输入长宽弄反,然后以为是卡精 C++ Code: ...
- Luogu-3829 [SHOI2012]信用卡凸包
这道题的转化很巧妙,可以把信用卡四个角的圆心看做平面上的点来做凸包,\(ans\)就是凸包周长加上一个圆的周长 // luogu-judger-enable-o2 #include<cmath& ...
- BZOJ 2829 信用卡凸包 ——计算几何
凸包裸题 #include <map> #include <cmath> #include <queue> #include <cstdio> #inc ...
- 【BZOJ 2829】 2829: 信用卡凸包 (凸包)
2829: 信用卡凸包 Description Input Output Sample Input 2 6.0 2.0 0.0 0.0 0.0 0.0 2.0 -2.0 1.5707963268 Sa ...
- bzoj 2829 信用卡凸包(凸包)
2829: 信用卡凸包 Time Limit: 10 Sec Memory Limit: 128 MB Submit: 1342 Solved: 577 [Submit][Status][Disc ...
随机推荐
- JS判断字符串是否为空或是否全为空格
var test = " "; //为空或全部为空格 if (test.match(/^[ ]*$/)) { console.log("all space or empt ...
- Spring Cloud(Dalston.SR5)--Eureka 常用配置
配置参数 默认值 说明 服务注册中心配置 Bean类:org.springframework.cloud.netflix.eureka.server.EurekaServerConfigBean eu ...
- BT.656 NTSC制式彩条生成模块(verilog)
BT.656 NTSC制式彩条生成模块(verilog) 1.知识储备 隔行扫描是将一副图像分成两场扫描,第一场扫描第1,2,5,7...等奇数行,第二场扫描2,4,6,8...等偶数行,并把扫奇数行 ...
- jQuery自定义扩展写法
jQuery自定义扩展有两种写法: $.extend({}) ; 调用方法: $.方法 $.fn.extend({}); 调用方法: ...
- linux svn客户端通过 https访问windows VisualSVN Server Manager
1)需求: 已经在阿里云windwos系统 下面安装了VisualSVN Server Manager 做为svn服务器: 现在要在腾讯云源码安装新版本客户端 2)开始源码编译安装TortoiseSV ...
- 【linux】常用命令-端口
端口操作 手动更改配置文件开放端口 vim /etc/sysconfig/iptables -A INPUT -p tcp -m state --state NEW -m tcp --dport 81 ...
- window下sh文件在linux转码
本文转载自:https://blog.csdn.net/omiconlee/article/details/53407438 1)在Windows下转换: 利用一些编辑器如UltraEdit或Edit ...
- 【问题解决】docker中创建volume时,无访问权限
挂载宿主机已存在目录后,在容器内对其进行操作,报"Permission denied". 可通过两种方式解决: 1> 关闭selinux. 临时关闭:# setenforce ...
- Linux环境下配置maven环境
1.下载安装包并解压 安装包直接去apache官网下载:https://maven.apache.org/download.cgi 将安装包放在自己指定的目录:/home/software/apach ...
- file /usr/lib64/mysql/plugin/dialog.so from install of Percona-Server-server-56-5.6.24-rel72.2.el6.x86_64 conflicts with file from package mariadb-libs-1:5.5.60-1.el7_5.x86_64
!!!点下面!!! https://www.cnblogs.com/chuijingjing/p/10005922.html