이전 기사는 간격 소수의 합을 찾는 방법에 대해 썼고, 간격 소수를 찾는 것에 대해 또 다른 기사를 썼습니다. 사실, 소수를 찾기 위해 min25 체의 일부만 사용되며 변화는 크지 않습니다. 구체적인 아이디어는 내 min25 체 연구 노트를 참조하십시오.
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod=998244353;
const int maxn = 1e6 + 10;
bool valid[maxn];
ll prime[maxn];
ll sum[maxn];
int tot;
void init()
{
tot = 0;
valid[1] = true;
for(int i = 2; i < maxn; i++)
{
if(!valid[i])
{
prime[++tot] = i;
sum[tot] = sum[tot - 1] + ll(1);
}
for(int j = 1; j <= tot && i * prime[j] < maxn; j++)
{
valid[i * prime[j]] = true;
if(i % prime[j] == 0) break;
}
}
}
ll id1[maxn], id2[maxn];
ll w[maxn], g[maxn];
ll solve(ll n)
{
ll sqr = sqrt(n + 0.5);
int ntot = 0;
for(ll i = 1; i <= tot; i++)
{
if(prime[i] <= sqr) ntot++;
else break;
}
int m = 0;
for(ll i = 1; i <= n; i = n / (n / i) + 1)
{
w[++m] = n / i;
if(w[m] <= sqr) id1[w[m]] = m;
else id2[n / w[m]] = m;
g[m] = w[m] - 1;
}
for(int j = 1; j <= ntot; j++)
{
for(int i = 1; i <= m && (ll)prime[j] * (ll)prime[j] <= w[i]; i++)
{
ll k = (w[i] / prime[j] <= sqr) ? id1[w[i] / prime[j]] : id2[n / (w[i] / prime[j])];
g[i] = g[i] - (g[k] - sum[j - 1]);
}
}
return g[id2[1]];
}
int main()
{
init();
ll n;
scanf("%lld", &n);
printf("%lld\n", solve(n));
return 0;
}