素数与筛法

素数判别

1.O（x）  [根本不用]

 

2.sqrt判别　　O（√N）

如果x可以表示为两个因子相乘

x=a*b    假设a<=b

那么x>=a*a

a<=√x

只需要枚举a<=√x就可以了

 

 

3.Miller-Rabin 素性测试

 

素性测试：

 

一个素数n，除了2，n-1　一定为偶数

 

a 属于（1，n-1 ），n是素数

 

条件1， 2 至少有一成立

 

Pｓ：s.t.使得

 

 

Longlong范围内测试这8个质数一般不会错

 

代码见下：

 

#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<string>
#include<cstring> #include<cstdlib> using namespace std; int gg[8] = {2,3,5,7,13,29,37,89}; int quickpow(int a,int b,int n) { int ans=1; while(b) { if(b%2==1) ans=ans*a%n; a=a*a%n; b/=2; } return ans; } bool miller_rabin(int a,int n)　　／／判断能否通过 { int d=n-1,r=0; while (d%2==0) d/=2,r++; int x = quickpow(a,d,n); if (x==1) return true; for (int i=0;i<r;i++) { if (x==n-1) return true; x=(long long)x*x%n; } return false; } bool is_prime(int n)　　　　／／Is-prime判断是否为素数

{
    if (n<=1) return false;　　／／特判
    for (int a=0;a<8;a++)
        if (n==gg[a]) return true;
    for (int a=0;a<8;a++)
        if (!miller_rabin(gg[a],n)) return false; return true; } int main() { long long N; cin>>N; if(is_prime(N)) cout<<"Prime"<<endl; else cout<<"Not Prime"<<endl; return 0; }

 

 

 

 

 

 

素数筛法

求1—n中所有质数

可以带入判素数   O（n√n）   但是会炸

1.未知名筛法

   先假定所有数字都是质数

   一个数的倍数一定不是质数，标记为合数

   最后没有被标记的，就是质数，因为不存在一个数i使得数i的倍数为这个数

 

算法复杂度达不到n²

b的枚举次数：

 

 

2.埃式筛法  O（nloglogn） 筛法

   合数的倍数一定会在筛素数倍数时候被筛掉

   所以只筛素数就好，只把质数的素数筛掉

   就是找到一个质数，把它的倍数全部标记为合数（但是你会发现有的数字会被标记多次，比如 6     被     2，3都标记，这样会浪费时间。。）

 

Ps:  1+1/2+1/3+1/4+......+1/n （调和级数） 这样的时间复杂度是 log n

 

 

 

 3.线性筛法

   P3383 【模板】线性筛素数

 

 

 

   线性筛法保证了每个数只会被他的最小质因子标记    复杂度降O（n）

  举个栗子：

   i = 5^3 * 7^2

   i2=5 * 5^3 * 7^2

    i3=3 * i2

     i4=2 * i3

从最小质因子开始，一直×到最小的2

   

 

3        i是质数，把它加入质数表

4        枚举质数表的数

5        枚举i的质数倍

6        i的质数倍不是质数

7        如果i是j的倍数，跳出

Why?

此时素数表里有 2 3 5 7 11 13

那就不会继续晒  15*3   15*7   15*11  15*13

 

所以就保证每个数只会被最小质因子筛掉

 

转载于:https://www.cnblogs.com/xiaoyezi-wink/p/10657881.html