#include <iostream>
#include <iomanip>
#include <cmath>

using lint = long long;
using ldouble = long double;

void solve() {
    lint a, b, c, d;
    std::cin >> a >> b >> c >> d;

    lint p = a - c, q = b - d;
    lint det = p * p - q * 8;

    if (det <= 0) {
        std::cout << (det < 0 ? "No" : "Yes") << "\n";
        return;
    }

    ldouble x1 = (-p + std::sqrt(ldouble(det))) / 4,
            x2 = (-p - std::sqrt(ldouble(det))) / 4;

    ldouble y1 = x1 * x1 + a * x1 + b,
            y2 = x2 * x2 + a * x2 + b;

    ldouble m = (y2 - y1) / (x2 - x1);
    ldouble n = y1 - m * x1;

    std::cout << std::fixed << std::setprecision(10)
              << m << " " << n << "\n";
}

int main() {
    std::cin.tie(nullptr);
    std::ios::sync_with_stdio(false);

    solve();

    return 0;
}