UVa10652

题目链接

简介：凸包

分析：简介就已经邪泄露了正解

tip

这道题求坐标是一个难点
有一个很好的方法
我们知道点+向量=点
所以我们可以先算出中点到达每块木板的四个点的方向向量，点的坐标直接用中点坐标加上方向向量即可

a=-torad(j);                                  //变成弧度制 node o=node(x,y);                             //顺时针转动，点+向量=点 po[++cnt]=o+rotate(node(-w/2,-h/2),a);                 //左下 po[++cnt]=o+rotate(node(-w/2,h/2),a);                  //左上 po[++cnt]=o+rotate(node(w/2,h/2),a);                   //右上 po[++cnt]=o+rotate(node(w/2,-h/2),a);                  //右下 area1+=w*h;

在计算面积的时候，一定要把坐标写准了

double SS(int n){    double ans=0;    for (int i=2;i<=n;i++)        ans+=Cross(po[sta[i]]-po[sta[1]],po[sta[i+1]]-po[sta[1]]);    return fabs(ans)/2;}

//这里写代码片#include<cstdio>#include<cstring>#include<iostream>#include<algorithm>#include<cmath>using namespace std;const double Pi=acos(-1.0);const double eps=1e-8;const int N=3000;int n,sta[N];struct node{    double x,y;    node (double xx=0,double yy=0)    {        x=xx;y=yy;    }};node po[N];node operator + (const node &a,const node &b){return node(a.x+b.x,a.y+b.y);}node operator - (const node &a,const node &b){return node(a.x-b.x,a.y-b.y);}node operator * (const node &a,const double &b){return node(a.x*b,a.y*b);}node operator / (const node &a,const double &b){return node(a.x/b,a.y/b);}bool operator < (const node &a,const node &b){return a.x<b.x||(a.x==b.x&&a.y<b.y);}bool operator == (const node &a,const node &b){return a.x==b.x&&a.y==b.y;}int dcmp(double x){    if (fabs(x)<eps) return 0;    else if (x>0) return 1;    else return -1;}double torad(double n){    return n*Pi/180;}double Cross(node x,node y){return x.x*y.y-x.y*y.x;}double Dot(node x,node y){return x.x*y.x+x.y*y.y;}node rotate(node x,double a){    return node(x.x*cos(a)-x.y*sin(a),x.x*sin(a)+x.y*cos(a));}int TuB(int n){    int top=0;    sort(po+1,po+1+n);    for (int i=1;i<=n;i++)    {        while (top>1&&dcmp( Cross( po[sta[top]]-po[sta[top-1]] , po[i]-po[sta[top-1]] ) )<=0) top--;     //顺时针转动不合法         sta[++top]=i;    }    int k=top;    for (int i=n-1;i>=1;i--)    {        while (top>k&&dcmp( Cross( po[sta[top]]-po[sta[top-1]] , po[i]-po[sta[top-1]] ))<=0) top--;        sta[++top]=i;    }    if (n>1) top--;    return top;}double SS(int n){    double ans=0;    for (int i=2;i<=n;i++)        ans+=Cross(po[sta[i]]-po[sta[1]],po[sta[i+1]]-po[sta[1]]);    return fabs(ans)/2;}int main(){    int T;    scanf("%d",&T);    while (T--)    {        double x,y,w,h,j,a;        double area1=0;        scanf("%d",&n);        int cnt=0;        for (int i=1;i<=n;i++)        {            scanf("%lf%lf%lf%lf%lf",&x,&y,&w,&h,&j);            a=-torad(j);                                  //变成弧度制             node o=node(x,y);                             //顺时针转动，点+向量=点             po[++cnt]=o+rotate(node(-w/2,-h/2),a);                 //左下             po[++cnt]=o+rotate(node(-w/2,h/2),a);                  //左上             po[++cnt]=o+rotate(node(w/2,h/2),a);                   //右上             po[++cnt]=o+rotate(node(w/2,-h/2),a);                  //右下             area1+=w*h;        }             int m=TuB(cnt);        double area2=SS(m);        printf("%.1lf %%\n",area1*100/area2);    }}