Friday 20 January 2017

Chapter 22 Exercise 7, Introduction to Java Programming, Tenth Edition Y. Daniel LiangY.

22.7 (Closest pair of points) Section  22.8 introduced an algorithm for finding the closest pair of points using a divide-and-conquer approach. Implement the
algorithm to meet the following requirements:

■Define the classes Point and CompareY in the same way as in Programming
Exercise 20.4.
■Define a class named Pair with the data fields p1 and p2 to represent two
points, and a method named getDistance() that returns the distance
between the two points.
■Implement the following methods:

/** Return the distance of the closest pair of points */
public static Pair getClosestPair(double[][] points)
/** Return the distance of the closest pair of points */
public static Pair getClosestPair(Point[] points)
/** Return the distance of the closest pair of points
* in pointsOrderedOnX[low..high]. This is a recursive
* method. pointsOrderedOnX and pointsOrderedOnY are
* not changed in the subsequent recursive calls.
*/
public static Pair distance(Point[] pointsOrderedOnX,
int low, int high, Point[] pointsOrderedOnY)

/** Compute the distance between two points p1 and p2 */
public static double distance(Point p1, Point p2)
/** Compute the distance between points (x1, y1) and (x2, y2) */
public static double distance(double x1, double y1,
double x2, double y2)


import java.util.Comparator;
class CompareY implements Comparator<Point> {
 
 @Override
 public int compare(Point o1, Point o2) {
  if (o1.getY() > o2.getY()) {
   return 1;
  } else if (o1.getY() < o2.getY()) {
   return -1;
  } else {
   if (o1.getX() > o2.getX()) {
    return 1;
   } else if (o1.getX() < o2.getX()) {
    return -1;
   } else {
    return 0;
   }
  }
 } 

}

class Point implements Comparable<Point> {
 private double x;
 private double y;  
 
 public Point() {
  this(0, 0);
 } 
 
 public Point(double x, double y) {
  super();
  this.x = x;
  this.y = y;
 }
 
 public double getX() {
  return x;
 }
 public void setX(double x) {
  this.x = x;
 }
 public double getY() {
  return y;
 }
 public void setY(double y) {
  this.y = y;
 }

 @Override
 public int compareTo(Point o) {
  if (x > o.x) {
   return 1;
  } else if (x < o.x) {
   return -1;
  } else {
   if (y > o.y) {
    return 1;
   } else if (y < o.y) {
    return -1;
   } else {
    return 0;
   }
  }
 }  
 
 @Override
 public String toString() {
  return "x: " + x + ",\ty: " + y;
 }
}

import java.util.ArrayList;
class Pair {
 private Point p1;
 private Point p2;
 public Point getP1() {
  return p1;
 }
 public void setP1(Point p1) {
  this.p1 = p1;
 }
 public Point getP2() {
  return p2;
 }
 public void setP2(Point p2) {
  this.p2 = p2;
 }
 public Pair(Point p1, Point p2) {
  super();
  this.p1 = p1;
  this.p2 = p2;
 }
 public double getDistance() {
  return distance(p1, p2);
 }
 
 public static double distance(Point p1, Point p2) {
  return distance(p1.getX(), p1.getY(), p2.getX(), p2.getY());
 }
  
 public static double distance(double x1, double y1, double x2, double y2) {
  return Math.sqrt((x1 - x2) * (x1 - x2) + (y1 - y2) * (y1 - y2));
 }
 
 public static Pair distance(Point[] pointsOrderedOnX, int low, int high, Point[] pointsOrderedOnY) {
  if(low >= high) {
   return null;
  } else if (low + 1 == high) {
   return new Pair(pointsOrderedOnX[low], pointsOrderedOnX[high]);
  }
  
  int halfSize = (low + high) / 2;
  Pair p1 = distance(pointsOrderedOnX, low, halfSize, pointsOrderedOnY);
  Pair p2 = distance(pointsOrderedOnX, halfSize + 1, high, pointsOrderedOnY);
  
  double distance = 0;
  Pair p = null;

  if (p1 == null && p2 == null) {
   distance = Double.MAX_VALUE;
  } else if (p1 == null) {
   distance = p2.getDistance();
   p = p2;
  } else if (p2 == null) {
   distance = p1.getDistance();
   p = p1;
  } else {
   distance = Math.min(p1.getDistance(), p2.getDistance());
   p = ((p1.getDistance() <= p2.getDistance())? p1 : p2);
  }

  ArrayList<Point> stripL = new ArrayList<Point>();
  ArrayList<Point> stripR = new ArrayList<Point>();
  for (int i = 0; i < pointsOrderedOnY.length; i++) {
   if ((pointsOrderedOnY[i].getX() <= pointsOrderedOnX[halfSize].getX())&&
     (pointsOrderedOnY[i].getX() >= pointsOrderedOnX[halfSize].getX() - distance)) {
    stripL.add(pointsOrderedOnY[i]);
   } else {
    stripR.add(pointsOrderedOnY[i]);
   }
  }

  
  double d3 = distance;
  int r = 0;
  for (int i = 0; i < stripL.size(); i++) {
   while (r < stripR.size() && stripL.get(i).getY() > stripR.get(r).getY() + distance) {
    r++;
   }

   int r1 = r;
   while (r1 < stripR.size() && stripR.get(r1).getY() <= stripL.get(i).getY() + distance) {
    if (d3 > distance(stripL.get(i), stripR.get(r1))) {
     d3 = distance(stripL.get(i), stripR.get(r1));
     p.p1 = stripL.get(i);
     p.p2 = stripR.get(r1);
    }
    r1++;
   }
  }

  return p;
 }
 
 public static Pair getClosestPair(Point[] points) {
  java.util.Arrays.sort(points);
  Point[] pointsOrderedOnY = points.clone();
  java.util.Arrays.sort(pointsOrderedOnY, new CompareY());
  return distance(points, 0, points.length - 1, pointsOrderedOnY);
 }
 
 public static Pair getClosestPair(double[][] points) {
  Point[] points2 = new Point[points.length];
  for (int i = 0; i < points.length; i++) {
   points2[i] = new Point(points[i][0], points[i][1]);
  }
  return getClosestPair(points2);
 }
}

public class Exercise07 {
 public static void main(String[] args) {
  Point[] points = new Point[500];
  for (int i = 0; i < points.length; i++) {
   points[i] = new Point((int)(Math.random() * 1000000), (int)(Math.random() * 1000000));
  }
  Pair pair = Pair.getClosestPair(points);
  System.out.println(pair.getP1());
  System.out.println(pair.getP2());
  System.out.println(pair.getDistance());
 }
}

No comments :

Post a Comment