凸包例题-UVa-10652

题意：有n块矩形木板，你的任务是用一个面积尽量小的凸多边形把它们包起来，并计算出木板占整个包装面积的百分比。

                                                                                                            

第一行输入的是数据的组数，（T<50)。每组数据第一行为木板个数，（N<600)。以下n行每行五个实数，x,y,w,h,j。j是顺时针旋转的角度。

Sample Input

1

4

4 7.5 6 3 0

8 11.5 6 3 0

9.5 6 6 3 90

4.5 3 4.4721 2.2361 26.565

Sample Output

64.3 %

思路：这就是一个简单的凸包问题，输入的所有点为凸包中的点，求包含这些点的凸包。

 

代码：

#include <iostream>
#include <algorithm>
#include <cmath>
#include <cstring>
#include <cstdio>

using namespace std;
const int N=2500;
#define eps 1e-8

//点的定义 
struct Point{
	double x,y;
	Point(double x=0,double y=0):x(x),y(y){};
	
	void input() {
		scanf("%lf%lf",&x,&y);
	}
}; 

//向量 
typedef Point Vector;


//判断相等 
int dcmp(double x)
{
	if(x<-eps) return -1;
	if(x>eps) return 1;
	return 0;
}

template <class T> T sqr(T x){
	return x*x;
}

//向量的运算 
Vector operator + (Vector A,Vector B)
{
	return Vector(A.x+B.x,A.y+B.y);
}

Vector operator -(Vector A,Vector B)
{
	return Vector(A.x-B.x,A.y-B.y);
}

Vector operator *(Vector A,double p)
{
	return Vector(A.x*p,A.y*p);
}

Vector operator /(Vector A,double p)
{
	return Vector(A.x/p,A.y/p);
}
double Dot(Vector A,Vector B)
{
	return A.x*B.x+A.y*B.y;
 } 
 
 bool operator < (const Point &a ,const Point &b)
{
	return a.x<b.x || (a.x==b.x && a.y<b.y);
}

bool operator >=(const Point &a,const Point &b)
{
	return a.x>=b.x && a.y>=b.y;
}

bool operator <= (const Point &a,const Point &b)
{
	return a.x<=b.x && a.y<=b.y;
}

bool operator == (const Point &a,const Point &b)
{
	return dcmp(a.x-b.x)==0 && dcmp(a.y-b.y==0);
}

//向量的长度 
double Length(Vector A)
{
	return sqrt(Dot(A,A));
}

//求俩向量的夹角 
double Angle(Vector A,Vector B)
{
	return acos(Dot(A,B))/Length(A)/Length(B);
}

//求差积 
/* 
在二维空间中，叉乘还有另外一个几何意义就是：
aXb等于由向量a和向量b构成的平行四边形的面积。
*/
double Cross(Vector A,Vector B)
{
	return A.x*B.y-A.y*B.x;
}

//角度转为弧度
double torad(double deg)
{
   	return  deg*acos(-1)/180;
}

//求凸包 
int ConvexHull(Point *p,int n,Point *ch)
{
	sort(p,p+n);
	int m=0;
	for(int i=0;i<n;i++)
	{
		while(m>1 && Cross(ch[m-1]-ch[m-2],p[i]-ch[m-2])<=0) m--;
		ch[m++]=p[i];
	}
	
	int k=m;
	for(int i=n-2;i>=0;i--)
	{
		while(m>k && Cross(ch[m-1]-ch[m-2],p[i]-ch[m-2])<=0) m--;
		ch[m++]=p[i];
	}
	
	if(n>1) m--;
	return m;
	
}

//向量旋转 
Vector Rotate(Vector A,double rad)
{
	return Vector(A.x*cos(rad)-A.y*sin(rad),A.x*sin(rad)+A.y*cos(rad));
}

//求多边形面积 
double PolygonArea(Point *p,int n)
{
	double area=0;
	for(int i=1;i<n-1;i++)
		area+=Cross(p[i]-p[0],p[i+1]-p[0]);
		
	return area/2;
 } 

int main()
{
	Point P[N],ch[N]; 
	int t;
	scanf("%d",&t);
	while(t--)
	{
		int n,pc=0;
		scanf("%d",&n);
		double area1=0;
		for(int i=0;i<n;i++)
		{
			double x,y,w,h,j;
			scanf("%lf%lf%lf%lf%lf",&x,&y,&w,&h,&j);
			Point o(x,y);
			double ang=torad(j);
			P[pc++]=o+Rotate(Vector(-w/2,-h/2),ang);
			P[pc++]=o+Rotate(Vector(w/2,-h/2),ang);
			P[pc++]=o+Rotate(Vector(-w/2,h/2),ang);
			P[pc++]=o+Rotate(Vector(w/2,h/2),ang);
			area1+=w*h;	
		}
		int m=ConvexHull(P,pc,ch);
		double area2=PolygonArea(ch,m);
		
		printf("%.1lf %%\n",area1*100/area2);
	}
	
	return 0;
}