#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; typedef long long int lint; #define rep(i, n) for (lint i = 0; i < n; i++) #define sort(v) sort((v).begin(), (v).end()) #define reverse(v) reverse((v).begin(), (v).end()) #define upper(v,hoge) upper_bound(v.begin(),v.end(),hoge) #define lower(v,hoge) lower_bound(v.begin(),v.end(),hoge) #define llower(v,hoge) *lower_bound(v.begin(), v.end(), hoge) #define lupper(v,hoge) *upper_bound(v.begin(), v.end(), hoge) #define mp make_pair #define IP pair #define enld endl bool isOK(lint N, lint P, lint Q) { if (N*N>P+Q*N*log2(double(N))) return true; else return false; } // 汎用的な二分探索のテンプレ lint binary_search(lint P, lint Q) { lint left = -1; //「index = 0」が条件を満たすこともあるので、初期値は -1 lint right = 1000000000; // 「index = a.size()-1」が条件を満たさないこともあるので、初期値は a.size() /* どんな二分探索でもここの書き方を変えずにできる! */ while (right - left > 1) { lint mid = left + (right - left) / 2; if (isOK(mid, P,Q)) right = mid; else left = mid; } /* left は条件を満たさない最大の値、right は条件を満たす最小の値になっている */ return left; } int main() { double P, Q; cin >> P >> Q; lint now1 = binary_search(P, Q); double now = double(now1); double d = 0.00001; while (1) { if (1-Q/now*log2(now)>P/now/now) { cout << fixed << setprecision(10) << now - d << endl; return 0; } else { now += d; } } }