素数:
素数指的是:因子只有 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;
}
<