[LOJ6202][Min_25]叶氏筛法

LOJ6202

Min_25模板题
Min_25是一种能快速求出积性函数前缀和的东西
首先定义$f(i)$表示${\sum }_{j=1}^{i}f(j)$，P为质数集合，pri[i]为第i个质数,mnpri[i]为i的最小质因子
$g(i,j)$表示$\sum f(k)[k\in P||mnpri[k]>pri[j]]$
其实就是所有质数以及埃氏筛j次之后未被筛掉的数的f的和
考虑递推g，如果当前质数$pri[i{]}^2>n$那不会被筛掉任何数，直接递推
否则就是有质因子pri[i]的被筛去，而且除去pri[i]的最小质因子大于pri[i]
所以减去$f(pri[i])\ast g(\frac{n}{pri[i]},i-1)$
这时候我们把pri[1]-pri[i-1]这些质数的f减去了，所以再加上去
然后g就可以递推了
我们要求的是$g(n,|P|)$，所以用滚动数组优化空间
然后因为大于$\sqrt{n}$的质因子只有一个，所以可以开两个数组表示$g(i,j)$和$g(\frac{n}{i},j)$

我们还可以处理合数的答案，这个在另外的题里会体现

Code：

#include<bits/stdc++.h>
#define ll long long
#define db double
using namespace std;
inline ll read(){
	ll res=0,f=1;char ch=getchar();
	while(!isdigit(ch)) {if(ch=='-') f=-f;ch=getchar();}
	while(isdigit(ch)) {res=(res<<1)+(res<<3)+(ch^48);ch=getchar();}
	return res*f; 
}
const int N=340005;
int pt[N],pri[N>>2],tot=0;
ll sum[N];
inline void init(){
	for(int i=2;i<=N-5;i++){
		sum[i]=sum[i-1];
		if(!pt[i]) pri[++tot]=i,sum[i]+=i;
		for(int j=1;j<=tot && i*pri[j]<=N-5;j++){
			pt[i*pri[j]]=1;
			if(i%pri[j]==0) break;
		}
	} 
}
inline db S(ll x){return (db)x*(x+1)/2;}
int cnt,t1[N],t2[N];
ll last[N<<1];
db h[N<<1];
inline db calc(ll x){
	cnt=0,memset(t1,0,sizeof(t1));memset(t2,0,sizeof(t2));memset(h,0,sizeof(h));memset(last,0,sizeof(last)); 
	int id;ll l,r,t;
	for(l=1;l<=x;l=r+1){
		t=(x/l);r=x/t;last[++cnt]=t;h[cnt]=S(t)-1;
		(t>=N-5?t2[x/t]:t1[t])=cnt; 
	}
	for(int j=1;j<=tot;j++){
		ll pw=(ll)pri[j]*pri[j];
		if(pw>x) break;
		for(int i=1;i<=cnt && last[i]>=pw;i++){
			t=last[i]/pri[j],id=t>=N-5?t2[x/t]:t1[t];
			h[i]-=(db)pri[j]*(h[id]-sum[pri[j]-1]);
		}
	}
	return h[x>=N-5?t2[1]:t1[x]];
}
int main(){
	init();
	ll l=read(),r=read();
	printf("%.f",calc(r)-calc(l-1));
	return 0;
}