【UVA10652】Board Wrapping(凸包+坐标旋转+多边形面积)
提交地址:https://vjudge.net/problem/UVA-10652
题目:
pdf:https://uva.onlinejudge.org/external/106/10652.pdf
解题思路:
凸包模版题,求完凸包再求这个凸包(多边形)的面积。
有两个需要注意的小点:
(1)转角的范围: 且这个木板是以木板中心(x,y)为旋转中心旋转的,若要求出旋转后的四个点坐标,坐标中心应该是(x,y),不能先求出平放时的四点坐标再旋转,因为那样旋转中心是(x,y)对应的(0,0)。
(2)Rotate(Vecotr p, double rad)中的参数rad是逆时针旋转时的弧度数,要将输入的度数转化一下。
ac代码:
#include <bits/stdc++.h>
using namespace std;
const int maxn = 4*610;
const double pi = acos(-1.0);
const double eps = 1e-10;
struct Point
{
double x,y;
Point(double x=0,double y=0):x(x),y(y){}
friend bool operator < (Point a, Point b)
{
return a.x == b.x ? a.y < b.y : a.x < b.x;
}
};
typedef Point Vector;
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);
}
int dcmp(double x)//精度三态函数(>0,<0,=0)
{
if (fabs(x) < eps)return 0; //等于
else return x < 0 ? -1 : 1;//小于,大于
}
double Cross(Vector a, Vector b)//外积
{
return a.x * b.y - a.y * b.x;
}
Vector Rotate(Vector a, double rad)//逆时针旋转rad弧度
{
return Vector(a.x*cos(rad) - a.y*sin(rad), a.x*sin(rad) + a.y*cos(rad));
}
int ConvexHull(Point *p, int n, Point* ch)//凸包
{
sort(p, p+n);//先x后y
int m = 0;//ch存最终的凸包顶点,下标从0开始
for(int i = 0; i < n; i++)
{
//m是待确定的
while(m > 1 && dcmp(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 && dcmp(Cross(ch[m-1]-ch[m-2], p[i]-ch[m-2])) <= 0) m--;//同样的方向
ch[m++] = p[i];
}
if(n>1) m--;//m是凸包上的点数+1,相当于n,ch[m]和ch[0]相同,遍历的时候<m即可
return m;
}
double PolygonArea(Point* p, int n)//多边形面积
{
double area = 0;
for(int i = 1; i < n - 1; i++)//下标从0开始存坐标,共n个点
area += Cross(p[i] - p[0], p[i + 1] - p[0]);
return area/2.0;
}
int main()
{
//freopen("/Users/zhangkanqi/Desktop/11.txt","r",stdin);
int t, n;
double x, y, w, h, ang;
Point allp[maxn], ch[maxn];//p存所有点
scanf("%d", &t);
while(t--)
{
double barea = 0;
int cnt = 0;
scanf("%d", &n);
for(int i = 0; i < n; i++)
{
scanf("%lf %lf %lf %lf %lf", &x, &y, &w, &h, &ang);
double hw = w/2.0, hh = h/2.0;
if(ang >= 0) ang = 2*pi-(pi*ang)/180.0;
else ang = -pi*ang/180.0;//本身就是逆时针的度数
Point c = Point(x, y);
//绕着中心旋转,中心c不动,而不是水平时的四个坐标旋转(此时的旋转中心是(0,0),中心的坐标也会变)
allp[cnt++] = c + Rotate(Point(-hw, -hh), ang);
allp[cnt++] = c + Rotate(Point(hw, -hh), ang);
allp[cnt++] = c + Rotate(Point(-hw, hh), ang);
allp[cnt++] = c + Rotate(Point(hw, hh), ang);
barea += w*h;
}
int m = ConvexHull(allp, cnt, ch);
double tarea = PolygonArea(ch, m);
printf("%.1lf %%\n", barea*100/tarea);
}
return 0;
}