poj 1061 青蛙的约会

最新推荐文章于 2021-10-06 17:26:28 发布

转载最新推荐文章于 2021-10-06 17:26:28 发布 · 46 阅读

0 ·

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/overflow/p/3187292.html

本文通过一个具体的POJ题目实例，详细介绍了如何运用扩展欧几里得算法解决特定类型的数学问题。文章首先给出了完整的代码实现，并对其背后的数学原理进行了深入浅出的解释。

题目：http://poj.org/problem?id=1061

思路：扩展欧几里得

#include <iostream>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <cstdio>
using namespace std;
long long gcd(long long a,long long b)
{
    if(b==0)
        return a;
    return gcd(b,a%b);
}
long long exgcd(long long a,long long b,long long &x,long long &y)
{
    if(b==0)
    {
        x=1;
        y=0;
        return a;
    }
    else
    {
        long long ans=exgcd(b,a%b,x,y);
        long long t=x;
        x=y;
        y=t-a/b*y;
        return ans;
    }
}
int main()
{
    long long x,y,m,n,l;
    ///  (x+mt) - (y+nt) = kl
    // (x-y) + (m-n)t = kl
    //  (m-n)t - kl = y-x
    cin>>x>>y>>m>>n>>l;
    long long xx,yy;
    long long A=-(n-m),B=-l,C=-x+y;
    long long g=gcd(A,B);
    if(C%g!=0)
        cout<<"Impossible\n";
    else
    {
        A/=g;
        B/=g;
        C/=g;
        exgcd(A,B,xx,yy);
        xx*=C;
        xx=(xx%B+B)%B;
        cout<<xx<<endl;
    }
    return 0;
}