POJ 2115 C Looooops

~~~题目链接~~~

题目大意：根据题中的c语言for循环表达式， 判断这个循环要几次结束， 或永远不结束， 当内存达到上限时又从0开始。

思路： 与POJ 1061一样， 根据表达式(a+c*x)%1<<k = b%1<<k  推出  c*x+(1<<k)*y = b-a;

Ps:要用(__int64)1<<k ，估计默认情况下1为带符号的32位

code:

#include <stdio.h>
#include <math.h>
void gcd(__int64 s1, __int64 s2, __int64 &d, __int64 &x, __int64 &y)
{
    if(!s2){d = s1;x = 1;y = 0;}
    else {gcd(s2, s1%s2, d, y, x);y -= x*(s1/s2);}
}
int main()
{
    __int64  a = 0, b = 0, c = 0, d = 0, k = 0, s1 = 0, s2 = 0, x = 0, y = 0, r = 0;
    while(scanf("%I64d %I64d %I64d %I64d", &a, &b, &c, &k) , a+b+c+k != 0)
    {
        gcd(c, (__int64)1<<k, d, x, y);
        if((b-a)%d)
            printf("FOREVER\n");
        else
        {
            x *= (b-a)/d;
            r = ((__int64)1<<k)/d;
            printf("%I64d\n", (x%r+r)%r);
        }
    }
    return 0;
}