P1516 青蛙的约会

andyc_03

于 2020-10-30 23:57:16 发布

阅读量736

点赞数

CC 4.0 BY-SA版权

分类专栏：算法数论文章标签：算法

本文链接：https://blog.youkuaiyun.com/andyc_03/article/details/109395603

算法同时被 2 个专栏收录

218 篇文章

订阅专栏

数论

10 篇文章

订阅专栏

这篇博客介绍了如何利用扩展欧几里得算法(exgcd)解决形如(n-m)T+LP=X-Y的线性同余方程，并给出了一段C++代码实现。通过exgcd找到方程的最小正整数解，文章探讨了当n小于m时的处理方式以及解不存在的情况。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

(n-m)T+LP=X-Y

用exgcd去求解这个方程，然后找到最小正整数解即可

代码

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll x,y,m,n,l;
void exgcd(ll a,ll b, ll &d, ll &x,ll &y)
{
	if(!b)
	{
		x=1; y=0; d=a;
		return;
	}
	exgcd(b,a%b,d,x,y);
	ll t=x; x=y; y=t-a/b*y;
}
int main()
{
	freopen("frog.in","r",stdin);
	freopen("frog.out","w",stdout);
	ll a,b,d;
	while(scanf("%lld%lld%lld%lld%lld",&x,&y,&m,&n,&l)==5 && x!=EOF)
	{
		if(n<m) swap(n,m),swap(x,y);
		exgcd(n-m,l,d,a,b);
		if((x-y)%d!=0 || m==n)
			printf("Impossible\n");
		else printf("%lld\n",(a*(x-y)/d%(l/d)+(l/d))%(l/d));
	}
	
	return 0;
}