AtCoder Beginner Contest 356 G. Freestyle(凸包+二分)

news2024/11/9 0:48:15

题目

思路来源

quality代码

题解

对n个泳姿点(ai,bi)建凸包,实际上是一个上凸壳,

对于询问(ci,di)来说,抽象画一下这个图,箭头方向表示询问向量

按x轴排增序,并且使得后面的y不小于前面的y,因为总可以多耗费体力达到相同的米数

新增一个点(0,0),新增一个点(1e9+1,y),其纵坐标与最后一个点纵坐标相同

对于询问的向量,求它与上凸壳的交点,显然用这个方向的向量是最优的

1. 如果询问的向量在上凸壳第一个点逆时针方向,无解

2. 如果与上凸壳有交点,

显然这个封闭图形内的点都是可以用各种泳姿1s凑出来的,那么交点方向最优,

只需看(ci,di)中的di,是交点处的纵坐标的几倍即可,答案即为这个倍数

3. 如果向量与竖直的这条边有交点,那么最后相当于多耗费体力的情况下1s也可以跑这么多米

那么,两条直线还是有交点的,只是向量与水平直线的交点,

实际用的仍然是交点这个向量,只是实际的含义是y相同的情况下多耗费了体力

然后qls的二分写的就很妙,统一了这三种情况,以及交点恰好是凸包的一个点的情况

我写的就很不优雅,需要讨论三种情况,判断l=0,l=n-1

代码1(qls)

#include<bits/stdc++.h>
#include<iostream>
using namespace std;
typedef long long ll;
typedef long double db;
struct Point
{
    ll x,y;
    int id;
    Point() {}
    Point(ll _x,ll _y):x(_x),y(_y) {}
    Point operator - (const Point& t)const
    {
        return Point(x-t.x,y-t.y);
    }
    ll operator * (const Point& t)const
    {
        return x*t.y-y*t.x;
    }
};
vector<Point> Graham(vector<Point> p)
{
    p.insert(p.begin(),Point(0,0));
    sort(p.begin(),p.end(),[](const Point& lhs,const Point& rhs)
    {
        return lhs.x==rhs.x ? lhs.y<rhs.y : lhs.x<rhs.x;
    });
    for(size_t i=1; i<p.size(); i++)
        p[i].y=max(p[i].y,p[i-1].y);
    p.emplace_back(1000000001,p.back().y);
    vector<Point> res;
    for(auto& t : p)
    {
        while(res.size()>1 && (t-res.back())*(t-res[res.size()-2])<=0)res.pop_back();
        res.push_back(t);
    }
    res.erase(res.begin());
    return res;
}
int main()
{
    int n;
    scanf("%d",&n);
    vector<Point> p(n);
    for(int i=0; i<n; i++)
        scanf("%lld%lld",&p[i].x,&p[i].y);
    p=Graham(p);
    n=p.size();
    int q;
    scanf("%d",&q);
    while(q--)
    {
        Point que;
        scanf("%lld%lld",&que.x,&que.y);
        if(que*p[0]<0)
        {
            printf("-1\n");
            continue;
        }
        int l=0,r=n-2;//l<=r+1
        while(l<r)//>0 =0 <0 或>0 <0渐变
        {
            int m=(l+r+1)/2;
            if(que*p[m]>0)l=m;
            else r=m-1;
        }
        db k1=1.0L*(p[l+1].y-p[l].y)/(p[l+1].x-p[l].x);
        db b1=p[l].y-k1*p[l].x;
        db k2=1.0L*que.y/que.x;
        db x=b1/(k2-k1),y=k2*x;
        printf("%.18Lf\n",que.y/y);
    }
    return 0;
}

代码2(我的二分)

