#pragma GCC optimize("Ofast") #include #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; #define rep(i,n) for(int i=0;i<(n);++i) #define reps(i,n) for(int i=1;i<=(n);++i) #define all(x) (x).begin(),(x).end() #define Fixed fixed << setprecision(14) #define int int64_t using pii = pair; constexpr int INF = 0x3f3f3f3f; constexpr long long LINF = 0x3f3f3f3f3f3f3f3fLL; constexpr int mod1 = 1e9+7; constexpr int mod2 = 998244353; constexpr int dy[] = {0,1,0,-1,1,-1,-1,1}; constexpr int dx[] = {1,0,-1,0,1,1,-1,-1}; template inline bool chmax(A &a, const B &b) { return b > a && (a = b, true); } template inline bool chmin(A &a, const B &b) { return b < a && (a = b, true); } template inline A updiv(const A a,const B b){ return (a + b - 1) / b; } template inline A gcd(const A a,const B b){ return b ? gcd(b,a % b) : a; } template inline A lcm(const A a,const B b){ return a / gcd(a,b) * b; } template using min_heap = priority_queue,greater >; template using max_heap = priority_queue; template using umap = unordered_map; signed main(void){ cin.tie(nullptr); ios_base::sync_with_stdio(false); int p, q; cin >> p >> q; long double ok = 1, ng = 1e15; while(ng - ok > 0.00000001){ long double mid = (ng + ok) / 2; if(mid <= p / mid + q * log2(mid)) ok = mid; else ng = mid; } cout << Fixed << ok << '\n'; return 0; }