尝试理解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);
}
}