质数筛选
给定一个整数N,求出1到N之间的所有质数,称为质数的筛选。
质数筛法
埃氏筛法-Eratosthenes
基于任意质数m的倍数2m,3m,…都不是质数,实现算法。
时间复杂度:O(NloglogN)
//核心代码
void EraPrime(int n)
{
memset(st, true,sizeof st);
cnt = 0;
for(int i = 2; i <= n ; i ++)
{
if(!st[i]) continue;
else primes[cnt ++] = i;//记录质数
for(int j = i; j <= n /i; j ++)
st[i * j] =false;//i*j是合数,被筛掉
}
}
线性筛法
埃氏筛法会重复标记合数。
线性筛法用累积质因子的方式标记每个合数,标记每个合数只筛一次。
每个合数 i * p 只会被它的最小质因子p筛一次。
时间复杂度O(N)
//核心代码
void linePrime(int n)
{
for(int i = 2; i <=n; i ++)
{
if(!st[i]) prime[cnt ++] = i;
for(int j = 0; prime[j] * i <= n; j ++)
{
st[prime[j] * i] = true;//deleted
if(i % prime[j] == 0) break;// i是质数,i== prime[j]时结束循环
// i是合数,能整除时结束循环
}
}
}
猜想
哥德巴赫猜想的内容如下:
任意一个大于 4的偶数都可以拆成两个奇素数之和。
例如:
8=3+5
20=3+17=7+13
42=5+37=11+31=13+29=19+23
现在,你的任务是验证所有小于一百万的偶数能否满足哥德巴赫猜想。
输入格式
输入包含多组数据。
每组数据占一行,包含一个偶数 n。
输出格式
对于每组数据,输出形如 n = a + b,其中 a,b是奇素数。
若有多组满足条件的 a,b,输出 b−a最大的一组。
若无解,输出 Goldbach’s conjecture is not correct。
数据范围
6≤n<106
输入样例:
8
100
180
0
输出样例:
8 = 3 + 5
100 = 3 + 97
180 = 7 + 173
#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
using namespace std;
const int N = 1e6 + 7;
int primes[N], cnt;
bool st[N];
void linePrime(int n)
{
for(int i = 2; i <=n; i ++)
{
if(!st[i]) primes[cnt ++] = i;
for(int j = 0; primes[j] * i <= n; j ++)
{
st[primes[j] * i] = true;//deleted
if(i % primes[j] == 0) break;
}
}
}
void asPrime(int n)
{
memset(st, true,sizeof st);
cnt = 0;
for(int i = 2; i <= n ; i ++)
{
if(!st[i]) continue;
else primes[cnt ++] = i;
for(int j = i; j <= n /i; j ++)
st[i * j] =false;
}
}
int main()
{
linePrime(N - 1);
int n;
while(cin >> n && n)
{
for(int i = 1; ; i ++)
{
int x = primes[i];//from 3
int y = n - x;
if(!st[y])
{
printf("%d = %d + %d\n",n,x,y);
break;
}
}
}
return 0;
}