日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

bzoj2226-[Spoj5971]LCMSum【欧拉函数,GCD】

發布時間:2023/12/3 编程问答 36 豆豆
生活随笔 收集整理的這篇文章主要介紹了 bzoj2226-[Spoj5971]LCMSum【欧拉函数,GCD】 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

正題

題目鏈接:
https://www.lydsy.com/JudgeOnline/problem.php?id=2226
https://www.luogu.org/problem/SP5971


題目大意

∑i=1nlcm(n,i)\sum_{i=1}^n lcm(n,i)i=1n?lcm(n,i)


解題思路

∑i=1nlcm(n,i)\sum_{i=1}^n lcm(n,i)i=1n?lcm(n,i)
∑i=1nnigcd(n,i)\sum_{i=1}^n \frac{ni}{gcd(n,i)}i=1n?gcd(n,i)ni?
n∑d∣n∑i=1n[(n,i)==d]in\sum_{d\mid n}\sum_{i=1}^n[(n,i)==d]indn?i=1n?[(n,i)==d]i
n∑d∣n∑i=1n[(n/d,i/d)==1]in\sum_{d\mid n}\sum_{i=1}^n[(n/d,i/d)==1]indn?i=1n?[(n/d,i/d)==1]i
n∑d∣nd?φ(d)2n\sum_{d|n}\frac{d*\varphi(d)}{2}ndn?2d?φ(d)?

然后預處理O(nlog?n)O(n\log n)O(nlogn)然后O(1)O(1)O(1)回答即可

時間復雜度O(nlog?n+T)O(n\log n+T)O(nlogn+T)


codecodecode

#include<cstdio> #include<cstring> #include<algorithm> #define ll long long using namespace std; const ll N=1e6; ll T,phi[N+10],ans[N+10],n; int main() {for (ll i=2;i<=N;i++) phi[i]=i;for (ll i=2;i<=N;i++){bool flag=(phi[i]==i);for (ll j=i;j<=N;j+=i){if(flag)phi[j]=phi[j]/i*(i-1);ans[j]+=phi[i]*i/2;}}scanf("%lld",&T);while(T--){scanf("%lld",&n);printf("%lld\n",n*(ans[n]+1));} }

總結

以上是生活随笔為你收集整理的bzoj2226-[Spoj5971]LCMSum【欧拉函数,GCD】的全部內容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。