最近对问题的分治算法(C++)

 

#include < stdio.h >
#include
< sstream >
#include 
< map >
#include 
< iostream >
#include 
< math.h >
using   namespace  std;
#define  N 5;
#define  M 5000;
typedef  
struct  Point
{
 
int  x;
 
int  y;
} Point;
Point P[
50 ] ;
int  GetNearest( int  i, int  j)
{  
 
if (i == j - 1 )
 {
  
return  (P[i].x - P[j].x) * (P[i].x - P[j].x) + (P[i].y - P[j].y) * (P[i].y - P[j].y);
 }
 
else  
 {
  
int  mid = (i + j) / 2 ;
  
int  n1,n2;
  
if (i < mid)
  {
    n1
= GetNearest(i,mid);    
  }
  
else
  {
   n1
= M;
  }
  
if (j > mid + 1 )
  {
    n2
= GetNearest(mid + 1 ,j);    
  }
  
else
  {
   n2
= M;
  }
  
int  min = n1;
  
if (n2 < n1)
  {
   min
= n2;
  }
  
for ( int  k = i;k <= mid;k ++ )
  {
   
if (P[k].x >= P[mid].x - n1)
   {
    
for ( int  l = mid + 1 ;l <= j;l ++ )
    {
     
if (P[l].x <= P[mid].x + n1)
     {
      
int  dis = (P[i].x - P[j].x) * (P[i].x - P[j].x) + (P[i].y - P[j].y) * (P[i].y - P[j].y);
      
if (dis < min)
      {
       min
= dis;
      }
     }
    }
   }
  }
  
return  min;
 }
}
void  main()
{  
 
 
for ( int  i = 0 ;i < 50 ;i ++ )
 {
  P[i].x
= i;
  P[i].y
= 2 * i;
 }
 P[
5 ].x = 5 ;
 P[
5 ].y = 9 ;
 
int  mindis = GetNearest( 0 , 49 );

}
 
最近问题是指在平面上给定n个,求其中距离最近的两个之间的距离。 分治算法思路如下: 1. 将按照x坐标排序,递归地将集分成左右两部分。 2. 对于左右两个子集,分别递归求解最近问题。 3. 计算跨越左右两部分的最近对,这个最近对必然是左右两个子集中最近对的其中一个。 4. 返回三个最近对中的最小值。 代码如下: #include <iostream> #include <algorithm> #include <cmath> #include <vector> using namespace std; struct point { double x, y; }; bool cmp(point a, point b) { return a.x < b.x; } double dist(point a, point b) { double dx = a.x - b.x; double dy = a.y - b.y; return sqrt(dx * dx + dy * dy); } double brute_force(vector<point>& a, int left, int right) { double min_dist = 1e9; for (int i = left; i <= right; i++) { for (int j = i + 1; j <= right; j++) { min_dist = min(min_dist, dist(a[i], a[j])); } } return min_dist; } double closest_pair(vector<point>& a, int left, int right) { if (left == right) { return 1e9; } if (left + 1 == right) { return dist(a[left], a[right]); } int mid = (left + right) >> 1; double d = min(closest_pair(a, left, mid), closest_pair(a, mid + 1, right)); vector<point> strip; for (int i = left; i <= right; i++) { if (abs(a[i].x - a[mid].x) < d) { strip.push_back(a[i]); } } sort(strip.begin(), strip.end(), [](point a, point b) { return a.y < b.y; }); int n = strip.size(); for (int i = 0; i < n; i++) { for (int j = i + 1; j < n && strip[j].y - strip[i].y < d; j++) { d = min(d, dist(strip[i], strip[j])); } } return d; } int main() { int n; cin >> n; vector<point> a(n); for (int i = 0; i < n; i++) { cin >> a[i].x >> a[i].y; } sort(a.begin(), a.end(), cmp); cout << closest_pair(a, 0, n - 1) << endl; return 0; }
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值