[莫比乌斯反演 复杂度分析] 51Nod 1222 最小公倍数计数

这么反演一通

转化为abc<=n的解数 [暴力]
HDU 4473 Exam 

a 以立方根的复杂度枚举一通再讨论下就好了

既然糖老师强调复杂度分析 那就来分析分析

对于子问题  abc<=n

复杂度是n的2/3次的

再套一重外循环是n的5/6次？

肯定是我分析错了

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include <tr1/unordered_map>
typedef long long ll;
using namespace std;
using namespace std::tr1;

const int maxn=10000000;

int prime[1000000],num;
int vst[maxn+5],miu[maxn+5];

inline void Pre(){
miu[1]=1;
for (int i=2;i<=maxn;i++){
if (!vst[i]) prime[++num]=i,miu[i]=-1;
for (int j=1;j<=num && (ll)i*prime[j]<=maxn;j++){
vst[i*prime[j]]=1;
if (i%prime[j]==0){
miu[i*prime[j]]=0;
break;
}
miu[i*prime[j]]=miu[i]*miu[prime[j]];
}
}
}

inline ll calc(ll n){
ll ans=0;
for(ll a=1;a*a*a<=n;a++)
for(ll b=a;a*b*b<=n;b++){
ll c=n/(a*b);
if(c<b) break;
if(a==b)
ans+=(c-b)*3+1;
else
ans+=(c-b)*6+3;
}
return ans;
}

inline ll Solve(ll n){
if (!n) return 0;
ll ret=0;
for (ll k=1;k*k<=n;k++)
if (miu[k])
ret+=miu[k]*calc(n/k/k);
ret=(ret+n)/2;
return ret;
}

int main(){
ll l,r;
freopen("t.in","r",stdin);
freopen("t.out","w",stdout);
scanf("%lld%lld",&l,&r); Pre();
printf("%lld\n",Solve(r)-Solve(l-1));
return 0;
}