51nod1244 莫比乌斯函数之和(杜教筛)
题面
题解
我……我忘记把预处理的块的大小调成\(n^{\frac{2}{3}}\)了……(仰天)
首先\(\mu*1=e\)
然后杜教筛就行了
//minamoto
#include<bits/stdc++.h>
#define R register
#define ll long long
#define IT map<ll,ll>::iterator
#define fp(i,a,b) for(R int i=a,I=b+1;i<I;++i)
#define fd(i,a,b) for(R int i=a,I=b-1;i>I;--i)
#define go(u) for(int i=head[u],v=e[i].v;i;i=e[i].nx,v=e[i].v)
using namespace std;
const int N=1e7+5;
bitset<N>vis;int p[N],mul[N];map<ll,ll>mp;
ll l,r;int m,sqr;
void init(int n){
mul[1]=1;
fp(i,2,n){
if(!vis[i])p[++m]=i,mul[i]=-1;
for(R int j=1;j<=m&&1ll*i*p[j]<=n;++j){
vis[i*p[j]]=1;
if(i%p[j]==0)break;
mul[i*p[j]]=-mul[i];
}
}
fp(i,2,n)mul[i]+=mul[i-1];
}
ll Mul(ll n){
if(n<=sqr)return mul[n];
IT it=mp.find(n);
if(it!=mp.end())return it->second;
ll res=1;
for(ll i=2,j;i<=n;i=j+1)
j=n/(n/i),res-=Mul(n/i)*(j-i+1);
return mp[n]=res;
}
int main(){
// freopen("testdata.in","r",stdin);
scanf("%lld%lld",&l,&r),init(sqr=N-5);
printf("%lld\n",Mul(r)-Mul(l-1));
return 0;
}
深深地明白自己的弱小