1. 程式人生 > >[BZOJ1857][SCOI2010]傳送帶(三分套三分)

[BZOJ1857][SCOI2010]傳送帶(三分套三分)

題目:

我是超連結

題解:

先畫個圖
這裡寫圖片描述
我們可以發現,答案必然是從AB上的某一個點出發,經過平面,到達CD上的另一點然後到達D。即答案可以表示為
|AF|/P+|FE|/R+|ED|/Q,FABECD
我們要找到最優解,其實就是在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)); }