关于素数

埃氏筛法：

用筛法求素数的基本思想：把从1开始的、某一范围内的正整数从小到大顺序排列， 1不是素数，首先把它筛掉。剩下的数中选择最小的数是素数，然后去掉它的倍数。依次类推，直到筛子为空时结束。如有：

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20

3 5 7 9 11 13 15 17 19 （从2 开始，筛去2的倍数之后剩余的数）

2 3 5 7 11 13 17 19 （从3开始，筛去3的倍数之后）

所以即为

最后剩余的人素数的个数。

素数筛代码：

int prime[MAXN]；
int isprime[MAXN + 1];
int con(int n)
{
    int p = 0;
    for(int i = 0; i <=n; i++)
    {
        isprime[i] = 1;
    }
    isprime[0] = isprime[1] = 0;
    for(int i = 2; i<=n; i++)
    {
        if(isprime[i])
        {
            prime[p++] = i;
            
        }
        for(int j = 2*i; j <= n; j = j +i)
        {
            isprime[j] = 0;
        }
        
    }
    return p;
}

最直接的应用：sdnu 1522 

2编号1，3编号2，5编号3........任给一个数，求这个数的所有的质因子里面的最大编号是多少

注：0和1的编号都是0。

思路：利用素数筛思想，把存每个素数倍数的数组改成存该素数，则每次查找时就是最大的质因子，用另一个数组存该最大质因子的编号

第一次提交的时候把prime（）函数写在while里了，T，放在外边进行预处理只进行一次，A了。

#include<cstdio>

#include<iostream>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#define maxa 1000000
using namespace std;
int m[maxa];
int a[maxa];
int prime()
{
    int i,j;
    memset(m,0,sizeof(m));
    m[0] = m[1] = 1;
    int cnt = 0;
    for(i = 2; i <= maxa; i++)
    {
        if(!m[i])
        {
            cnt++;

            for(j = 1; j*i <= maxa; j++)
            {
                m[i*j] = i;

                a[i] = cnt;
            }
        }

    }
}

int main()
{
    int n;

prime（）;

   while(scanf("%d",&n)!=EOF)
    {
        if(n==0||n==1)
            printf("1\n");
        else
        {
        //    prime();

            printf("%d\n",a[m[n]]);
        }
    }

}

区间筛法：  a < b  < 10^12

                b  - a < 10^6

a,b范围太大，不能直接求int范围的素数。而区间间隔比较小，只有1e6，而且对于int范围内的合数来说，最小质因子必定小于2^16。所以可以求出[a,b]中合数，转而求出素数，然后暴力枚举所有素数对即可。

如何求区间[a,b]中的合数：合数的最小质因子小于2^16，即小于50000。所以先求出小于50000的所有素数。则区间[l,u]中的合数，必定可以表示为小于50000的素数的倍数。枚举所有小于50000的素数，再找出[a，b]中所有的合数。

typedef long long ll;
bool is_prime[1000000+5];
bool is_prime_small[1000000+5];
//对区间[a,b)内的整数执行筛法。is_prime[i-a]=true→i是素数
void segment_sieve(ll a,ll b)
{
    for(int i=0;(ll)i*i<b;i++) is_prime_small[i]=true;
    for(int i=0;i<b-a;i++) is_prime[i]=true;
    for(int i=2;(ll)i*i<b;i++){
        if(is_prime_small[i]){
            for(int j=2*i;(ll)j*j<b;j+=i) is_prime_small[i]=false;//筛[2,√b)
            for(ll j=max(2ll,(a+i-1)/i)*i;j<b;j+=i) is_prime[j-a]=false;//筛[a,b)
            //2LL是2的长整数形式
            //((a+i-1)/i)*i是满足>=a&&%i==0的离a最近的数
            //也可以写成(a%i==0)?a:(a/i+1)*i
        }
    }
}