POJ 2115

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
typedef long long LL;
void Exgcd(LL a, LL b, LL& d, LL& x, LL& y) {
    if(b == 0) {d = a, x = 1, y = 0;}
    else Exgcd(b,a%b,d,y,x),y -= x*(a/b);
}
int main() {
    LL A,B,C,K,D;
    while(cin >> A >> B >> C >> K && (A + B + C + K)) {
        LL X,Y;
        LL tmp = B - A;
        LL bb = 1LL << K;
        Exgcd(C,bb,D,X,Y);
        if(tmp % D) cout << "FOREVER\n";
        else {
            X = tmp / D * X;
            LL tt = bb / D;
            cout << (X%tt+ tt) % tt << endl;
        }
    }
}