「一本通 6.4 例 1」 luogu P1516 青蛙的约会

analysis

首先，模拟不能过
那么就是一道数学：
$设青蛙a从位置A开始，步长为m,青蛙b从位置B开始,步长为n$
$则a的位置:A+km,b的位置:B+kn,其中,k为跳跃步数$
$若a,b相遇,则有同余方程:A+km\equiv B+kn(\mathrm{m}\mathrm{o}\mathrm{d}L)$
$这个等价于A+km-(B+kn)=xL,x\in Z$
$整理方程得到:k(m-n)-xL=B-A$
$由于这是个丢番图方程(不定方程),x,k有多个解,那么不妨设x>0,则:$
$k(n-m)+xL=A-B$
$此时将(n-m)看成a,L看成b,A-B看成c$
$那么这就是解一个不定方程ak+by=c中k的最小正整数解了$
$当然,不定方程有解的判断条件:(a,b)|c$
$若(n-m,L)|(A-B),我们就可以用扩展欧几里得算出k(n-m)+xL=(n-m,L)的一组解了$
注意,这里的各项系数a,b要大于等于0,否则欧几里得定理就不成立了
若a小于0,反号就好,b本来就是可以通过x的取值做到大于0的
$求出一个特解k_0后,k(n-m)+xL=(n-m,L)的通解k_i=k_0-\phi \frac{b}{(a,b)},\phi \in Z$
(证明在文后)
$那么,对于原方程,通解的系数要乘上一个\frac{A-B}{(a,b)}$
$对于最小解,就是通解中最靠近0的那一个,处理方式见代码$

code

#include <bits/stdc++.h>
using namespace std;
#define loop(i, start, end) for (register int i = start; i <= end; ++i)
#define anti_loop(i, start, end) for (register int i = start; i >= end; --i)
#define clean(arry, num) memset(arry, num, sizeof(arry))
#define ll long long
#define isdegit(a) ((a >= '0' && a <= '9'))
template <typename T>
void read(T &x) {
    x = 0;
    char r = getchar();
    T neg = 1;
    while (!isdegit(r)) {
        if (r == '-')
            neg = -1;
        r = getchar();
    }
    while (isdegit(r)) {
        x = (x << 1) + (x << 3) + r - '0';
        r = getchar();
    }
    x *= neg;
}

ll exgcd(ll a, ll b, ll &x, ll &y) {
    if (b == 0) {
        x = 1;
        y = 0;
        return a;
    }
    ll res = exgcd(b, a % b, x, y);
    ll x1 = x;
    x = y;
    y = x1 - (a / b) * y;
    return res;
}

int main() {
#ifndef ONLINE_JUDGE
    freopen("datain.txt", "r", stdin);
#endif
    ll A, n, B, m, L;
    read(A);
    read(B);
    read(m);
    read(n);
    read(L);
    ll x = 0, y = 0;
    if (n < m)
        swap(n, m), swap(A, B);
    ll b = L, a = n - m, c = A - B;
    ll gcd = exgcd(a, b, x, y);
    if (c % gcd || m == n) {
        printf("Impossible\n");
    } else {
        printf("%lld\n", ((x * c / gcd) % (b / gcd) + b / gcd) % (b / gcd));
    }
    return 0;
}

相关证明

• 关于通解
  对于方程:$ax+by=(a,b)$
  $假设已经求出一组解x_0,y_0,且其他解设为x_i,y_i$
  $则a{x}_0+b{y}_0=(a,b)$
  $a{x}_i+b{y}_i=(a,b)$
  $则a{x}_0+b{y}_0=a{x}_i+b{y}_i$
  $同时除以(a,b)得:\frac{a}{(a,b)}(x_0-x_i)=\frac{b}{(a,b)}(y_i-y_0)$
  $即\frac{a}{(a,b)}|\frac{b}{(a,b)}(y_i-y_0)$
  $由(\frac{a}{(a,b)},\frac{b}{(a,b)})=1,有(x_0-x_1)|\frac{b}{(a,b)}$
  $等价变换后得x_i=x_0-\phi \frac{b}{(a,b)},\phi \in Z$