【BZOJ】1857: [Scoi2010]传送带（三分）

http://www.lydsy.com/JudgeOnline/problem.php?id=1857

好神奇的三分。。

第一次写三分啊sad。。看了题解啊题解QAQ

首先发现无论怎么走一定是在AB和CD上选了两个点然后走的（包括ABCD四个点），所以我们就是要找出这两个点就行了。

且AB上有且只有一个最优点，而每一个AB上的点也对应CD唯一一个最优点orz。

所以我们三分AB上的点（酷炫但是不知道为啥这是个单峰函数orz），然后对应三分CD上的点（这个比前面的好证。。这个初中就学了的吧。sad）然后就行了。。

距离和费用就是三段的和。

+ View Code?

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52

#include <cstdio> #include <cstring> #include <cmath> #include <string> #include <iostream> #include <algorithm> #include <queue> using namespace std; #define rep(i, n) for(int i=0; i<(n); ++i) #define for1(i,a,n) for(int i=(a);i<=(n);++i) #define for2(i,a,n) for(int i=(a);i<(n);++i) #define for3(i,a,n) for(int i=(a);i>=(n);--i) #define for4(i,a,n) for(int i=(a);i>(n);--i) #define CC(i,a) memset(i,a,sizeof(i)) #define read(a) a=getint() #define print(a) printf("%d", a) #define dbg(x) cout << (#x) << " = " << (x) << endl #define printarr2(a, b, c) for1(_, 1, b) { for1(__, 1, c) cout << a[_][__]; cout << endl; } #define printarr1(a, b) for1(_, 1, b) cout << a[_] << '\t'; cout << endl inline const int getint() { int r=0, k=1; char c=getchar(); for(; c<'0'||c>'9'; c=getchar()) if(c=='-') k=-1; for(; c>='0'&&c<='9'; c=getchar()) r=r*10+c-'0'; return k*r; } inline const int max(const int &a, const int &b) { return a>b?a:b; } inline const int min(const int &a, const int &b) { return a<b?a:b; }   const double eps=1e-3; int ax, ay, bx, by, cx, cy, dx, dy; double P, Q, R; double dis(double x1, double y1, double x2, double y2) { return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2)); } double cal(double x, double y) {     double x1, x2, y1, y2, lx=cx, ly=cy, rx=dx, ry=dy;     while(abs(rx-lx)>eps || abs(ry-ly)>eps) {         x1=lx+(rx-lx)/3;    y1=ly+(ry-ly)/3;         x2=lx+(rx-lx)/3*2;  y2=ly+(ry-ly)/3*2;         double d1=dis(x1, y1, x, y)/R+dis(x1, y1, dx, dy)/Q;         double d2=dis(x2, y2, x, y)/R+dis(x2, y2, dx, dy)/Q;         if(d1<d2) rx=x2, ry=y2;         else lx=x1, ly=y1;     }     return dis(x, y, ax, ay)/P+dis(lx, ly, x, y)/R+dis(lx, ly, dx, dy)/Q; } int main() {     cin >> ax >> ay >> bx >> by >> cx >> cy >> dx >> dy >> P >> Q >> R;     double x1, x2, y1, y2, lx=ax, ly=ay, rx=bx, ry=by;     while(abs(rx-lx)>eps || abs(ry-ly)>eps) {         x1=lx+(rx-lx)/3;    y1=ly+(ry-ly)/3;         x2=rx-(rx-lx)/3;    y2=ry-(ry-ly)/3;         double d1=cal(x1, y1), d2=cal(x2, y2);         if(d1<d2) rx=x2, ry=y2;         else lx=x1, ly=y1;     }     printf("%.2lf", cal(lx, ly));     return 0; }

 

 

---

 

 

Description

在一个2维平面上有两条传送带，每一条传送带可以看成是一条线段。两条传送带分别为线段AB和线段CD。lxhgww在AB上的移动速度为P，在CD上的移动速度为Q，在平面上的移动速度R。现在lxhgww想从A点走到D点，他想知道最少需要走多长时间

Input

输入数据第一行是4个整数，表示A和B的坐标，分别为Ax，Ay，Bx，By 第二行是4个整数，表示C和D的坐标，分别为Cx，Cy，Dx，Dy 第三行是3个整数，分别是P，Q，R

Output

输出数据为一行，表示lxhgww从A点走到D点的最短时间，保留到小数点后2位

Sample Input

0 0 0 100
100 0 100 100
2 2 1

Sample Output

136.60

HINT

对于100%的数据，1<= Ax，Ay，Bx，By，Cx，Cy，Dx，Dy<=1000
1<=P，Q，R<=10

Source

Day2