POJ2115解题报告【拓展欧几里得模板题】

题目地址：

　　http://poj.org/problem?id=2115

题目概述：

　　求解线性同余方程：

大致思路：

　　这是个模板题，用拓展欧几里得即可解决。

代码：

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <vector>
#include <ctime>
#include <map>
#include <stack>
#include <set>
#include <queue>
#include <cstring>
#include <algorithm>
using namespace std;

#define sacnf scanf
#define scnaf scanf
#define maxn 47780
#define maxm 35
#define inf 1061109567
#define Eps 0.000001
const double PI=acos(-1.0);
#define mod 1000000007
#define MAXNUM 10000
void Swap(int &a,int &b) {int t=a;a=b;b=t;}
int Abs(int x) {return (x<0)?-x:x;}
typedef long long ll;

ll P[maxm];

ll ex_gcd(ll a,ll b,ll &x,ll &y)
{
    if(b==0) {x=1;y=0;return a;}
    ll ans=ex_gcd(b,a%b,x,y);
    ll temp=x;x=y;y=temp-(a/b)*y;
    return ans;
}

int main()
{
    //freopen("data.in","r",stdin);
    //freopen("data.out","w",stdout);
    //clock_t st=clock();
    ll a,b,c,m,x,y;int k;
    P[0]=1;for(int i=1;i<=32;i++) P[i]=P[i-1]*2;
    while(~scanf("%lld%lld%lld",&a,&b,&c))
    {
        scanf("%d",&k);
        if(a==0&&b==0&&c==0&&c==0) break;
        m=P[k];b=(b-a+m)%m;ll p=ex_gcd(c,m,x,y);
        if(b%p!=0) printf("FOREVER
");
        else
        {
            ll q=m/p;
            ll ans=(b/p*x%q+q)%q;
            printf("%lld
",ans);
        }
    }
    //clock_t ed=clock();
    //printf("

Time Used : %.5lf Ms.
",(double)(ed-st)/CLOCKS_PER_SEC);
    return 0;
}