#include<cstdio>
#include<math.h>
#include<algorithm>
using namespace std;
long long l,x,x1s,x1j,x2s,x2j,A,B,C,k,auxp1,auxp2,auxq1,auxq2;
struct sol
{
	long long x,y,z,t;
} m[100000];
bool comp(const sol &A,const sol &B)
{
	if(A.x>B.x) return false;
	if(A.x<B.x) return true;
	if(A.y>B.y) return false;
	return true;
}
long long modul(long long x)
{
	return x<0?(-x):x;
}
void getsol(long long d)
{
	long long p1,p2,q1,q2;
	p1=d;
	q1=-x1s*d/x1j;
	p2=A/d;
	q2=-x2s*A/(d*x2j);
	if((p1*p2==A && q1*q2==C && p1*q2+p2*q1==B) || (p1*p2==-A && q1*q2==-C && p1*q2+p2*q1==-B))
	{
		bool ok=true;
		for(int i=1;i<=l;i++)
			if(m[i].x==p1 && m[i].y==q1 && m[i].z==p2 && m[i].t==q2)
			{
				ok=false;
				break;
			}
		if(ok==true)
		{
			++l;
			m[l].x=p1;
			m[l].y=q1;
			m[l].z=p2;
			m[l].t=q2;
		}
	}
	p1=d;
	q1=-x2s*d/x2j;
	p2=A/d;
	q2=-x1s*A/(d*x1j);
	if((p1*p2==A && q1*q2==C && p1*q2+p2*q1==B) || (p1*p2==-A && q1*q2==-C && p1*q2+p2*q1==-B))
	{
		bool ok=true;
		for(int i=1;i<=l;i++)
			if(m[i].x==p1 && m[i].y==q1 && m[i].z==p2 && m[i].t==q2)
			{
				ok=false;
				break;
			}
		if(ok==true)
		{
			++l;
			m[l].x=p1;
			m[l].y=q1;
			m[l].z=p2;
			m[l].t=q2;
		}
	}
}
long long getx1x2()
{
	long long delta=B*B-4*A*C;
	if((long long)sqrt(delta)*sqrt(delta)!=sqrt(delta)*sqrt(delta))
		return 0;
	x1s=-B-(long long)sqrt(delta);
	x1j=2*A;
	x2s=-B+(long long)sqrt(delta);
	x2j=2*A;
	return 1;
}
int main()
{
	freopen("ecuatie.in","r",stdin);
	freopen("ecuatie.out","w",stdout);
	scanf("%lld%lld%lld%lld",&A,&B,&C,&k);
	if(!getx1x2())
	{
		printf("-1");
		return 0;
	}
	x=(int)sqrt(modul(A));
	for(int i=1;i<=x;i++)
		if(A%i==0)
			getsol(i), getsol(-i), getsol(A/i), getsol(-A/i);
	sort(m+1,m+l+1,comp);
	if(!m[k].x) printf("-1");
	else
	{
		printf("(");
		if(m[k].x==-1) printf("-");
			else if(m[k].x!=1) printf("%lld",m[k].x);
		printf("x");
		if(m[k].y<0) printf("%lld",m[k].y);
			else printf("+%lld",m[k].y);
		printf(")(");
		if(m[k].z==-1) printf("-");
			else if(m[k].z!=1) printf("%lld",m[k].z);
		printf("x");
		if(m[k].t<0) printf("%lld",m[k].t);
			else printf("+%lld",m[k].t);
		printf(")");
	}
	return 0;
}