最一般的筛素数法

 1 int countPrimes(int n) {
 2         vector<bool> prime(n,true);
 3         prime[1]=prime[0]=false;
 4         for(int i=3;i<n;i++)
 5             if(i%2==0)
 6                 prime[i]=false;
 7
 8         for(int i=3;i<sqrt(n);i++)
 9             if(prime[i])
10                 for(int j=i*i;j<n;j+=i)
11                     prime[j]=false;
12
13         return count(prime.begin(),prime.end(),true);
14     }

时间: 2024-10-30 16:15:31

最一般的筛素数法的相关文章

筛素数法小结

筛选素数方法小结: 最简单的筛素数法方法就是从2开始,将所以2的倍数去掉,然后从3开始,将3的倍数去掉,依次进行下去即可.根据这样很容易写出代码,下面代码就是是筛素数法得到100以内的素数并保存到primes[]数组中. 1 const int MAXN = 100; 2 bool flag[MAXN]; 3 int primes[MAXN / 3], pi; 4 void GetPrime_1() 5 { 6 int i, j; 7 pi = 0; 8 memset(flag, false,

模板:筛素数法

参考:http://blog.csdn.net/liukehua123/article/details/5482854 1.开一个大的bool型数组prime[],大小就是n+1就可以了.先把所有的下标为奇数的标为true,下标为偶数的标为false. 2.然后: for( i=3; i<=sqrt(n); i+=2 ) {   if(prime[i]) for( j=i+i; j<=n; j+=i ) prime[j]=false; } 3.最后输出bool数组中的值为true的单元的下标,

Poj2689筛素数

题目大意: 给定一个区间l,r,求这个区间内相邻的质数中最近的两个和最远的两个.区间范围是1-2^31,区间的长度最多是10^6. 思路: 刚开始对筛选法的理解不深,不知道如何筛选任意一段区间的素数,看了题解恍然大悟,原来用的筛选法总是筛选从1-n的素数,对于为何这样筛选理解不深刻.说下1-n的筛选法,就是用一个数组is_prime[1..n]标记1-n中哪个是素数哪个不是,从2(第一个素数)开始扫描,所有2的倍数都不是素数,那么下一个素数是谁?就是我们向后从is_prime[]中找到的第一个是

线性筛素数详细整理

如果你在1个月前让我判断素数,我一定会猛敲出以下代码: bool check( int num ) { int tmp =sqrt( num); for(int i= 2;i <=tmp; i++) if(num %i== 0) return 0 ; return 1 ; //实在是太慢了! } $ $ 下面给大家带来3种筛选素数和一种直接判断素数 $ $ $ $ 什么是线性筛? 对于求多个质数时与其一个个判断不如用排除法,用空间换取大量时间. $ $ $ $ $ $ 一般筛法(埃拉托斯特尼筛法

Choose and divide UVA - 10375(筛素法+唯一分解定理的应用)

通过题目给的定义C(m,n)=m!/(n!(m-n)!),以及题目要求计算的C(p,q)/C(r,s)联立可得 p!s!(r-s)!/q!r!(p-q)! 看到这个式子,我们可以分析一下,我们可以将每个阶乘,都通过唯一分解定理将它们分解 (具体教程可见:https://blog.csdn.net/qq_39439314/article/details/78270905) 所以,首先我们应该先求出10000以内的所有素数,然后通过唯一分解定理将各个阶乘都分解,求出所有存在的可用质数以及对应的质数的

常见模板(欧拉筛素数,最小生成树,快排,并查集,单源最短路)

欧拉筛素数: #include<cstdio> #define maxn 10000000+10 using namespace std; int n,prime[5000001],num_prime=0,m; bool if_prime[maxn]; void euler(int limit) { for(int i=2;i<=limit;i++) { if(!if_prime[i]) prime[++num_prime]=i; for(int j=1;prime[j]*i<=l

洛谷 P3383 【模板】线性筛素数

P3383 [模板]线性筛素数 题目描述 如题,给定一个范围N,你需要处理M个某数字是否为质数的询问(每个数字均在范围1-N内) 输入输出格式 输入格式: 第一行包含两个正整数N.M,分别表示查询的范围和查询的个数. 接下来M行每行包含一个不小于1且不大于N的整数,即询问概数是否为质数. 输出格式: 输出包含M行,每行为Yes或No,即依次为每一个询问的结果. 输入输出样例 输入样例#1: 100 5 2 3 4 91 97 输出样例#1: Yes Yes No No Yes 说明 时空限制:5

筛素数

整理一下筛素数的方法 我在网上了解到两种筛素数的方法 一种是  1/3n*判断  的时间复杂度 一种是的时间复杂度应该是比这个低 先说一下第一种的思路 首先:一个数如果他除以一个素数除不尽,那么他除以该素数的倍数也除不尽 所以我们可以这么考虑 如果一个数是二或三的倍数 那么它一定不是素数 于是 对于  1 2 3 4 5 6 7 8 9 10 11 12…… 那么排除2和3的倍数 剩下的是 1 5 7 11 …… 对于六个数 6*n   6*n+1   6*n+2   6*n + 3   6*n

bjfu1211 推公式,筛素数

题目是求fun(n)的值 fun(n)= Gcd(3)+Gcd(4)+…+Gcd(i)+…+Gcd(n).Gcd(n)=gcd(C[n][1],C[n][2],……,C[n][n-1])C[n][k] means the number of way to choose k things from n things. n最大一百万,马上反映到可能是递推打表. 首先肯定是推公式了,fun(n)其实就是Gcd(n)的一个前n项和,没意义,直接看Gcd(n),把前几项列出来,发现公式是Gcd(n) =