Codeforces Round #357 (Div. 2) E. Runaway to a Shadow 计算几何

题意

有一个蟑螂，在x0,y0点，每秒移动，可以移动T秒，T秒后不在阴影中就会被拍死

阴影都是圆，现在给你圆心坐标和圆的半径。

这个蟑螂是随机选择一个方向走的

问你这个蟑螂活下来的概率是多少

题解：

蟑螂和圆都会有一个角度的区间，表示在T秒内能够到达这个阴影中

然后把所有区间拿出来，取并集

然后再除以2pi就好了

思路很简单。

取区间的这个东西，用简单的初中几何知识就能得到。

然后这道题就结束了。

#include<bits/stdc++.h>
using namespace std;
const double eps = 1e-6;
const double pi = acos(-1.0);
double sqr(double x)
{
    return x*x;
}
double dis(double x,double y,double x1,double y1)
{
    return sqrt(sqr(x-x1)+sqr(y-y1));
}
vector<pair<double,int> >a;
double x,y,v,t,r;
int n;
int main()
{

    scanf("%lf%lf%lf%lf",&x,&y,&v,&t);
    r=v*t;
    scanf("%d",&n);
    for(int i=0;i<n;i++)
    {
        double x0,y0,r0;
        scanf("%lf%lf%lf",&x0,&y0,&r0);
        double D=dis(x,y,x0,y0);
        if(D<=r0)
        {
            cout<<"1.000000000"<<endl;
            return 0;
        }
        if(r+r0+eps<D)continue;
        double angl,angr,ang;

        double angm=atan2(y0-y,x0-x);

        if(angm<0)angm+=2*pi;

        double len1 = sqrt(D*D-r0*r0);
        if(len1<r+eps){
            ang=asin(r0/D);
        }
        else{
            ang=acos((D*D+r*r-r0*r0)/(2.0*D*r));
        }

        angl=angm-ang;
        angr=angm+ang;

        if(angl<0){
            a.push_back(make_pair(angl+2*pi,1));
            a.push_back(make_pair(2*pi,-1));
            a.push_back(make_pair(0,1));
            a.push_back(make_pair(angr,-1));
        }
        else if(angr>2*pi){
            a.push_back(make_pair(angl,1));
            a.push_back(make_pair(2*pi,-1));
            a.push_back(make_pair(0,1));
            a.push_back(make_pair(angr-2*pi,-1));
        }
        else{
            a.push_back(make_pair(angl,1));
            a.push_back(make_pair(angr,-1));
        }
    }
    sort(a.begin(),a.end());
    double ans = 0;
    double last = 0;
    int now = 0;
    for(int i=0;i<a.size();i++)
    {
        if(now>0)
            ans+=a[i].first-last;
        last=a[i].first;
        now+=a[i].second;
    }
    printf("%.12f\n",ans/(2*pi));
}