CF 70D Professor's task题解

动态凸包。。只能照着模板敲了。。还不太懂

#include<stdio.h>
#include<math.h>
#include<vector>
#include<algorithm>
using namespace std;
const double pi=acos(-1.0);
double ax,ay,bx,by;
double dis(double x,double y)
{
    return sqrt(x*x+y*y);
}
double get(double a,double w)
{
    double v=0;
    int num=100;
    for(int i=0;i<=100;i++)
    {
        double x=ax+(bx-ax)*i/(double)num;
        double y=ay+(by-ay)*i/(double)num;
        double p=atan2(y,x)-a;
        if(p>2*pi) p-=2*pi;
        while(p<0) p+=2*pi;
        if(p>2*pi-p) p=2*pi-p;
        double t=((w<1e-5)?(i==0?0:1e10):p/w);
        v=max(v,hypot(x-ax,y-ay)/t);
        //printf("%.5lf\n",t);
    }
    return v;
}
double x[10010],y[10010],a[10010],w[10010];
int main()
{
    int n,k;
    while(scanf("%lf%lf%lf%lf",&ax,&ay,&bx,&by)!=EOF)
    {
        scanf("%d",&n);
        for(int i=0;i<n;i++)
            scanf("%lf%lf%lf%lf",&x[i],&y[i],&a[i],&w[i]);
        vector<double>v;
        v.clear();
        for(int i=0;i<n;i++)
        {
            ax-=x[i];ay-=y[i];bx-=x[i];by-=y[i];
            double vv=get(a[i],w[i]);
            printf("%.5lf\n",vv);
            v.push_back(vv);
            ax+=x[i];ay+=y[i];bx+=x[i];by+=y[i];
        }
        sort(v.begin(),v.end());
        reverse(v.begin(),v.end());//·´×ª
        for(int i=0;i<v.size();i++)
            printf("%.5lf\n",v[i]);
        scanf("%d",&k);
        if(k>=v.size()) printf("0.00000\n");
        else printf("%.5lf\n",v[k]);
    }
    return 0;
}