PTA数据结构7-2 求素数个数

pta7-2 求素数个数
求素数的个数。本题要求编写一个程序,求1~n的素数个数。 要求至少给出两种解法,对于相同的n,给出这两种解法的结果,通过相关数据进行测试,目的是通过对比同一问题不同解法的绝对执行时间体会如何设计“好”的算法。
输入格式:
输入在一行中给出1个整数n(<= 10 000 000)。
输出格式:
对每一组输入,在一行中输出1~n的素数个数。
输入样例1:
5
输出样例1:
3
输入样例2:
14
输出样例2:
6

#include <iostream>
#include <algorithm>
#include <cstring>//memset函数头文件
using namespace std;
const int MAX_N=10000100;
bool p[MAX_N];
int pr[MAX_N];
void Euler() {//欧拉筛法
    int cnt = 0;
    memset(p, true, sizeof(p));//将p内的数都变成true
    p[0] = p[1] = false;//让p[0].p[1]都等于false
    for (int i = 2; i < MAX_N; ++i) {
        if (p[i]==1) 
        pr[cnt++] = i;    
        for (int j = 0; j < cnt && pr[j] * i < MAX_N; ++j) {
            p[pr[j] * i] = false;    
            if (i % pr[j] == 0) 
            break;   
        }
    }
}
int main(){
	Euler();
	int n,count=0;
	cin >> n;
	for(int i = 0;i <= MAX_N;i++){
		
    	if(pr[i] <= n)//计算n以内的素数个数
    		count++;
    	else
    		break;
	}
	cout << count << endl;
	return 0;
}

memset函数
memset:char型初始化函数
头文件:<string.h> 或 <memory.h>
函数原型:void *memset(void *s , int ch , size_t n )
memset(结构体/数组名 , 用于替换的ASCII码对应字符 , 前n个字符 );
memset(结构体/数组名 , "用于替换的字符“ , 前n个字符 );
函数解释:将s中的前n个字节用ch替换并且返回s
函数作用:在一段内存块中填充某一个给定的值,常用于较大的对结构体和数组的清零操作。
memset()函数原型是extern void *memset(void *buffer, int c, int count) buffer:为指针或是数组,
c:是赋给buffer的值,
count:是buffer的长度.
这个函数在socket中多用于清空数组.
如:原型是memset(buffer,0,sizeof(buffer))
Memset 用来对一段内存空间全部设置为某个字符,一般用在对定义的字符串进行初始化为‘ ’或‘/0’;
在本题中,一开始首先想到的是C和C++中学到的求素数的方法,但是在当中发现运行超时,他人告知有求素数的另外方法,即欧拉筛法和埃氏筛法,但是埃氏筛法也是超时的
欧拉筛法
欧拉筛法的基本思想 :在埃氏筛法的基础上,让每个合数只被它的最小质因子筛选一次,以达到不重复的目的。
代码 :

int prime[maxn];
int visit[maxn];
void Prime(){
    mem(visit,0);
    mem(prime, 0);
    for (int i = 2;i <= maxn; i++) {
        cout<<" i = "<<i<<endl;
        if (!visit[i]) {//!取反,即visit[i]==1时执行
            prime[++prime[0]] = i; //纪录素数, 这个prime[0] 相当于cnt,用来计数
        }
        for (int j = 1; j <=prime[0] && i*prime[j] <= maxn; j++) {
// cout<<"  j = "<<j<<" prime["<<j<<"]"<<" = "<<prime[j]<<" i*prime[j] = "<<i*prime[j]<<endl;
            visit[i*prime[j]] = 1;
            if (i % prime[j] == 0) {
                break;
            }
        }
    }
}

1.对于visit[i*prime[j]] = 1 的解释: 这里不是用i的倍数来消去合数,而是把 prime里面纪录的素数,升序来当做要消去合数的最小素因子。
现i在消去合数中的作用是当做倍数的。
2.对于 i%prime[j] == 0 就break的解释 :当 i是prime[j]的倍数时,i = kprime[j],如果继续运算 j+1,i * prime[j+1] = prime[j] * k prime[j+1],这里prime[j]是最小的素因子,当i = k * prime[j+1]时会重复,所以才跳出循环。
举个例子 :i = 8 ,j = 1,prime[j] = 2,如果不跳出循环,prime[j+1] = 3,8 * 3 = 2 * 4 * 3 = 2 * 12,在i = 12时会计算。因为欧拉筛法的原理便是通过最小素因子来消除。

  • 3
    点赞
  • 17
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值