正題

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

---

題目大意

給出數(shù)列$A$求$$\sum _{i=1}^n\sum _{j=1}^{n}lcm(A_i,A_j)$$

---

解題思路

設(shè)$c_i$表示$A_j=i$的個數(shù)，然后答案就是（下面$n=5e4$）$$\sum _{i=1}^n\sum _{j=1}^{n}lcm(i,j)c_i{c}_j$$
$$\sum _{i=1}^n\sum _{j=1}^n\frac{ij}{gcd(i,j)}{c}_i{c}_j$$
$$\sum _{x=1}^{n}x\sum _{i=1}^{?\frac{n}{x}?}\sum _{j=1}^{?\frac{n}{x}?}[gcd(i,j)==1]ij{c}_{i?x}{c}_{j?x}$$
然后反演一下
$$\sum _{x=1}^n\sum _{x|d}\mu (\frac{d}{x})d\sum _{i=1}^{?\frac{n}{d}?}\sum _{j=1}^{?\frac{n}{d}?}ij{c}_{i?d}{c}_{j?d}$$
$$\sum _{d=1}^{n}d\sum _{x|d}\mu (\frac{d}{x})\sum _{i=1}^{?\frac{n}{d}?}\sum _{j=1}^{?\frac{n}{d}?}ij{c}_{i?d}{c}_{j?d}$$
然后兩部分分開處理就好了，時間復(fù)雜度$O(n\log n)$

---

$code$

#include<cstdio> #include<cstring> #include<algorithm> #define ll long long using namespace std; const ll N=5e4+10; ll n,cnt,pri[N],mu[N],c[N],g[N],f[N],ans; bool v[N]; void prime(){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]*mu[pri[j]];}}for(ll i=1;i<N;i++)for(ll j=i;j<N;j+=i)g[j]+=mu[i]*i;return; } int main() {prime();scanf("%lld",&n);for(ll i=1;i<=n;i++){ll x;scanf("%lld",&x);c[x]++;}for(ll i=1;i<N;i++)for(ll j=1;i*j<N;j++)f[i]+=c[i*j]*j;for(ll i=1;i<N;i++)ans+=f[i]*f[i]*i*g[i];printf("%lld\n",ans);return 0; }

總結(jié)

以上是生活随笔為你收集整理的P3911-最小公倍数之和【莫比乌斯反演】的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

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