UVA 1421 Archery(二分+维护区间)

二分的难点还是在判断上。

写这题的时候需要注意精度和L，R分别为角度。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
const int maxn=5000+10;
const double eps=1e-6;
struct node{
    double d,l,r;
}line[maxn];
int n;
double w;
int cmp(const node a,const node b) {return a.d<b.d;}
int check(double mid)
{
    double L=atan2(line[0].d,line[0].r-mid);//角度是逆时针转的
    double R=atan2(line[0].d,line[0].l-mid);
    for(int i=1;i<n;i++)
    {
        double l=atan2(line[i].d,line[i].r-mid);
        double r=atan2(line[i].d,line[i].l-mid);
        if(L-r>eps) return -1;
        if(l-R>eps) return 1;
        L=max(L,l);
        R=min(R,r);
    }
    return 0;
}
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%lf%d",&w,&n);
        for(int i=0;i<n;i++) scanf("%lf%lf%lf",&line[i].d,&line[i].l,&line[i].r);
        sort(line,line+n,cmp);
        double l=0,r=w,mid;
        bool ok=false;
        while(r-l>eps)
        {
            mid=(l+r)/2;
            int k=check(mid);
            if(k==0) {ok=true;break;}
            else if(k==-1) r=mid;
            else l=mid;
        }
        if(ok) printf("YES\n");
        else printf("NO\n");
    }
    return 0;
}