【扩展GCD】

最新推荐文章于 2022-10-19 23:23:23 发布

dawxy

最新推荐文章于 2022-10-19 23:23:23 发布

阅读量487

点赞数

分类专栏：数论模版文章标签：扩展GCD

本文链接：https://blog.csdn.net/u012760629/article/details/38040123

版权

模版同时被 2 个专栏收录

42 篇文章 0 订阅

订阅专栏

数论

8 篇文章 0 订阅

订阅专栏

题目链接 HDU2669

【题意】求出X*a + Y*b = 1的最小的非负整数x时的一个解 X，Y；如果不存在解输出sorry；

【分析】扩展GCD模版，紫书上的蛮简洁的。

【AC CODE】0ms

#include <cstdio>
#include <cstring>
#include <cstdlib>
typedef long long LL;

void e_gcd(LL a, LL b, LL& d, LL& x, LL& y)
{
	if (!b) d = a, x = 1, y = 0;
	else e_gcd(b, a%b, d, y, x), y -= x*(a / b);
}
int main()
{
#ifdef SHY
	freopen("e:\\1.txt", "r", stdin);
#endif
	LL a, b, d, x, y;
	while (~scanf("%I64d %I64d%*c", &a, &b))
	{
		e_gcd(a, b, d, x, y);
		LL c = 1;
		if (0 == c%d)
		{
			x *= c / d; //求出x0后要乘上c/d倍,y也一样:y *= c/d;
			/*此时的x,y才是其中一个解,下面可以用通项公式推出所有解*/
			LL gx = abs(b/d);//通项公式:(x0+k*gx,y0-k*gy)；gx = b/gcd(a,b),gy = a/gcd(a,b)互素,k为任意整数
			x = (x%gx+gx)%gx;//求出>=0的最小的x,如果要求>=0的最小的y则只需要y(y%gy+gy)%gy;
			/*也可以下面这样循环求，效率低
			if(gx < 0) gx = -gx, gy = - gy;//防止负数，但是这里都是整数其实可以去掉
			while (x > 0) x -= gx, y += gy;
			while (x < 0) x += gx, y -= gy;
			printf("%I64d %I64d\n", x, y);
			*/
			printf("%I64d %I64d\n", x, (c-a*x)/b);
		}
		else puts("sorry");
	}
	return 0;
}