loj #6235. 区间素数个数

#6235. 区间素数个数

题目描述

求 1∼n 1\sim n1∼n 之间素数个数。

输入格式

一行一个数 n nn 。

输出格式

一行一个数，表示答案。

样例

样例输入

10

样例输出

4

样例解释 1

2,3,5,72,3,5,72,3,5,7

数据范围与提示

对于 100% 100\%100% 的数据，2≤n≤1011 2 \leq n \leq 10^{11}2≤n≤10​11​​。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define maxn 316228
using namespace std;
long long n;
bool vis[maxn];
int lim,p[],sum[maxn],last[maxn*],cnt;
long long val[maxn*],f[maxn*];
void prepare(){
    for(int i=;i<=lim;i++){
        if(!vis[i])p[++p[]]=i;
        sum[i]=sum[i-]+!vis[i];
        for(int j=;j<=p[]&&i*p[j]<=lim;j++){
            vis[i*p[j]]=;
            if(i%p[j]==)break;
        }
    }
}
int main(){
    cin>>n;lim=sqrt(n);
    prepare();
    for(long long i=;i<=n;i=n/(n/i)+)
        val[++cnt]=n/i;
    reverse(val+,val+cnt+);
    copy(val+,val+cnt+,f+);
    for(int i=;i<=p[];i++){
        for(int j=cnt;j>=;j--){
            long long k=val[j]/p[i];
            long long pos=k<=lim?k:cnt+-n/k;
            if(k<p[i])break;
            f[j]-=f[pos]+last[pos]-i+;
            last[j]=i;
        }
    }
    printf("%lld",sum[lim]+f[cnt]-);
    return ;
}