#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #define _overload3(_1,_2,_3,name,...) name #define _rep(i,n) repi(i,0,n) #define repi(i,a,b) for(ll i=ll(a);ibool chmax(T &a, const T &b) { if (abool chmin(T &a, const T &b) { if (b void cins(itr first,itr last){ for (auto i = first;i != last;i++){ cin >> (*i); } } template void array_output(itr start,itr goal){ string ans = "",k = " "; for (auto i = start;i != goal;i++) ans += to_string(*i)+k; if (!ans.empty()) ans.pop_back(); PRINT(ans); } ll gcd(ll a, ll b) { return a ? gcd(b%a,a) : b; } const ll INF = 1e18; const ll MOD = 1000000007; const ll MOD2 = 998244353; typedef pair P; typedef pair point; const ll MAX = 200000; constexpr ll nx[8] = {1,0,-1,0,1,-1,1,-1}; constexpr ll ny[8] = {0,1,0,-1,1,1,-1,-1}; int main(){ cin.tie(0); ios::sync_with_stdio(false); ll a,b,c,d; cin >> a >> b >> c >> d; ll D = (a-c)*(a-c)-4*2*(b-d); if (D < 0){ PRINT("No"); return 0; } else if (D == 0){ PRINT("Yes"); return 0; } double m = a-c,n = b-d; double p = (-m+sqrt((double)D))/4.0,q = (-m-sqrt((double)D))/4.0; double py = p*p+a*p+b,qy = q*q+a*q+b,t = (py-qy)/(p-q); printf("%.10lf %.10lf\n",t,py-t*p); }