POJ 1873 The Fortified Forest 计算凸包

题目描述：http://poj.org/problem?id=1873

题目大意：

一片森林中有n棵树，现在需要砍掉其中的一些树做成围栏把剩下的树围起来。给定每棵树的(x,y)坐标、价值和砍掉之后能做成围栏的长度。找到最优的一种可行方案使得砍掉的树的价值最小。如果两种方案砍掉树的价值相同，则选择砍掉数目数量少的哪一个方案。

多组数据，2 <= n <= 15。

解题思路：

看到n的范围马上想到枚举。枚举每棵树砍掉或者不砍，算出砍掉的树能做成围栏的长度length 和 剩下的树的凸包长度fence，如果length>=fence那么可行，然后在所有可行解中取最优解。

知识点：

枚举、计算凸包。

总结：

思路很简单，但是WA了几次。原因是计算凸包写错了。

我用的是对xy坐标排序的graham-scan算法，在排序时只对x坐标排序，对y不排序，然后在构造过程中判断是否左转的时候把共线也判断成左转。这样写的话计算出的凸包没有问题，凸包的边上三点共线时会把中间的点踢掉。

还有一种情况是不把共线判断为左转，这样的话计算出的凸包就会有共线的点。但是要求排序时对xy双关键字排序，不然会出错。

所以以后还是老老实实写双关键字的排序吧。如果需要构造有三点共线的凸包，那么判断左转时不加=0。如果构造含无三点共线的凸包，判断左转时加上=0。

/*
 *
 *	Problem: 1873		
 *	Memory: 192K		Time: 547MS
 *	Language: C++		Result: Accepted
 *
 */

#include "stdio.h"
#include "math.h"
#define MAXN 17
#define EPS 1e-8
#define INFV 99999999

struct Point {
	double x, y;
}p[MAXN];

struct Tree {
	struct Point p;
	int v;
	int l;
}t[MAXN];


bool cut[MAXN];
int minValue, top, n;
double extra;
int ansList[MAXN], ansn;

void sort(int l, int r) {
	int i=l, j=r;
	double x = p[(l+r)>>1].x, y = p[(l+r)>>1].y;
	do {	
		while (p[i].x<x || (p[i].x==x && p[i].y<y)) i++;
		while (p[j].x>x || (p[j].x==x && p[j].y>y)) j--;
		if (i<=j) {
			Point t = p[i];
			p[i] = p[j];
			p[j] = t;
			i++; j--;
		}
	} while (i<=j);
	if (i<r) sort(i, r);
	if (l<j) sort(l, j);
}

void init() {
	minValue = INFV;
	ansn     = n+1;
	for (int i=1;i<=n;i++) {
		scanf("%lf %lf %d %d", &t[i].p.x, &t[i].p.y, &t[i].v, &t[i].l);
	}
}

bool turnLeft(struct Point a, struct Point b, struct Point c) {
	double x1 = b.x-a.x, y1 = b.y-a.y;
	double x2 = c.x-a.x, y2 = c.y-a.y;
	return x1*y2-x2*y1>=0;   
}

double dis(struct Point a, struct Point b) {
	return sqrt(pow(a.x-b.x, 2) + pow(a.y-b.y, 2));
}

double convexHull() {
	double ans = 0;
	struct Point q[MAXN];
	int    n1 = 0;

	sort(1, top);
	
	q[1] = p[1]; n1 = 1;
	for (int i=2; i<=top; i++) {
		while (n1>=2 && turnLeft(q[n1-1], q[n1], p[i])) n1--;
		q[++n1] = p[i];
	}
	for (int i=1; i<n1; i++) ans += dis(q[i], q[i+1]);

	q[1] = p[top]; n1 = 1;
	for (int i=top-1;i>=1;i--) {
		while (n1>=2 && turnLeft(q[n1-1], q[n1], p[i])) n1--;
		q[++n1] = p[i];
	}
	for (int i=1; i<n1; i++) ans += dis(q[i], q[i+1]);

	return ans;
}

void dfs(int dp) {
	if (dp == n+1) {
		double length = 0;
		int cnt   = 0; 
		int value = 0;
		top = 0;
		for (int i=1;i<=n;i++) 
		if (cut[i]) {
			length  += t[i].l;
			value   += t[i].v;
			cnt     ++;
		} else {
			p[++top] = t[i].p;
		}
		double fence = convexHull(); // 剩下的树构成的凸包长度
		if (fence-length<=EPS && (value<minValue || (value==minValue && cnt<ansn))) {
			minValue = value;
			extra    = length - fence;
			ansn     = 0;
			for (int i=1;i<=n;i++) 
			if (cut[i]) ansList[++ansn] = i;
		}

	} else {
		cut[dp] = true;  dfs(dp+1);
		cut[dp] = false; dfs(dp+1);
	}
}

int main() {
	
	scanf("%d", &n);

	int t = 0;

	while (n) {

		t ++;

		init();

		// 习惯性写深搜了，也可以按位枚举
		dfs(1);

		printf("Forest %d\n", t);
		printf("Cut these trees: ");
		for (int i=1;i<=ansn;i++) printf("%d ", ansList[i]);
		printf("\n");
		printf("Extra wood: %.2lf\n", extra);
		printf("\n");

		scanf("%d", &n);

	}
	return 0;
}