各种不会,在网上搜了不少的资料才得以通过
#include <iostream>
#include <fstream>
#include <cstdio>
using namespace std;
long long x,y,q;
long long extend_Eulid(long long a,long long b){
if(b == 0){
x = 1;y = 0;q = a;
}else{
extend_Eulid(b,a%b);
long long temp = x;x = y;y = temp - a/b*y;
}
}
int main(){
// ifstream cin("in.txt");
// ofstream out("out.txt");
long long A,B,C,k;
while(cin>>A>>B>>C>>k)
{
if(A==0&&B==0&&C==0&&k==0)
break;
long long a=C,n=1ll<<k,b=B-A;
extend_Eulid(a,n);
//printf("%d=(%d)*%d+(%d)*%d\n",q,x,a,y,b);
if(b%q==0)
{
// int e=(x%(n/q)+n/q);//我的代码出错了
long long e=(x*(b/q))%n;//太神奇了，为什么啊因为有的结果为负，所以 先加再
//取模
e=(e%(n/q)+n/q)%(n/q);
cout<<e<<endl;
// out<<e<<endl;
}
else
cout<<"FOREVER"<<endl;
// out<<"FOREVER"<<endl;
}
// out.close();
return 0;
}