[FJSC2014]圈地

浏览数：15 / 时间：2015年06月09日

【题目描述】

　　2维平面上有n个木桩，黄学长有一次圈地的机会并得到圈到的土地，为了体现他的高风亮节，他要使他圈到的土地面积尽量小。圈地需要圈一个至少3个点的多边形，多边形的顶点就是一个木桩，圈得的土地就是这个多边形内部的土地。（因为黄学长非常的神，所以他允许圈出的第n点共线，那样面积算0）

【输入格式】

　　第一行一个整数n，表示木桩个数。
　　接下来n行，每行2个整数表示一个木桩的坐标，坐标两两不同。

【输出格式】

　　仅一行，表示最小圈得的土地面积，保留2位小数。

【样例输入】

样例1：
3

0 0

0 1

1 0

样例2：
4

0 0

0 1

0 2

1 1

【样例输出】

样例1：
0.50
样例2：
0.00
【数据范围】

　　对于10%的数据，n<=100；
　　对于30%的数据，n<=300；
　　对于50%的数据，n<=500；
　　对于100%的数据，n<=1000。

Solution

　　这题被我乱搞搞掉了= =

　　暴力n^3不多说。但是有很多时候的情况是没有用的。于是我们把这些点分成sqrt(n)块，块内暴力，轻松愉快。

　　这样做不靠谱，所以我们可以随机旋转坐标系，rand个四五十次就可以把这题水过了。

出题人题解：

　　显然，这时候暴力枚举会T。于是我们转变一下思路，如果我们确定了2个点以后，第三个点有必要去盲目的枚举吗？答案是否定的。实际上我们把经过这两点的线看成一个斜率，把他当成y轴你会发现第三个点明显是在坐标系左右找一个离”y轴”最近的点来算面积更新答案。然后我们可以继续思考，发现我们可以把点按照某个斜率当成”y轴”进行“从左到右”的排序，这样当2点共线的时候，用这两个点的左右2个点去更新答案就好了。也就是说我们采用旋转坐标系的方法，一开始按x坐标排好序，认为直接用竖着的那条斜率，然后维护的话每次其实当两点共线后只要交换他们就能得到斜率转过该事件点的序列。所以我们可以预处理出所有可行的斜率，当成事件点，不断转动坐标系更新答案就好。这样复杂度只有n^2，期望得分100.(这确实只是个暴力的优化啊。。。不要砸我T_T)

（好像很复杂的样子。。。）

 1 #include<cstdio>
 2 #include<ctime>
 3 #include<cstdlib>
 4 #include<cmath>
 5 #include<algorithm>
 6 const double pi=3.1415926535897932384626;
 7 struct P{double x,y;}a[1010],z,b[1010];int n,k;double ans=1e10,tmp,t;
 8 bool operator<(P const&a,P const&b){return a.x<b.x||(a.x==b.x&&a.y<b.y);}
 9 P operator+(P const&a,P const&b){return (P){a.x+b.x,a.y+b.y};}
10 P operator-(P const&a,P const&b){return (P){a.x-b.x,a.y-b.y};}
11 P operator*(P const&a,double p){return (P){a.x*p,a.y*p};}
12 P operator*(P const&a,P const&b){return (P){a.x*b.x-a.y*b.y,a.x*b.y+a.y*b.x};}
13 double cross(P const&a,P const&b){return a.x*b.y-a.y*b.x;}
14 double dot(P const&a,P const&b){return a.x*b.x+a.y*b.y;}
15 int main()
16 {
17     scanf("%d",&n);int i,j;k=(int)sqrt(n)+10;//srand(time(0));
18     for(i=1;i<=n;i++)scanf("%lf%lf",&a[i].x,&a[i].y);
19     if(n<=300)
20     {
21         for(i=1;i<n;i++)
22         for(j=i+1;j<n;j++)
23         for(k=j+1;k<=n;k++)
24         {
25             tmp=std::abs(cross(a[k]-a[i],a[j]-a[i]));
26             if(ans>tmp)ans=tmp;
27         }
28         printf("%.2lf\n",ans/2.0);
29     }
30     else{int p1=31,p2=23;
31         for(int tt=1;tt<=p1;tt++)
32         {
33             z=(P){cos(pi/p2*tt),sin(pi/p2*tt)};
34             for(i=1;i<=n;i++)b[i]=a[i]*z;
35             std::sort(b+1,b+1+n);
36             for(int kk=1;(kk-1)*k<=n;kk++)
37             for(i=(k*(kk-1)+1);i<n&&i<k*kk;i++)
38             for(j=i+1;j<n&j<k*kk;j++)
39             for(int l=j+1;l<=n&&l<k*kk;l++)
40             {
41                 tmp=std::abs(cross(b[l]-b[i],b[j]-b[i]));
42                 if(ans>tmp)ans=tmp;
43             }
44         }
45         printf("%.2lf\n",ans/2.0);
46     }
47 }