[BZOJ1857][SCOI2010]傳送帶(三分套三分)
阿新 • • 發佈:2019-02-20
題目:
題解:
先畫個圖
我們可以發現,答案必然是從AB上的某一個點出發,經過平面,到達CD上的另一點然後到達D。即答案可以表示為
我們要找到最優解,其實就是在AB上確定一個點F,然後在CD上確定一個E
根據我的胡亂猜測,這個點的位置對於答案的大小是單峰的,單峰就用三分啊
那麼我們三分F點,然後三分E點,也就是三分套三分,可以解決此題
題解部分引自這個up,侵刪
程式碼:
#include <cstdio>
#include <cmath>
using namespace std;
const double eps=1e-4;
int 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 f(double x,double y)
{
double lx=cx,ly=cy,rx=dx,ry=dy,x1,y1,x2,y2;
while (fabs(rx-lx)>= eps || fabs(ry-ly)>=eps)
{
x1=lx+(rx-lx)/3; y1=ly+(ry-ly)/3;
x2=rx-(rx-lx)/3; y2=ry-(ry-ly)/3;
double f1=dis(x,y,x1,y1)/(r*1.0)+dis(x1,y1,dx,dy)/(q*1.0);
double f2=dis(x,y,x2,y2)/(r*1.0)+dis(x2,y2,dx,dy)/(q*1.0);
if (f1<f2) rx=x2,ry=y2;
else lx=x1,ly=y1;
}
return dis(x,y,lx,ly)/(r*1.0)+dis(lx,ly,dx,dy)/(q*1.0)+dis(x,y,ax,ay)/(p*1.0);
}
int main()
{
scanf("%d%d%d%d",&ax,&ay,&bx,&by);
scanf("%d%d%d%d",&cx,&cy,&dx,&dy);
scanf("%d%d%d",&p,&q,&r);
double lx=ax,ly=ay,rx=bx,ry=by,x1,y1,x2,y2;
while (fabs(rx-lx)>=eps || fabs(ry-ly)>=eps)
{
x1=lx+(rx-lx)/3; y1=ly+(ry-ly)/3;
x2=rx-(rx-lx)/3; y2=ry-(ry-ly)/3;
if (f(x1,y1)<f(x2,y2)) rx=x2,ry=y2;
else lx=x1,ly=y1;
}
printf("%.2lf",f(lx,ly));
}