问题描述:
1023 The Best Polygon (35 分)
An n-gon is a polygon with n sides. For example, a triangle is a 3-gon. Now you are asked to find the best n-gon in a given convex N-gon. The vertices of the n-gon are selected from vertices of the N-gon. The n-gon you are supposed to find must have the largest area among all possible n-gons, which means that it approximates the N-gon the best. The figure below shows the best 6-gon (the shaded part) in a 10-gon.
Input Specification:
Each input file contains one test case. For each case, the first line gives two positive integers N ( 6≤N≤300 ) which is the total number of vertices of the convex N-gon, and n ( 6≤n≤min(10,N) ) which is the total number of vertices of the approximating convex n-gon. Then N lines follow, the i-th line gives the 2-D coordinates (x, y) of the i-th vertex ( i=0,⋯,N−1 ). The x and y coordinates in a line are real numbers with their absolute values no more than 1000, and they are separated by a space.
Output Specification:
Print in a line all the vertex indices of the best n-gon in descending order. All the numbers are separated by a space and there must be no extra space at the beginning or the end of the line. It is guaranteed that the solution is unique.
Sample Input:
10 6
133.0 1.0
544.0 71.0
558.0 206.0
536.0 338.0
463.0 436.0
330.0 503.0
188.0 499.0
305.0 2.0
55.0 410.0
2.0 140.0
Sample Output:
9 8 5 3 1 0
时隔一年之后,博主终于在这题上拿到了AC。。。
博主第一次看到这题大约是2018年3月下旬。那时候,博主当然用了自己唯一掌握的的dfs程序,经过了数天的努力,终于在2,4,5号测试点获得了运行超时 。显然dfs的时间复杂度太高。。。身心俱疲的博主抱着试一试的心理,在自己学校的评优网站上晒出了此题。我校的ACM大佬迅速地给出了动态规划的算法。因为博主那时候比较忙,就拖着没有去仔细地阅读她的程序。。。
(图中的关键信息已经被隐去。。。)
这一拖就拖了一年。。。一年之后的今天,再次闲下来的博主再一次拿出了此题。花了两天,写出了自己的动态规划程序。。。
这一题总共分为三步:1.利用极角排序算法把题目中给的散乱的点列以顺/逆时针方向排序。2.用动态规划的方法计算出最大面积。3.回溯动态规划过程,找到这个过程所选取的点。
关于1:对于给出的所有点,先计算出它们的几何中心。调用math库中的atan2()函数把所有点对于几何中心的极角计算出来并依次排序;
关于2:开三维数组vdp[i][j][k](i<j,k<j-i),表示以第i个顶点为起始,以第j个顶点为结尾的子图中所有包含i-j边的k边形的最大面积。状态转移方程如下:vdp[i][j][k]=max[vdp[i][jj][k-1]+S(i,jj,j)] 其中jj在(i,j)内变动,S(i,jj,j)表示三角形(i,jj,j)的面积的两倍。
关于3:vdp[i][j][k]不仅保存一个表示最大面积的double变量s,也保存每一次状态转移int型的jj。这样,最后在vdp[m][][]里找到最大面积时,可以通过jj找到上一个状态的位置。
关于优化:在初始化vdp数组,即k=3时,可以通过旋转卡壳法求出面积最大的三角形。
(p.s.当年ACM大佬的动态规划在2号测试点没能通过的原因是:博主给她的极角排序算法存在隐藏的bug。。。)
久违的AC代码:
#include<bits/stdc++.h>
using namespace std;
double x,y;
int n,m;
const double eps=1e-8;
inline int sgn(double d)
{return d<-eps?-1:d>eps;}
struct point
{
double x;
double y;
int id;
} po;
struct cmp1
{
bool operator() (const point& p1,const point& p2) const
{
return atan2((p1.y-y),(p1.x-x))<atan2((p2.y-y),(p2.x-x));
}
};
vector<point> v;
vector<double> va;
inline double area(int n1,int n2)
{
if(!sgn(va[n1*n+n2]))
{
va[n1*n+n2]=v[n1].x*v[n2].y-v[n1].y*v[n2].x;
va[n2*n+n1]=-va[n1*n+n2];
}
return va[n1*n+n2];
}
int main()
{
// freopen("data1.txt","r",stdin);
scanf("%d %d",&n,&m);
va.resize(n*n,0.0);
x=0.0;
y=0.0;
for(int i=0;i<n;i++)
{
scanf("%lf %lf",&po.x,&po.y);
x+=po.x;
y+=po.y;
po.id=i;
v.emplace_back(po);
}
x=x/n;
y=y/n;
sort(v.begin(),v.end(),cmp1());
vector<pair<double,int> > vdp((n-2)*(n-2)*(m-2),make_pair(0.0,0));
for(int i=0;i<n-2;i++)
{
for(int j=i+2;j<n;j++)
{
double maxs=0.0;
int k=i+1;
for(;k<j;)
{
double s=area(i,k)+area(k,j)+area(j,i);
if(sgn(s-maxs)>0)
{
k++;
maxs=s;
}
else
break;
}
vdp[i*(n-2)+j-2].first=maxs;
if(k<j)
vdp[i*(n-2)+j-2].second=k-1;
else
vdp[i*(n-2)+j-2].second=j-1;
}
}
for(int l=1;l<m-2;l++)
{
for(int i=0;i<n-l-2;i++)
{
double maxs=0.0;
int kk=0;
for(int j=i+l+2;j<n;j++)
{
for(int jj=i+l+1;jj<j;jj++)
if(sgn(maxs-vdp[(l-1)*(n-2)*(n-2)+i*(n-2)+jj-2].first-area(jj,j)-area(j,i)-area(i,jj))<0)
{
maxs=vdp[(l-1)*(n-2)*(n-2)+i*(n-2)+jj-2].first+area(jj,j)+area(j,i)+area(i,jj);
kk=jj;
}
vdp[l*(n-2)*(n-2)+i*(n-2)+j-2].first=maxs;
vdp[l*(n-2)*(n-2)+i*(n-2)+j-2].second=kk;
}
}
}
double maxs=0.0;
int ii=0,jj=0;
for(int i=0;i<(n-m+1);i++)
{
for(int j=i+m-3;j<n-2;j++)
{
if(sgn(maxs-vdp[(m-3)*(n-2)*(n-2)+i*(n-2)+j].first)<0)
{
maxs=vdp[(m-3)*(n-2)*(n-2)+i*(n-2)+j].first;
ii=i;
jj=j;
}
}
}
vector<int> vr;
vr.emplace_back(v[ii].id);
vr.emplace_back(v[jj+2].id);
for(int i=m-3;i>=0;i--)
{
vr.emplace_back(v[vdp[i*(n-2)*(n-2)+ii*(n-2)+jj].second].id);
jj=vdp[i*(n-2)*(n-2)+ii*(n-2)+jj].second-2;
}
sort(vr.begin(),vr.end(),greater<int>());
printf("%d",vr[0]);
for(int i=1;i<m;i++)
printf(" %d",vr[i]);
return 0;
}