【前言】话说好久没有写题解了。到暑假了反而忙。o(╯□╰)o

【原题】

2458: [BeiJing2011]最小三角形

Time Limit: 10 Sec  Memory Limit: 128 MB

Submit: 574  Solved: 177

[

id=2458" style="color:blue; text-decoration:none">Submit][Status]

Description

Xaviera如今遇到了一个有趣的问题。

平面上有N个点。Xaviera想找出周长最小的三角形。

因为点许多。分布也很乱,所以Xaviera想请你来解决问题。

为了减小问题的难度,这里的三角形也包含共线的三点。

Input

第一行包括一个整数N表示点的个数。

接下来N行每行有两个整数,表示这个点的坐标。

Output

输出仅仅有一行。包括一个6位小数,为周长最短的三角形的周长(四舍五入)。

Sample Input

4

1 1

2 3

3 3

3 4


Sample Output

3.414214




HINT

100%的数据中N≤200000。

Source

【分析】今天新学了解决这类问题的方法——分治。

没错,就是分治。

先讲一下n是10^5级别的平面近期点对吧(CF 245 DIV 2 D)。非常easy懂。

具体的原理能够參考这个博客。讲的非常具体。

(非常多时候仅仅要感性认识原理就可以)以下讲一下具体做法。

①对于平面上的点。按x坐标排序(这是永久排序)。

②每次递归(l。r)。函数的返回值是第l个到第r个之间的全部点的近期点对。

③假设l=r。那么返回无穷大;假设l+1=r。就直接返回两个点的距离。

④每次先递归(l。mid)和(mid+1,r)。显然,这两个会有两个返回值。最好还是设为d1和d2。

首先我们设D=MIN(d1,d2)。即当前的最优值临时是D。

⑤显然,另一种情况。左边那块的某个点和右边那块的某个点产生关系。那么,我们能够从mid这个位置向左跑到mid-D,向右跑到mid+D,然后把这一段中的点都拎出来——由于仅仅有这两段中的点才有可能产生小于D的贡献。

⑥这时候我们要意识到潜在复杂度的保证(实际原理也不难懂呵)。

首先,假设直接枚举两两点要N^2。我们先把拎出来的点按y排序。

(NlogN)然后看似也是N^2的枚举,仅仅是加了一个优化(从底下開始枚举i。假设Y[J]-Y[I]>D就直接break)——这样可证明差点儿是线性。

总复杂度N*logN*logN。

再讲一下本题。也是几乎相同道理。

由于是三角形,我们把一些细节改一下就可以。

【代码】

#include<cstdio>
#include<cmath>
#include<algorithm>
#define N 200005
#define INF 210000000000.0
using namespace std;
struct arr{int x,y;}a[N],num[N];int n,i,Test;
inline bool cmpx(const arr &a,const arr &b){return a.x<b.x;}
inline bool cmpy(const arr &a,const arr &b){return a.y<b.y;}
inline double dis(const arr &a,const arr &b){return sqrt((a.x-b.x)*1.*(a.x-b.x)+(a.y-b.y)*1.*(a.y-b.y));}
inline double work(int l,int r)
{
if (l==r) return INF;
if (l+1==r) return INF;
if (l+2==r) return dis(a[l],a[l+1])+dis(a[l+1],a[r])+dis(a[l],a[r]);
int mid=(l+r)>>1;
double d1=work(l,mid),d2=work(mid+1,r);
double D=min(d1,d2),ans=D,DD=D/2.0;int cnt=0;
for (int i=l;i<=r;i++)
if (fabs(a[mid].x-a[i].x)<=DD) num[++cnt]=a[i];
sort(num+1,num+cnt+1,cmpy);
for (int i=1;i<cnt-1;i++)
for (int j=i+1;j<cnt;j++)
{
if (num[j].y-num[i].y>DD) break;
for (int k=j+1;k<=cnt;k++)
{
if (num[k].y-num[i].y>DD) break;
double temp=dis(num[i],num[j])+dis(num[i],num[k])+dis(num[j],num[k]);
if (temp<ans) ans=temp;
}
}
return ans;
}
int main()
{
read(n);//读入优化就不贴了。 for (i=1;i<=n;i++)
read(a[i].x),read(a[i].y);
sort(a+1,a+n+1,cmpx);
printf("%.6lf",work(1,n));
return 0;
}

bzoj 2458: [BeiJing2011]最小三角形 题解的更多相关文章

  1. BZOJ 2458: [BeiJing2011]最小三角形 | 平面分治

    题目: 给出若干个点 求三个点构成的周长最小的三角形的周长(我们认为共线的三点也算三角形) 题解: 可以参考平面最近点对的做法 只不过合并的时候改成枚举三个点更新周长最小值,其他的和最近点对大同小异 ...

  2. [BZOJ]2458: [BeiJing2011]最小三角形

    题目大意:给出平面上n个点,求最小的由这些点组成的三角形的周长.(N<=200,000) 思路:点按x坐标排序后分治,每次取出与排在中间的点的横坐标相差不超当前答案一半的点,按y坐标排序后再暴力 ...

  3. BZOJ 2458: [BeiJing2011]最小三角形 (分治)

    分治就是了. 类似于分治找最近/远点对. CODE #include <bits/stdc++.h> using namespace std; const double eps = 1e- ...

  4. bzoj-2458 2458: [BeiJing2011]最小三角形(计算几何+分治)

    题目链接: 2458: [BeiJing2011]最小三角形 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 1101  Solved: 380 Des ...

  5. 分治 - 计算几何 - BZOJ2458,[BeiJing2011]最小三角形

    http://www.lydsy.com/JudgeOnline/problem.php?id=2458 [BeiJing2011]最小三角形 描述 Frisk现在遇到了一个有趣的问题. 平面上有N个 ...

  6. bzoj2458: [BeiJing2011]最小三角形(分治+几何)

    题目链接:bzoj2458: [BeiJing2011]最小三角形 学习推荐博客:分治法编程问题之最接近点对问题的算法分析 题解:先将所有点按x值排列,然后每次将当前区间[l,r]分成左右两半递归求解 ...

  7. BZOJ2458:[BJOI2011]最小三角形——题解

    http://www.lydsy.com/JudgeOnline/problem.php?id=2458 Description Xaviera现在遇到了一个有趣的问题. 平面上有N个点,Xavier ...

  8. BZOJ2458 Beijing2011最小三角形(分治)

    类似于平面最近点对,考虑分治,即分别计算分割线两侧的最小三角形再考虑跨过线的三角形. 复杂度证明也是类似的,对于某一个点,在另一侧可能与其构成最小三角形的点在一个d*d/2的矩形内(两边之和大于第三边 ...

  9. [BZOJ2458][BeiJing2011]最小三角形(分治)

    求平面上n个点组成的周长最小的三角形. 回忆平面最近点对的做法,找到横坐标的中点mid分治到两边,合并时考虑离mid横坐标不超过当前最小值d的所有点,按y排序后暴力更新答案. 这个题也一样,先分治到两 ...

随机推荐

  1. 在C#中将数字转换成中文

    上篇我们讲了在MSSQL中将数字转换成中文,这篇我们讲讲在C#中将数字转换成中文 下篇将讲一下如何将金额转换成中文金额,废话不多说,具体代码如下: /// <summary> /// 数字 ...

  2. 同步mysql数据到ElasticSearch的最佳实践

    Elasticsearch是一个实时的分布式搜索和分析引擎.它可以帮助你用前所未有的速度去处理大规模数据.ElasticSearch是一个基于Lucene的搜索服务器.它提供了一个分布式多用户能力的全 ...

  3. Future使用场景与分析

    前面分享了CountDownLatch的用法,但是由于分享过程中,发现有些朋友,问我Future与CountDownLatch的有什么区别? 答案:只是concurrent包下的并发帮助工具类,两者并 ...

  4. 关于Date数据类型格式化的转换

    例如:   jsp页面读取数据库中日期格式的列时可能显示为1988-05-03 00:00:00 格式,但是我们不想要后面的00:00:00时间,只想要前面的年月日,那么该如何做出修改呢? 方法一:我 ...

  5. Python与数据结构[0] -> 链表/LinkedList[2] -> 链表有环与链表相交判断的 Python 实现

    链表有环与链表相交判断的 Python 实现 目录 有环链表 相交链表 1 有环链表 判断链表是否有环可以参考链接, 有环链表主要包括以下几个问题(C语言描述): 判断环是否存在: 可以使用追赶方法, ...

  6. thymeleaf初步使用

    thymeleaf模板引擎初步使用 #thymelea模板配置 spring.thymeleaf.prefix=classpath:/templates/ spring.thymeleaf.suffi ...

  7. [POI2014]Salad Bar

    题目大意: 一个长度为$n(n\leq10^6)$的字符串,每一位只会是$p$或$j$.你需要取出一个子串$S$(从左到右或从右到左一个一个取出),使得不管是从左往右还是从右往左取,都保证每时每刻已取 ...

  8. NOIP 2017 赛后反思 [补档]

    首先写一下比赛的情况: D1: T1: 之前做过类似的题目, 因而知道大致的结论, 迅速完成. T2: 貌似直接模拟就可以了, 涉及到字符串信息提取, 比较麻烦, 因而想放到最后做. T3: 非常简洁 ...

  9. 四. Java继承和多态4. 多态和动态绑定

    在Java中,父类的变量可以引用父类的实例,也可以引用子类的实例. 请读者先看一段代码: public class Demo { public static void main(String[] ar ...

  10. pr_debug、dev_dbg等动态调试一

    内核版本:Linux-3.14 作者:彭东林 邮箱:pengdonglin137@163.com pr_debug: #if defined(CONFIG_DYNAMIC_DEBUG) /* dyna ...