The 2018 ACM-ICPC Chinese Collegiate Programming Contest B. Rolling The Polygon

最新推荐文章于 2019-08-31 22:19:10 发布

倚剑笑紅尘

最新推荐文章于 2019-08-31 22:19:10 发布

阅读量271

点赞数

本文链接：https://blog.csdn.net/qq_36782366/article/details/81124819

版权

题目链接：

题目大意：按照逆时针绕向给出一个凸多边形的 n 个顶点 P 0 ,P 1 ,··· ,P n−1 ，再给出凸多边形内部（含边界）
一点 Q。现在要将这个凸多边形在地上无滑动地滚动一周，初始时 P 0 P 1 边与地面接触，假设当前是P i P (i+1) mod n 边与地面接触，那么滚动一下之后则是 P (i+1) mod n P (i+2) mod n 边与地面接触。不难发现，从初始状态滚动 n 下之后 P 0 P 1 边再一次与地面接触，这时认为凸多边形已经滚动了一周。现在你需要求出凸多边形滚动一周之后点 Q 经过的轨迹长度。

题目思路：每次以pi~pi+1滚动，Q滚动的距离即绕pi+1这个点，半径为Q到pi+1的距离，旋转角度为向量pi~pi+1和向量pi+1~pi+2之间的夹角。然后依次累加即可

#include<bits/stdc++.h>
using namespace std;
typedef long double db;
const int MAXN=55;
const db PI=acos(-1.0);
struct Point
{
    db x,y;
    Point(){}
    Point(db _x,db _y):x(_x),y(_y){}
    Point operator - (const Point &t)const
    {
        return Point(x-t.x,y-t.y);
    }
    db operator ^ (const Point &t)const
    {
        return x*t.x+y*t.y;
    }
    db len()const
    {
        return sqrt(x*x+y*y);
    }
}p[MAXN];
db ang(Point a,Point b,Point c)
{
    db t=((b-a)^(c-a))/(b-a).len()/(c-a).len();
    return acos(max(-1.0L,min(1.0L,t)));//这步关键啊  我WA了无数次
}
int main()
{
    int T;
    scanf("%d",&T);
    for(int ca=1;ca<=T;ca++)
    {
        int n;
        scanf("%d",&n);
        for(int i=0;i<=n;i++)
        {
            double x,y;
            scanf("%lf%lf",&x,&y);
            p[i]=Point(x,y);
        }
        db res=0;
        for(int i=0;i<n;i++)
            res+=(p[n]-p[i]).len()*(PI-ang(p[i],p[(i-1+n)%n],p[(i+1)%n]));
        printf("Case #%d: %.3f\n",ca,(double)res);
    }
    return 0;
}