题目大意:
N个点的半径为M的圆,求最小周长的方案使他们都被包含在所画的图形中(可以相隔)。
因为N<=9,显然可以想到状压dp。一开始我死活想不出来求圆弧的面积(我是撒币吗),最后经dalao指点才想到圆弧的总长就是一个半径为M的圆的周长(显然得证法证明)。
那么可以得到方程f[i]=min(f[j]+f[i^j]),j为i的所有子集。先处理g[i]表示i组成一个图形的周长,周长是i中所有点的凸包总长+2πM。那么问题就迎刃而解了。
附代码(小学dalao好厉害):
#include<bits/stdc++.h>
using namespace std;
struct Node{int x,y;}a[15],s[15],x[15];
inline bool cmp(Node a,Node b){return a.x<b.x;}
int n,ls,lx,m;
inline double dis(Node l,Node r){return sqrt((l.x-r.x)*(l.x-r.x)+(l.y-r.y)*(l.y-r.y));}
double f[2049];
double yz;
int main(){
int tot=0;
f[0]=0;
while(scanf("%d%d",&n,&m),n+m){
for(int i=0;i<1024;i++) f[i]=1008610086;
yz=2*m*acos(-1);
for(int i=0;i<n;i++) scanf("%d%d",&a[i].x,&a[i].y);
sort(a,a+n,cmp);
for(int j=(1<<n)-1;j>=1;j--){
ls=lx=0;
for(int i=0;i<n;i++)
if((1<<i)&j){
if(!ls) s[ls++]=x[lx++]=a[i];
else{
while(ls>1&&(s[ls-1].y-s[ls-2].y)*(a[i].x-s[ls-1].x)<(a[i].y-s[ls-1].y)*(s[ls-1].x-s[ls-2].x)) ls--;
s[ls++]=a[i];
while(lx>1&&(x[lx-1].y-x[lx-2].y)*(a[i].x-x[lx-1].x)>(a[i].y-x[lx-1].y)*(x[lx-1].x-x[lx-2].x)) lx--;
x[lx++]=a[i];
}
}
f[j]=yz;
for(int i=1;i<ls;i++) f[j]+=dis(s[i],s[i-1]);
for(int i=1;i<lx;i++) f[j]+=dis(x[i],x[i-1]);
}
for(int i=1;i<(1<<n);i++)
for(int j=1;j<i;j++)
if(i|j==i) f[i]=min(f[i],f[i^j]+f[j]);
tot++;
printf("Case %d: length = %.2lf\n",tot,f[(1<<n)-1]);
}
return 0;
}