POJ 2115 C Looooops (扩展欧几里得算法)

#include<iostream>
#include<cstring>
#include<string>
#include<cstdio>
#include<algorithm>
using namespace std;
typedef long long LL;
LL gcdEx(LL a,LL b,LL& x,LL& y)
{   //求ax+by=gcd(a,b) 的整数解 返回gcd(a,b); if(b==0) {x=1; y=0; return a;}LL x1,y1;int gcd=gcdEx(b,a%b,x1,y1);x=y1;y=x1-a/b*y1;return gcd;
}
int main()
{LL A,B,C,K,x,y,d,k,s,ans;while(cin>>A>>B>>C>>K){if(A==0&&B==0&&C==0&&K==0) break;d=gcdEx(C,(1LL<<K),x,y);  //要把1转化为long long 类型，不然会溢出 if((B-A)%d){cout<<"FOREVER"<<endl; continue;}k=(B-A)/d;ans=k*x; s=(1LL<<K)/d;cout<<(ans%s+s)%s<<endl;}return 0;
}