C++ :

#include<bits/stdc++.h>
using namespace std;
int a,b,c,d,s[10005];
int tf(int x,int y)
{
	int a1=x,b1=y,c1=a1%b1;
	while(c1!=0)
	{
		a1=b1;
		b1=c1;
		c1=a1%b1;
	}
	return x*y/b1;
}
void yf(int x,int y)
{
	for(int i=2;i<=10005;i++)
		if(s[i]==0)
			for(int j=i+i;j<=10005;j+=i)
				s[j]=1;
	for(int i=2;i<=min(x,y);)
	{
		if(x%i==0 && y%i==0) x/=i,y/=i;
		else i++;
		if(s[x]==0 && s[y]==0 ) break;
	}
	cout<<x<<" "<<y;
}
int main()
{
    cin>>a>>b>>c>>d;
    int fm=tf(b,d),fz=(fm/b)*a-(fm/d)*c;
    yf(fz,fm);
}