码迷,mamicode.com
首页 > 其他好文 > 详细

poj 1061 扩展欧几里德同余方程求解

时间:2015-10-06 11:39:21      阅读:113      评论:0      收藏:0      [点我收藏+]

标签:

摘要写在一瞪眼。

#include<iostream>
using namespace std;
long long exgcd(long long a,long long b,long long &k,long long &t)
{
    if (b==0)
    {
        k=1;
        t=0;
        return a;

    }
    else
    {
        long long tp_gcd;
        tp_gcd=exgcd(b,a%b,k,t);
        long long temp;
        temp=k;
        k=t;
        t=temp-(a/b)*t;
        return tp_gcd;
    }
}
// (n-m)*t-k*l=x-y
// --> a*x-b*y=c     a=n-m,b=l;
// -->使用扩展欧几里得能够求得a*x0-b*y0=gcd(a,b),d=gcd(a,b) 的解x0,y0
// x0*(l/d)就得到了第2行中,x的一个特解x‘
// x=x‘+k*(l/d)
int main()
{
    long long x,y,n,m,l;
    long long x0,y0;
    while(cin>>x>>y>>m>>n>>l)
    {
        long long d=exgcd(n-m,l,x0,y0); //d是gcd(a,b)
        if((x-y)%d!=0) cout<<"Impossible"<<endl;
        else
        {
            x0*=(x-y)/d;
            x0=(x0%(l/d)+(l/d))%(l/d);
            cout<<x0<<endl;
        }
    }
    return 0;
}


版权声明:本文博主原创文章,博客,未经同意不得转载。

poj 1061 扩展欧几里德同余方程求解

标签:

原文地址:http://www.cnblogs.com/lcchuguo/p/4856800.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!