可爱的(爆炸的)caioj.cn
题面描述
思路
简要题意,给定一些树,砍掉一些有价值的树后,得到一些木材,这些木材的长度,要不小于剩下的树围成的凸包的周长,现在要求满足最小价值时,最小砍伐树木的方案。
把状态压缩为二进制
v a l val val表示费用, s t a t e state state表示压缩的状态, n u m num num表示砍掉的树的数目, l e n len len为砍伐后得到的木材长度, u s e d used used表示凸包周长。
每一次更新一下就好了。
有一个剪枝,若当前状态价值大于 a n s v a l ans_{val} ansval,就直接跳过。
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<cstdlib>
using namespace std;
const int N=17;
const double inf=1e9;
struct node{double x,y,v,l;/*val,len*/}a[N],p[N],h[N];
bool cmp(node a,node b){return a.x==b.x?a.y<b.y:a.x<b.x;}
int tp,sta[N],n,m;
bool v[N];
double dis(node a,node b){return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));}
double mul(node p1,node p2,node p0)
{
double x1=p1.x-p0.x,y1=p1.y-p0.y;
double x2=p2.x-p0.x,y2=p2.y-p0.y;
return x1*y2-x2*y1;
}
void andrew()
{
memset(v,false,sizeof(v));
tp=0;sort(p+1,p+m+1,cmp);sta[++tp]=1;
for(int i=2;i<=m;i++)
{
while(tp>1&&mul(p[sta[tp]],p[i],p[sta[tp-1]])<=0)v[sta[tp--]]=0;
v[i]=1;sta[++tp]=i;
}
int top=tp;
for(int i=m-1;i;i--)
if(!v[i])
{
while(tp>top&&mul(p[sta[tp]],p[i],p[sta[tp-1]])<=0)--tp;
sta[++tp]=i;
}
for(int i=1;i<=tp;i++)h[i]=p[sta[i]];
}
int main()
{
int t=0;
while(scanf("%d",&n)!=EOF&&n)
{
for(int i=1;i<=n;i++)
scanf("%lf%lf%lf%lf",&a[i].x,&a[i].y,&a[i].v,&a[i].l);
double ans_val=inf,ans_left;int ans_num,ans_state;
for(int i=1;i<(1<<n)-1;i++)
{
m=0;
double val=0,len=0,left=0;
for(int j=1;j<=n;j++)
{
if(!(i>>(j-1)&1))
p[++m]=a[j];
else
{
val+=a[j].v;
len+=a[j].l;
}
}
if(val>ans_val)continue;
double used=0;
andrew();
for(int j=2;j<=tp;j++)used+=dis(h[j-1],h[j]);
left=len-used;
if(left>=0)
{
if(val<ans_val)
{
ans_state=i;
ans_num=n-m;
ans_val=val;
ans_left=left;
}
else if(val==ans_val&&n-m<ans_num)
{
ans_state=i;
ans_num=n-m;
ans_left=left;
}
}
}
printf("Forest %d\nCut these trees:",++t);
for(int i=0;i<n;i++)
if((1<<i)&ans_state)
printf(" %d",i+1);
printf("\nExtra wood: %.2f\n\n",ans_left);
}
return 0;
}