半平面交+数组模拟双端队列

人生第一次代码过两百行啊...加油加油

#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<cstdio>
#include<iomanip> #define ld long double
using namespace std; const int maxn = 1e3;
const ld pi = acosl(-);
const ld eps = 1e-;//控制精度 int n;
int head, tail;
ld r, cost1, cost2; typedef struct Grid
{//棋盘
ld x,y;
Grid(double a = , double b = ){x = a, y = b;}//构造函数,结构体对象创建时执行
} point , Vector ;//point:类型为Grid的数据类型;Vector与point结构相同,数据类型名称不同 typedef struct Net
{
point a,b;
Net() {}//重载构造函数,目的是以形如line x;的方式定义直线,避免定义时确定参数的复杂性
Net(point i, point j){a = i, b = j;}
} line ; point node[maxn], p[maxn];
line l[maxn], L[maxn], q[maxn]; Vector operator - (point a, point b) // V需大写,避免与vector混淆
{
return Vector(b.x - a.x, b.y - a.y);
}//矢量定义 double operator * (Vector a, Vector b)
{
return (a.x * b.y - a.y * b.x);
}//叉乘 ld getangle(Vector t)
{
return atan2(t.y, t.x);//返回一个(-pi,pi]的角度
}//获得某个矢量的极角 ld getangle(line t)
{
return atan2(t.b.y - t.a.y, t.b.x - t.a.x);
}//获得某条直线的极角 bool cmp(line a, line b)
{
Vector v1 = (a.b - a.a), v2 = (b.b - b.a);
ld A = getangle(v1), B = getangle(v2);
if(fabs(A - B) < eps){
return (v1 * (b.b - a.a) >= );//靠左边的丢在后面,方便双端队列去重
}
return A < B;
}//决定根据极角排序的直线的优先级 bool check()
{
ld ans = ;
for(int i = ; i < n- ; i++){
ans += (node[i] - node[]) * (node[i+] - node[]);
}
if(ans < ) return false;
return true;
}//检查point输入顺序进行调整 true为逆时针 point getnode(line a, line b)//求两直线(线段)交点,一般式三参数暴力求解
{
ld a1 = a.b.y - a.a.y, b1 = a.a.x - a.b.x, c1 = a.b.x * a.a.y - a.a.x * a.b.y;
ld a2 = b.b.y - b.a.y, b2 = b.a.x - b.b.x, c2 = b.b.x * b.a.y - b.a.x * b.b.y;
ld D = a1 * b2 - a2 * b1;
//Dx/D,Dy/D,注意c1、c2移至右侧变号
return point((c2 * b1 - b2 * c1) / D, (c1 * a2 - a1 * c2) / D);
} bool on_right(line a, line b, line c)//观察a、b交点与c的关系
{
point n = getnode(a, b);
if((c.b - c.a) * (n - c.a) < ){//交点在右侧/下方
return true;
}
return false;
} line narrow_line(line x)//实现直线平移,勾画出圆心轨迹
{
ld dx = x.b.x - x.a.x, dy = x.b.y - x.a.y;
Vector turnv;//旋转矢量
turnv.x = -dy, turnv.y = dx; ld m = sqrtl(dx * dx + dy * dy);
turnv.x = turnv.x / m * r;
turnv.y = turnv.y / m * r;
//求对应坐标
x.b.x += turnv.x;
x.b.y += turnv.y;
x.a.x += turnv.x;
x.a.y += turnv.y;
return x;
} void narrow_polygon()//模拟圆心轨迹
{
for(int i = ; i <= n- ; i++){
L[i] = narrow_line(l[i]);
}
} bool HPI()
{
sort(L,L+n,cmp);//按极角排序
head = , tail = ;
int cnt = ;
for(int i = ; i < n- ; i++){
Vector v1 = L[i].b - L[i].a, v2 = L[i+].b - L[i+].a;//?
if(fabs(getangle(v1) - getangle(v2)) < eps){
continue;
}
L[cnt++] = L[i];
}
L[cnt++] = L[n-]; /////正片开始
for(int i = ; i < cnt ; i++){//注意下标
while(tail - head > && on_right(q[tail-], q[tail-], L[i])) tail--;//符合onright规则,删去上一条直线
while(tail - head > && on_right(q[head], q[head+], L[i])) head++;//头部同理
q[tail++] = L[i];
}
//换个角度看世界,最后检验q[head]以及q[tail-1],q[tail]实际为空
while(tail - head > && on_right(q[tail-], q[tail-], q[head])) tail--;
while(tail - head > && on_right(q[head], q[head+], q[tail-])) head++; if((tail - head) >= ) return true;
return false;
} ld get_S()//求初始面积
{
ld ans = ;
for(int i = ; i < n ; i++){
ans += fabs((node[i] - node[]) * (node[(i+)%n] - node[]));
}
return ans / ;
} ld dis(point a, point b)
{
ld dx = b.x - a.x, dy = b.y - a.y;
return sqrtl(dx * dx + dy * dy);
}//两点间距离 ld get_inS()//内部半平面交 面积
{
ld res = ;
int tot = ;
for(int i = head ; i < tail- ; i++){
p[tot++] = getnode(q[i], q[i+]);
}
p[tot++] = getnode(q[tail-], q[head]); for(int i = ; i < tot- ; i++){
res += fabs((p[i] - p[]) * (p[i+] - p[]));
}
res /= ; for(int i = ; i < tot- ; i++){
res += dis(p[i], p[i+]) * r;
}// c * r
res += dis(p[], p[tot - ]) * r; return res + r * r * pi;
} int main(){
//freopen("mow.in","r",stdin);
int T;cin >> T;
while( T-- )
{
cin >> n >> r >> cost1 >> cost2;
for(int i = n- ; i >= ; i--){
cin >> node[i].x >> node[i].y;
}//本题中为顺时针方向,因而逆序存点,以确保直线左侧为内侧 if(check()){//逆
for(int i = ; i < n- ; i++){//
l[i] = line(node[i], node[i+]);
}
l[n-] = line(node[n-], node[]);
}else{//顺
for(int i = ; i < n- ; i++){
l[i] = line(node[i+],node[i]);
}
l[n-] = line(node[], node[n-]);
}//0 ~ n-1 narrow_polygon(); ld s = get_S();//草坪面积
ld ans = s * cost1;//手动割草费用 if(HPI()){
ld mowable = get_inS();//机器割草面积
ans = min(ans, mowable * cost2 + (s - mowable) * cost1);
}
cout << fixed << setprecision() << ans << "\n"; } return ;
}
//实现直线平移,勾画出圆心轨迹

