P1516 青蛙的约会

题目传送门

思路

因为两个青蛙同时跳到同一个点上才算碰面,设 t t t 为跳的次数, p p p 为两个青蛙跳的圈数之差,有如下式子:

( x + m × t ) − ( y + n × t ) = p × L (x+m \times t ) - ( y+n \times t ) = p \times L (x+m×t)(y+n×t)=p×L

整理得:

( n − m ) × t + L × p = x − y (n-m) \times t + L \times p = x - y (nm)×t+L×p=xy

首先,要判断 gcd ⁡ ( n − m , L ) ∤ x − y \gcd ( n-m , L ) \nmid x-y gcd(nm,L)xy 的情况。

然后,直接套拓展欧几里得算法,求出 ( n − m ) × t + L × p = gcd ⁡ ( n − m , L ) (n-m) \times t + L \times p = \gcd ( n-m , L ) (nm)×t+L×p=gcd(nm,L) 一组解 t 0  ,  p 0  ,  g c d = gcd ⁡ ( n − m , L ) t_0 \text{ , } p_0 \text{ , } gcd= \gcd ( n-m , L ) t0 , p0 , gcd=gcd(nm,L)

特解:求 x x x 的最小正整数解,即 t = b / gcd ⁡ ( a , b )  ,  x = ( x   m o d   t + t )   m o d   t t = b / \gcd(a,b) \text{ , } x= ( x \bmod t + t ) \bmod t t=b/gcd(a,b) , x=(xmodt+t)modt

带入特解(本题求的是 t t t 的最小解),求出 ( n − m ) × t + L × p = gcd ⁡ ( n − m , L ) (n-m) \times t + L \times p = \gcd ( n-m , L ) (nm)×t+L×p=gcd(nm,L) t 0 t_0 t0 的最小解 t ′ t' t

t ′ = ( t 0   m o d   L g c d + L g c d )   m o d   L g c d t'= ( t_0 \bmod \dfrac{L}{gcd} + \dfrac{L}{gcd} ) \bmod \dfrac{L}{gcd} t=(t0modgcdL+gcdL)modgcdL

t = t ′ × x − y g c d t = t' \times \dfrac{x - y}{gcd} t=t×gcdxy

所以

t = ( x − y g c d × t 0   m o d   L g c d + L g c d )   m o d   L g c d t=(\dfrac{x - y}{gcd} \times t_0 \bmod \dfrac{L}{gcd} + \dfrac{L}{gcd} ) \bmod \dfrac{L}{gcd} t=(gcdxy×t0modgcdL+gcdL)modgcdL

C o d e Code Code

#define int long long
signed main()
{
	int x,y,m,n,l;
	int X,Y,gcd;
	cin>>x>>y>>m>>n>>l;
	int A=x-y,B=n-m;
	if(B<0)//注意取反!!!
	{
		B=-B;
		A=-A;
	}
	if(A%(gcd=exgcd(B,l,X,Y))!=0)
		cout<<"Impossible"<<endl;
	else
		cout<<(  A/gcd * X % (l/gcd) + (l/gcd)  ) % (l/gcd)<<endl;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值