题解：洛谷P1891 疯狂LCM-CSDN博客

原题链接

题目描述

描述：

众所周知，czmppppp是数学大神犇。一天，他给众蒟蒻们出了一道数论题，蒟蒻们都惊呆了。。。

给定正整数N，求LCM(1,N)+LCM(2,N)+...+LCM(N,N)。

输入输出格式

输入格式：

第一行一个数T，表示有T组数据。

对于每组数据，一行，一个正整数N。

输出格式：

T行，每行为对应答案。

输入输出样例

输入样例#1： 复制

输出样例#1： 复制

1
4
55

说明

对于30%的数据，1≤T≤5，1≤N≤100000

对于100%的数据，1≤T≤300000，1≤N≤1000000

题解

挺妙的一道题。

要求的是\(\sum_{i=1}^{n}lcm(i, n)\)，一看好像没什么思路（可能是本人菜）。~~当然只会想到暴力打表和看题解了……?~~

首先这步应该~~地球人~~都会去试：
\[\sum_{i=1}^{n}lcm(i, n) = n \sum_{i=1}^{n}\frac{i}{\gcd(i, n)}\]
然后很关键的一步是像证明莫反一样添一个\(\sum\)（说一句废话），式子就变成了：
\[n\sum_{d\mid n}\sum_{i=1}^{n}\frac{i}{d}[d=\gcd(i, n)]\]
中括号里是一个类似布尔型的变量。

对于中括号里的式子，两边同除以\(d\)：
\[n\sum_{d\mid n}\sum_{i=1}^{n}\frac{i}{d}[\gcd(\frac{i}{d}, \frac{n}{d}) = 1]\]
把\(\frac{i}{d}\)变成\(i\)，得到：
\[n\sum_{d\mid n}\sum_{i=1}^{\frac{n}{d}}i[\gcd(i, \frac{n}{d}) = 1]\]
把\(d\)倒着枚举，我们就得到了酱紫的东西：
\[n\sum_{d\mid n}\sum_{i=1}^{d}i[\gcd(i,d) = 1]\]
然后我们来看\(\sum_{i=1}^{d}i[\gcd(i,d) = 1]\)：与\(d\)互质的所有数之和。