扩展欧几里得:
1.给任意a,b,c,,,求ax+by=c,这里的c必须是gcd(x,y)的整数倍才会又整数解,否则没有整数解.
(这里的a,b,c可以是个小式子,可能为负数)
2.利用扩展欧几里得求ax+by=gcd(a,b)
(这里的a和b可以是一个小式子,可能为负数)
可以得到一组(x0,y0)解,那么就可以得到所有的解
x=x0+k(b/gcd(a,b))
y=y0-k(a/gcd(a,b)) 这里的k是任意正整数。
(若是ax-by,那么这里的y=y0+k(a/gcd(a,b))
3.一般用来求ax≡p( mod b )等价于ax+by=p,如果p不是gcd(a,b)的倍数,那么无解,否则,求出ax+by=gcd(a,b)然后两边同时乘以p/gcd(a,b)即可。(这里的a,b,p可以是个小式子,可能为负数)
题目:https://www.acwing.com/problem/content/224/
题意:俩个青蛙在一条长l的环上跳动,有一个0点坐标,然后a,b青蛙分别有个a,b初始坐标,然后a每次跳的距离是m,b每次跳的距离是n,俩只青蛙每次跳动的时间一样,问:俩只青蛙能否跳到同一个点,若能跳到同一个点,那么最少跳的次数。
题解:设跳动次数x次,那么可以化简得出公式(m-n)x+ly=b-a,这里的b-a,m-n,gcd(m-n,l)都可能为负数,然后套用扩展欧几里得
若不太理解可以先做稍微简单的一道:https://blog.csdn.net/Are_you_ready/article/details/112006948
#include <iostream>
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e6+9;
ll exgcd(ll a,ll b,ll &x,ll &y)
{
if(!b)
{
x=1,y=0;
return a;
}
ll d=exgcd(b,a%b,y,x);
y-=a/b*x;
return d;
}
int main()
{
ll a,b,m,n,l;
cin>>a>>b>>m>>n>>l;
ll x,y;
ll d=exgcd(m-n,l,x,y); //套用扩展欧几里得
if((b-a)%d)printf("Impossible\n");//如果不是倍数,那么没有正整数解
else
{
x*=(b-a)/d;//套用完以后,先让x乘以倍数,然后再进行求其他坐标
ll t=abs(l/d);//这里l/d可能为负数,取个绝对值,然后求最小的正整数x
printf("%lld\n",(x%t+t)%t);
}
return 0;
}