參考：https://www.cnblogs.com/SilverNebula/p/7045199.html
所是反演其實(shí)反演作用不大，又是一道做起來(lái)感覺(jué)詭異的題
轉(zhuǎn)成前綴和相減的形式
\[ \sum_{i=1}^{n}\sum_{j=1}^{n}[\frac{i*j}{gcd(i,j)}\leq n] \]
\[ \sum_{d=1}^{n}\sum_{i=1}^{\left \lfloor \frac{n}ozvdkddzhkzd\right \rfloor}\sum_{j=1}^{\left \lfloor \frac{n}ozvdkddzhkzd\right \rfloor}[gcd(i,j)==1][i*j\leq\left \lfloor \frac{n}ozvdkddzhkzd \right \rfloor] \]
\[ \sum_{k=1}^{n} \mu(k)\sum_{d=1}^{\left \lfloor \frac{n}{k} \right \rfloor}\sum_{i=1}^{\left \lfloor \frac{n}{dk} \right \rfloor}\sum_{j=1}^{\left \lfloor \frac{n}{dk} \right \rfloor}[i*j*d\leq\left \lfloor \frac{n}{k^2} \right \rfloor] \]
然后是非常神奇的縮小范圍……
\[ \sum_{k=1}^{\sqrt{n}}\mu(k)\sum_{d=1}^{\left \lfloor \frac{n}{k^2} \right \rfloor}\sum_{i=1}^{\left \lfloor \frac{n}{dk^2} \right \rfloor}\sum_{j=1}^{\left \lfloor \frac{n}{dk^2} \right \rfloor}[i*j*d\leq\left \lfloor \frac{n}{k^2} \right \rfloor] \]
然后對(duì)于這個(gè)友好的范圍直接枚舉就可以了。

#include<iostream> #include<cstdio> #include<cmath> using namespace std; const int N=1000005,m=1000000; int q[N],mb[N],tot; long long a,b; bool v[N]; long long wk(long long n) {if(!n)return 0;long long re=0ll,tmp=0ll,a=sqrt(n);for(long long k=1;k<=a;k++)if(mb[k]){tmp=0;long long b=n/k/k;for(long long i=1;i*i*i<=b;i++){for(long long j=i+1;j*j*i<=b;j++)tmp+=(b/(i*j)-j)*6+3;tmp+=(b/(i*i)-i)*3;tmp++;}re+=mb[k]*tmp;}return (re+n)/2; } int main() {mb[1]=1;for(int i=2;i<=m;i++){if(!v[i]){q[++tot]=i;mb[i]=-1;}for(int j=1;j<=tot&&i*q[j]<=m;j++){int k=i*q[j];v[k]=1;if(i%q[j]==0){mb[k]=0;break;}mb[k]=-mb[i];}}scanf("%lld%lld",&a,&b);printf("%lld\n",wk(b)-wk(a-1));return 0; }

轉(zhuǎn)載于:https://www.cnblogs.com/lokiii/p/8324776.html

與50位技術(shù)專家面對(duì)面20年技術(shù)見(jiàn)證，附贈(zèng)技術(shù)全景圖

總結(jié)

以上是生活随笔為你收集整理的51nod 1222 最小公倍数计数【莫比乌斯反演】的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。