poj 2115 扩展欧几里德

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

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

 void init(){
     for(int i=;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=; y=; 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!=){
                 LL d,x,y;
         if(a==b){
             puts(""); 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\n",x);
         }
     }
 }