下面是我做这道题时的一些分析 大部分是错的 大家参考参考 我也是写下自己错的好好反思一下以后不犯这种错误
题意:这道题找的是花费最小的价值把树木围起来 花多少树木都没事 但是最后至少能剩下一颗 剩下两颗也不一定能用前面的所有树围起来 所以必须枚举到只剩一颗的情况
然后问题就来了:
1:怎么枚举所有情况?
2:怎么确定花费的最小价值?
3:枚举到只剩一棵树的情况 则找出价值最大的那一颗数就行了
4:枚举到只剩两颗棵树的时候也不用求凸包 计算剩余点的距离(两倍)和价值 找到一个价值最大的输出
5:其他情况 计算凸包的周长 和凸包包围树木的价值(标记数组看看树木是不是在凸包中 枚举的时候赋值 每次枚举都要初始化) 找到能满足篱笆条件的最大的一个组合且价值最大的情况 保存当前最佳状态 (思考一下能不能不用每次都计算凸包的周长和价值 找一找有什么剪枝优化没 )
6:把价值从小到大排 价值相等高度从小到大排
这样排序后从头到尾枚举 每种情况只要找到满足条件就可以保证这是当前枚举情况的最优解(比如枚举的的是砍掉3颗树的情况
只要有一个满足就可以输出了)
枚举N个当中选一个的反向就是N个中选8个的选法 所以枚举操作可以大大简化
最后做的时候意识到既然我都枚举了所有情况那就都算一遍得了 反正数据量小也不会超时 而且不用判断特殊情况 因为所有情况在计算凸包的时候都计算了
刚开始做的时候自己写递归求组合弄了一上午没弄出来 主要是应该返回到什么位置没写出来要不都用自己写的了 后来就去看了别人的 代码用的是别人写的递归求的组合
要是只剩两颗树就是两棵树的距离的2倍 一棵树的情况也一样 没有漏判的情况
题意:这道题找的是花费最小的价值把树木围起来 花多少树木都没事 但是最后至少能剩下一颗 剩下两颗也不一定能用前面的所有树围起来 所以必须枚举到只剩一颗的情况
然后问题就来了:
1:怎么枚举所有情况?
2:怎么确定花费的最小价值?
3:枚举到只剩一棵树的情况 则找出价值最大的那一颗数就行了
4:枚举到只剩两颗棵树的时候也不用求凸包 计算剩余点的距离(两倍)和价值 找到一个价值最大的输出
5:其他情况 计算凸包的周长 和凸包包围树木的价值(标记数组看看树木是不是在凸包中 枚举的时候赋值 每次枚举都要初始化) 找到能满足篱笆条件的最大的一个组合且价值最大的情况 保存当前最佳状态 (思考一下能不能不用每次都计算凸包的周长和价值 找一找有什么剪枝优化没 )
6:把价值从小到大排 价值相等高度从小到大排
这样排序后从头到尾枚举 每种情况只要找到满足条件就可以保证这是当前枚举情况的最优解(比如枚举的的是砍掉3颗树的情况
只要有一个满足就可以输出了)
枚举N个当中选一个的反向就是N个中选8个的选法 所以枚举操作可以大大简化
最后做的时候意识到既然我都枚举了所有情况那就都算一遍得了 反正数据量小也不会超时 而且不用判断特殊情况 因为所有情况在计算凸包的时候都计算了
刚开始做的时候自己写递归求组合弄了一上午没弄出来 主要是应该返回到什么位置没写出来要不都用自己写的了 后来就去看了别人的 代码用的是别人写的递归求的组合
要是只剩两颗树就是两棵树的距离的2倍 一棵树的情况也一样 没有漏判的情况
下面看看代码 有注释
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define eps 1e-8
#define maxn 1000000
using namespace std;
struct Tree
{
int x,y,v,l,n,ok;
}t[20],num[20],s[20];
struct H //保存最优解
{
int v,number[16],len;
double l;
}h;
int T,e=1;
typedef struct Tree tree;
int cmp(tree a,tree b)//排序参数 求凸包用
{
if(a.x<b.x)return 1;
else if(a.x==b.x&&a.y<b.y)return 1;
return 0;
}
double Dis(tree a, tree b)//计算距离
{
return sqrt((double)((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y)));
}
//极角比较, 返回-1: p0p1在p0p2的右侧,返回0:p0,p1,p2共线
int Cross(tree p1, tree p2, tree pb)
{
int delta=(p1.x-pb.x)*(p2.y-pb.y)-(p2.x-pb.x)*(p1.y-pb.y);
if (delta<0) return 1;
else if (delta==0) return 0;
else return -1;
}
int ConvexHull(tree *p,int n,tree *ch)//Andrew法寻找凸包 必须先排序
{
int m = 0;
for(int i = 0; i < n; i++)//寻找下凸包
{
if(p[i].ok)continue;
while(m > 1 && Cross(ch[m-1],p[i],ch[m-2]) <= 0 )m--;
ch[m++] = p[i];
}
int k = m;
for(int i = n-2; i >= 0; i--)//寻找上凸包
{
if(p[i].ok)continue;
while(m > k && Cross(ch[m-1],p[i],ch[m-2]) <= 0 )m--;
ch[m++] = p[i];
}
if(n > 1)m--;
return m;
}
void Solve(int m,int v,double l)
{
double sum = 0;
for(int i = 0; i < m; i++)//计算凸包的周长
sum +=Dis(s[i],s[i+1]);
if(v<h.v&&l-sum>0||fabs(l-sum)<eps)//价值不比原来的小就退出当前组合继续计算下一种组合
{
h.len = 0;
h.v = v;//保存花费的最小价值
h.l = l-sum;//保存剩余的树木长度
for(int i = 0; i < T; i++)//更新砍掉的树木,最后排序从小到大输出
if(t[i].ok)h.number[h.len++]=t[i].n;
}
for(int i = 0; i < T; i++)//这一定要初始化
t[i].ok = 0;
}
void combine(int n, int m,const int M )
{
for(int i=n; i>=m; i--) // 注意这里的循环范围
{
num[m-1].n = i-1;
if (m > 1)
combine(i-1,m-1,M);
else // m == 1, 输出一个组合
{
int v = 0;
double l = 0;
for(int j=M-1; j>=0; j--)
{
t[num[j].n].ok=1;//标记此树已被砍掉
v +=t[num[j].n].v;//计算当前组合砍掉的树的价值;
l +=(double)t[num[j].n].l;//计算当前组合砍掉的树的长度
}
//其实这里num这个定义成数组就行了 因为就只是递归记录下标而已
int m = ConvexHull(t,T,s);//找出凸包
Solve(m,v,l);
}
}
}
int main()
{
#ifdef LOCAL
freopen("in.txt","r",stdin);
freopen("o.txt","w",stdin);
#endif // LOCAL
while((scanf("%d",&T),T)!=0)
{
for(int i = 0; i < T; i++)
{
scanf("%d%d%d%d",&t[i].x,&t[i].y,&t[i].v,&t[i].l);
t[i].n = i+1;//标号
}
sort(t,t+T,cmp);//排序求凸包的必要步骤
h.v = maxn;
for(int i = 1; i < T; i++)
combine(T,i,i);
/*枚举砍掉的树的所有组合找出凸包*/
sort(h.number,h.number+h.len);
printf("Forest %d\nCut these trees:",e++);
for(int i = 0; i < h.len; i++)
printf(" %d",h.number[i]);
printf("\nExtra wood: %.2f\n\n",h.l);
}
return 0;
}