【拓展欧几里得】机器人的相遇

#include<iostream>
#include<bits/stdc++.h>
using namespace std;

long long gcd(long long a, long long b)
{
	if(b==0)
		return a;
	return gcd(b,a%b);
}

long long gcd1(long long a,long long b,long long &x,long long &y)
{
	if(b==0)
	{
        x=1;
        y=0;
        return a;
    }
	long long r=gcd1(b,a%b,x,y);
	long long t=x;
	x=y;
	y=t-a/b*y;
	return r;
}

int main()
{
	long long x,y,m,n,L;
	cin>>x>>y>>m>>n>>L;
	long long a=n-m;
	long long b=L;
	long long c=x-y;
	long long d=gcd(a,b);
	if(c%d!=0)
	{
		cout<<"Impossible";
		return 0;
	}
//	a/=d;b/=d;
	int beishu=c/d;
	long long t,k; //t为步数，k为圈数，分别为a b的系数。 
	gcd1(a,b,t,k);
	t=t*beishu;    //因为其实是满足 ta+kb=d的，我需要=c的t，所以t乘上c/d
//	t=((t%b)+b)%b;
//	if(!t) t+=b;
    t=(t%b+b)%b;
    cout<<t;
	return 0;
}

需要注意的点：

一、欧几里得其实都会写，但是拓展欧几里得就要注意了：一个是形参的系数x y要用引用类型，二个是注意递归到b==0处要设置x=1,y=0；

二、关于本题思路上的突破（神奇的欧几里得）：

其中，a=n-m，b=l，然后x-y也是已知量，想要求的未知量是系数t和k （答案要的是步数t），其实就是拓展欧几里得的公式摆在眼前了。目前就要抓住拓展欧几里得的条件：ax+by=d，这个d必须是a和b的最大公约数gcd。那么，在本式中，只要x-y是n-m和l的最大公约数的倍数即可，因为这样所想要的系数也只是乘了这个倍数而已。

三、关于最后的答案操作。

    t=(t%b+b)%b;

相当于固定操作了，可以满足t得到的是正数。记忆：“模加模b”。