poj 2115 C Looooops

来源：互联网发布：win7指纹加密软件编辑：程序博客网时间：2024/05/19 23:53

// 推广的欧几里德算法

#include <iostream>
using namespace std;

__int64 exgcd(__int64 a,__int64 b,__int64 &x,__int64 &y)
{
if (b==0)
{
x=1;y=0;
return a;
}
__int64 d=exgcd(b,a%b,x,y);
__int64 temp=x;
x=y;
y=temp-(a/b)*y;
return d;
}

int main()
{
__int64 a,b,c,k;
while (scanf("%I64d%I64d%I64d%I64d",&a,&b,&c,&k) )
{

  if (!a && !b && !c && !k) break;
  if (a==b)
  {
   printf("0/n");
   continue;
  }
  if (c==0)
  {
   printf("FOREVER/n");
   continue;
  }
  __int64 B=c;
  __int64 A=(__int64)1<<k;
  __int64 C=b-a;
  __int64 x,y;
  __int64 d=exgcd(A,B,x,y);
  if (C%d)   //无解的情况（可以查看算法导论）
  {
   printf("FOREVER/n");
   continue;
  }

y=(y*(C/d))%A; //得到第一个解

  __int64 t=A/d;
  y=(y%t+t)%t;   //C++的%运算符好像对负数的模运算会先算绝对值再把负号加上去
                 //用这个式子可以保证得到最小的正余数

printf("%I64d/n",y);
}

return 0;
}