正題

題目鏈接:https://www.luogu.com.cn/problem/P1829

---

題目大意

給出$n,m$求$$\sum _{i=1}^n\sum _{j=1}^{m}lcm(i,j)$$

---

解題思路

$$\sum _{i=1}^n\sum _{j=1}^m\frac{ij}{gcd(i,j)}$$
$$\sum _{x=1}^{n}x\times (\sum _{i=1}^{?\frac{n}{x}?}\sum _{j=1}^{?\frac{m}{x}?}ij\times [gcd(i,j)==1])$$
$$\sum _{x=1}^{n}x\times \sum _{x|d}\mu (\frac{d}{x})d^2?\sum _{i=1}^{?\frac{n}{x}?}\sum _{j=1}^{?\frac{m}{x}?}ij$$
單獨拎出后面的來看
$$\sum _{x|d}\mu (\frac{d}{x})d^2?\sum _{i=1}^{?\frac{n}{x}?}\sum _{j=1}^{?\frac{m}{x}?}ij$$
$$\sum _{x|d}\mu (\frac{d}{x})d^2?\sum _{i=1}^{?\frac{n}{d}?}i?\sum _{j=1}^{?\frac{m}{d}?}j$$
然后整除分塊處理這一個
之后再在外面套一個整除分塊即可。

時間復雜度$O(n+m)$

---

$code$

#include<cstdio> #include<cstring> #include<algorithm> #define ll long long using namespace std; const ll N=1e7+10,XJQ=20101009; ll n,m,mu[N],pri[N],g[N],cnt,ans; bool v[N]; void ycl(ll n){mu[1]=1;for(ll i=2;i<=n;i++){if(!v[i])pri[++cnt]=i,mu[i]=-1;for(ll j=1;j<=cnt&&i*pri[j]<=n;j++){v[i*pri[j]]=1;if(i%pri[j]==0)break;mu[i*pri[j]]=-mu[i];}}for(ll i=1;i<=n;i++)mu[i]=(mu[i-1]+mu[i]*i*i%XJQ)%XJQ;return; } ll solve(ll n,ll m){ll ans=0;for(ll l=1,r;l<=n;l=r+1){r=min(n/(n/l),m/(m/l));ll q=n/l,p=m/l;ans=(ans+(q*(q+1)/2%XJQ)*(p*(p+1)/2%XJQ)%XJQ*(mu[r]-mu[l-1]+XJQ)%XJQ)%XJQ;}return ans; } int main() {scanf("%lld%lld",&n,&m);if(n>m)swap(n,m);ycl(m);for(ll l=1,r;l<=n;l=r+1){r=min(n/(n/l),m/(m/l));(ans+=(r+l)*(r-l+1)/2%XJQ*solve(n/l,m/l)%XJQ)%=XJQ;}printf("%lld\n",ans); }

總結(jié)

以上是生活随笔為你收集整理的P1829-[国家集训队]Crash的数字表格/JZPTAB【莫比乌斯反演】的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

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