素数筛时间空间优化

素数筛时间空间优化

时间优化

普通素数筛原理 ：素数的倍数一定不是素数；据此筛选素数。

void work1()
{
    int cnt=0;
    for(int i=2; i<=n; i++)
    {
        if(!flag[i])
        {
            prime[cnt++]=i;
            /**在此处可有个小优化 j=i*i ，但是i*i可能爆int */
            /**因为>=2 && <i 那部分在 i=2~i-1 中已经标记了 */
            for(int j=i*2; j<=n; j+=i)
                flag[j]=1;
        }

    }
}

可以发现，在标记某些数字的时候，出现了重复：比如 12，在 i=2 时访问到以后，在 i=3 使也能访问到。

那么就进行改进 ，保证每个数字只被访问一次。

时间优化原理 ：每个合数必有一个最小素因子，根据每个最小素因子去访问合数就能防止重复访问。

void work2()
{
    int cnt=0;
    memset(flag,0,sizeof(flag));
    for(int i=2; i<=n; i++)
    {
        if(!flag[i])
            prime[cnt++]=i;
    
        for(int j=0; (i*prime[j]<=n) && (j<cnt) ; j++)
        {
            flag[i*prime[j]]=1;///素数的倍数
            if(i%prime[j]==0)/**保证每个数只被标记一次*/
                break;
        }
    }
}

理解 break 语句，当 i = 4 时 ，只标记8，然后跳出；其实可标记的还有：12，20，28……但是，因为 4%2==0 ,所以 i = 6时标记12，i = 10 时标记 20，i = 14 时标记28；

它们均可由 2 * i 得到。

空间优化

空间优化原理 ：因为 int 型数字有4字节 = 32位，我们可用每一位上的 0 1 来作为标记。

PS：这个方法给我的感觉是操作系统中的分页管理存储方式

void work3()
{
    int cnt=0;
    for(int i=2;i<=n;i++)
    {
        /**相当于分页 ，第 i/32 页第 (i%32) 位*/
        if(!(f[i/32] &(1<<(i%32))))
            prime[cnt++]=i;

        for(int j=0;(j<cnt)&&(i*prime[j]<=n);j++ )
        {
            int t=i*prime[j];
            f[t/32] |= (1<<(t%32));
            if(t==0)
                break;
        }
    }
}

撸代码：

#include<stdio.h>
#include<string.h>
#define N 10000000
const int n=1e2;
bool flag[N];
int prime[5900000];
/**
素数的倍数一定不是素数
*/
void work1()
{
    int cnt=0;
    for(int i=2; i<=n; i++)
    {
        if(!flag[i])
        {
            prime[cnt++]=i;
            for(int j=i*2; j<=n; j+=i)
                flag[j]=1;
        }

    }
}

/**每个合数必有一个最小素因子，
根据每个最小素因子去访问合数就能防止重复访问*/
void work2()
{
    int cnt=0;
    memset(flag,0,sizeof(flag));
    for(int i=2; i<=n; i++)
    {
        //printf("\n %d: *********\n",i);
        if(!flag[i])
            prime[cnt++]=i;
        for(int j=0; (i*prime[j]<=n) && (j<cnt) ; j++)
        {
            flag[i*prime[j]]=1;///素数相乘的倍数
          //  printf("%d ",i*prime[j]);
            if(i%prime[j]==0)/**保证每个数只被标记一次*/
                break;
        }
    }
}
/**空间优化*/
/***/
int f[N/32];
void work3()
{
    int cnt=0;
    for(int i=2;i<=n;i++)
    {
        /**相当于分页 ，第 i/32 页第 (i%32) 位*/
        if(!(f[i/32] &(1<<(i%32))))
            prime[cnt++]=i;

        for(int j=0;(j<cnt)&&(i*prime[j]<=n);j++ )
        {
            int t=i*prime[j];
            f[t/32] |= (1<<(t%32));
            if(t==0)
                break;
        }
    }
}
int main()
{
    work3();
    for(int i=0; i<100; i++)
        printf("%d ",prime[i]);
    return 0;
}