码迷,mamicode.com
首页 > 移动开发 > 详细

【UVA10652】Board Wrapping(求凸包面积)

时间:2018-12-10 19:37:34      阅读:214      评论:0      收藏:0      [点我收藏+]

标签:大致   tps   name   告诉   pre   mes   ping   space   put   

点此看题面

大致题意: 告诉你若干个矩形的重心坐标相对\(y\)轴的偏转角度,求矩形面积和能围住这些矩形的最小凸包面积之比。


矩形面积和

这应该是比较好求的吧。

已经给了你长和宽,直接乘起来累加即可。


最小凸包面积

这道题关键还是在于求凸包面积。

首先,我们要注意将题目中给出的角度转换成弧度,还要记得取相反数,不然调死你。

这段代码可以与旋转函数放在一起:

inline Point Rotate(Vector A,double deg)
{
    register double rad=deg*acos(-1)/180;dcmp(rad)<0&&(rad+=2*acos(-1)),rad=-rad;//将角度转换成弧度,并取相反数
    return Vector(A.x*cos(rad)-A.y*sin(rad),A.x*sin(rad)+A.y*cos(rad));//返回旋转后的点的坐标
} 

则我们就可以很方便地求出每个矩形四个顶点的坐标了。

枚举四个方向,分别将对应的向量(\(±\frac{s1}2\)\(±\frac{s2}2\))进行旋转,然后加上原先点的坐标即可(\(s1,s2\)分别为长和宽)。

这段代码如下:

A=Point(x,y),
P.p[++P.n]=A+Rotate(Point(s1/2,s2/2),deg),
P.p[++P.n]=A+Rotate(Point(s1/2,-s2/2),deg),
P.p[++P.n]=A+Rotate(Point(-s1/2,s2/2),deg),
P.p[++P.n]=A+Rotate(Point(-s1/2,-s2/2),deg);

这样,我们就可以得到\(4n\)个点,然后将这些点求一个凸包,然后作出凸包面积即可。

这段过程比较板子,可以参考这篇博客:初学计算几何(四)——初识凸包

具体实现详见代码吧。


代码

#include<bits/stdc++.h>
#define N 600
using namespace std;
int n;
namespace ComputationGeometry
{
    #define eps 1e-10
    inline int dcmp(double x) {return fabs(x)<eps?0:(x>0?1:-1);}
    struct Point
    {
        double x,y;
        Point(double nx=0,double ny=0):x(nx),y(ny){}
    };
    typedef Point Vector;
    inline Vector operator + (Point A,Point B) {return Vector(A.x+B.x,A.y+B.y);}
    inline Vector operator - (Point A,Point B) {return Vector(A.x-B.x,A.y-B.y);}
    inline Vector operator * (Vector A,double x) {return Vector(A.x*x,A.y*x);}
    inline Vector operator / (Vector A,double x) {return Vector(A.x/x,A.y/x);}
    inline bool operator < (Vector A,Vector B) {return fabs(A.x-B.x)>eps?A.x<B.x:A.y<B.y;}
    inline double Cro(Vector A,Vector B) {return A.x*B.y-A.y*B.x;}
    inline Point Rotate(Vector A,double deg)//旋转
    {
        register double rad=deg*acos(-1)/180;dcmp(rad)<0&&(rad+=2*acos(-1)),rad=-rad;
        return Vector(A.x*cos(rad)-A.y*sin(rad),A.x*sin(rad)+A.y*cos(rad));
    } 
    struct Polygon
    {
        int n;Point p[(N<<2)+5];
        Polygon() {n=0;}
    };
    typedef Polygon ConvexHull;
    inline double PolygonArea(Polygon P)//求多边形(凸包)面积
    {
        register int i;register double res=0;
        for(i=2;i<P.n;++i) res+=Cro(P.p[i]-P.p[1],P.p[i+1]-P.p[1])/2;
        return res;
    }
    inline ConvexHull GetConvexHull(Polygon P)//求凸包
    {
        register int i,t;register Polygon S=P;register ConvexHull res;
        for(sort(S.p+1,S.p+S.n+1),i=1;i<=S.n;++i)
        {
            while(res.n>1&&dcmp(Cro(res.p[res.n]-res.p[res.n-1],S.p[i]-res.p[res.n-1]))<=0) --res.n;
            res.p[++res.n]=S.p[i];
        }
        for(t=res.n,i=S.n-1;i;--i)
        {
            while(res.n>t&&dcmp(Cro(res.p[res.n]-res.p[res.n-1],S.p[i]-res.p[res.n-1]))<=0) --res.n;
            res.p[++res.n]=S.p[i];
        }
        return res;
    }
};
using namespace ComputationGeometry;
int main()
{
    freopen("a.in","r",stdin);
    register int test_tot,i;register double x,y,s1,s2,deg,sum;register Point A;register Polygon P;register ConvexHull H;scanf("%d",&test_tot);
    while(test_tot--)
    {
        for(scanf("%d",&n),P.n=sum=0,i=1;i<=n;++i) 
        {
            scanf("%lf%lf%lf%lf%lf",&x,&y,&s1,&s2,&deg),A=Point(x,y),sum+=s1*s2,//同时统计矩形面积和
            P.p[++P.n]=A+Rotate(Point(s1/2,s2/2),deg),P.p[++P.n]=A+Rotate(Point(s1/2,-s2/2),deg),P.p[++P.n]=A+Rotate(Point(-s1/2,s2/2),deg),P.p[++P.n]=A+Rotate(Point(-s1/2,-s2/2),deg);//存储点
        }
        printf("%.1lf %%\n",100*sum/PolygonArea(GetConvexHull(P)));//计算答案
    }
    return 0;
}

【UVA10652】Board Wrapping(求凸包面积)

标签:大致   tps   name   告诉   pre   mes   ping   space   put   

原文地址:https://www.cnblogs.com/chenxiaoran666/p/UVA10652.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!