Goldbach`s Conjecture LightOJ - 1259 欧拉筛(数论专题)

n=a+b
满足a,b(a<=b)都是质数,a+b=n方案个数
Input
输入一个T,代表样例数量 T(≤ 300)。
每一行输入一个n (4 ≤ n ≤ 107)。
奥奥 n是个偶数。
Output
每一行输出样例的行号以及可以表示n的两个素数之和的方式数。
Sample Input
2

6

4

Sample Output
Case 1: 1

Case 2: 1
这道题爆了好多次,不是RE就是TLE,甚至还有MLE,最后在筛法中欧拉筛是最快的,这里我放一个TLE的,先对每一个数字进行标记数据范围内所有的质数,然后根据所给条件判断。
错误代码

#include<stdio.h>
#include<string.h>
#include<math.h>
#include<map>
#include<stack>
#include<set>
#include<queue>
#include<vector>
#include<stdlib.h>

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<cstdlib>

#define inf 0x3f3f3f3f
using namespace std;
typedef long long ll;
const int maxn=1e7+1;
const int mod=1000;
ll t,n,m;
bool vis[maxn];
void init(int n){
    for(int i=2;i<=n;i++){
        if(vis[i])continue;
        for(int j=i*2;j<=n;j+=i)vis[j]=1;
    }
}
int main()
{
    init(10000000);
    int t;
    cin>>t;
    int cas=0;
    while(t--)
    {
        ll n;
        int sum=0;
        cin>>n;
        for(int i=2;i<=n/2;i++){
           int j=n-i;
           if(vis[i]==0&&vis[j]==0){
               sum++;
               cout<<sum<<endl;
           }
        }
        printf("Case %d: %d\n",++cas,sum);
    }
	system("pause");
    return 0;
}

正确是使用欧拉筛,对1e7的数据范围内记录质数,因为1e7中质数个数并没有太多所以记得要对记录质数的数据的范围开小到1e6,不然会爆

#include<map>
#include<stack>
#include<queue>
#include<string>
#include<math.h>
#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
//#include<bits/stdc++.h>
using namespace std;
const int maxn=1e7+1;
typedef long long ll;
const int inf=0x3f3f3f3f;
const int minn=0xc0c0c0c0;
bool vis[maxn];
ll n,t,ans,cnt,prime[1000000];
void eul(ll n)
{
	cnt=0;
	for(ll i=2;i<=n;i++)
	{
		if(!vis[i])
			prime[++cnt]=i;
		for(ll j=1;j<=cnt&&i*prime[j]<=n;j++)
		{
			vis[i*prime[j]]=1;
			if(i%prime[j]==0)
				break;
		}
	}
}
int main()
{
	ios::sync_with_stdio(false);
	cin.tie(0);cout.tie(0);
	eul(10000000);
	scanf("%lld",&t);
	for(int ca=1;ca<=t;ca++)
	{
		ans=0;
		scanf("%lld",&n);
		for(int i=1;i<=cnt;i++)
		{
			if(!vis[n-prime[i]]&&prime[i]<=n-prime[i])
				ans++;
			if(prime[i]>n)
				break;
		}
		printf("Case %d: %lld\n",ca,ans);
	}
	return 0;
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值