不图啥,只是为了我敲代码时,敲的更快一点 heihei
质数是指在大于1的自然数中,除了1和它本身以外不再有其他因数的自然数。
-----百度百科
1,判定质数
先上我们暴力的算法,叫做试除法,
bool prime(int n){
if(n<2)return false;
for(int i=2;i<n;i++)
if(n%i==0)return false;
return true;
}
时间复杂度 O(n)
假设我们判断12是否为质数,3是12的约数,4也是12的约束,又比如说6,2是6的约数,3也是6的约数,如此我们发现,n的约数总是成对出现的,所以我们只需要扫描比较小的约数,
bool prime(int n){
if(n<2)return false;
for(int i=2;i<n/i;i++)
if(n%i==0)return false;
return true;
}
2,分解质因数
暴力
void divide(int n){
for(int i=2;i<=n;i++){
if(n%i==0){
int cnt=0;
while(n%i==0){
n/=i;
cnt++;
}
printf("%d %d",i,cnt)
}
}
}
还是用我们上边的思想
void divide(int n){
for(int i=2;i<=n/i;i++){
if(n%i==0){
int cnt=0;
while(n%i==0){
n/=i;
cnt++;
}
printf("%d %d\n",i,cnt);
}
}
if(n>1)printf("%d %d\n",n,1);
}
3,筛质数
我们从最小的2开始,依次从前向后的去掉他倍数的数,如果F(n)没有被去掉那就是质数,我们设置了st[]数组来标记是否被去掉,如果没有去掉那就让cnt加1
void get_primes(int n){
for(int i=2;i<=n;i++){
if(!st[i]){
cnt++;
}
for(int j=i+i;j<=n;j+=i)st[j]=true;
}
}
印式筛选
上面我们不管这个数是不是质数都给他进行了一次洗牌的过程,其实我们不需要筛掉所有数的,只需要筛掉质数的倍数的数就可以了
void get_primes(int n){
for(int i=2;i<=n;i++){
if(!st[i]){
cnt++;
for(int j=i+i;j<=n;j+=i)st[j]=true;
}
}
}
所以我们把筛除的操作放到if 条件里面。
线性筛选
线性筛选每次删选的是他们的最小质因子,我们从最小的质数筛选,能删掉i的一定是他的最小质因子
埃氏筛选法的效率已经不低了,但是会有很多的重复计算,即一个数会被多次踢出,浪费了计算时间。对于任何一个合数z,一定存在质数x,使得:
z = x * y (y 不为 1)
但是这样的x不唯一,比如6=23=32,这里x = 2或者3。所以我们现在要找一个办法,让每一个合数只能被踢出一次而不是多次。这时注意到x的最小值是一定的,所以我们可以用每一个合数的最小质因子来踢出。
void get_primes(int n){
for(int i=2;i<=n;i++){
if(!st[i])primes[cnt++]=i;
for(int j=0;primes[j]<=n/i;j++)
{
st[primes[j]*i]=true;
if(i%primes[j]==0)break;
}
}
}