The Fortified Forest[POJ1873]

可爱的(爆炸的)caioj.cn

题面描述

传送门

思路

简要题意，给定一些树，砍掉一些有价值的树后，得到一些木材，这些木材的长度，要不小于剩下的树围成的凸包的周长，现在要求满足最小价值时，最小砍伐树木的方案。

把状态压缩为二进制

$val$表示费用，$state$表示压缩的状态，$num$表示砍掉的树的数目，$len$为砍伐后得到的木材长度，$used$表示凸包周长。

每一次更新一下就好了。

有一个剪枝，若当前状态价值大于$an{s}_{val}$，就直接跳过。

#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;
}