今天,我们要认识一个新....(朋友)
它是————唯一分解定理
在认识它之前,我们要请出它的一位搭档---------筛法(欧拉筛,埃氏筛...)[可自行选择观看,或跳过]
【筛法求质数/约数之和/最小因子和....】
先说一下>>
埃氏筛法:
它的基本思想就是建一个数组,用来标记它有没有被筛掉,如果没有,那就从它的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的因子,这个时候,我们在做一下判断就好了。
今天就到这里,拜拜!