POJ ~ 1061 ~ 青蛙的约会（扩展欧几里得）

最新推荐文章于 2022-04-14 00:39:57 发布

张松超

最新推荐文章于 2022-04-14 00:39:57 发布

阅读量243

点赞数

分类专栏：【素数/同余/不定方程/乘性函数】

本文链接：https://blog.csdn.net/ZscDst/article/details/82917939

版权

【素数/同余/不定方程/乘性函数】专栏收录该内容

10 篇文章 0 订阅

订阅专栏

题解

假设答案为a，其实就是求解： $x+ma \equiv y+na (mod\quad L)$ ，化为 $(m-n)a + Lk = y-x$ 。

对应到 $ax+by=c$ 中，a = m-n，b = L， c = y-x。x为a，y为k。要求最小的非负整数x。

假设 $ax+by=c$ 的一组解为(x0，y0)，那么通解为 $(x0+b't,y0-a't)\qquad a'=a/\gcd \quad b'=b/\gcd$

所以最小非负解为 $(x_0\%b'+b')\%b'$ 。

//#include<bits/stdc++.h>
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
using namespace std;
typedef long long LL;
void exgcd(LL a, LL b, LL& d, LL& x, LL& y)
{
    if (!b) { d=a; x=1; y=0; }
    else { exgcd(b, a%b, d, y, x); y -= x*(a/b); }
}
LL cal(LL a, LL b, LL c)
{
    LL GCD, x0, y0;
    exgcd(a, b, GCD, x0, y0);
    if (c%GCD) return -1;
    LL bb = b/GCD;
    bb = llabs(bb);//GCD可能为负数，需要把bb变为正
    x0 *= c/GCD; y0 *= c/GCD;
    return (x0%bb+bb)%bb;
}
int main()
{
    LL x, y, m, n, L; scanf("%lld%lld%lld%lld%lld", &x, &y, &m, &n, &L);
    LL ans = cal(m-n, L, y-x); /// x+ma ≡ y+na (mod L) => (m-n)a + L*k = y-x;
    if (ans != -1) printf("%lld\n", ans);
    else printf("Impossible\n");
    return 0;
}
/*
1 2 3 4 5
*/