poj-2115

各种不会,在网上搜了不少的资料才得以通过

补充:欧几里得算法的核心是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;
}