POJ 1873 The Fortified Forest（凸包）题解

题意：二维平面有一堆点，每个点有价值v和删掉这个点能得到的长度l，问你删掉最少的价值能把剩余点围起来，价值一样求删掉的点最少

思路：n<=15，那么直接遍历2^15，判断每种情况。这里要优化一下，如果价值比当前最优大了continue。POJ的G++输出要用%f...orz，还是乖乖用C++...

代码：

#include<set>
#include<map>
#include<stack>
#include<cmath>
#include<queue>
#include<vector>
#include<string>
#include<cstdio>
#include<cstring>
#include<sstream>
#include<iostream>
#include<algorithm>
typedef long long ll;
using namespace std;
const int maxn =  + ;
const int MOD = 1e9 + ;
const int INF = 0x3f3f3f3f;
struct node{
    double x, y, l;
    int v, id;
}p[maxn], s[maxn], q[maxn];
int n, top;
double dis(node a, node b){
    return sqrt((a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y));
}
bool cmp(node a, node b){
    double A = atan2((a.y - p[].y), (a.x - p[].x));
    double B = atan2((b.y - p[].y), (b.x - p[].x));
    if(A != B) return A < B;
    else{
        return dis(a, p[]) < dis(b, p[]);
    }
}
double cross(node a, node b, node c){   //(a->b)X(a->c)
    return (b.x - a.x) * (c.y - a.y) - (c.x - a.x) * (b.y - a.y);
}
void solve(){
    int pos = ;
    for(int i = ; i <= n; i++){
        if(p[i].y < p[pos].y || (p[i].y == p[pos].y && p[i].x < p[pos].x)){
            pos = i;
        }
    }
    swap(p[], p[pos]);
    sort(p + , p + n + , cmp);
    s[] = p[], s[] = p[];
    top = ;
    for(int i = ; i <= n; i++){
        while(top >=  && cross(s[top - ], p[i], s[top]) >= ){
            top--;
        }
        s[++top] = p[i];
    }
}
double need(double len){
    if(n <= ) return len;
    if(n == ){
        return len - 2.0 * dis(p[], p[]);
    }
    solve();
    double Need = ;
    for(int i = ; i < top; i++){
        Need += dis(s[i], s[i + ]);
    }
    Need += dis(s[top], s[]);
    return len - Need;
}
int main(){
    int N, ca = ;
    while(~scanf("%d", &N) && N){
        for(int i = ; i < N; i++){
            scanf("%lf%lf%d%lf", &q[i].x, &q[i].y, &q[i].v, &q[i].l);
            q[i].id = i + ;
        }
        int sz = , que[], tempQue[];
        int MinValue = INF;
        double ans = ;
        for(int i = ; i < ( << N); i++){
            int num = , value = ;
            double len = ;
            n = ;
            for(int j = ; j < N; j++){
                if(i & ( << j)){
                    p[++n] = q[j];
                }
                else{
                    value += q[j].v;
                    len += q[j].l;
                    tempQue[num] = q[j].id;
                    num++;
                }
            }
            if(value > MinValue) continue;
            double tmp = need(len);
            if(tmp < ) continue;
            if(value < MinValue || (value == MinValue && num < sz)){
                MinValue = value;
                ans = tmp;
                sz = num;
                for(int j = ; j < num; j++){
                    que[j] = tempQue[j];
                }
            }
        }
        if(ca != ) printf("\n");
        printf("Forest %d\n", ca++);
        printf("Cut these trees: ");
        for(int i = ; i < sz; i++){
            printf("%d ", que[i]);
        }
        printf("\n");
        printf("Extra wood: %.2lf\n", ans);
    }
    return ;
}