标签：传送带 double rx cx ay ax NC20276 lx SCOI2010

题目

• 原题地址：[SCOI2010]传送带
• 题目编号：NC20276
• 题目类型：三分
• 时间限制：C/C++ 1秒，其他语言2秒
• 空间限制：C/C++ 262144K，其他语言524288K

1.题目大意

• 平面两条线段AB、CD，想从A走到D，在两条线段和平面上走分别有各自的速度，求最短时间

2.题目分析

• 在AB上从A走到X，再在平面中从X走到线段CD上的Y点，再从Y点走到D
• 学习一下三分的板子

3.题目代码

#include<bits/stdc++.h>

using namespace std;

const double eps=1e-4;//锁精度

double ax,ay,bx,by,cx,cy,dx,dy;
double P,Q,R;
double a1,a2,b1,b2;

inline void swap(double &x,double &y){
    double c;
    c=x,x=y,y=c;
}
inline double f(double x,double a,double b){//已知x以及两系数,求y 
    return a*x+b;
}
inline double far(double x,double y,double m,double n){//求两点之间的距离
    return sqrt((x-m)*(x-m)+(y-n)*(y-n));
}
inline double fi(double x,double y,double a,double b){
    double tim = far(x,y,a,b)/R;
    tim += far(x,y,dx,dy)/Q;
    return tim;
}
inline double suan(double x,double y){
    double tim = far(x,y,ax,ay)/P;//A到(x,y)的时间
    double l,r,ans,L,R;
    if(cx==dx){
        l = cy,r = dy;
        if(l>r){
            swap(l,r);
        }
        while(r-l>=eps){
            double lx=l+(r-l)/3,rx=r-(r-l)/3;
            L = fi(cx,lx,x,y), R = fi(cx,rx,x,y);
            if(L<=R){
                ans = L;
                r = rx;
                continue;
            } else {
                ans = R;
                l = lx;
            }
        }
    } else {
        l = cx,r = dx;
        if(l>r){
            swap(l,r);
        }
        while(r-l>=eps){//三分由(x,y)到CD上的某点 
            double lx = l + (r-l)/3,rx = r - (r-l)/3;
            L = fi(lx,f(lx,a2,b2),x,y), R = fi(rx,f(rx,a2,b2),x,y);
            if(L<=R){
                ans = L;
                r = rx;
                continue;
            } else {
                ans = R;
                l = lx;
            }
        }
    }
    return tim + ans;
}
int main(){
    cin >> ax >> ay >> bx >> by >> cx >> cy >> dx >> dy;
    cin >> P >> Q >> R;
    if(ax!=bx){
        a1 = (ay-by) / (ax-bx);
        b1 = (ay-a1*ax);
    }
    if(cx!=dx){
        a2 = (cy-dy) / (cx-dx);
        b2 = (cy-a2*cx);
    }
    double l,r,ans,L,R;
    if(ax==bx){
        l = ay,r = by;
        if(l>r){
            swap(l,r);
        } 
        while(l<=r){
            double lx = l+(r-l)/3,rx = r-(r-l)/3;
            L = suan(ax,lx),R = suan(ax,rx);
            if(L<=R){
                ans = L;
                r = rx - eps;
                continue;
            } else {
                ans = R;
                l = lx + eps;
            }
        }
    } else {
        l = ax,r = bx;
        if(l>r){
            swap(l,r);
        }
        while(r-l>=eps){//三分由A到AB中的一点(不走AB即为从A到A) 
            double lx=l+(r-l)/3,rx=r-(r-l)/3;
            L=suan(lx,f(lx,a1,b1)),R=suan(rx,f(rx,a1,b1));
            if(L<=R){
                r = rx;
                ans = L;
                continue;
            } else {
                l = lx;
                ans = R;
            } 
        }
    } 
    printf("%.2lf",ans);
    return 0;
}

• 参考题解

标签：传送带,double,rx,cx,ay,ax,NC20276,lx,SCOI2010
来源： https://www.cnblogs.com/zhangyi101/p/16612444.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。