Uoj #188. 【UR #13】Sanrd(min_25筛)
题目大意
现在有一函数 f ( x ) f(x) f(x)其函数值为x的次大素因数,如为素数则为0,现求 ∑ i = l r f ( i ) \sum_{i=l}^rf(i) ∑i=lrf(i)
解题思路
注意到min_25筛的更新中有一步为
g
(
n
,
m
)
=
∑
m
<
p
≤
n
,
p
e
≤
n
,
e
≥
1
,
p
∈
P
f
(
p
e
)
(
[
e
>
1
]
+
g
(
⌊
n
p
e
,
p
⌋
)
)
+
h
(
n
)
−
h
(
m
)
g(n,m)=\sum_{m<p\le \sqrt{n},p^e\le n,e\ge1,p\in P}f(p^e)\left([e>1]+g(\lfloor\frac{n}{p^e},p\rfloor)\right)+h(n)-h(m)
g(n,m)=m<p≤n,pe≤n,e≥1,p∈P∑f(pe)([e>1]+g(⌊pen,p⌋))+h(n)−h(m)
其中h函数设为到因变量为止的素数的数量,那么
h
(
n
)
−
h
(
m
)
h(n)-h(m)
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));
}