HDU - 2115 - Loooooops （扩展欧几里得定理）

题目链接：

http://poj.org/problem?id=2115

题解：

设kk = 2^k;
本题公式为（A + CX） % KK = B，求X，可以转化为
C*X - KK*Y = B - A，设gcd（C , KK）= D，求出D后，若（B - A） % D = 0即为有解的情况， 此时求出来的X需要进行转化， x= x* (b-a) / ans , r = kk / ans , x = (x%r + r) % r, 这样x即为满足条件的最小正整数。

AC代码：

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
using namespace std;
__int64 exgcd(__int64 a,__int64 b,__int64 &x,__int64 &y)
{
    if(b == 0)
    {
        x = 1,y = 0;
        return a;
    }
    __int64 re = exgcd(b, a%b, x, y);
    __int64 tmp = x;
    x = y;
    y = tmp - a/b*y;
    return re;
}
int main()
{
    __int64 a,b,c,x,y,k,kk,ans;
    while(cin >> a >> b >> c >> k)
    {
        if(a == 0 && b == 0 && c == 0 && k == 0)
            break;
        kk=(__int64)1<<k;
        ans = exgcd(c,kk,x,y);
        if((b-a) % ans != 0)// gcd(a,b)是ax + by 的最小正元素
            printf("FOREVER\n");
        else
        {
            x = x* (b-a) / ans;
            __int64 r = kk / ans;
            x = (x%r + r) % r;
            cout << x << endl;
        }
    }
    return 0;
}