POJ - 2115 C Looooops （扩展欧几里得）

拓展欧几里德  算是套公式的题目吧，反正开始很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;
}