POJ 2115

思路：构造模线性方程，然后用扩欧求解。

方程：（x*c+a）%2^k=b 所以 x*c-y*2^k=b-a ，套ax+by=c，然后扩欧。

注意：我开始写b=-2^k然后WA成狗，后来改成2^k AC了 。b没必要是负的，反正正负a和b的线性组合集都一样，况且此题不需要y。因为是求最小正整数解，所以是%t/d。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<string>
#include<vector>
#include<stack>
#include<bitset>
#include<cstdlib>
#include<cmath>
#include<set>
#include<list>
#include<deque>
#include<map>
#include<queue>
using namespace std;
#define inf 0x3f3f3f3f
#define ll long long
#define fo freopen("in.txt","r",stdin)
#define fc fclose(stdin)
#define fu0(i,n) for(i=0;i<n;i++)
#define fu1(i,n) for(i=1;i<=n;i++)
#define fd0(i,n) for(i=n-1;i>=0;i--)
#define fd1(i,n) for(i=n;i>0;i--)
#define mst(a,b) memset(a,b,sizeof(a))
#define sd(n) scanf("%d",&n)
#define sdd(n,m) scanf("%d %d",&n,&m)
#define ss(s) scanf("%s",s)
#define sddd(n,m,k) scanf("%d %d %d",&n,&m,&k)
#define pans(ans) prllf("%d\n",ans)
#define all(a) a.begin(),a.end()
#define sc(c) scanf("%c",&c)
#define we(a) while(scanf("%d",&a)!=EOF)
const ll maxn=200005;
const double eps=1e-8;
const double PI = acos(-1.0);
ll gcd(ll a,ll b)
{
    return b==0?a:gcd(b,a%b);
}
void kgcd(ll a,ll b,ll &d,ll &x,ll &y)
{
    if(!b)
    {
        d=a;
        x=1;
        y=0;
    }
    else
    {
        kgcd(b,a%b,d,y,x);
        y-=x*(a/b);
    }
}
int main()
{
    ll a,b,c,k;
    while(cin>>a>>b>>c>>k)
    {
        if(!a&&!b&&!c&&!k)
            break;
        ll x,y,t=(1LL<<k);
        ll d;
        kgcd(c,t,d,x,y);
        if(d==0)
        {
            cout<<0<<endl;
            continue;
        }
        ll mod =t/d;
        if((b-a)%d)
            cout<<"FOREVER"<<endl;
        else
        {
            x=(((b-a)/d)*x%(mod)+(mod))%mod;
            cout<<x<<endl;
        }

    }
    return 0;
}