#include"stdafx.h"
#include<iostream>
#include<cmath>
#define TRUE 1
#define FALSE 0
using namespace std;
typedef struct Node//坐标点
{
 double x;
 double y;
}Node;  
typedef struct List
{
 Node* data;      //点
 int count;      //点的个数
}List;
typedef struct CloseNode
{
 Node a;
 Node b;     //计算距离的两个点
 double space;     //距离平方
}CloseNode;
int n;     //点的数目
//输入各点到List中
void create(List &L)
{
 cout << "请输入平面上点的数目:\n";
 cin >> n;
 L.count = n;
 L.data = new Node[L.count];      //动态空间分配
 cout << "输入各点坐标 :x_y):" << endl;
 for (int i = 0; i<L.count; ++i)
  cin >> L.data[i].x >> L.data[i].y;
}
//求距离的平方
double square(Node a, Node b)
{
 return ((a.x - b.x)*(a.x - b.x)) + ((a.y - b.y)*(a.y - b.y));
}
//冒泡排序
void BubbleSort(Node r[], int length)
{
 int change, n;
 n = length; change = TRUE;
 double b, c;
 for (int i = 0; i<n - 1 && change; ++i)
 {
  change = FALSE;
  for (int j = 0; j<n - i - 1; ++j)
  {
   if (r[j].x>r[j + 1].x)
   {
    b = r[j].x; c = r[j].y;
    r[j].x = r[j + 1].x; r[j].y = r[j + 1].y;
    r[j + 1].x = b; r[j + 1].y = c;
    change = TRUE;
   }
  }
 }
}
//分治法中先将坐标按X轴从小到大的顺序排列
void paixu(List L)
{
 BubbleSort(L.data, L.count);   //调用冒泡排序
}
//左右各距中线d的区域的最近对算法
void middle(const List & L, CloseNode &cnode, int mid, double midX)
{
 int i, j;    //分别表示中线左边,右边的点
 double d = sqrt(cnode.space);
 i = mid;
 while (i >= 0 && L.data[i].x >= (midX - d))    //在左边的d区域内
 {
  j = mid;
  while (L.data[++j].x <= (midX + d) && j <= L.count)    //在右边的d区域内
  {
   if (L.data[j].y<(L.data[i].y - d) || L.data[j].y>(L.data[i].y + d))   //判断纵坐标是否在左边某固定点的2d区域内
    continue;
   double space = square(L.data[i], L.data[j]);
   if (cnode.space>space)    //在满足条件的区域内依次判断
   {
    cnode.a = L.data[i];
    cnode.b = L.data[j];
    cnode.space = space;
   }
  }
  --i;
 }
}
//分治法求最近对
void DivideConquer(const List &L, CloseNode &closenode, int begin, int end)
{
 if (begin != end)
 {
  int mid = (begin + end) / 2;     //排列后的中间的那个点
  double midX = L.data[mid].x;
  DivideConquer(L, closenode, begin, mid);      //继续在左半边用分治法求最近对
  DivideConquer(L, closenode, mid + 1, end);      //继续在右半边用分治法求最近对
  middle(L, closenode, mid, midX);               //判断左右各距中线d的区域,是否有最近对
 }
}
void main()
{
 //初始化
 List list;
 CloseNode closenode;
 closenode.space = 10000;    
 create(list);    
 cout << "各点坐标为:" << endl;
 for (int i = 0; i<list.count; ++i)
  cout << "X=" << list.data[i].x << "   Y=" << list.data[i].y << "\n";
 cout << "用分治法求最近对:" << endl;
 paixu(list);
 cout << "经过排序后的各点:" << endl;
 for (int j = 0; j<list.count; ++j)
  cout << "X=" << list.data[j].x << "   Y=" << list.data[j].y << "\n";
 DivideConquer(list, closenode, 0, list.count - 1);
 cout << "最近对为点 (" << closenode.a.x << "," << closenode.a.y << ")和点(" << closenode.b.x << "," << closenode.b.y << ")\n" << "最近距离为: " << sqrt(closenode.space) << endl;
}
 

