各种不会,在网上搜了不少的资料才得以通过
补充:欧几里得算法的核心是gcd(a,b)=gcd(b,a%b)我们可以通过递归或者非递归写出代码来.关于它的证明在这里
扩展的欧几里得算法是若c=gcd(a,b),则必有一组整数想x,y使得c=ax+by,下面程序的6~13行就是该算法的代码实现,函数中的x,y就是我们求的一组x,y,q就是gcd(a,b),都可以直接引用,应用时还应该注意其他限制!
本题的数据范围比较大,要注意!
#include <iostream>
#include <fstream>
#include <cstdio>
using namespace std;
long long x,y,q;
long long extend_Eulid(long long a,long long b){//扩展欧几里得算法
if(b == 0){
x = 1;y = 0;q = a;
}else{
extend_Eulid(b,a%b);
long long temp = x;x = y;y = temp - a/b*y;
}
}
int main(){
long long A,B,C,k;
while(cin>>A>>B>>C>>k)
{
if(A==0&&B==0&&C==0&&k==0)
break;
long long a=C,n=1ll<<k,b=B-A;
extend_Eulid(a,n);
if(b%q==0)
{
// int e=(x%(n/q)+n/q);//我的代码出错了
long long e=(x*(b/q))%n;//太神奇了,为什么啊因为有的结果为负,所以 先加再
//取模
e=(e%(n/q)+n/q)%(n/q);
cout<<e<<endl;
// out<<e<<endl;
}
else
cout<<"FOREVER"<<endl;
}
return 0;
}