POJ1113 Wall

题目：http://acm.pku.edu.cn/JudgeOnline/problem?id=1113

思路：求出凸包加上一个圆的周长，第一道凸包题

 

 

1. import java.io.BufferedReader;
2. import java.io.IOException;
3. import java.io.InputStreamReader;
5. public class Main {
7.     class Point {
8.         int x;
9.         int y;
11.         public Point(int x, int y) {
12.             this.x = x;
13.             this.y = y;
14.         }
15.     }
17.     public static void main(String[] args) throws IOException {
18.         Main main = new Main();
19.         BufferedReader read = new BufferedReader(new InputStreamReader(
20.                 System.in));
21.         String[] s;
22.         s = read.readLine().split(" ");
23.         int n = Integer.parseInt(s[0]);
24.         int l = Integer.parseInt(s[1]);
25.         int x, y;
26.         Point[] p = new Point[n];
27.         for (int i = 0; i < n; i++) {
28.             s = read.readLine().split(" ");
29.             x = Integer.parseInt(s[0]);
30.             y = Integer.parseInt(s[1]);
31.             p[i] = main.new Point(x, y);
32.         }
33.         Point[] ch = new Point[n];
34.         int len = 0;
35.         if (n >= 3) {
36.             len = Graham_scan(p, ch, n);
37.         }
38.         double sum = 0;
39.         for (int i = 0; i < len - 1; i++) {
40.             sum += distance(ch[i], ch[i + 1]);
41.         }
42.         if (len > 1) {
43.             sum += distance(ch[len - 1], ch[0]);
44.         }
45.         sum += 2 * l * Math.PI;
46.         System.out.println(Math.round(sum));
48.     }
50.     public static double multiply(Point p1, Point p2, Point p0) {
51.         return ((p1.x - p0.x) * (p2.y - p0.y) - (p2.x - p0.x) * (p1.y - p0.y));
52.     }
54.     public static double distance(Point p1, Point p2) {
55.         return (Math.sqrt((p1.x - p2.x) * (p1.x - p2.x) + (p1.y - p2.y)
56.                 * (p1.y - p2.y)));
57.     }
59.     public static int Graham_scan(Point[] PointSet, Point[] ch, int n) {
60.         int i, j, k = 0, top = 2;
61.         Point tmp;
62.         for (i = 1; i < n; i++)
63.             if ((PointSet[i].y < PointSet[k].y)
64.                     || ((PointSet[i].y == PointSet[k].y) && (PointSet[i].x < PointSet[k].x)))
65.                 k = i;
66.         tmp = PointSet[0];
67.         PointSet[0] = PointSet[k];
68.         PointSet[k] = tmp;
69.         for (i = 1; i < n - 1; i++) {
70.             k = i;
71.             for (j = i + 1; j < n; j++)
72.                 if ((multiply(PointSet[j], PointSet[k], PointSet[0]) > 0)
73.                         || ((multiply(PointSet[j], PointSet[k], PointSet[0]) == 0) && (distance(
74.                                 PointSet[0], PointSet[j]) < distance(
75.                                 PointSet[0], PointSet[k]))))
76.                     k = j;
77.             tmp = PointSet[i];
78.             PointSet[i] = PointSet[k];
79.             PointSet[k] = tmp;
80.         }
81.         ch[0] = PointSet[0];
82.         ch[1] = PointSet[1];
83.         ch[2] = PointSet[2];
84.         for (i = 3; i < n; i++) {
85.             while (top > 0 && multiply(PointSet[i], ch[top], ch[top - 1]) >= 0)
86.                 top--;
87.             ch[++top] = PointSet[i];
88.         }
89.         return top + 1;
90.     }
91. }