#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
__int64 euclid(__int64 a,__int64 b,__int64 &x,__int64 &y)
{
	if(b==0)
	{
		x=1;
		y=0;
		return a;
	}
	__int64 t=euclid(b,a%b,x,y);
	__int64 temp=x;
	x=y;
	y=temp-a/b*y;
    return t;
}
int main()
{
	__int64 A,B,C,k,d,a,b,c,x,y;
	while(scanf("%I64d%I64d%I64d%I64d",&A,&B,&C,&k))
	{
		if(A+B+C+k==0)
			break;
		a=C;
		b=1<<k;
		c=B-A;
		d=euclid(a,b,x,y);
		if(c%d)
		{
			cout<<"FOREVER"<<endl;
			continue;
		}
		x=((c/d)*x);
		y=c/d*y;
		x=(x%(b/d)+(b/d))%(b/d);
		printf("%I64d\n",x);
	}
	return 0;
}