poj 2115 扩展欧几里德

#include<stdio.h>
#include<string.h>
#define max 32

typedef long long LL;
LL pow2[max+5];

void init(){
    for(int i=1;i<=max;i++){
        pow2[i]=1LL<<i;
    }
}

LL a,b,c,k;

void  gcd(LL a,LL b,LL& d,LL& x,LL& y){
    if(!b){
        d=a; x=1; y=0; return;
    }
    gcd(b,a%b,d,y,x); y-=(a/b)*x;
}

int main(){
    init();
    while(scanf("%lld%lld%lld%lld",&a,&b,&c,&k)&&a+b+c+k!=0){
                LL d,x,y;
        if(a==b){
            puts("0"); continue;
        }
        gcd(c,pow2[k],d,x,y);
        if((b-a)%d){
            puts("FOREVER");
        }
        else{
            x=x*(b-a)/d;
            x=(x%(pow2[k]/d)+(pow2[k]/d))%(pow2[k]/d);
            printf("%lld
",x);
        }
    }
}