7-5 最强素数 线性筛+前缀和+暴力

小李在你帮助之下轻松战胜了他的同学们,于是满怀恶意的同学出了一个题目来为难小李,作为小李神一样的队友,你又要出力了。 素数41能写成连续6个素数之和:41=2+3+5+7+11+13。 现在要求n以内的素数中,能表示为最多连续素数之和的那个数,如果有多个答案,请输出最大的那个素数。
输入格式:

仅一行,一个整数n。
输出格式:

输出就一个整数,为所求的能表示为最多连续素数和的那个素数。
输入样例:

100

输出样例:

41

【样例说明】 41=2+3+5+7+11+13 【数据规模】 30%的数据,1<=n<=1000 60%的数据,1<=n<=10000 80%的数据,1<=n<=100000 100%的数据,1<=n<=1000000

线性筛

void get_prime() {
	int N = n + 3;
	vis[0] = vis[1] = true;
	for(int i=2; i<N; i++) {
		if(!vis[i]) prime[++pn] = i;
		for(int k=1; k<=pn && prime[k]*i<N; k++) {
			vis[prime[k]*i] = true;
			if(i % prime[k] == 0) break;
		}
	}
}
#ifdef debug
#include <time.h>
#include "/home/majiao/mb.h"
#endif

#include <iostream>
#include <algorithm>
#include <vector>
#include <string.h>
#include <map>
#include <set>
#include <stack>
#include <queue>
#include <math.h>

#define MAXN ((int)1e6+7)
#define ll long long int
#define INF (0x7f7f7f7f)
#define fori(lef, rig) for(int i=lef; i<=rig; i++)
#define forj(lef, rig) for(int j=lef; j<=rig; j++)
#define fork(lef, rig) for(int k=lef; k<=rig; k++)
#define QAQ (0)

using namespace std;

#define show(x...) \ do { \ cout << "\033[31;1m " << #x << " -> "; \ err(x); \ } while (0)

void err() { cout << "\033[39;0m" << endl; }
template<typename T, typename... A>
void err(T a, A... x) { cout << a << ' '; err(x...); }

namespace FastIO{

	char print_f[105];
	void read() {}
	void print() { putchar('\n'); }

	template <typename T, typename... T2>
		inline void read(T &x, T2 &... oth) {
			x = 0;
			char ch = getchar();
			ll f = 1;
			while (!isdigit(ch)) {
				if (ch == '-') f *= -1; 
				ch = getchar();
			}
			while (isdigit(ch)) {
				x = x * 10 + ch - 48;
				ch = getchar();
			}
			x *= f;
			read(oth...);
		}
	template <typename T, typename... T2>
		inline void print(T x, T2... oth) {
			ll p3=-1;
			if(x<0) putchar('-'), x=-x;
			do{
				print_f[++p3] = x%10 + 48;
			} while(x/=10);
			while(p3>=0) putchar(print_f[p3--]);
			putchar(' ');
			print(oth...);
		}
} // namespace FastIO
using FastIO::print;
using FastIO::read;

int n, m, Q, K, prime[MAXN], pn, sum[MAXN];
bool vis[MAXN];

void get_prime() {
	int N = n + 3;
	vis[0] = vis[1] = true;
	for(int i=2; i<N; i++) {
		if(!vis[i]) prime[++pn] = i;
		for(int k=1; k<=pn && prime[k]*i<N; k++) {
			vis[prime[k]*i] = true;
			if(i % prime[k] == 0) break;
		}
	}
}

#define ISPRIME(x) (!vis[x])

int main() {
#ifdef debug
	freopen("test", "r", stdin);
	clock_t stime = clock();
#endif
	cin >> n;
	get_prime();
	int anslen = 0, ans = 0;
	for(int i=1; i<=pn; i++) sum[i] = sum[i-1] + prime[i];
	for(int i=1; i<=pn; i++) { //暴力
		for(int j=i; j<=pn; j++) {
			int s = sum[j] - sum[i-1];
			if(s > n) break;
			int len = j - i + 1;
			if(ISPRIME(s)) {
				if(anslen < len) anslen = len, ans = s;
				else if(anslen == len) ans = max(ans, s);
			}
		}
	}
	printf("%d\n", ans);




#ifdef debug
	clock_t etime = clock();
	printf("rum time: %lf 秒\n",(double) (etime-stime)/CLOCKS_PER_SEC);
#endif 
	return 0;
}




全部评论

相关推荐

点赞 收藏 评论
分享
牛客网
牛客企业服务