Description
求
∑i=lrμ(i)
l,r<=10^10
Solution
设
M(n)=∑i=1nμ(i)
我们知道,
∑d|nμ(d)=[n=1]
那么
1=∑i=1n∑d|iμ(d)
=∑T=1n∑d|Tμ(d)
=∑i=1n∑d=1⌊ni⌋μ(d)
∑i=1nM(⌊ni⌋)
于是,
M(n)=1−∑i=2nM(⌊ni⌋)
后面的东西可以用分块来加速。
然后打上记忆化标记。
或者可以先预处理出一段的前缀和,这样会快一点。
Code
#include<cstdio>
#include<cstring>
#include<algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define rep(i,a) for(int i=last[a];i;i=next[i])
#define N 5000000
using namespace std;
typedef long long ll;
const int mo=6666666;
int mu[N+5],p[N+5],cnt;
int last[mo],next[mo],v[mo];
ll t[mo],l,r;
bool bz[N+5];
void add(int x,ll y,int z) {
t[++cnt]=y;v[cnt]=z;next[cnt]=last[x];last[x]=cnt;
}
int calc(ll x) {
if (x<=N) return mu[x];int ans=1,k=x%mo;
rep(i,k) if (t[i]==x) return v[i];
for(ll l=2,r;l<=x;l=r+1) r=x/(x/l),ans-=calc(x/l)*(r-l+1);
add(k,x,ans);
return ans;
}
int main() {
fo(i,2,N) {
if (!bz[i]) p[++p[0]]=i,mu[i]=-1;
fo(j,1,p[0]) {
int k=i*p[j];if (k>N) break;
bz[k]=1;if (!(i%p[j])) break;
mu[k]=-mu[i];
}
}
mu[1]=1;fo(i,1,N) mu[i]+=mu[i-1];
scanf("%lld%lld",&l,&r);printf("%d",calc(r)-calc(l-1));
}