算法:欧几里得和拓展欧几里得(gcd和exgcd)以及一些应用

研究这种数论的可能只有偏信安方向的(虽然本人不是)吧。

欧几里得算法-即使用辗转相除法计算最大公因数的方法。gcd意思是最大公因数。算法来源于经过数学推理得到的这么一些结论:

1. gcd(a,b)=gcd(b,a%b)

2.0和一个其他的数的最大公因数就是那个其他的数

#include<iostream>
using namespace std;

//递归算法
int gcd(int a, int b)
{
	if (b == 0)return a;
	else return gcd(b, a%b);
}

//迭代算法
int gcd_recursion(int a, int b)
{
	while (b != 0)
	{
		int t = b;
		b = a%b;
		a = t;
	}
	return a;
}
int main()
{
	int a, b;
	while (cin >> a >> b)
	{
		cout << gcd(a,b)<< endl;
		cout << gcd_recursion(a, b) << endl;
	}
	return 0;
}

拓展欧几里得-即exgcd。来源于这么一个结论:不都是0的非负数a和b,和gcd(a,b)满足一个式子:存在整数xy,ax+by=gcd(a,b),这个算法研究的是:给出a和b,同时计算gcd(a,b)、x、 y这3个数的问题。

推理过程:

如果其中一个数是0 比如b=0 gcd(a,b)=a,x=1,y=0

如果都不是0:

利用之前的式子:gcd(a,b)=gcd(b,a%b) 再带入分析这2对儿不同的ab的xy的关系

ax1+by1=bx2+(a%b)y2

              =bx2+(a-(a/b)*b)y2;

ax1+by1=ay2-bx2-((a/b)*b) y2;  //这步整理有点跳跃

              =ay2=b(x2-(a/b)*y2)

这两对xy的关系为 

x1=y2;

y1=x2-(a/b)y2; 

这个代码大家可以跑一下试试。

#include<iostream>
using namespace std;

int exgcd(int a, int b, int& x, int& y)
{
	int t,gcd;
	if (b == 0)
	{
		x = 1;y = 0;return a;
	}
	gcd=exgcd(b, a%b, x, y);
	t = x;
	x = y;
	y = t - a/b*y;
	return gcd;
}

int main()
{
	int a, b,x=0,y=0;//因为没初始化会报错
	while (cin >> a >> b)
	{
		cout << exgcd(a, b, x, y) << endl;
		cout<<" "<< x << " " << y << endl;
	}
	return 0;
}

 

这个算法应用(全是信安的数学问题):

1解方程 pa+qb=c

若c mod gcd(p,q)=0   则解为:之前算法算出的解  乘上c/gcd(p,q)。否则没有解。求出一组解p0,q0 后,其他整数解满足

 p=p0+b/gcd(p,q)*t

q=q0-a/gcd(p,q)*t , t是正整数

2 解方程 a x ≡ b mod n

即 求 ax-ny =b  的解 ,其中已知a b n 计算 x y 转化为了上面的问题  有解时解的数量为gcd(a,n)

void linear_mod_equation (int a, int b, int n, int *sol)
{
    int d, x, y;
    gcd (a, n, d, x, y );
    if (b % d) d = 0;
    else
    {
        sol [0] = x * ( b / d ) % n ;
        for (int i = 1; i < d; ++i)
            sol[i] = (sol[i - 1] + n / d) % n ;
    }
}

3 求逆元 即 把上一个问题的b变成1, 解只有一个。 答案就是直接运行exgcd的代码里的算法解出的结果x






     







  


  • 0
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值