forked from ngthanhtrung23/CompetitiveProgramming
-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathC.cpp
More file actions
85 lines (75 loc) · 2.81 KB
/
C.cpp
File metadata and controls
85 lines (75 loc) · 2.81 KB
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
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
#include <bits/stdc++.h>
using namespace std;
#define fi first
#define se second
#define pb push_back
#define mp make_pair
#define Fit(i,c) for (__typeof__((c).begin()) i = (c).begin(); i != (c).end(); i++)
#define inf 1000000005
#define all(a) (a).begin(), (a).end()
#define ms(a,x) memset(a, x, sizeof(a))
#define mod 1000000007
#define sz(a) ((int)(a).size())
template<class T> int getbit(T s, int i) { return (s >> i) & 1; }
template<class T> T onbit(T s, int i) { return s | (T(1) << i); }
template<class T> T offbit(T s, int i) { return s & (~(T(1) << i)); }
template<class T> int cntbit(T s) { return __builtin_popcount(s);}
#define Rep(i,n) for(int i = 0; i < (n); ++i)
#define Repd(i,n) for(int i = (n)-1; i >= 0; --i)
#define For(i,a,b) for(int i = (a); i <= (b); ++i)
#define Ford(i,a,b) for(int i = (a); i >= (b); --i)
typedef unsigned long long ull;
typedef long long ll;
typedef double ld;
#define eps 1e-12
typedef pair<int, int> II;
typedef pair<ll, ll> LL;
template<class T> T gcd(T a, T b){ T r; while (b != 0) { r = a % b; a = b; b = r; } return a;}
template<class T> T lcm(T a, T b) { return a / gcd(a, b) * b; }
#define PI (2 * acos(0))
#define linf (1ll << 60)
#define maxn 100005
ld vw, vc;
int n;
int x[maxn], y[maxn];
ld f[maxn];
ld sqr(ld A){
return A * A;
}
int main(){
// freopen("in.txt", "r", stdin);
freopen("c.in", "r", stdin);
freopen("c.out", "w", stdout);
ios_base::sync_with_stdio(false);
cin.tie(NULL);
cin >> n >> vw >> vc;
for(int i = 0; i <= n; i++) f[i] = (1ll << 60);
For(i, 1, n) cin >> x[i] >> y[i];
f[1] = 0;
for(int i = 1; i <= n; i++){
if(i != 1){
if(y[i] < y[i - 1]){
for(int j = i - 1; j > 1; j--) if(y[j] > y[i] && y[j - 1] <= y[i]){
ld X = (x[j] * (ld)(1.0) * (y[i] - y[j - 1]) + x[j - 1] * (ld)(1.0) * (y[j] - y[i])) / (y[j] - y[j - 1]);
ld d = sqrt(sqr(X - x[j - 1]) + sqr(y[i] - y[j - 1]));
f[i] = min(f[j - 1] + d / vw + (x[i] - X) / vc, f[i]);
break;
}
}
}
if(i < n){
ld d = sqrt(sqr(x[i] - x[i + 1]) + sqr(y[i] - y[i + 1]));
f[i + 1] = min(f[i] + d / vw, f[i + 1]);
if(y[i] < y[i + 1]){
for(int j = i + 1; j < n; j++) if(y[j] > y[i] && y[j + 1] <= y[i]){
ld X = (x[j] * (ld)(1.0) *(y[i] - y[j + 1]) + x[j + 1] * (ld)(1.0) * (y[j] - y[i])) / (y[j] - y[j + 1]);
ld d = sqrt(sqr(X - x[j + 1]) + sqr(y[i] - y[j + 1]));
f[j + 1] = min(f[i] + d / vw + (X - x[i]) / vc, f[j + 1]);
break;
}
}
}
}
cout << fixed << setprecision(6) << f[n] << endl;
return 0;
}