唯一分解定理

今天，我们要认识一个新....（朋友）

它是————唯一分解定理

在认识它之前，我们要请出它的一位搭档---------筛法（欧拉筛，埃氏筛...）[可自行选择观看，或跳过]

【筛法求质数/约数之和/最小因子和....】

先说一下>>

埃氏筛法：

它的基本思想就是建一个数组，用来标记它有没有被筛掉，如果没有，那就从它的2倍开始循环，一直到n（结束点），将它的倍数全部标记上（如2：就标记4、6、8、10.....）循环一遍之后，在

看下一个数有没有被标记，如果没有重复上面操作，否则，跳过。

具体代码如下：

#include <iostream>
using namespace std;
int book[100010];
int main(){
	int i,j,n;
	cin >> n;
	for(i = 2;i <= n;i++){
		if(book[i] == 0){
			for(j = i + i;j <= n;j += i){
				book[j] = 1;
			}
		}
	}
	for(i = 2;i <= n;i++){
		if(book[i] == 0)cout << i << " ";
	}
	return 0;
} 

但是，这种代码在计算一些大数据时，可能会超时（如：10的100次方），同时，有一部分数字，它会被筛掉很多次（不止一次），比如6（被2,3筛掉），21（被3,7筛掉）......

因此，又有了一个优化版——————欧拉筛

look 一下代码

#include <iostream>
using namespace std;
int prime[100010],cnt = 0,i,j;
bool st[100010];
void ola_prime(int n){
	for(i = 2;i <= n;i++){
		if(st[i] == 0)prime[++cnt] = i;
		for(j = 1;j <= cnt && prime[j] * i <= n;j++){
			st[i * prime[j]] = 1;
			if(i % prime[j] == 0)break;
		}
	}
	return ;
}
int main(){
	int k;
	cin >> k;
	ola_prime(k);
	for(i = 1;i <= cnt;i++)cout << prime[i] << " "; 
	return 0;
}

欧拉筛确保了一个数字只会被它最小的素因子筛掉，避免了重复去筛。

若大家有什么不懂的地方，可自行搜索，或去看CSDN上其他大佬的详解

啰嗦了这么久，开始介绍今天的主角————唯一分解定理

乘积为n的素数

时间限制：1秒        内存限制：128M

题目描述

给你一个正整数n(n>=2&&n<=1000000)根据唯一分解定理输出素数x1、x2、x3…xk，x1<=x2<=x3<=..<=xk，并且x1*x2*x3*…*xk等于n。

输入描述

输入一个正整数n。

输出描述

第一行输出k，
第二行输出素数x1，x2…xk使它们的乘积为n。

样例

输入

10

输出

2
2 5

​看不懂，没关系，先上代码

#include <iostream>
#include <cmath>
using namespace std;
int A[100010],S[100010],t = 0,i;
void ans(int h){
	int p = sqrt(h);
	for(i = 2;i <= p;i++){
		if(h % i == 0){
			A[++t] = i;
			while(h % i == 0){
				h /= i;
				S[t]++;
			}	
		}
	}
	if(h > 1)A[++t] = h,S[t] = 1; 
	return ;
}
int main(){
	int n,j,m = 0;
	cin >> n;
	ans(n);
	for(i = 1;i <= t;i++)m += S[i];
	cout << m << endl;
	for(j = 1;j <= m;j++){
		for(i = 1;i <= S[j];i++)cout << A[j] << " ";
	}
	return 0;
}

其中，ans函数就是唯一分解定理的模板代码。

首先，建两个变量，用来存n的因子和幂，再写个for循环，从2（最小质数）循环到根号n，如果i是n的因数， 就再加一个while循环，从n中除去i，并将i的幂加1。还有一种特殊情况，那就是n最多有1个大于根号n的因子，这个时候，我们在做一下判断就好了。

今天就到这里，拜拜！