高程期末复习 3 素数判定-线性筛和埃氏筛考试前冲刺简约版

一个大于2的偶数总可以分解成两个素数的和

 拆分素数和问题，你看看你不查看任何资料的情况下能不能自己打出来，用传统的方法  进行遍历

#include <iostream>
#include <math.h>

using namespace std;
int prime(int a)
{
    int i;
    for(i=2;i<=(int)sqrt(a);i++)
    {
        if(a%i==0)
        return 0;
    }
    return 1;
}
int main()
{
    int n,ans=0,i;
    cin>>n;
    for(i=2;i<n/2;i++)
    {
        if(prime(i)==1&&prime(n-i)==1)
        ans++;
    }
    cout<<ans<<endl;
    return 0;
}

我自己打了一遍：我的严重问题是：在拆分素数主函数for循环逐一判断时，循环的终止条件我居然写的是sqrt(n),没有认真去想边界条件，应该是n/2,这个错了基本运行不出来了，sqrt（n）范围太小了，第二个应该注意的是素数函数中sqrt前面最好加(int )进行强制类型转换，还有素数的判断条件是%之后得0.

注意：才拆分两个不同的素数和有两种方法，一种是在if条件中加i!=(n-i)   第二种是判断边界退出的时候去掉等号。

素数定理：

 

 计算一个数有多少位   （int)log10(n)+1

 

#include <iostream>
#include <math.h>
using namespace std;

int main()
{
   double n,m;
   while(cin>>n)
   {
      m=n-log10(n)-log10(log(10))+1;
      cout<<int(m)<<endl;
   }

   return 0;
}

这道题结合素数定理和求位数的方法很容易求解，可知这样的题肯定不是考试的重点，我们先练练手，重点在下面

埃氏筛

 筛 ：怎么实现?  (c语言中借助整型数组，0是素数，1不是)程序开始时都设成1，默认都是素数

利用埃氏筛

#include <iostream>

using namespace std;

#define M 101
int b[M]={0};
int prime[M];
int num=0;
int find_prinme()
{
    int i,j;
    b[0]=b[1]=0;
    for(i=2;i<M;i++)
    {
        if(b[i]==0)
        {
            prime[++num]=i;
            for(j=i+i;j<M;j+=i)
            {
                b[j]=1;
            }
        }
    }
}
int main()
{
    int i;
    find_prinme();
    for(i=1;i<=num;i++)
    {
        cout<<prime[i]<<endl;
    }
    return 0;
}

自己打代码出现的问题是num已经全局变量定义了，就不要在找素数的函数中重新定义，一定要注意，这个不容易发现错误，并且最后结果什么也不输出。

线性筛

埃氏筛有的数被重复的删掉，所以我们引入了线性筛。

求100以内所有素数用线性筛的做法

#include <iostream>

using namespace std;
#define M 101
int b[M]={1,1};
int prime[M],num;
int find_prime()
{
    int i,j;
    for(i=2;i<M;i++)
    {
            if(b[i]==0)
            prime[++num]=i;
            for(j=1;j<=num&&i*prime[j]<M;j++)
            {
                b[i*prime[j]]=1;
                if(i%prime[j]==0)
                break;
            }

    }
}
int main()
{
    int i;
    find_prime();
    for(i=1;i<=num;i++)
    cout<<prime[i]<<endl;
    return 0;
}