【数论】素数

素数：

素数指的是：因子只有  1  和它本身。

素数的分布：

 

素数分布的应用

例题：素数个数的位数（nefu 117)

题目链接

素数个数的位数

Description

小明是一个聪明的孩子，对数论有着很浓烈的兴趣。他发现求1到正整数10n 之间有多少个素数是一个很难的问题，该问题的难以决定于n 值的大小。现在的问题是，告诉你n的值，让你帮助小明计算小于10n的素数的个数值共有多少位？

Input

输入数据有若干组，每组数据包含1个整数n（1 < n < 1000000000），若遇到EOF则处理结束。

Output

对应每组数据，将小于10n 的素数的个数值的位数在一行内输出，格式见样本输出。同组数据的输出，其每个尾数之间空一格，行末没有空格。

Sample Input

3

7

Sample Output

3

6

#include<cstdio>
#include<cmath>
#include<iostream>
using namespace std;
int main()
{
    double n;
    double t;
    while(cin>>n){
        t=double(n-log10(n)-log10(log(10)));
        cout<<int(t)+1<<endl;
    }
    return 0;
}

 

---

例题：哥德巴赫猜想（nefu 2）：

题目链接

Description

哥德巴赫（Goldbach ]C.，1690.3.18~1764.11.20）是德国数学家；出生于格奥尼格斯别尔格（现名加里宁城）；曾在英国牛津大学学习；原学法学，由于在欧洲各国访问期间结识了贝努利家族,所以对数学研究产生了兴趣；曾担任中学教师。1725年，到了俄国，同年被选为彼得堡科学院院士；1725年~1740年担任彼得堡科学院会议秘书；1742年，移居莫斯科，并在俄国外交部任职。

1742年，哥德巴赫在教学中发现，每个不小于6的偶数都是两个素数（只能被1和它本身整除的数）之和。如6＝3＋3，14＝3＋11等等。公元1742年6月7日哥德巴赫写信给当时的大数学家欧拉，欧拉在6月30日给他的回信中说，他相信这个猜想是正确的，但他不能证明。叙述如此简单的问题，连欧拉这样首屈一指的数学家都不能证明，这个猜想便引起了许多数学家的注意。从哥德巴赫提出这个猜想至今，许多数学家都不断努力想攻克它，但都没有成功。

我们不需要你去证明哥德巴赫猜想。

如果哥德巴赫猜想是正确的，一个（不小于6的）偶数，都是两个素数之和。那么这个偶数能被至少一个素数对表示，如14，即可以表示为14＝3＋11，也可以表示为14=7+7。不同的偶数对应的素数对的数目是不一样的，如偶数6，就只能表示为6=3+3。对于每个给定的偶数，我们希望知道有多少素数对的和等于该偶数。

Input

有多组测试数据。每组测试数据占一行，包含唯一的一个正偶数n.(6 <= n <= 2^24,)。 输出以EOF结束。

Output

对于每个输入的偶数，输出一行包含唯一的一个整数：表示有多少个素数对的和是输入的偶数。

Sample Input

6

14

Sample Output

1

2

---

 

【题解】：

直接暴力即可，首先我们需要打表出来，自愧不如，原来自己素数打表现在还是出毛病。

书本上标准的答案:

#include<stdio.h>
#include<math.h>
#include<string.h>
using namespace std;
const int N=16777226;
bool x[N];
int main()
{
    memset(x,true,sizeof(x));
    int e=static_cast < int > (sqrt((float)N));
    x[0]=x[1]=false;
    for(int i=4;i<=N;i+=2){
        x[i]=false;
    }
    for(int i=3;i<=e;i+=2){
        if(x[i]){
            for(int j=i*i;j<=N;j+=i){
                x[j]=false;
            }
        }
    }
    int n;
    while(scanf("%d",&n)!=EOF&&n){
        int ans=0;
        for(int i=2;i<n/2+1;i++){
            if(x[i]&&x[n-i]){
                ans++;
            }
        }
        printf("%d\n",ans);
    }
    return 0;
}

 

埃拉托色尼筛法：

#include<stdio.h>
#include<math.h>
#include<string.h>
using namespace std;
const int N=16777226;
typedef long long ll;
bool x[N];
int main()
{
    memset(x,true,sizeof(x));
    x[0]=x[1]=false;
    for(ll i=2;i<N;i++){
        if(x[i]){
            for(ll j=i*i;j<N;j+=i){
                x[j]=false;
            }
        }
    }
    int n;
    while(scanf("%d",&n)!=EOF&&n){
        int ans=0;
        for(int i=2;i<n/2+1;i++){
            if(x[i]&&x[n-i]){
                ans++;
            }
        }
        printf("%d\n",ans);
    }
    return 0;
}

 

6N±1    筛选法：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e6;
int prime[N],cnt=0;
inline bool Isprime(int x){
    if(x==2)return true;
    if(x%2==0){
        return false;
    }
    for(ll i=3;i*i<=x;i+=2){
        if(x%i==0){
            return false;
        }
    }
    return true;
}
void is_prime(){
    int t=0;
    prime[t++]=2;
    prime[t++]=3;
    for( ll i=6;i<=N;i+=6){
        for( ll j = -1; j <=1 ; j +=2){
            if(Isprime(i+j)){
                prime[t++]=(i+j);
            }
        }
    }

    printf("%d\n",t);
    cnt=t;
}
int main()
{
    int t1,t2;
    t1=clock();
    is_prime();
    t2=clock();
    printf("%d\n",t2-t1);
    /*for(int i=0;i<100;i++){
        printf("  ：%d\n",prime[i]);
    }*/
    return 0;
}
<