BZOJ 1857 [Scoi2010]传送带

三分。

三分AB上的点,再三分CD上的点即可(目测但是我太弱不会证,以后会证了再补证明过程)。

要注意线段退化成点的情况,所以最好不要一边三分一边统计答案(在这里被坑QAQ)

#include<cmath>
#include<cstdio>
#include<algorithm>
using namespace std;
const double EPS=1e-3;
double ax,ay,bx,by,cx,cy,dx,dy,P,Q,R;
double dis(double x1,double y1,double x2,double y2)
{
    return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
}
double tri(double x, double y)
{
    double lx=cx,ly=cy,rx=dx,ry=dy,mx1,my1,mx2,my2,r1,r2;
    while(abs(lx-rx)>EPS || abs(ly-ry)>EPS)
    {
        mx1=lx+(rx-lx)/3;mx2=lx+2*(rx-lx)/3;
        my1=ly+(ry-ly)/3;my2=ly+2*(ry-ly)/3;
        r1=dis(x,y,mx1,my1)/R+dis(mx1,my1,dx,dy)/Q;
        r2=dis(x,y,mx2,my2)/R+dis(mx2,my2,dx,dy)/Q;
        if(r1>r2)lx=mx1,ly=my1;
        else rx=mx2,ry=my2;
    }
    return dis(x,y,lx,ly)/R+dis(lx,ly,dx,dy)/Q;
}
int main()
{
    scanf("%lf%lf%lf%lf%lf%lf%lf%lf%lf%lf%lf",&ax,&ay,&bx,&by,&cx,&cy,&dx,&dy,&P,&Q,&R);
    double lx=ax,ly=ay,rx=bx,ry=by,mx1,my1,mx2,my2,r1,r2;
    while(abs(lx-rx)>EPS || abs(ly-ry)>EPS)
    {
        mx1=lx+(rx-lx)/3;mx2=lx+2*(rx-lx)/3;
        my1=ly+(ry-ly)/3;my2=ly+2*(ry-ly)/3;
        r1=tri(mx1,my1)+dis(ax,ay,mx1,my1)/P;
        r2=tri(mx2,my2)+dis(ax,ay,mx2,my2)/P;
        if(r1>r2)lx=mx1,ly=my1;
        else rx=mx2,ry=my2;
    }
    printf("%.2lf\n",tri(lx,ly)+dis(ax,ay,lx,ly)/P);
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值