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;}         elsecout<<"FOREVER"<<endl;}    return 0;}