【解题报告】【HDOJ1392】【Graham凸包】Surround the Trees

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=1392

典型凸包

1.当只有一个点时 输出时0.00

2.当只有两个点时 输出时两个点的距离

3.当n>2 且共线时，输出需要*2,（这个是题目的奇怪之处，当只有2点时不用）

 

#include<stdio.h>
#include<string.h>
#include<math.h>
#define max(a,b) (a>b?b:a)

struct node{
    double x,y;
    int i;
};

struct node stack[101];
long int point1,point2;
double  x[101];
double  y[101];

void instack(int x,int y,int i);//入栈
void outstack();//出栈

int main()
{
    long int n,i,j,t1;
    double sum,a,b,x1,x2,y1,y2,t,p;
    //freopen("1.txt","r",stdin);
    
    while(scanf("%ld",&n)!=EOF&&n)
    {

        t1=0;//代表着y0
        for(i=0;i<n;i++)
        {
            scanf("%lf%lf",&x[i],&y[i]);
            if((y[i]<y[t1])||(x[i]>x[t1]&&y[i]==y[t1]))//寻找Y最小的点，若同，则寻X最大
                t1=i;
        }

        if(n==1) { printf("0.00\n"); continue;}//当只有一个点时
        if(n==2) //当只有两个点时
        {
            a=(x[1]-x[0])*(x[1]-x[0]);
            b=(y[1]-y[0])*(y[1]-y[0]);
            printf("%.2lf\n",sqrt(a+b));
            continue;
        }

        //将y0点放在x[0],y[0];
                p=x[t1];
                x[t1]=x[0];
                x[0]=p;
                p=y[t1];
                y[t1]=y[0];
                y[0]=p;    
        
        //将x[0],y[0]作为标准
        for(i=1;i<n;i++)
        {
            x[i]=x[i]-x[0];
            y[i]=y[i]-y[0];
        }
        x[0]=0;y[0]=0;
    
        //按角度从小到大排序
        for(i=1;i<n;i++)
        {
            for(j=i;j<n;j++)
                if(x[i]*y[j]-x[j]*y[i]<0)
                {
                    t=x[i];
                    x[i]=x[j];
                    x[j]=t;
                    t=y[i];
                    y[i]=y[j];
                    y[j]=t;
                }
        }

        //将前三点入栈
        stack[0].x=x[0];stack[0].y=y[0];stack[0].i=0;
        stack[1].x=x[1];stack[1].y=y[1];stack[1].i=1;
        stack[2].x=x[2];stack[2].y=y[2];stack[2].i=2;        
        point1=1;
        point2=2;
        //利用拐向来判断是否凸包内的
        for(i=3;i<n;)
        {
            
            x1=x[i]-stack[point1].x;
            y1=y[i]-stack[point1].y;
            
            x2=stack[point2].x-stack[point1].x;
            y2=stack[point2].y-stack[point1].y;
            
            if(x1*y2-x2*y1>=0)
            {
                outstack();
                continue;
            }
            else
            {
                instack(x[i],y[i],i);
            }
            i++;
        }
        
        //计算路径长度
        sum=0;
        for(i=1;i<=point2;i++)
        {
            
            a=(stack[i].x-stack[i-1].x)*(stack[i].x-stack[i-1].x);
            b=(stack[i].y-stack[i-1].y)*(stack[i].y-stack[i-1].y);
            sum+=sqrt(a+b);
        }

        //输出
        if(point2==1) printf("%.2lf\n",sum*2);//当n>2 且全部共线
        else//正常情况的输出
        {
            a=stack[point2].x*stack[point2].x;
            b=stack[point2].y*stack[point2].y;
            sum+=sqrt(a+b);
            printf("%.2lf\n",sum);
        }
        
    }
    return 0;
}

//***********下面为子函数***************
void instack(int x,int y,int i)//入栈
{
    point1++;
    point2++;
    stack[point2].x=x;
    stack[point2].y=y;
    stack[point2].i=i;
}

void outstack()//出栈
{
    point1--;
    point2--;
}