Balloons in box

data:
3
0 0
10 10
3 5
2 1
7 7

1. #include<iostream>
2. #include<cmath>
3. #include<string>
4. #include<vector>
5. using namespace std;
7. const double PI = 3.14;
8. int N;
9. double lu[2];
10. double rd[2];
12. vector<vector<double> > data;
13. vector<vector<double> > dist;
15. //recursive search the area
16. double getArea(int cur,double* r,bool* visited,int* next,double area)
17. {
18.     double w = min(abs(lu[0] - data[cur][0]),abs(rd[0] - data[cur][0]));
19.     double h =  min(abs(lu[1] - data[cur][1]),abs(rd[1] - data[cur][1]));
20.     double radius = min(w,h);
22.     for(int i = 0;i<N;i++)
23.     {
24.         if(i == cur)
25.             continue;
26.         double d = dist[cur][i];
27.         if(visited[i])
28.         {
29.             radius = min(radius,d - r[i]);
30.         }
31.         else
32.         {
33.             radius = min(radius,d);
34.         }
35.     }
37.     area += PI * radius * radius;
38.     visited[cur]  = true;
39.     r[cur] = radius;
41.     double subarea = 0;
42.     int index = -1;
43.     for(int i = 0; i<N ;i++)
44.     {
45.         if(visited[i])
46.             continue;
48.         double aa = getArea(i,r, visited,next,area);
50.         if(aa > area)
51.         {
52.             area = aa;
53.             index = i;
54.         }
55.     }
56.     next[cur] = index;
57.     visited[cur] = false;
59.     return area;
60. }
62. int main()
63. {
64.     // read the input
65.     cin>>N;
66.     cin>>lu[0]>>lu[1];
67.     cin>>rd[0]>>rd[1];
69.     data.resize(N);
70.     for(int i = 0;i<N;i++)
71.     {
72.         data[i].resize(2);
73.         for(int j= 0;j<2;j++)
74.         {
75.             cin>>data[i][j];
76.         }
77.     }
79.     /*
80.      * calculate
81.      */
82.     
83.     //all distance
84.     dist.resize(N);
85.     for(int i = 0;i<N;i++)
86.     {
87.         dist[i].resize(N);
88.         for(int j  = 0;j<i;j++)
89.         {
90.             dist[i][j] = sqrt((data[i][0]- data[j][0])* (data[i][0]- data[j][0]) + (data[i][1]- data[j][1]) * (data[i][1]- data[j][1]) );
91.             dist[j][i] = dist[i][j];
92.         }
93.     }
95.     double* r = new double[N];
96.     bool* visited = new bool[N];
97.     int* next = new int[N];
99.     double maxArea = 0;
100.     int* order = new int[N];
101.     int index = -1;
102.     for(int i = 0;i< N;i++)
103.     {
104.         for(int j = 0;j<N;j++)
105.         {
106.             r[j]= 0;
107.             visited[j] = false;
108.             next[j] = -1;
109.         }
110.     // begin the search from the point i
111.         double area = getArea(i,r,visited,next,0);
112.         if(area >maxArea)
113.         {
114.             maxArea = area;
115.             index = i;
116.             for(int j = 0;j<N;j++)
117.             {
118.                 order[j] = next[j];
119.             }
120.         }
121.     }
123.     while(index!= -1)
124.     {
125.         cout<<index<<endl;
126.         index = order[index];
127.     }
129.     delete[] visited;
130.     delete[] r;
131.     delete[] next;
132. }
133.