hdu6755 Mow的更多相关文章

  1. HDU 6762 Mow (2020 Multi-University Training Contest 1 1012) 半平面交

    Mow 题目链接 分析 将多边形的边向内部缩 r 个单位长度,然后这些边所围成的内部区域,就是圆心的合法范围,该范围也是一个多边形,假设面积是\(a\),周长是\(b\),那么可以知道圆可以覆盖的面积 ...

  2. java中获取路径的几种基本的方法

    package com.ygh.blog.realpath; import java.io.File; import java.io.IOException; import java.io.Input ...

  3. TPC-H生成.tbl文件导入postgresql数据库的坑

    数据库project好好的不用主流的MySQL和Microsoft server而要求用听都没听过的postgresql (当然,可能你三个都没听过) 这里的坑主要是把生成的那八张.tbl的表导入pg ...

  4. js实现多张图片每隔一秒换一张图片

    aaarticlea/png;base64,iVBORw0KGgoAAAANSUhEUgAAAlkAAAHWCAIAAADLlAuAAAAgAElEQVR4nOzd5XNc157w+/l7bt2n6t

  5. mysql触发器,答题记录表同步教学跟踪(用户列表)

    aaarticlea/png;base64,iVBORw0KGgoAAAANSUhEUgAABVQAAAOOCAIAAABgEw4AAAAgAElEQVR4nOy92VcT27r/zX+xLtflvt

  6. java调用sqlldr导入csv文件数据到临时表

    package cn.com.file;import java.io.BufferedReader;import java.io.BufferedWriter;import java.io.File; ...

  7. CSS布局(二)

    本节内容:position.float.clear.浮动布局例子.百分比宽度 position CSS中的position属性设置元素的位置.属性值:static.relative.fixed.abs ...

  8. 这些HTML、CSS知识点,面试和平时开发都需要 No5-No7

    系列知识点汇总 这些HTML.CSS知识点,面试和平时开发都需要 No1-No4(知识点:HTML.CSS.盒子模型.内容布局) 这些HTML.CSS知识点,面试和平时开发都需要 No5-No7(知识 ...

  9. C#基础整理参数

    基本概念 把数据传入方法中,可以使方法有多个返回值. 参数的传递 值参数,通过将实参的值复制到形参的方式传递数据.值参数的实参可以是变量或者是表达式

随机推荐

  1. 关于SPSS Modeler18 提示:用于定义的观测值的字段的值无效

    今天在做实验的时候,按照实验步骤严格设置了参数,当运行节点的时候,一直提示:用于定义的观测值的字段的值无效,如下图 我把我的流文件发给同学,同学的机器上是可以运行的,但是我的不行,不知道什么原因,有知 ...

  2. pikachu靶场-CSRF

    xss和csrf区别: CSRF是借用户的权限完成攻击,攻击者并没有拿到用户的权限,而XSS是直接盗取到了用户的权限,然后实施破坏. PS: 由于之前将php5升级到php7,导致靶场环境出现以下问题 ...

  3. Apache Dubbo Provider默认反序列漏洞复现(CVE-2020-1948)

    Apache Dubbo Provider默认反序列漏洞(CVE-2020-1948) 0x01 搭建漏洞环境 漏洞介绍 2020年06月23日, 360CERT监测发现Apache Dubbo 官方 ...

  4. Spring 获取单例流程(三)

    读完这篇文章你将会收获到 Spring 何时将 bean 加入到第三级缓存和第一级缓存中 Spring 何时回调各种 Aware 接口.BeanPostProcessor .InitializingB ...

  5. 循环&&数组&&方法&&面向对象

    day03 数值的默认值 类型 初始化的值 byte,short,int,long 0 float,double  0.0 char 空格 boolean false 引用类型 null JVM的内存 ...

  6. Oracle Online Patching报错"This is not a RAC setup. OPatch cannot determine the local node name"

    Oracle Online Patching报错"This is not a RAC setup. OPatch cannot determine the local node name&q ...

  7. css实现1px 像素线条_解决移动端1px线条的显示方式

    使用CSS 绘制出 1px 的边框,在移动端上渲染的效果会出现不同,部分手机发现1px 线条变胖了,这篇文章整理2种方式实现1px 像素线条. 1.利用box-shadow + transform & ...

  8. js 字符串转方法,this域绑定

    闲着没事,开发一个列表页面配置的功能,其中涉及到了按钮点击事件,在页面进行编辑,保存到数据库中.写好的js脚本,function是字符串格式,所以要让生成的脚本生效,还要做一些操作. 1.首先保存在数 ...

  9. Spring Boot 2 实战:常用读取配置的方式

    1. 前言 在Spring Boot项目中我们经常需要读取application.yml配置文件的自定义配置,今天就来罗列一下从yaml读取配置文件的一些常用手段和方法. 2. 使用@Value注解 ...

  10. Tomcat双击startup.bat闪退的原因及解决方式

    很久不碰Tomcat了,最近因为种种原因需要重新投入到Java Web的怀抱,所以又重新接触了Tomcat 我下载了tomcat的压缩包将其解压缩到某个位置,我这里是D盘下的tomcat文件夹中,但是 ...