Advertisement
Guest User

Untitled

a guest
Dec 6th, 2019
111
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C++ 1.44 KB | None | 0 0
  1. #include <bits/stdc++.h>
  2. using namespace std;
  3.  
  4. #define FOR(i, k, l) for(int i = k; i < l; i++)
  5. #define DFOR(i, k, l) for(int i = k; i > l; i--)
  6. #define FA(i, k) for(int i = 0; i < int(k.size()); i++)
  7. #define DFA(i, k) for(int i = int(k.size()) - 1; i > -1; i--)
  8. #define ASKS(i) for(cin >> (i); (i)--;)
  9.  
  10. #define MAXINT 2147483647
  11. #define endl '\n'
  12. #define all(x) x.begin(), x.end()
  13. #define fi first
  14. #define se second
  15.  
  16. typedef long long ll;
  17. typedef long double ld;
  18. typedef pair<int, int> pii;
  19. typedef pair<ll, ll> pll;
  20. typedef vector<int> vi;
  21. typedef vector<ll> vll;
  22. typedef set<int> si;
  23. typedef set<ll> sll;
  24. typedef priority_queue<int> pqi;
  25. typedef priority_queue<ll> pqll;
  26. typedef priority_queue<int, vector<int>, greater<int>> pqig;
  27. typedef priority_queue<ll, vector<ll>, greater<ll>> pqllg;
  28.  
  29. ld g(ld p, ld f, ld a, ld b){
  30.     ld tf = sqrt(a * a + (1 - b) * (1 - b)) / f;
  31.     ld tp = sqrt((1 - a) * (1 - a) + b * b) / p;
  32.     return tp + tf;
  33. }
  34.  
  35. int main(){
  36.     ios::sync_with_stdio(0);
  37.     cin.tie(0);
  38.     freopen("forest.in", "r", stdin);
  39.     freopen("forest.out", "w", stdout);
  40.     cout << fixed << setprecision(12);
  41.  
  42.     ld p, f, a, l = 0, r = 1, m1, m2;
  43.     cin >> p >> f >> a;
  44.     while(r - l > 0.00000000000001){
  45.         m1 = l + (r - l) / 3;
  46.         m2 = r - (r - l) / 3;
  47.         if(g(p, f, a, m1) > g(p, f, a, m2))
  48.             l = m1;
  49.         else
  50.             r = m2;
  51.     }
  52.     cout << l;
  53.  
  54.  
  55.     return 0;
  56. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement