拓展欧几里德 算是套公式的题目吧,反正开始很sb,没用脑子
下面的代码中 前面很好理解(前提是看懂了拓展欧几里得) 最后面 求 x ,得是 最小正整数解, 那两行也是套路吧,不是很理解,好几个提都得这样
ps:欧几里得不是很懂的 建议先搞明白,然后附带这逆元也会了
/*
首先反思,垃圾水题没思考就做
正确题意是 扩展欧几里得
由原式:(b-a) + i * (1<<k) % c == 0 推出
本结论是 (b-a) | gcd(c, 1<<k) ,,, 若存在,公式中最后的x 就是 最小的i
*/
#include<iostream>
#include<cstdio>
#include<cstring>
#include<math.h>
using namespace std;
const int maxn = 100000 + 10;
const int INF = 1e9 +7;
typedef long long LL;
LL a, b, c, k;
void find_(LL a, LL b, LL &d, LL &x, LL &y) {
if(!b) { d = a; x = 1; y = 0; return; }
find_(b, a%b, d, x, y);
LL t = x;
x = y;
y = t - (a/b) * y; // 详细写出精华部分 == 深刻理解 + 记忆
}
void solve() {
LL b_ = (1LL << k);
LL d;
LL x, y;
//find_(6, 15, d, x, y);
//cout << x << "===" << y << endl;
find_( c, b_, d, x, y);
if((b-a) % d) { cout << "FOREVER" << endl; return; }
else {
x *= ((b-a) / d);
x = (x%(b_/d) + b_/d) % (b_/d); // x%(b_/d) + b_/d ==>> 达到最小公倍数需要的 c 的个数
cout << x << endl;
return;
}
}
int main() {
while(~scanf("%lld%lld%lld%lld", &a, &b, &c, &k)) {
if(!a && !b && !c && !k) break;
solve();
}
return 0;
}