#include <iostream>
#include <math.h>
using namespace std;
int main ()
{
double a,xz,eps;
int p;
cin>>a>>p>>xz>>eps;
double x=pow(a,1.0/p);
double xn=xz,x_prev=-10000000;
int i=0;
while (fabs(xn-x_prev)>=eps)
{
x_prev=xn;
xn=(x_prev/(p*p))*((p*p-1)+0.5*(p+1)*a/(pow(x_prev,p))-0.5*(p-1)*pow(x_prev,p)/(a));
i++;
}
cout<<i<<' '<<xn<<' '<<x<<endl;
return 0;
}
I2luY2x1ZGUgPGlvc3RyZWFtPgojaW5jbHVkZSA8bWF0aC5oPgp1c2luZyBuYW1lc3BhY2Ugc3RkOwoKaW50IG1haW4gKCkgCnsKICAgIGRvdWJsZSBhLHh6LGVwczsKICAgIGludCBwOwogICAgY2luPj5hPj5wPj54ej4+ZXBzOwogICAgZG91YmxlIHg9cG93KGEsMS4wL3ApOwogICAgZG91YmxlIHhuPXh6LHhfcHJldj0tMTAwMDAwMDA7CiAgICBpbnQgaT0wOwogICAgd2hpbGUgKGZhYnMoeG4teF9wcmV2KT49ZXBzKQogICAgeyAgCiAgICAgICAgeF9wcmV2PXhuOwogICAgICAgIHhuPSh4X3ByZXYvKHAqcCkpKigocCpwLTEpKzAuNSoocCsxKSphLyhwb3coeF9wcmV2LHApKS0wLjUqKHAtMSkqcG93KHhfcHJldixwKS8oYSkpOwogICAgICAgIGkrKzsKICAgICAgIAogICAgfQogICAgY291dDw8aTw8JyAnPDx4bjw8JyAnPDx4PDxlbmRsOyAKICAgIHJldHVybiAwOwp9