poj2115(扩展欧几里得)

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
typedef __int64 ll;
ll exgcd(ll a,ll b,ll &x,ll &y)
{
ll d;
if(b==0)
{
x=1;y=0;
d=a;
}
else
{
d=exgcd(b,a%b,y,x);
y-=(a/b)*x;
}
return d;
}
ll Abs(ll a)
{
if(a<0) return -a;
return a;
}
int main()
{
//freopen("in.txt","r",stdin);
ll a,b,c,k;
while(scanf("%I64d%I64d%I64d%I64d",&a,&b,&c,&k),a||b||c||k)
{
ll x,y;
k=(1ll<<k);
ll d=exgcd(k,c,x,y);   //这里的(x,y)为kx+cy=gcd(k,c)的解
if((b-a)%d==0)
{
ll e=(b-a)/d*y%k+k;
cout<<e%(k/d)<<endl;
}
else
cout<<"FOREVER"<<endl;
}
return 0;
}