浅谈唯一分解定理

最新推荐文章于 2022-07-20 10:55:27 发布

原创最新推荐文章于 2022-07-20 10:55:27 发布 · 2.5k 阅读

12 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#gcd #算法导论 #唯一分解定义 #算法

本文深入解析唯一分解定理，即算术基本定理，阐述一个大于一的正整数可唯一分解为质数乘积的形式。文章涵盖定理定义、证明方法，以及在求因子个数、因子之和、最大公约数和最小公倍数中的应用。

原文链接

唯一分解定理是数论里面非常重要的一部分，必须要熟练掌握💪

一、定义

唯一分解定理又称算术基本定理，指：一个大于一的正整数N都可以唯一分解成有限个质数的乘积， N=p1^a1 * p2^a2 * p3^a3 * … * pn^an，这里p1< p2< p3 <…< pn均为质数，ai均为正整数.这样的式子成为N的标准分解式

二、证明

一个数肯定是由合数和质数构成的，合数又可以分解成质数和合数，最后递归下去就会变成质数的乘积
例如： 364 * 9||6 * 62 * 3 * 2 * 3==2² * 3²
详细可以看具体严谨证明

三、应用

1.求出数n的因子个数

个数==（1+a1）* (1+a2) * (1+a3) * (1+a4)…*(1+an) 这里的ai表示N的标准分解式的幂次数
因为当我的a1=3时那我n的因子肯定会有 p1⁰ p1¹ p1² p1³ 这四个数
然后再和p2的个数搭配起来就是两个数的因子数相乘了 p1^x 可以与 p2^y 随意搭配，所以进行乘法
就比如上面的36：
在这里插入图片描述
图上很容易看出来一个数的因子一定可以表示成其质因子的乘积形式，所以就有了第二个应用

2.求所有的因子之和

还是上面的图，所有因子和就是上面的9个数的和，而这9个数是由左面的两竖列两两相乘得到的，所以就有：
1+3+9+2+6+18+4+12+36==(1+2+4) * (1+3+9)==(2⁰+2¹+2³) * (3⁰+3¹+3²)
公式:（q1⁰+q1¹+q1²…q1^a1）* (q2⁰+q2¹+q2²…q2^a2) * …*（qn⁰+qnⁿ+qn²…qn^an）

3.求gcd和lcm

给你两个数如何用唯一分解定理求其gcd和lcm，假设两个数是a和b则：
ap1^a1 * p2^a2 * p3^a3 * … * pn^an
bp1^b1 * p2^b2 * p3^b3 * … * pn^bn
两者
gcdp1^min(a1,b1) * p2^min(a2,b2) * p3^min(a3,b3) * … * pn^min(an,bn)
反之求lcm就是把上面得min全部换成max
例如：
42^2
6==2*3
gcd(4,6)==2^min(2,1) * 3^min(0,1)==2
lcm(4,6)==2^max(2,1) * 3^max(0,1)==12

补充

一个数的标准分解式不会太长，因为1，2，3，5，7，11，13，17，19，23就这几个数乘积就超过1e8了

代码

#include<bits/stdc++.h>
using namespace std;
const int MAXN=1e5;
int isprime(int n){
	if(n==1) return 0;
	for(int i=2;i*i<=n;i++){
		if(n%i==0) return 0;
	}
	return 1;
}
int pr[50],cnt;
int main()
{
	int n;
	cin>>n;
	for(int i=2;i<=50;i++){ //筛选出50以内的所有指数
		if(isprime(i)){
			pr[++cnt]=i;
		}
	}
	for(int i=1;i<=cnt&&pr[i]*pr[i]<=n;i++){//注意这里有优化pr[i]*pr[i]<=n
		if(n%pr[i]==0){
			int e=0;
			while(n%pr[i]==0){
				e++;
				n/=pr[i];
			}
			cout<<pr[i]<<"^"<<e<<endl;
		}
	}
	if(n>1) cout<<n<<"^"<<1<<endl;//当n>1时n一定也是质数且其一定大于sqrt(n)所以其指数一定为1
	return 0;
}