分治法求解最近对问题(c++)的更多相关文章

  1. [C++] 分治法之棋盘覆盖、循环赛日程表

    一.分治的基本思想 将一个难以直接解决的大问题,分割成一些规模较小的相同问题,以便各个击破,分而治之. 对于一个规模为 n 的问题,若问题可以容易地解决,则直接解决,否则将其分解为 k 个规模较小的子 ...

  2. Java算法——分治法

         一.基本概念 在计算机科学中,分治法是一种很重要的算法.字面上的解释是“分而治之”,就是把一个复杂的问题分成两个或更多的相同或相似的子问题,再把子问题分成更小的子问题……直到最后子问题可以简 ...

  3. 分治法避免定义多个递归函数,应该使用ResultType

    总结:对二叉树应用分治法时,应避免定义多个递归函数,当出现需要递归求解多种的结果时,尽量使用ResultType来让一次递归返回多种结果. 题目:Binary Tree Maximum Path Su ...

  4. 分治法(一)(zt)

    这篇文章将讨论: 1) 分治策略的思想和理论 2) 几个分治策略的例子:合并排序,快速排序,折半查找,二叉遍历树及其相关特性. 说明:这几个例子在前面都写过了,这里又拿出来,从算法设计的策略的角度把它 ...

  5. p1257 平面上最接近点对---(分治法)

    首先就是一维最接近点的情况... #include<iostream> #include<cstdio> #include<cstring> #include< ...

  6. 分治法及其python实现例子

    在前面的排序算法学习中,归并排序和快速排序就是用的分治法,分治法作为三大算法之一的,有非常多的应用例子. 分治法概念 将一个复杂的问题分成两个或更多的相同或相似的子问题,再把子问题分成更小的子问题-- ...

  7. 分治法 - Divide and Conquer

    在计算机科学中,分治法是一种很重要的算法.分治法即『分而治之』,把一个复杂的问题分成两个或更多的相同或相似的子问题,再把子问题分成更小的子问题……直到最后子问题可以简单的直接求解,原问题的解即子问题的 ...

  8. poj 3714 Raid【(暴力+剪枝) || (分治法+剪枝)】

    题目:  http://poj.org/problem?id=3714 http://acm.hust.edu.cn/vjudge/contest/view.action?cid=27048#prob ...

  9. python编写PAT 1007 Maximum Subsequence Sum(暴力 分治法 动态规划)

    python编写PAT甲级 1007 Maximum Subsequence Sum wenzongxiao1996 2019.4.3 题目 Given a sequence of K integer ...

随机推荐

  1. 消息队列——RabbitMQ学习笔记

    消息队列--RabbitMQ学习笔记 1. 写在前面 昨天简单学习了一个消息队列项目--RabbitMQ,今天趁热打铁,将学到的东西记录下来. 学习的资料主要是官网给出的6个基本的消息发送/接收模型, ...

  2. C语言 · 时间转换

    问题描述 给定一个以秒为单位的时间t,要求用"<H>:<M>:<S>"的格式来表示这个时间.<H>表示时间,<M>表示分 ...

  3. 在线浏览PDF之PDF.JS (附demo)

    平台之大势何人能挡? 带着你的Net飞奔吧!:http://www.cnblogs.com/dunitian/p/4822808.html#skill 下载地址:http://mozilla.gith ...

  4. Sublime配置python开发环境

    Package Control Package Control 是Sublime 里直接安装附加插件的包管理器.可以通过以下步骤手动安装: 1.点击Preferences > Browse Pa ...

  5. 微软新神器-Power BI横空出世,一个简单易用,还用得起的BI产品,你还在等什么???

    在当前互联网,由于大数据研究热潮,以及数据挖掘,机器学习等技术的改进,各种数据可视化图表层出不穷,如何让大数据生动呈现,也成了一个具有挑战性的可能,随之也出现了大量的商业化软件.今天就给大家介绍一款逆 ...

  6. 数据的双向绑定 Angular JS

    接触AngularJS许了,时常问自己一些问题,如果是我实现它,会在哪些方面选择跟它相同的道路,哪些方面不同.为此,记录了一些思考,给自己回顾,也供他人参考. 初步大致有以下几个方面: 数据双向绑定 ...

  7. linux上使用google身份验证器(简版)

    系统:centos6.6 下载google身份验证包google-authenticator-master(其实只是一个.zip文件,在windwos下解压,然后传进linux) #cd /data/ ...

  8. 萌新笔记——vim命令“=”、“d”、“y”的用法(结合光标移动命令,一些场合会非常方便)

    vim有许多命令,网上搜有一堆贴子.文章列举出各种功能的命令. 对于"="."d"."y",我在无意中发现了它们所具有的相同的一些用法,先举 ...

  9. Linux下用netstat查看网络状态、端口状态(转)

    转:http://blog.csdn.net/guodongdongnumber1/article/details/11383019 在linux一般使用netstat 来查看系统端口使用情况步.  ...

  10. 机器学习之sklearn——EM

    GMM计算更新∑k时,转置符号T应该放在倒数第二项(这样计算出来结果才是一个协方差矩阵) from sklearn.mixture import GMM    GMM中score_samples函数第 ...