拓展欧几里得-青蛙得约会

最新推荐文章于 2022-11-26 11:15:29 发布

原创最新推荐文章于 2022-11-26 11:15:29 发布 · 465 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#拓展欧几里得

acm 同时被 2 个专栏收录

178 篇文章

订阅专栏

基本算法

166 篇文章

订阅专栏

#include <cstdio>
#include <iostream>
#include <cstring>
using namespace std;
long long x,y,m,n,l,a,b,p;

long long exgcd(long long a,long long b,long long &xx,long long &yy){
    if(b == 0){
        xx = 1;
        yy = 0;
        return a;
    }
    long long ans = exgcd(b,a%b,xx,yy);
    long long tmp = xx;
    xx = yy; yy = tmp - a/b*yy;
    return ans;
}

void init(){
    a = m-n;
    b = l;
    p = y-x;
    long long xx,yy;
    long long ans = exgcd(a,b,xx,yy);
    long long kb = p/ans;
    if(p % ans == 0){
       xx *= kb;
       b /= ans;
       if(b < 0) b = -b;
       long long nn = xx % b;
       if(nn <= 0)  nn+= b;
       printf("%lld\n",nn);
    }
    else
        printf("Impossible\n");
}

int main(){
    while(~scanf("%lld%lld%lld%lld%lld",&x,&y,&m,&n,&l)){
        init();
    }

}