POJ1873

題意:給你n個樹,每個樹都有自己的座標和價值、長度。現在要砍掉一些樹,來把剩下來的樹給圍起來。

優先砍掉價值便宜的樹,如果有幾種砍掉樹方案切價值一樣則選擇砍掉數目最少的方案。然後價值和數量一樣隨便輸出一種就行了。

思路:由於n比較小,就可以用dfs枚舉出每一顆樹的狀態。在用剩下樹 使用凸包求出需要圍城的周長。

這道花了一個早上,煩死了,才學習了凸包。好多東西使用的時候出現了問題。比如求凸包中的while寫成了if,莫名其妙就

原數組給替換了,而且題目給的數據還非常的水,根本檢測不出來我的問題,盯着自己的代碼,不斷找錯,不過最後終於A了。


#include<stdio.h>
#include<math.h>
#include<string.h>
#include<algorithm>
using namespace std;
const int maxn=1<<16;
const int maxm=1<<29;
int n,e,k;
double sum1=0,sum2=0,sum3=0;
double maxn1=0,maxn2=0,maxn3=0;
int ans[maxn],kj[maxn];
struct node
{
    double x,y;
    double val,l;
} st[maxn],num[maxn],num2[maxn];
int cnt=0;
double dis(node p1,node p2)
{
    return (p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y);
}
double multi(node p1,node p2,node p3)
{
    return (p1.x-p3.x)*(p2.y-p3.y)-(p2.x-p3.x)*(p1.y-p3.y);
}
bool cmp(node p1,node p2)
{
    if(multi(p1,p2,num2[0])>0) return true;
    if(multi(p1,p2,num2[0])==0&&dis(p1,num2[0])<dis(p2,num2[0]))return true;
    return false;
}
void dfs(int cur)
{
    if(cur==n)
    {
        sum1=0,sum2=0,sum3=0;
        cnt=0;
        k=0;
        for(int i=0; i<n; i++)
        {
            if(ans[i]==0)//0代表沒有砍掉,保存在另一個數組中
            {
                num2[cnt++]=num[i];
            }
            else
            {
                sum1+=num[i].val;
                sum2+=num[i].l;
            }
        }
        sum3=n-cnt;
        if(cnt==0||cnt==n) return;
        if(cnt==1)//當只剩下來一個樹的時候,不用花費長度來圍它。下面可以和下面的合併,懶得合併了
        {
            if(maxn1>sum1)
            {
                maxn1=sum1;
                maxn2=sum2;
                maxn3=sum3;
                for(int i=0; i<n; i++)
                {
                    kj[i]=ans[i];
                }

            }
            else if(maxn1==sum1)
            {
                if(sum3<maxn3)
                {
                    maxn1=sum1;
                    maxn2=sum2;
                    maxn3=sum3;
                    for(int i=0; i<n; i++)
                        kj[i]=ans[i];
                }
            }
            return ;
        }
        double sum=0;
        if(cnt==2)//當剩下來兩顆樹的時候,要花費2倍的他們的距離長度。這要分開求,不然會有錯誤。
        {
            sum=dis(num2[0],num2[1]);
            sum=2*sqrt(sum);

        }
        else
        {
            for(int i=0; i<cnt; i++)
            {
                if(num2[i].x<num2[k].x||(num2[i].x==num2[k].x&&num2[i].y<num2[k].y)) k=i;
            }
            swap(num2[0],num2[k]);
            sort(num2+1,num2+cnt,cmp);
            st[0]=num2[0];
            st[1]=num2[1];
            st[2]=num2[2];
            e=2;
            for(int i=3; i<cnt; i++)
            {
                while(e>1&&multi(num2[i],st[e],st[e-1])>=0) e--;
                st[++e]=num2[i];
            }

            for(int i=0; i<=e; i++)
                sum+=sqrt((double)dis(st[i],st[i==e?0:i+1]));
        }
        if(sum<=sum2)
        {
            if(maxn1>sum1)
            {
                /*
                for(int i=0; i<n; i++)
                    printf("%d ",ans[i]);
                puts("");
                printf("%lf %lf %lf\n",maxn1,sum,sum2);
                */
                maxn1=sum1;
                maxn2=sum2-sum;
                maxn3=sum3;
                for(int i=0; i<n; i++)
                    kj[i]=ans[i];
            }
            else if(maxn1==sum1)
            {
                if(sum3<maxn3)
                {
                    maxn1=sum1;
                    maxn2=sum2-sum;
                    maxn3=sum3;
                    for(int i=0; i<n; i++)
                        kj[i]=ans[i];
                }
            }
        }
        return ;
    }
    ans[cur]=0;
    dfs(cur+1);
    ans[cur]=1;
    dfs(cur+1);
}
int main()
{
    int tt=1;
    while(scanf("%d",&n)!=EOF)
    {
        if(n==0) break;
        maxn1=maxn,maxn2=maxm;
        for(int i=0; i<n; i++)
            scanf("%lf %lf %lf %lf",&num[i].x,&num[i].y,&num[i].val,&num[i].l);
        dfs(0);
        printf("Forest %d\n",tt++);
        printf("Cut these trees:");
        for(int i=0; i<n; i++)
            if(kj[i])
                printf(" %d",i+1);
        puts("");
        printf("Extra wood: %.2lf\n\n",maxn2);
    }
    return 0;
}






發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章