LOJ 6235 洲阁筛

尝试理解Min_25筛失败，绝望地用洲阁筛切了个板题。
LOJ 6235
求小于等于n的素数个数。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int LIM=400000,P=30000;
ll sqr;
bool vis[LIM];
int prime[P],sum[LIM],last[LIM<<1],cnt=0;
ll val[LIM*2],f[LIM*2];
int tot;
void init(){
    tot=0;
    for(int i=2;i<=sqr;i++){
        if(!vis[i]){
            prime[++tot]=i;
            for(int j=2;j*i<=sqr;j++){
                vis[i*j]=true;
            }
        }
    }
}
ll n;
int main(){
    scanf("%lld",&n);
    sqr=sqrt(n);
    init();
    for(ll i=1;i<=n;i=n/(n/i)+1){
        val[++cnt]=n/i;
    }
    reverse(val+1,val+1+cnt);
    for(int i=1;i<=cnt;i++)f[i]=val[i];
    for(register int i=1;i<=tot;i++){
        for(register int j=cnt;j>0;j--){
            ll k=val[j]/prime[i],pos=k<=sqr?k:cnt+1-n/k;
            if(k<prime[i])break;
            f[j]-=f[pos]+last[pos]-i+1;
            last[j]=i;
        }
    }
    printf("%lld\n",tot+f[cnt]-1);
}

求小于等于n的素数和。

#include <bits/stdc++.h>
using namespace std;
const int mod=998244353;
const int inv2=(mod+1)/2;
typedef long long ll;
const int LIM=400000,P=30000;
ll sqr;
bool vis[LIM];
int prime[P],sum[LIM],last[LIM<<1],cnt;
ll val[LIM*2],f[LIM*2];
int tot;
void init(){
    tot=cnt=0;
    memset(vis,false,sizeof(vis));
    for(int i=2;i<=sqr;i++){
        if(!vis[i]){
            prime[++tot]=i;
            sum[tot]=(sum[tot-1]+i)%mod;
            for(int j=2;j*i<=sqr;j++){
                vis[i*j]=true;
            }
        }
    }
    memset(last,0,sizeof(last));
}
ll n;
int main(){
    while(~scanf("%lld",&n)){
        sqr=sqrt(n);
        init();
        for(ll i=1;i<=n;i=n/(n/i)+1){
            val[++cnt]=n/i;
        }
        reverse(val+1,val+1+cnt);
        for(int i=1;i<=cnt;i++){
            f[i]=val[i]%mod*((val[i]+1)%mod)%mod*inv2%mod;
        }
        for(int i=1;i<=tot;i++){
            for(int j=cnt;j>0;j--){
                ll k=val[j]/prime[i],pos=k<=sqr?k:cnt+1-n/k;
                if(k<prime[i])break;
                f[j]-=prime[i]*(f[pos]+sum[last[pos]]-sum[i-1]);
                f[j]%=mod;
                if(f[j]<0)f[j]+=mod;
                last[j]=i;
            }
        }
        printf("%lld\n",(sum[tot]+f[cnt]-1+mod)%mod);
    }
}