前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >【CodeForces 613A】Peter and Snow Blower

【CodeForces 613A】Peter and Snow Blower

作者头像
饶文津
发布2020-06-02 14:37:29
3280
发布2020-06-02 14:37:29
举报
文章被收录于专栏:饶文津的专栏饶文津的专栏

题意

给出原点(不是(0,0)那个原点)的坐标和一个多边形的顶点坐标,求多边形绕原点转一圈扫过的面积(每个顶点到原点距离保持不变)。

分析

多边形到原点的最小距离和最大距离构成的两个圆之间的圆环就是所求面积。

判断最大距离一定在顶点上,最小距离可能在点上也可能在边上。

如果原点到一个顶点的连线和它所在的边构成钝角,那么最小距离在点上,否则最小距离就是顶点和该边构成三角形的原点所在的高,可以用海伦公式求三角形面积,然后求高。

不过我用的方法是点到直线的距离公式。

然后π用 acos(-1.0) 表示。

代码

代码语言:javascript
复制
#include <stdio.h>
#include<iostream>
#include <cmath>
#define dd double
using namespace std;

int n;
dd far,nea=1e19,t,P=acos(-1.0);

struct po
{
    double x,y;
} p[100005];

dd P_P(po a,po b)
{
    return (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);
}

dd P_L(po a,po b,po c)//point c --> line ab
{
    dd l1=P_P(a,b);
    dd l2=P_P(a,c);
    dd l3=P_P(b,c);

    if(l2*l2 >= (l1*l1+l3*l3))
        return l3;

    if(l3*l3>=(l2*l2+l1*l1))
        return l2;

    dd A=a.y-b.y;
    dd B=b.x-a.x;
    dd C=a.x*b.y - a.y*b.x;
//line ab: Ax+By+C=0;
    dd D=A*c.x + B*c.y + C;

    return fabs(D*D/(A*A+B*B));
}

int main()
{
    scanf("%d",&n);
    for(int i=0;i<=n;i++)
        scanf("%lf%lf",&p[i].x,&p[i].y);

    for(int i=1;i<=n;i++)
    {
        if (i<n) nea = min( P_L(p[i], p[i+1], p[0]), nea);
        else nea = min( P_L(p[n], p[1], p[0]), nea);

        far=max(far,P_P(p[0],p[i]));
    }
    printf("%.18lf\n",P*(far-nea));

    return 0;
}
本文参与 腾讯云自媒体同步曝光计划,分享自作者个人站点/博客。
原始发表:2016-02-17 ,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 作者个人站点/博客 前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体同步曝光计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
目录
  • 分析
  • 代码
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档