#include<bits/stdc++.h>
#include<iostream>
using namespace std;
typedef long long ll;
typedef long double db;
struct Point
{
    ll x,y;
    int id;
    Point() {}
    Point(ll _x,ll _y):x(_x),y(_y) {}
    Point operator - (const Point& t)const
    {
        return Point(x-t.x,y-t.y);
    }
    ll operator * (const Point& t)const
    {
        return x*t.y-y*t.x;
    }
};
vector<Point> Graham(vector<Point> p)
{
    p.insert(p.begin(),Point(0,0));
    sort(p.begin(),p.end(),[](const Point& lhs,const Point& rhs)
    {
        return lhs.x==rhs.x ? lhs.y<rhs.y : lhs.x<rhs.x;
    });
    for(size_t i=1; i<p.size(); i++)
        p[i].y=max(p[i].y,p[i-1].y);
    p.emplace_back(1000000001,p.back().y);
    vector<Point> res;
    for(auto& t : p)
    {
        while(res.size()>1 && (t-res.back())*(t-res[res.size()-2])<=0)res.pop_back();
        res.push_back(t);
    }
    res.erase(res.begin());
    return res;
}
int main()
{
    int n;
    scanf("%d",&n);
    vector<Point> p(n);
    for(int i=0; i<n; i++)
        scanf("%lld%lld",&p[i].x,&p[i].y);
    p=Graham(p);
    n=p.size();
    int q;
    scanf("%d",&q);
    while(q--)
    {
        Point que;
        scanf("%lld%lld",&que.x,&que.y);
        if(que*p[0]<0)
        {
            printf("-1\n");
            continue;
        }
        int l=0,r=n-1;
        while(l<=r){
            int m=(l+r)/2;
            if(que*p[m]>0)l=m+1;
            else r=m-1;
        }
        l--;
        if(l<0)l=0;
        if(l==n-1)l--;
        //if(l==n-1)l--;
        db k1=1.0L*(p[l+1].y-p[l].y)/(p[l+1].x-p[l].x);
        //k1=0 b1=5 k2=8/6 x=5/(8/6)=30/8 y=30/6=5 
        db b1=p[l].y-k1*p[l].x;
        db k2=1.0L*que.y/que.x;
        db x=b1/(k2-k1),y=k2*x;
        printf("%.18Lf\n",que.y/y);
    }
    return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/1808017.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

Docker高级篇之Docker-compose容器编排

文章目录 1. Docker-compse介绍2. Docker-compse下载3. Docker-compse核心概念4. Docker-compse使用案例 1. Docker-compse介绍 Docker-compose时Docker官方的一个开源的项目&#xff0c;负责对Docker容器集群的快速编排。Docker-compose可以管理多个Docker容器组成一个应用&a…

【单片机毕业设计9-基于stm32c8t6的酒窖监测系统】

【单片机毕业设计9-基于stm32c8t6的酒窖监测系统】 前言一、功能介绍二、硬件部分三、软件部分总结 前言 &#x1f525;这里是小殷学长&#xff0c;单片机毕业设计篇9基于stm32的酒窖监测系统 &#x1f9ff;创作不易&#xff0c;拒绝白嫖可私 一、功能介绍 -------------------…

Docker高级篇之轻量化可视化工具Portainer

文章目录 1. 简介2. Portainer安装 1. 简介 Portianer是一款轻量级的应用&#xff0c;它提供了图形化界面&#xff0c;用于方便管理Docker环境&#xff0c;包括单机环境和集成环境。 2. Portainer安装 官网&#xff1a;https://www.portainer.io 这里我们使用docker命令安装&…

8.让画面动起来

一、Unity Shader中的内置变量&#xff08;时间篇&#xff09; 动画效果往往都是把时间添加到一些变量的计算中&#xff0c;以便在时间变化的同时也可以随之变化。Unity shader提供了一系列关于时间的内置变量来允许我们方便地在Shader中访问运行时间&#xff0c;实现各种动画…

STM32 | 独立看门狗 | RTC(实时时钟)

01、独立看门狗概述 在由单片机构成的微型计算机系统中,由于单片机的工作常常会受到来自外界电磁场的干扰,造成程序的跑飞,而陷入死循环,程序的正常运行被打断,由单片机控制的系统无法继续工作,会造成整个系统的陷入停滞状态,发生不可预料的后果,所以出于对单片机运行状…

Elasticsearch:Open Crawler 发布技术预览版

作者&#xff1a;来自 Elastic Navarone Feekery 多年来&#xff0c;Elastic 已经经历了几次 Crawler 迭代。最初是 Swiftype 的 Site Search&#xff0c;后来发展成为 App Search Crawler&#xff0c;最近又发展成为 Elastic Crawler。这些 Crawler 功能丰富&#xff0c;允许以…

基于Java+SpringBoot制作一个景区导览小程序

基于Java+SpringBoot制作一个景区导览小程序。其中系统前端功能包括注册登录、景区采风、旅游导览、地图导航、发布采风、门票预订、修改个人信息;系统后台功能包括用户管理、景区管理、采风管理等模块。 摘要一、小程序1. 创建小程序2. 首页3. 景区采风页4. 旅游导览页5. 发布…

人工智能_机器学习097_PCA数据降维算法_数据去中心化_特征值_特征向量_EVD特征值分解原来和代码实现---人工智能工作笔记0222

降维算法的原理,一会咱们再看,现在先看一下,算法 可以看到PCA算法的,原理和过程,我们先看一下代码 为了说明PCA原理,这里,我们,先来计算一下X的方差,可以看到 先把数据进行去中心化,也就是用数据,减去数据的平均值. B = X-X.mean(axis=0) 这段代码是用于计算矩阵X的每一列减去该…

【Web世界探险家】3. CSS美学(二)文本样式

&#x1f4da;博客主页&#xff1a;爱敲代码的小杨. ✨专栏&#xff1a;《Java SE语法》 | 《数据结构与算法》 | 《C生万物》 |《MySQL探索之旅》 |《Web世界探险家》 ❤️感谢大家点赞&#x1f44d;&#x1f3fb;收藏⭐评论✍&#x1f3fb;&#xff0c;您的三连就是我持续更…

⌈ 传知代码 ⌋ 基于曲率的图重新布线

&#x1f49b;前情提要&#x1f49b; 本文是传知代码平台中的相关前沿知识与技术的分享~ 接下来我们即将进入一个全新的空间&#xff0c;对技术有一个全新的视角~ 本文所涉及所有资源均在传知代码平台可获取 以下的内容一定会让你对AI 赋能时代有一个颠覆性的认识哦&#x…

SQL159 每个创作者每月的涨粉率及截止当前的总粉丝量

描述 用户-视频互动表tb_user_video_log iduidvideo_idstart_timeend_timeif_followif_likeif_retweetcomment_id110120012021-09-01 10:00:002021-09-01 10:00:20011NULL210520022021-09-10 11:00:002021-09-10 11:00:30101NULL310120012021-10-01 10:00:002021-10-01 10:00…

如何获取MySQL中表的大小?(官方校正版)

与大多数关系数据库一样&#xff0c;MySQL 提供了有关数据库本身的有用元数据。虽然大多数其他数据库将此信息称为 catalog&#xff0c; 但MySQL 官方文档INFORMATION_SCHEMA 将元数据 称为 tables。 目录 1 列出单个数据库中的单表大小 2 列出所有数据库中的所有表大小 以下…

从年金理论到杠杆效应,再到财务报表与投资评估指标

一、解释普通年金终值和普通年金现值的概念。 普通年金终值&#xff1a;以利率为1%&#xff0c;每期收款100元&#xff0c;5期为例&#xff0c;普通年金终值的折算过程如图&#xff1a; 普通年金现值&#xff1a;以利率为1%&#xff0c;每期收款100元&#xff0c;5期为例&am…

【C++修行之道】类和对象(六)再谈构造函数(初始化列表)| explicit关键字 | static成员 | 友元|匿名对象|拷贝时一些编译器优化

目录 一、再谈构造函数 1.1 构造函数体赋值 1.2 初始化列表 1. 所有的成员,既可以在初始化列表初始化,也可以在函数体内初始化 2. 每个成员变量在初始化列表中只能出现一次(初始化只能初始化一次) 3. 类中包含以下成员&#xff0c;必须放在初始化列表位置进行初始化&…

【设计模式】行为型设计模式之 备忘录模式(快照模式)

介绍 备忘录应用场景明确并且有限&#xff0c;一般用来数据的防丢失、撤销和恢复。对大对象的备份和恢复&#xff0c;备忘录模式能有效的节省时间和空间开销。 定义 备忘录模式&#xff1a;也称为快照模式&#xff0c;在不违背封装原则的前提下&#xff0c;捕获一个对象的内…

20240610 基于QGIS生成地区示意图的地图shp文件

目录 本文目标前置条件具体步骤1. 创建Project2. 插入世界地图3. 对地区示意图进行地理匹配4. 创建shp文件&#xff0c;勾画轨迹 注意事项 本文目标 基于QGIS生成地区示意图的地图shp文件&#xff0c;此shp文件可以用来学习&#xff0c;但是未经审批不可用于发表。 前置条件 …

水滴型锤片粉碎机:多功能粉碎利器

在现代工业生产中&#xff0c;粉碎机作为一种重要的机械设备&#xff0c;广泛应用于饲料、化工、木材等多个领域。其中&#xff0c;水滴型锤片粉碎机凭借其设计和粉碎能力&#xff0c;成为市场上的热门产品。 水滴型锤片粉碎机其设计灵感来源于水滴的形态。这种设计使得机器在…

【全开源】房屋出租出售预约系统(FastAdmin+ThinkPHP+Uniapp)

房屋出租出售预约系统&#xff1a;一站式解决房产交易难题 一款基于FastAdminThinkPHPUniapp开发的房屋出租出售预约系统&#xff0c;支持小程序、H5、APP&#xff0c;包含房客、房东(高级授权)、经纪人(高级授权)三种身份。核心功能有&#xff1a;新盘销售、房屋租赁、地图找…

Android 事件分发机制详解(上)

前言 Android事件分发机制是Android开发者必须了解的基础。 目录 一. 基础认知 1.1 事件分发的由来 安卓的View是树形结构的&#xff0c;View可能会重叠在一起&#xff0c;当我们点击的地方有多个View都可以响应的时候&#xff0c;这个点击事件应该给谁呢&#xff1f;为了解…

cefsharp124.x升级125.x(cef125.0.21/Chromium 125.0.6422.142)

一、版本说明 1.1 依赖关系变化 依赖移除:cef.redist.x64,cef.redist.x86增加新支持chromiumembeddedframework.runtime 旧版本需要移除依赖cef.redist.x64和cef.redist.x86否则会初始化异常。 自版本121.*以后common依赖关系变化 chromiumembeddedframework.runtime.win-x6…