首页 > 其他分享 >洛谷 P4423 [BJWC2011]最小三角形 题解

洛谷 P4423 [BJWC2011]最小三角形 题解

时间:2022-11-10 19:35:45浏览次数:40  
标签:y2 洛谷 题解 BJWC2011 x2 double y1 theta x1

求平面最近点对的时候有这样一种思路:

将所有点全部绕原点旋转同一个角度,然后按横坐标排序。

根据数学直觉,在随机旋转后,答案中的两个点在数组中肯定不会离得太远。

把这种思路运用到这题当中,我们旋转随机角度将点打乱之后枚举三角形的三个顶点即可。

坐标旋转公式:

\[x2=x1\cos\theta-y1\sin\theta\\ y2=x1\sin\theta+y1\cos\theta \]

注意角度制和弧度制的转化即可。

/*
 * Title: P4423 [BJWC2011]最小三角形
 * Source: 洛谷
 * URL: https://www.luogu.com.cn/problem/P4423
 * Author: Steven_lzx
 * Command: -std=c++23 -O2 -Wall -fno-ms-extensions
 * Date: 2022.11.10
 */
#pragma GCC optimize(2)
#include <bits/stdc++.h>
using namespace std;
const double PI=acos(-1);
struct Point{double x,y;}p[200005];
int n;
double ans=1e15;
double dis(double x1,double y1,double x2,double y2){return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));}
void rotate(double theta)
{
    double nx,ny,tx,ty;
    theta=theta/180*PI;
    for(int i=1;i<=n;i++)
    {
        nx=p[i].x;
        ny=p[i].y;
        tx=nx*cos(theta)-ny*sin(theta);
        ty=nx*sin(theta)+ny*cos(theta);
        p[i]=(Point){tx,ty};
    }
    sort(p+1,p+n+1,[](Point a,Point b){return a.x<b.x;});
    for(int i=1;i<n-1;i++)//暴力枚举三个点
        for(int j=i+1;j<min(i+20,n);j++)
            for(int k=j+1;k<=min(i+20,n);k++)
                ans=min(ans,dis(p[i].x,p[i].y,p[j].x,p[j].y)+dis(p[j].x,p[j].y,p[k].x,p[k].y)+dis(p[i].x,p[i].y,p[k].x,p[k].y));
    return;
}
int main()
{
    mt19937 rnd(20080422);
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
        scanf("%lf%lf",&p[i].x,&p[i].y);
    rotate(rnd()%360);
    rotate(rnd()%360);
    printf("%.6lf\n",ans);
    return 0;
}

标签:y2,洛谷,题解,BJWC2011,x2,double,y1,theta,x1
From: https://www.cnblogs.com/2020gyk080/p/16878119.html

相关文章

  • P4407 [JSOI2009] 电子字典 题解
    题目:P4407这题差不多就是P1688的改版。参考一下我在P1688的做法,我们继续使用Hash,然后只要考虑如何去重就好了。于是就有了这个暴力的想法:#代表修改,@代表添加,$代......
  • [AGC040F] Two Pieces 题解
    linkSolution这个题真的挺难的。/kk看了一个下午的题解才搞懂。/fn我们发现我们如果设状态\((x,d)\)表示前面的一个点在\(x\),另一个在\(x-d\),那么三种操作相当于:......
  • 题解 P3974【[TJOI2015]组合数学】
    postedon2022-10-2814:11:41|under题解|sourceproblem给出一个网格图,其中某些格子有财宝,每次从左上角出发,只能向下或右走。问至少走多少次才能将财宝捡完。此对......
  • 【题解】【切开字符串】
    P8631[蓝桥杯2015国AC]切开字符串Sol首先问题可以转化为对每个前缀求出本质不同奇回文子串数,和对每个后缀求出本质不同子串数和本质不同奇回文子串数。本质不同子......
  • simpread-(127 条消息) fastAPI 中的跨域问题解决_小童同学_的博客 - CSDN 博客_fasta
    本文由简悦SimpRead转码,原文地址blog.csdn.netCrossOrigin前言前端采用Vue,后端采用fastAPI的CV项目在开发时遇到跨域问题,记录学习过程与解决方案。概念C......
  • 体验 Python 剪辑视频以及相关问题解决, 一劳永逸!
    前言对于使用Python对视频进行剪辑我们最常用的就是Moviepy,我之前也写过一篇​​《必杀技--使用FFmpeg命令快速精准剪切视频》​​,这篇文章单纯使用的是FFmpeg,他是通过......
  • 题解 P8827 [传智杯 #3 初赛] 森林
    本题解提供两种做法。做法一为了叙述方便,先引入\(n\)级母树的概念。定义\(1\)级母树即为该子树被删去前,其所在的原来的完整的树。如下图,以\(5\)为根的一级母树......
  • node -v显示信息 npm -v 不显示信息问题解决
    两个方法:第一个,1.将打开nodejs文件夹(如果你是安装到D盘,就打开D盘!就是nodejs文件所在目录)2.分别右击该文件,点击列表属性,选择安全,编辑,勾选写入,应用,确定。(目的是为了一会要......
  • 2022 icpc 沈阳站 记录(非题解)
    赛前大概是赛前三周才突然知道拥有了比赛机会。赛前训练和vp频率很高,有一段时间cf上都是绿的。比赛的那一周只有一天没在vp,到了周六热身赛我人都有点麻木。(可能正赛也是......
  • 洛谷刷题_ISBN 号码
    P1055[NOIP2008普及组]ISBN号码题目链接:https://www.luogu.com.cn/problem/P1055这道题从题意上来说还是比较简单的,刚开始想用整形直接输入一个一个数字,没有想到scan......