哥德巴赫猜想

质数筛选

给定一个整数N,求出1到N之间的所有质数,称为质数的筛选。

质数筛法

埃氏筛法-Eratosthenes

基于任意质数m的倍数2m,3m,…都不是质数,实现算法。
时间复杂度:O(NloglogN)

//核心代码
voidraPrime(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;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值