poj 1113 Wall

题意：给你城堡的坐标，建造城堡外面的墙，城堡到墙有一定的距离，该距离给定，求墙的最小面积

分析：凸包求最外围的点，周长和就是凸包周长+2*PI*L(运用外角和等于360)

#include<stdio.h>
#include<algorithm>
#include<math.h>
#include<iostream>
using namespace std;
const double PI=(acos(-1.0));
const int MN=1500;

struct POINT 
{
    int x,y;
    int flag;
};
POINT list[MN];
int top,stack[MN];

void swap(POINT &a,POINT &b)//一定要传地址，不然形参改，实参不改
{
	POINT t;
	t=a;a=b;b=t;  
}

int CrossProd(POINT p0,POINT p1,POINT p2)
{
	return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y); 
}

int dist(POINT p1,POINT p2)
{
    return (p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y);
}

bool cmp(POINT p1,POINT p2)  
{  
    return CrossProd(list[0],p1,p2)>0;   
}  

int select(int n)
{
    int i,j,num;
    for (i=1;i<n;i++)
    {
        if(list[i].flag)
        {
            for (j=i+1;j<n;j++)
            {
                if(list[j].flag)
                {//判断是否共线
                    int t=CrossProd(list[0],list[i],list[j]);
                    if(t==0)
                    {
                        if(dist(list[0],list[i])>dist(list[0],list[j]))
							list[j].flag=0;
                        else list[i].flag=0; 
                    }
                }
            }
        }
    }
    int cas=1;
    for (j=1;j<n;j++)
    {
        if(list[j].flag)
			list[cas++]=list[j];
    }
    num=cas-1;
    sort(list+1,list+num+1,cmp);//sort前闭后开
    return num;
}

void graham(int n,int L)
{
    double sum=0;
    int i;
    for (i=0;i<=2;i++) stack[i]=i;
    top=2;
    for(i=3;i<=n;i++)
    {
        while(CrossProd(list[stack[top-1]],list[stack[top]],list[i])<0) top--;//
        top++;
        stack[top]=i;
    }
	for(i=0;i<=top;i++)
	{
		sum+=sqrt((double)dist(list[stack[i]],list[stack[(i+1)%(top+1)]]));
	}
	printf("%.0lf\n",sum+2*PI*L);// 
	
}

int main()
{
    int N,L,num,i;
    while(scanf("%d%d",&N,&L)!=EOF)
    {
		for (i=0;i<N;i++)
		{
			scanf("%d%d",&list[i].x,&list[i].y);
			list[i].flag=1;
			//找到y最小x最小的点
			if(list[i].y<list[0].y || (list[i].x<list[0].x && list[i].y==list[0].y))
				swap(list[i],list[0]);         
		}
		//排除同一直线的点
		num=select(N);
		graham(num,L);  
	}
	return 0;
}