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;
}