扩展欧几里得入门

最新推荐文章于 2021-03-14 06:43:19 发布

林木森.

最新推荐文章于 2021-03-14 06:43:19 发布

阅读量115

点赞数

分类专栏：基础数论

本文链接：https://blog.csdn.net/qq_42937891/article/details/100925622

版权

基础数论专栏收录该内容

5 篇文章 0 订阅

订阅专栏

青蛙的约会

经典扩欧板子题，稍微注意一点就可以A了，之所以要写呢，是因为发现这题在洛谷上的数据比较强，让人自闭。
方程大家应该都找的到吧，我也就不提了，直接说坑点，就是如果gcd之中如果有一个参数为负数的话，最后的结果符号会出现问题，不确定是否为正还是负，这一题里面很明显的出现了，如果（m-n）为负数的话，那么gcd可能为负数，那么ans= (ans % f + f) % f;（f=l/gcd；一定为负数，明显错了，那么x,y也受到了影响，会出现错误，所以我们在求之前先保证（m-n）为正数就可以了！

#include<string>
#include<iostream>
#include<algorithm>
#include<map>
#include<stdio.h>
using namespace std;
typedef long long ll;
ll x, y;
ll exgcd(ll a, ll b)
{
	if (b == 0)
	{
		x = 1;
		y = 0;
		return a;
	}
	ll g = exgcd(b, a % b);
	ll temp = x;
	x = y;
	y = temp -  (a / b) * y;
	return g;
}
int main()
{
	ll a , b , n , m , l;
	cin >> a >> b >> n >> m >> l;
	ll p = (a - b), q = (m - n);
	if (q < 0) q = -q, p = -p;//转化一下成正数
	ll gcd = exgcd(q, l);
	if (p % gcd != 0)
	{
		printf("Impossible\n");
	}
	else
	{
		ll f = l / gcd, ans = ((x * p / gcd));
		ans= (ans % f + f) % f;
		printf("%lld\n", ans);
	}
}