模线性方程,水题
#include <iostream>
#include <stdio.h>
using namespace std;
typedef long long ll;
ll exgcd(ll a,ll b,ll &x,ll &y)
{
if(b==0){x=1;y=0;return a;}
ll d=exgcd(b,a%b,x,y);
ll t=x;
x=y;
y=t-(a/b)*y;
return d;
}
void modeq(ll a,ll b,ll n)
{
ll e,i,d,x,y;
d=exgcd(a,n,x,y);
if(b%d>0) printf("FOREVER\n");
else
{
e=b/d*x;
n=n/d;
e=(e%n+n)%n;//防止负数
printf("%lld\n",e);
}
}
int main()
{
ll a,b,c,k,s;
freopen("in.txt","r",stdin);
while(scanf("%lld%lld%lld%lld",&a,&b,&c,&k)!=EOF)
{
if(a==0 && b==0 && c==0 && k==0)
break;
s=1;
b=b-a;
for(int i=1;i<=k;i++)
s*=2;
if(c<0){c=-c,b=-b;}
if(b<0) {b%=s,b+=s;}
//cout<<c<<" "<<b<<" "<<s<<endl;
modeq(c,b,s);
}
return 0;
}