题意:
对于C的for(i=A ; i!=B ;i +=C)循环语句,问在k位存储系统中循环几次才会结束。
若在有限次内结束,则输出循环次数。
否则输出死循环。
(k位==mod $2^{k}$)
列出方程:$A+Cx\equiv B(mode\quad 2^{k})$
转换一下:$Cx+ky=B-A$
用exgcd解出 $Cx+ky=gcd(C,k)$
然后把求出的$x*(B-A)/gcd(C,k)$
再$\% (k/gcd(C,k))$求个最小正整数解
end.
#include<iostream>
#include<cstdio>
#include<cstring>
#define re register
using namespace std;
typedef long long ll;
ll A,B,C,k,g,x0,y0,q;
void exgcd(ll a,ll b,ll &x,ll &y){
if(!b) g=a,x=,y=;
else exgcd(b,a%b,y,x),y-=x*(a/b);
}
int main(){
while(cin>>A>>B>>C>>k){
if(!A&&!B&&!C&&!k) break;
k=1ll<<k; //注意long long 用位运算 要 1ll
exgcd(C,k,x0,y0); q=k/g;
if((B-A)%g) cout<<"FOREVER"<<endl;
else{
x0=(x0*(B-A)/g%q+q)%q;
cout<<x0<<endl;
}
}return ;
}