poj 2115

// 题意：for (variable = A; variable != B; variable += C)
//给你一个for循环，问当循环几次时，可以跳出循环，输出循环次数
//否则，如果出现死循环，这输出FOREVER
//因此 我们可以得到 当在有限的循环中： ((B-A)+2^k)%2^k/C==x(x代表循环的次数)
//公式转化： ((B-A)+2^k)%2^k==C*X;=》C*X=(B-A)(mod 2^k); 即：我们可以得到 一个模线性方程
//令：a=C;b=B-A;n=2^k;d=gcd(a,n);b%d==0时，才有解
//1.x=(x*(b/d))%n//求出ax=b(mod n)的最小解
//2.x=(x*(b/d))%n;x=(x%(n/d)+n/d)%(n/d);求出ax=b(mod n)的最小整数解
//1.2公式：不知道是怎么推出来的 只是记住了。
/*注意：
计算n=2^k时，用位运算是最快的，1<<k （1左移k位）就是2^k
但是使用long long的同学要注意格式， 1LL<<k
使用__int64的同学要强制类型转换 (__int64)1<<k
*/
//d=ax+by=gcd(a,b);
#include<stdio.h>
__int64 ext_euc(__int64 a,__int64 b,__int64 &x,__int64 &y)
{
    if(b==0)
    {
        x=1;
        y=0;
        return a;
    }
    __int64 d=ext_euc(b,a%b,x,y);
    __int64 tem=x;
    x=y;
    y=tem-a/b*y;
    return d;
}
int main()
{
    __int64 A,B,C,K;
    while(scanf("%I64d%I64d%I64d%I64d",&A,&B,&C,&K),A+B+C+K)
    {
        __int64 x,y;
        __int64 a=C;
        __int64 b=B-A;
        __int64 n=(__int64)1<<K;
        __int64 d=ext_euc(a,n,x,y);//ax=b(mod n)
        if(b%d)//判断是不是有解。
            printf("FOREVER\n");
        else
        {
            x=(x*(b/d))%n;//求出最小解
            x=(x%(n/d)+n/d)%(n/d);//求出最小整数解。
            printf("%I64d\n",x);
        }
    }
    return 0;
}