扩展欧几里得 poj1061 青蛙的约会

中文题 题意好理解


两只青蛙跳了t步,A的坐标为x+mt  B的坐标为y+nt 相遇的充要条件
mt-y-nt=Lp
方程可列
(n-m)t+Lp=x-y
满足的最小t即为解

 

#include<iostream>
using namespace std;
typedef  __int64 ll;

ll exgcd(ll m ,ll & x ,ll n ,ll &y)
{
	ll x1,y1,x0,y0;
	x0=1;y0=0;
	x1=0;y1=1;
	ll r=(m%n+n)%n;
	ll q=(m-r)/n;
	x=0;y=1;
	while(r)
	{
		x=x0-q*x1;
		y=y0-q*y1;
		x0=x1;
		y0=y1;
		x1=x;
		y1=y;
		m=n;
		n=r;
		r=m%n;
		q=(m-r)/n;
	}
	return n;
}
/*
两只青蛙跳了t步,A的坐标为x+mt  B的坐标为y+nt 相遇的充要条件
mt-y-nt=Lp 
方程可列
(n-m)t+Lp=x-y 
满足的最小t即为解 
*/
int main()
{
	ll t,p;
	ll x,y,m,n,l;
	while(scanf("%I64d%I64d%I64d%I64d%I64d",&x,&y,&m,&n,&l)!=EOF)
	{
		ll M=exgcd(n-m,t,l,p);
		/*如果当前x-y是句初始点坐标差不能被gcd(n-m,L)的最大公约数整除
		证明除不尽则会一直在追赶,或者m==n那么就会一直在同样的距离 */
		if((x-y)%M || m==n)
		{
			printf("Impossible\n");
		}
		else
		{
			ll s=l/M;
			t=t*((x-y)/M);
			//+s在mod个s避免在第一圈的情况就遇到了,实际上是遇不到的 
			t=(t%s+s)%s;
			printf("%I64d\n",t);
		
		}
	}
	
	return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值