SP5971 LCMSUM - LCM Sum

SP5971 LCMSUM - LCM Sum

思路

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

=>∑i=1ningcd(i,n)=> \sum_{i = 1}^{n}\frac{i n}{gcd(i, n)}=>i=1ngcd(i,n)in

=>n∑i=1nigcd(i,n)=> n\sum_{i = 1}^{n}\frac{i}{gcd(i, n)}=>ni=1ngcd(i,n)i

我们按照P2303 [SDOI2012] Longge的思路枚举gcd(i,n)gcd(i, n)gcd(i,n)

=>n∑d∣n∑i=1nid(gcd(i,n)==1)=>n\sum_{d\mid n} \sum_{i = 1}^{n} \frac{i}{d}(gcd(i, n) == 1)=>ndni=1ndi(gcd(i,n)==1)

=>n∑d∣n∑i=1ndi(gcd(i,nd)==1)=>n\sum_{d\mid n}\sum_{i = 1}^{\frac{n}{d}} i(gcd(i, \frac{n}{d}) == 1)=>ndni=1dni(gcd(i,dn)==1)

这个式子就熟悉了∑i=1ni(gcd(i,n)==1)=nϕ(n)2\sum_{i = 1} ^{n} i(gcd(i, n) == 1) = \frac{n\phi(n)}{2}i=1ni(gcd(i,n)==1)=2nϕ(n),给定一个整数nnn小于nnn的数并且与nnn互质的数的和是这个式子,所以上式变成

=>=>n∑d∣ndϕ(d)2=>=>n\sum_{d\mid n}\frac{d\phi(d)}{2}=>=>ndn2dϕ(d)

所以我们只要枚举nnn的约数就行了。

这里还存在一个问题,对于gcd=d=1gcd = d = 1gcd=d=1的时候,我们统计的答案是没用贡献的,所以我们还要加上lcm(1,n)=nlcm(1, n) = nlcm(1,n)=n得到我们最后的答案,

这里我整体复杂度是O(n)+TnO(n) + T\sqrt nO(n)+Tn,先用素数筛筛选出所有的ϕ(i)\phi(i)ϕ(i),然后再每次计算答案。

第一道自己推出来的数学题,难得啊!!!

代码

