考试题T3

题意分析

这题一看没有什么思路

幸好我们机房的红太阳\(ghj1222\)切了这道题

首先我们考虑风跑一个来回之后人怎么样

就是跑了一个区间

也就是风跑了若干个来回之后 人跑了若干个区间

所以我们考虑求出距离最小的那个区间

距离是一个单峰函数 所以我们用三分求

然后的话 问题就是求两条有向线段之间的最小距离

这也是一个单峰函数 所以我们用三分求

---

我们需要处理出来风从一头走到另一头的人的走的向量

\[\frac{\overrightarrow{v}}{|\overrightarrow{v}|}\frac{v_w}{v_t}dis\]

$dis: $风走的距离

CODE:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<cstdlib>
#include<string>
#include<queue>
#include<map>
#include<stack>
#include<list>
#include<set>
#include<deque>
#include<vector>
#include<ctime>
#define ll long long
#define inf 0x7fffffff
#define N 500008
#define IL inline
#define M 1008611
#define D double
#define eps 1e-7
#define R register
using namespace std;
template<typename T>IL void read(T &_)
{
    T __=0,___=1;char ____=getchar();
    while(!isdigit(____)) {if(____=='-') ___=0;____=getchar();}
    while(isdigit(____)) {__=(__<<1)+(__<<3)+____-'0';____=getchar();}
    _=___ ? __:-__;
}
/*-------------OI使我快乐-------------*/
struct Node
{
    D xx,yy;
    friend Node operator +(const Node &A,const Node &B)
    {return (Node){A.xx+B.xx,A.yy+B.yy};}
    friend Node operator -(const Node &A,const Node &B) 
    {return (Node){A.xx-B.xx,A.yy-B.yy};}
    friend Node operator *(const Node &A,const D &B)
    {return (Node){A.xx*B,A.yy*B};}
    friend Node operator /(const Node &A,const D &B)
    {return (Node){A.xx/B,A.yy/B};}
}fro_x,to_x,fro_y,to_y,tmp;
D v_x,v_y,d_x,d_y;
IL D qury_dis(Node now)
{return sqrt((now.xx*now.xx)+(now.yy*now.yy));}
IL D qury_ans(const Node &sx,const Node &tx,const Node &sy,const Node &ty)
{//三分最小的距离
    Node tmpx=tx-sx,tmpy=ty-sy;    
    D le=0,ri=1,ans=0;
    while((ri-le)>eps)
    {
        D mid=(ri-le)/3.0+le,mmid=ri-(ri-le)/3.0;
        Node cdy,wzy,zjz,ghj;
        cdy=sx+tmpx*mid;wzy=sy+tmpy*mid;
        zjz=sx+tmpx*mmid;ghj=sy+tmpy*mmid;
//这里要注意的是 二者的进行是等比例的      
        D dis_cdy=qury_dis(wzy-cdy),dis_wzy=qury_dis(ghj-zjz);
        if(dis_cdy>dis_wzy) le=mid,ans=dis_wzy;
        else ri=mmid,ans=dis_cdy;
    }
    return ans;
}
IL D work(ll now)
{
    Node cdy,wzy,zjz;
    cdy=fro_x+tmp*now*2;wzy=cdy+tmp;zjz=wzy+tmp;
    return min(qury_ans(cdy,wzy,fro_y,to_y),qury_ans(wzy,zjz,to_y,fro_y));
}
int main()
{
    freopen("chaser.in","r",stdin);
    freopen("chaser.out","w",stdout);
    while(scanf("%lf%lf",&fro_x.xx,&fro_x.yy)!=EOF)
    {
        scanf("%lf%lf%lf",&to_x.xx,&to_x.yy,&v_x);
        scanf("%lf%lf%lf%lf%lf",&fro_y.xx,&fro_y.yy,&to_y.xx,&to_y.yy,&v_y);
        tmp=((to_x-fro_x)/qury_dis(to_x-fro_x)*v_x/v_y*qury_dis(to_y-fro_y));
//相当于是风跑一个来回之后    
//人在其方向上走的向量    
//      printf("now is %.2f %.2f\n",tmp.xx,tmp.yy);
        scanf("%lf%lf",&d_x,&d_y);
//三分区间对应的位置       
        ll le=0,ri=1e13;
        while(ri-le>=5)
        {
//ghj1222说 整数三分卡边界比较恶心 所以设置一个误差 然后暴力求就可以了
            ll mid=(ri-le)/3+le,mmid=ri-(ri-le)/3;
            D resa=work(mid),resb=work(mmid);
            if(resa>resb) le=mid;
            else ri=mmid;
        }
        D ans=1e18;
        for(R ll i=le;i<=ri;++i)
        ans=min(ans,work(i));
        if(ans<d_x) puts("Dangerous");
        else if(ans>=d_x&&ans<=d_y) puts("Perfect");
        else puts("Miss");
    }
    fclose(stdin);
    fclose(stdout);
    return 0;
}

转载于:https://www.cnblogs.com/LovToLZX/p/10702745.html