#include<iostream>
#include<iomanip>
#include<cmath>
using namespace std;
long a,b,c,d;
main()
{
	cin>>a>>b>>c>>d;
	long D=(a-c)*(a-c)-8*(b-d);
	if(D==0)
	{
		cout<<"Yes"<<endl;
	}
	else if(D<0)
	{
		cout<<"No"<<endl;
	}
	else
	{
		double ds=sqrt(D);
		double x1,x2;
		if(c-a>=0)
		{
			x2=(c-a+ds)/4;
			x1=(b-d)/x2/2;
		}
		else
		{
			x1=(c-a-ds)/4;
			x2=(b-d)/x1/2;
		}
		double y1=x1*x1+a*x1+b,y2=x2*x2+a*x2+b;
		cout<<fixed<<setprecision(16)<<(y2-y1)/(x2-x1)<<" "<<y1-(y2-y1)/(x2-x1)*x1<<endl;
	}
}