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;
}