广东工业大学22级集训专题4.G青蛙的约会

题目

思路

设speed（速度差），t（步数），l（圈长），n（圈数）,dis（初始距离）

由情景可以列出公式t*speed+dis=n*l ;

变换得到t*speed-n*l=-dis；此时，不妨设speed=n-m，dis=-dis=x-y；

我们就需要求解公式 at+bn=c，由裴蜀定理，该公式c是gcd（a，b）的倍数时一定有解，对于c不是这种情况时，尝试进行变化，即令g=gcd（a，b）；at+bn=g有解，

那么at*c/g+bn*c/g=c同样有解

这时能用扩展欧几里得算法得到一个解，但是要注意得到的这个解不是我们真正想要的，我们要的是t*c/g;

那么如何由这个解得到最小正整数解呢

我们可以证明a*c/g*(t+kb*c/g)+b*c/g*(n-ka*c/g)与上式等价；换句话说，我们得到的解可以减去若干个数，即t%s，要使t最小，则s应当尽量小，可见当s=b/g时有最小解；

那么我们要得到最小正整数解就应当（ans%b/g+b/g）%b/g；

AC代码

#include <iostream>
using namespace std;
long long gcd;
long long arr[2];
void exgcd(long long a,long long b)
{
    if(b==0)
    {
        arr[0]=1,arr[1]=0,gcd=a;
        return ;
    }
    exgcd(b,a%b);
    long long tem=arr[0];
    arr[0]=arr[1];
    arr[1]=tem-(a/b)*arr[1];
}

int main()
{
    long long x,y,m,n,l;
    scanf("%lld%lld%lld%lld%lld",&x,&y,&m,&n,&l);
    long long speed=n-m,dis=x-y;
    if(speed<0) dis=-dis,speed=-speed;
    exgcd(speed,l);
    if(dis%gcd) printf("Impossible");
    else
    {
        printf("%lld",((arr[0]*dis/gcd)%(l/gcd)+l/gcd)%(l/gcd));
    }
    return 0;
}

我认为这是道不错的数学题，很多细节很考验思维