exgcd

#include<cmath>
using namespace std;
typedef long long ll;
int (int a,int b,int &x,int &y){
    if(b==0) {x=1;y=0;return a;}
    int c=exgcd(b,a%b,x,y);
    int z=x;x=y;y=z-y*(a/b);
    return c;
}
int main(){
    int a,b,c;
    
    while(~scanf("%d%d%d",&a,&b,&c)){
        if(a==b&&b==c&&a==0) break;
        int x,y;
        int gcd=exgcd(a,b,x,y);
        // printf("x=%d y=%dn",x,y);
        int mod=b/gcd;
        int t=c/gcd;
        int x1=(t*x%mod+mod)%mod;
        int y1=(c-a*x1)/b;
        int ans1=abs(x1)+abs(y1);
        // printf("mod=%d x1=%d y1=%dn",mod,x1,y1);
        mod=a/gcd;
        int y2=(t*y%mod+mod)%mod;
        int x2=(c-y2*b)/a;
        int ans2=abs(x2)+abs(y2);
        // printf("mod=%d x2=%d y2=%dn",mod,x2,y2);
        if(ans1>ans2){
            printf("%d %dn",abs(x2),abs(y2));
        }
        else printf("%d %dn",abs(x1),abs(y1));
    }
}

//#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll (ll a,ll b,ll &x,ll &y){
    if(b==0) {x=1,y=0;return a;}
    ll d=exgcd(b,a%b,x,y);
    ll z=x;x=y;y=z-y*(a/b);
    return d;
}
int main(){
    ll a,b,c,k;
    while(~scanf("%lld%lld%lld%lld",&a,&b,&c,&k)){
        if(!a&&!b&&!c&&!k) break;
        ll x,y;
        ll l=1;
        l<<=k;
        ll gcd=exgcd(c,l,x,y);
        // printf("gcd=%lldn",gcd);
        if((b-a)%gcd==0){
            ll t=(b-a)/gcd;
            l/=gcd;
            printf("%lldn",(t*x%l+l)%l);
        }
        else printf("FOREVERn");
    }
    //system("pause");
}