UVA10652【凸包计算】

https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=1593

给一些矩形（给定中心（x,y）长w宽h,和旋转角j）

求一个最小多边形把所有的矩形包围，输出矩形所占多边形的面积比。

可以想到是把矩形拆成4个点然后跑凸包。

由于第一次做凸包不是很会，然后看了模板YY一下，就是sort之后每次当前点在凸包最末的那一条边的左边，才能加进来，否则弹掉最后一个点，直到可以加入。

这里还有矩形是旋转的，所以搞点的时候卡了一阵子。

可以看成是每个端点先相对于中心点作原点的向量旋转theta角度之后，把原点平移。

(然而我以为三角函数传入的是360°的角度，没想到竟然是弧度，然后就WA了- -||)

然后处理一下点，就当成裸的凸包吧：

.

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<queue>
#include<cmath>
#define eps 1e-10
const int maxn=620;
const double pi=3.1415926;
using namespace std;
int dcmp(double x)
{
	if(fabs(x)<eps)return 0;
	else return x<0?-1:1;
}
struct point
{
	double x,y;
	point(){}
	point(double _x,double _y)
	{
		x=_x;
		y=_y;
	}
	point operator +(const point &b)
	{
		return point(x+b.x,y+b.y);
	}
	point operator -(const point &b)
	{
		return point(x-b.x,y-b.y);
	}
	point operator *(const double &b)
	{
		return point(x*b,y*b);
	}
}p[4*maxn];
double cross(point a,point b)
{
	return a.x*b.y-a.y*b.x;
}
double dot(point a,point b)
{
	return a.x*b.x+a.y*b.y;
}
int n;
int cnt;
point rotate(point a,double rad)
{
	//cout<<a.x*cos(rad)<<" "<<a.y*sin(rad)<<" "<<endl;
	return point(a.x*cos(rad)+a.y*sin(rad),a.y*cos(rad)-a.x*sin(rad));
}
bool cmp(point x,point y)
{
	if(x.x!=y.x)return x.x<y.x;
	else return x.y<y.y; 
}
point q[4*maxn];
double ans1,ans2;
void solve()
{
	sort(p+1,p+cnt+1,cmp);
	int m=0;
	for(int i=1;i<=cnt;i++)
	{
		while(m>1&&dcmp(cross(q[m-1]-q[m-2],p[i]-q[m-2]))<=0)m--;
		q[m++]=p[i];
	}
	int k=m;
	for(int i=cnt-1;i>=1;i--)
	{
		while(m>k&&dcmp(cross(q[m-1]-q[m-2],p[i]-q[m-2]))<=0)m--;
		q[m++]=p[i];
	}
	if(cnt>1)m--;
	
	for(int i=0;i<m-1;i++)ans2+=cross(q[i]-q[0],q[i+1]-q[0]);
	ans2=0.5*fabs(ans2);
	printf("%.1f %\n",ans1/ans2);
}
int main()
{
	int T;
	scanf("%d",&T);
	while(T--)
	{
		cnt=0;
		ans1=ans2=0;
		scanf("%d",&n);
		for(int i=1;i<=n;i++)
		{
			double x,y,w,h,j;
			scanf("%lf%lf%lf%lf%lf",&x,&y,&w,&h,&j);
			ans1+=w*h*100;
			j=j/180*pi;
			//cout<<rotate(point(-w/2,h/2),j).x<<" "<<rotate(point(-w/2,h/2),j).y<<endl;
			p[++cnt]=rotate(point(-w/2,h/2),j)+point(x,y);
			p[++cnt]=rotate(point(w/2,h/2),j)+point(x,y);
			p[++cnt]=rotate(point(-w/2,-h/2),j)+point(x,y);
			p[++cnt]=rotate(point(w/2,-h/2),j)+point(x,y);
		}
		solve();
	}
	return 0;
}