POJ 2115 C Looooops（exgcd）

嗯...

 

题目链接：http://poj.org/problem?id=2115

 

(A+s*C)%2^k=B

(A+s*C)≡B(mod 2^k)

s*C-m*2^k=B-A

ax+by=c

有一个问题，b没必要是负的，反正正负a和b的线性组合集都一样，况且此题不需要y

 

 

AC代码：

#include<cstdio>
#include<iostream>

using namespace std;

inline void exgcd(long long a, long long b, long long &g, long long &x, long long &y){
    if(!b) { g = a; x = 1; y = 0;}
    else { exgcd(b, a % b, g, y, x); y -= x * (a / b);}
}

int main(){
    long long a, b, c, A, B, C, k, g, x, y;
    while(~scanf("%lld%lld%lld%lld", &A, &B, &C, &k)){
        if(!A && !B && !C && !k) break;
        a = C; b = 1ll << k; c = B - A;
        exgcd(a, b, g, x, y);
        if(c % g) printf("FOREVER\n");
        else{
            c /= g; b /= g;
            printf("%lld\n", (x % b * c % b + b) % b);
        }
    }
    return 0;
}