线性筛求素数线性筛求因子个数线性筛求因子合

最新推荐文章于 2021-07-02 22:26:27 发布

never疯

最新推荐文章于 2021-07-02 22:26:27 发布

阅读量871

点赞数

分类专栏： C语言算法

本文链接：https://blog.csdn.net/qq_40574571/article/details/85854444

版权

C语言算法专栏收录该内容

10 篇文章 1 订阅

订阅专栏

线性筛求素数线性筛求因子个数线性筛求因子合

普通线性筛筛素数

#include<stdio.h>
#define max_n 10005
int num[max_n]={0}
int main(){
	for(int i=2;i<max_n;i++){
		if(num[i]==0){
	  	num[++num[0]]=i;
		}
		for(int j = 1; j <= num[0] ; j++){
				if(i * num[j] > max_n) break;
				num[i * num[j]]=1;
				if(i % num[j]==0)break;	
		}	
	}
return 0;
}

线性筛求因子个数

#include<stdio.h>
#define max_n 10005
int num[max_n]={0};
int prime_num[max_n]={0};
int snum[max_n]={0};
int main(){
	for(int i=2;i<max_n;i++){
		if(num[i]==0){
	  		num[++num[0]]=i;
	  		snum[i] = 2;
	  		prime_num[i] = 1;
		}
		for(int j = 1; j < =num[0] ; j++){
				if(i * num[j] > max_n) break;
				num[i * num[j]]=1;
				if(i % num[j]==0){
					prime_num[i * num[j]]=prime_num[i] + 1;
					snum[i * num[j]] = snum[i] / (prime_num[i] + 1) * (prime_num[i] + 2);
					break;
				}else{
					snum[ i * num[j] ] = snum[i] * snum[num[j]];
					prime_num[ i * num[j] ] = 1;
				}	
		}	
	}
return 0;
}

线性筛求因子合

#include<stdio.h>
#define max_n 10005
int num[max_n]={0};
int fnum[max_n + 5]={0};
int min_num[max_n+5]={0};
void init(){
	for(int i=2;i<max_n;i++){
		if(num[i]==0){
	  	num[++num[0]]=i;
	  	min_num[i]=1;
	  	fnum[i]= 1 + i;
}
		for(int j = 1; j <= num[0] ; j++){
				if(i * num[j] > max_n) break;
				num[i * num[j]]=1;
				if(i % num[j]==0){
						min_num[num[j] * i]=min_num[i] + 1;
						fnum[num[j] * i] = fnum[i] / (pow(num[j],min_num[i]+1)-1)*( pow(num[j],min_num[i]+2)-1);
				}else{
				min_num[num[j]*i] = 1;
				fnum[ i * num [ j ]] =fnum[num[j]]*fnum[i]; 
				}	
		}	
	}
return ;
}