POJ 1061 青蛙的约会 (欧几里德+同余方程)

最新推荐文章于 2019-08-01 11:45:13 发布

原创最新推荐文章于 2019-08-01 11:45:13 发布 · 183 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

收录于

POJ 同时被 2 个专栏收录

23 篇文章

订阅专栏

数学基础

10 篇文章

订阅专栏

本文深入探讨了扩展欧几里得算法的原理及应用，通过实例代码讲解了如何求解最大公约数和贝祖等式的解。适用于数学算法、计算机科学领域的学习者。

#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<queue>
#include<math.h>
#include<iostream>
using namespace std;
void Ex_gcd(long long a, long long b, long long &x, long long &y)
{
    if(b == 0)//µÝ¹é³ö¿Ú
    {
        x = 1;
        y = 0;
        return;
    }
    long long x1, y1;
    Ex_gcd(b, a%b, x1, y1);
    x = y1;
    y = x1-(a/b)*y1;
}
long long gcd(long long a,long long b)
{
    return b?gcd(b,a%b):a;
}
void jiaohuan(long long &a,long long &b)
{
    long long c=a;
    a=b;
    b=c;
}
int main()
{

  long long x,y,m,n,L,a,b,c,g;
  scanf("%lld%lld%lld%lld%lld",&x,&y,&m,&n,&L);
  if(m<n)
  {
      jiaohuan(x,y);
      jiaohuan(m,n);
  }
  a=m-n;
  b=L;
  c=y-x;///是y-x  不是x-y
  if(a==0)
  {
      printf("Impossible\n");
      return 0;
  }
  g=gcd(a,b);
 // printf("%lld\n",g);
   if(c%g!=0)
   {
       printf("Impossible\n");
   }
   else
   {
      c=c/g;
      x=y=0;
       Ex_gcd(a,b,x,y);
       b=b/g;
        x=x*c;
        x=(x%b+b)%b;
        printf("%lld\n",x);
   }

    return 0;
}