青蛙的约会 POJ - 1061 扩展欧几里得解二元一次不定方程

扩展欧几里得解二元一次不定方程

两只青蛙，起点分别为x,y ,每次可以跳m,n长度。（同方向跳）在长为 l 的圆表面，最少跳多少次相遇。

设跳k次，则(k*m + x) - (k*n + y) = p*l 假设前面的多跳 p圈。 

变为ax + by = c的形式套用扩展欧几里得模版。这里是输出x的最小值。

得到答案 *c/p得 原方程特解之一，以x为例， 如果是非负数 直接%(b/gcd(a,b)) ,负数就先%(b/gcd(a,b)) ，再加%(b/gcd(a,b))

#include <iostream>
#define ll long long
using namespace std;
ll exgcd(ll a,ll b,ll &x,ll &y) {
	if(!b) {//解方程求 x,y a*x+b*y = gcd(a,b)
		x=1,y=0;
		return a;
	}
	ll re=exgcd(b,a%b,x,y),tmp=x;
	x=y,y=tmp-(a/b)*y;
	return re;
}
int main() {
	ll x,a,b,y,c,m,n,l,p;
	cin>>x>>y>>m>>n>>l;
	a = n-m;
	b = l;
	c = x-y;
	//求 ax + by = c
	p = exgcd(a,b,x,y);// 先求 ax + by = gcd(a,b)
	if(c%p!=0) {// 无解 
		printf("Impossible\n");
	} else {
		x = x*c/p; // 这里输出 x, 解*c/p 为通解
		if(x>=0) { //求最小解
			x%=b/p;
		} else {
			x = x%b/p + b/p;
		}
		printf("%lld\n",x);
	}

	return 0;
}