最近对问题的分治算法（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);

}