Uoj #188. 【UR #13】Sanrd（min_25筛）

Uoj #188. 【UR #13】Sanrd（min_25筛）

题目大意

现在有一函数$f(x)$其函数值为x的次大素因数，如为素数则为0，现求${\sum }_{i=l}^{r}f(i)$

解题思路

注意到min_25筛的更新中有一步为
$$g(n,m)=\sum _{m<p\le \sqrt{n},p^e\le n,e\ge 1,p\in P}f(p^e)\left([e>1]+g(\lfloor \frac{n}{p^e},p\rfloor )\right)+h(n)-h(m)$$
其中h函数设为到因变量为止的素数的数量，那么$h(n)-h(m)$就是n和m之间素数的数量，n与m之间素数的数量乘上上一层递归中的素数即为这一段的答案。对l-1和r分别求前缀和，做差即为答案

AC代码

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int size=2e6+5;
LL n;int s;
int tol,tot;
bool prime[size];int p[size];
LL num[size];
LL h[size];
inline int ID(LL x){return x<=s?x:tot-n/x+1;}
LL get_g(LL x,int k)
{
	if(p[k]>x||x<=1) return 0;
	LL ans=h[ID(x)]-(k-1);
	ans=ans*p[k-1];
	while(1LL*p[k]*p[k]<=x&&k<=tol)
	{
		LL t=p[k];
		for(int e=1;t*p[k]<=x;t=t*p[k],e++)
		{
			ans+=get_g(x/t,k+1)+p[k];
		}
		k++;
	}
	return ans;
}
LL get(LL sce)
{
	n=sce;
	s=sqrt(n);
	tol=0;tot=0;
	while(1LL*s*s<=n) s++;
	while(1LL*s*s>n) s--;
	for(int i=1;i<=s;i++) prime[i]=true;
	for(int i=2;i<=s;i++)
	{
		if(prime[i]) p[++tol]=i;
		for(int j=1;j<=tol&&i*p[j]<=s;j++)
		{
			prime[i*p[j]]=false;
			if(i%p[j]==0) break;
		}
	}
	p[tol+1]=0;
	for(int i=1;i<=s;i++) num[++tot]=i;
	for(int i=s;i>=1;i--) if(n/i>s) num[++tot]=n/i;
	for(int i=1;i<=tot;i++) h[i]=num[i]-1;
	int x=1;
	for(int i=1;i<=tol;i++)
	{
		while(num[x]<1LL*p[i]*p[i]) x++;
		for(int j=tot;j>=x;j--)
		{
			h[j]-=h[ID(num[j]/p[i])]-(i-1);
		}
	}
	return get_g(n,1);
}
int main()
{
	LL l,r;
	scanf("%lld%lld",&l,&r);
	printf("%lld\n",get(r)-get(l-1));
}