/*Author : lifehappy
*/
#pragma GCC optimize(2)
#pragma GCC optimize(3)
// #include <iostream>
// #include <algorithm>
// #include <stdlib.h>
// #include <cmath>
// #include <vector>
// #include <cstdio>
#include <bits/stdc++.h>
#define mp make_pair
#define pb push_back
#define endl '\n'using namespace std;typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int> pii;const double pi = acos(-1.0);
const double eps = 1e-7;
const int inf = 0x3f3f3f3f;inline ll read() {ll f = 1, x = 0;char c = getchar();while(c < '0' || c > '9') {if(c == '-') f = -1;c = getchar();}while(c >= '0' && c <= '9') {x = (x << 1) + (x << 3) + (c ^ 48);c = getchar();}return f * x;
}void print(ll x) {if(x < 10) {putchar(x + 48);return ;}print(x / 10);putchar(x % 10 + 48);
}const int N = 1e6 + 10;int eular[N], n;vector<int> prime;bool st[N];void init() {st[0] = st[1] = 1;eular[1] = 1;for(int i = 2; i < N; i++) {if(!st[i]) {prime.pb(i);eular[i] = i - 1;}for(int j = 0; j < prime.size() && i * prime[j] < N; j++) {st[i * prime[j]] = 1;if(i % prime[j]) {eular[i * prime[j]] = eular[i] * (prime[j] - 1);}else {eular[i * prime[j]] = eular[i] * prime[j];break;}}}
}int main() {// freopen("in.txt", "r", stdin);// freopen("out.txt", "w", stdout);// ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);init();int T = read();while(T--) {int n = read();ll ans = 0;for(int i = 1; i * i <= n; i++) {if(n % i == 0) {ans += 1ll * i * eular[i] / 2;if(i * i != n) {ans += 1ll * n / i * eular[n / i] / 2;}}cout << ans << endl;}printf("%lld\n", 1ll * n * ans + n);}return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/314433.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

程序员修神之路--用NOSql给高并发系统加速

领取福利记得长按&#xff0c;领取技术书籍哦随着互联网大潮的到来&#xff0c;越来越多网站&#xff0c;应用系统需要海量数据的支撑&#xff0c;高并发、低延迟、高可用、高扩展等要求在传统的关系型数据库中已经得不到满足&#xff0c;或者说关系型数据库应对这些需求已经显…

欧拉函数的性质及其证明

欧拉函数 ppp是素数&#xff0c;则有ϕ(p)p−1\phi(p) p - 1ϕ(p)p−1 证明&#xff1a;显然。 ppp是素数&#xff0c;npkn p ^ knpk&#xff0c;则ϕ(n)pk−pk−1\phi(n) p ^ k - p ^ {k - 1}ϕ(n)pk−pk−1 证明&#xff1a; [1,n][1, n][1,n]内&#xff0c;ppp的约数有p…

AtCoder Regular Contest 061

文章目录C - Many FormulasD - Snukes ColoringE - Snukes Subway TripF - Card Game for Three传送门 C - Many Formulas Score : 300300300 points 爆搜 直接dfsdfsdfs爆搜即可 复杂度O(2n)O(2^n)O(2n) 代码 D - Snuke’s Coloring Score : 400400400 points 考虑有标…

限时团购,6.5折:《C# 7.0 核心技术指南》

大家好&#xff0c;经过近两年的翻译&#xff0c;《C# 7.0 核心技术指南》终于和大家见面了。全书由 ThoughtWorks 高级咨询师&#xff0c;资深 .NET 专家刘夏翻译。作为一本第七次再版的图书&#xff0c;此次翻译对书中的字句进行了重新整理。期间和图书的原作者 Joe Albahari…

小A的最短路

小A的最短路 思路 树上问题求两个点的最短距离&#xff0c;显然能用lcalcalca来进行lognlog_nlogn​的查询&#xff0c;引入了两个无边权的点&#xff0c;所以我们的路劲就可以规划成三种x−>y&#xff0c;x−>u−>v−>y&#xff0c;x−>v−>u>−yx -&g…

Codeforces Beta Round #11 B. Jumping Jack 思维

文章目录题意&#xff1a;思路&#xff1a;传送门 题意&#xff1a; 初始在000点&#xff0c;第iii次跳跃可以跳iii的距离&#xff0c;可以选择向右或者向左跳&#xff0c;求最少多少次可以跳到xxx点。 −1e9≤x≤1e9-1e9\le x\le 1e9−1e9≤x≤1e9 思路&#xff1a; 我们贪…

Azure 命令行工具大混战,都是什么,该选哪个?

点击上方蓝字关注“汪宇杰博客”导语最近在学习 Azure 的命令行玩法&#xff0c;发现官方有不止一种命令行工具&#xff0c;容易对新手产生混淆&#xff0c;本文将介绍各种工具都是干啥的&#xff0c;以及如何选择。目前&#xff0c;微软官方有3个Azure命令行工具&#xff0c;分…

HDU 2588 GCD(欧拉函数)

GCD 思路 题目要求&#xff0c;对于给定的n,mn, mn,m要求有多少数∑i1ngcd(i,n)>m\sum _{i 1} ^{n} gcd(i, n) > m∑i1n​gcd(i,n)>m 我们可以对这个式子进行化简&#xff0c;通过枚举dgcd(i,n)d gcd(i, n)dgcd(i,n)有 ∑d∣n∑i1ngcd(i,d)d\sum _{d \mid n} \su…

2021 ccpc 哈尔滨 G. Damaged Bicycle 状压 + 期望dp

文章目录题意&#xff1a;思路&#xff1a;传送门 题意&#xff1a; 你需要从111走到nnn&#xff0c;初始速度是ttt&#xff0c;某些地方有自行车&#xff0c;每个位置自行车有pip_ipi​的概率是坏掉的&#xff0c;如果自行车没坏可以骑上自行车&#xff0c;速度是rrr&#x…

HDU 4983 Goffi and GCD(欧拉函数)

Goffi and GCD 思路 题目要求∑i1n∑j1ngcd(n−i,n)gcd(n−j,n)nk\sum_{i 1} ^{n} \sum _{j 1} ^{n} gcd(n - i, n)gcd(n - j, n) n ^ {k}∑i1n​∑j1n​gcd(n−i,n)gcd(n−j,n)nk 显然有gcd(n−i,n)<ngcd(n - i, n) < ngcd(n−i,n)<n对于k>3k > 3k>3直…

揭秘鸿蒙生态背后的DevOps实践

&#xff08;图片来源于网络&#xff09;8月9日&#xff0c;华为发布了鸿蒙操作系统&#xff0c;在发布会上我们看到了鸿蒙系统的研发历程&#xff1a;2017年&#xff0c;鸿蒙内核1.0完成技术验证&#xff1b;2018年&#xff0c;鸿蒙内核2.0用于终端TEE&#xff1b;2019年&…

AT3557 Four Coloring 切比雪夫距离 + 四色构造

传送门 由于曼哈顿距离在图上显示的是一个棱形&#xff0c;并不是很好看&#xff0c;所以我们将其旋转454545&#xff0c;转换成切比雪夫距离&#xff0c;这样就变成了一个正方形&#xff0c;正方形内部的点距离都不超过ddd&#xff0c;此时可以将正方形内部的点看成一个整体&a…

51nod 1188 最大公约数之和 V2(欧拉函数)

1188 最大公约数之和 V2 思路 用欧拉函数可以化简式子如下 ∑i1n∑j1i−1gcd(i,j)\sum_{i 1} ^{n} \sum _{j 1} ^{i - 1} gcd(i, j)i1∑n​j1∑i−1​gcd(i,j) ∑i1n∑j1igcd⁡(i,j)−(n1)(n)2 \sum_{i 1} ^{n} \sum_{j 1} ^{i} \gcd(i, j) - \frac{(n 1)(n)}{2}i1∑n​…

SonarQube系列一、Linux安装与部署

来源&#xff1a;https://www.cnblogs.com/7tiny/p/11269774.html【前言】随着项目团队规模日益壮大&#xff0c;项目代码量也越来越多。且不说团队成员编码水平层次不齐&#xff0c;即便是老手&#xff0c;也难免因为代码量的增加和任务的繁重而忽略代码的质量&#xff0c;最终…

P2906 [USACO08OPEN]Cow Neighborhoods G 切比雪夫距离 + 并查集 + set

传送门 考虑将曼哈顿距离转换成切比雪夫距离&#xff0c;这样问题就变成了max(∣x1−x2∣,∣y1−y2∣)≤dmax(|x_1-x_2|,|y_1-y_2|)\le dmax(∣x1​−x2​∣,∣y1​−y2​∣)≤d&#xff0c;这个式子就很好看了&#xff0c;我们首先按照(x,y)(x,y)(x,y)排序&#xff0c;让后我…

2018-2019 ACM-ICPC, Asia Shenyang Regional Contest E. The Kouga Ninja Scrolls 切比雪夫距离 +线段树

传送门 将曼哈顿距离转换成切比雪夫距离&#xff0c;现在就是求max(∣x1−x2∣,∣y1−y2∣)max(|x_1-x_2|,|y_1-y_2|)max(∣x1​−x2​∣,∣y1​−y2​∣)&#xff0c;显然我们可以将x,yx,yx,y分开考虑&#xff0c;下面以xxx为例。 考虑一段区间内不同门派的最大值和最小值&am…

ASP.NET Core 框架本质学习

本文作为学习过程中的一个记录。学习文章地址&#xff1a;https://www.cnblogs.com/artech/p/inside-asp-net-core-framework.html一. ASP.NET Core 框架上的 Hello World程序public class Program{public static void Main()> new WebHostBuilder() .UseKestrel() …

牛客小白月赛12:月月给华华出题(欧拉函数)

月月给华华出题 思路 ∑i1nigcd(i,n)\sum_{i 1} ^{n} \frac{i}{gcd(i, n)}i1∑n​gcd(i,n)i​ ∑d∣n∑i1nid(gcd(i,d)d) \sum _{d \mid n} \sum_{i 1} ^{n} \frac{i}{d} (gcd(i, d) d)d∣n∑​i1∑n​di​(gcd(i,d)d) ∑d∣n∑i1ndi(gcd(i,d)1) \sum_{d\mid n} \sum_{i 1…

AtCoder Regular Contest 064

文章目录C - Boxes and CandiesD - An Ordinary GameE - Cosmic RaysF - Rotated PalindromesC - Boxes and Candies Score : 300300300 points 贪心 每次比较相邻两个&#xff0c;贪心的给最后一个加即可。 代码 D - An Ordinary Game Score : 500500500 points 博弈 结…

基于C#实现的轻量级多线程队列

工作中我们经常会遇到一些一些功能需要实现造作日志&#xff0c;数据修改日志&#xff0c;对于这种业务需求如果我们以同步的方式实现&#xff0c;难免会影响到系统的性能。如下我列出集中解决方案。使用Thread异步处理。使用线程池或Task异步处理。以上两种方案确实能解决我们…