【题解】SDOI-2012 Longge的问题

转载注：
转载自PinkRabbit的博客园
原文地址
以下是转载内容：

题目传送门：链接。

能自己推出正确的式子的感觉真的很好！

题意简述：

求$\sum_{i=1}^n gcd(i,n)$。$n≤2^{32}$。

题解：

我们开始化简式子：

$\sum_{i=1}^ngcd(i,n)$

$=\sum_{j=1}^n (j×\sum_{i=1}^n[gcd(i,n)=j])$

$=\sum_{j=1}^n(j×\sum_{i=1}^n[gcd(i/j,n/j)=1](j|i,j|n))$

$=\sum_{j=1}^n(j×\phi(n/j)(j|n))$

$=\sum_{j|n}(j×\phi(n/j))$

到这里就可以直接计算了。

但是还可以进一步化简！（以下的p为质数）

$\sum_{j|n}(j×\phi(n/j))$

$=\sum_{j|n}(n/j×\phi(j))$

$=\sum{j|n}(n/j×(j⋅\prod_{p|j} \frac {p−1} p))$

$=\sum{j|n}(n⋅\prod_{p|j} \frac {p−1} p)$

$=n×\sum_{j|n}\prod_{p|j} \frac {p−1} p$

接下来我们令$n=p_1^{c_1}p_2^{c_2}p_3^{c_3}⋯p_k^{c_k}$，并定义$f_i=\frac {p_i−1} {p_i}$。

那么n的因子j可以表示为：$j=p_1^{c_1}p_2^{c_2}p_3^{c_3}⋯p_k^{c_k}$，满足$0≤c_i≤b_i$。

那么$\prod_{p|j}\frac {p−1} p=\prod_{i=1}^kfi[ci>0]$。

我们观察一类$\prod_{i=1}^kf_i[c_i>0]$相等的j，它们必要满足在$i$相等的情况下，$c_i$同时大于$0$或$c_i$同时等于$0$。

那么这一类的$j$有多少个呢？如果这类j有质因子$p_{q_1},p_{q_2},p_{q_3},⋯,p_{q_g}$。

那么这类$j$的答案为$\prod_{i=1}^gf_{q_i}$，而个数为$\prod_{i=1}^gb_{q_i}$。

$b_i$就是原来n的质因数分解的指数。

那么对答案的贡献为：$\prod_{i=1}^g\chi_{q_i}$。这里$\chi_i=f_i⋅b_i$。

发现每一个质因子的贡献都是独立的，那么最后我们枚举n的每一个质因子取不取，得到最后的答案：$n⋅\prod _{i=1}^k(\chi_i+1)$。

举个例子：如果n只有3个质因子，那么答案为$n⋅(1+\chi_1+\chi_2+\chi_3+\chi_1\chi_2+\chi_1\chi3+\chi_2\chi_3+\chi_1\chi_2\chi_3)$。

显然可以化简为：$n⋅(\chi_1+1)⋅(\chi_2+1)⋅(\chi_3+1)$。

当然可以类比到质因数更多的情况。

总之，答案就是：$n⋅\prod_{i=1}^k \frac {b_ip_i−bi+pi} {p_i}$。

代码：

#include<cstdio>
long long n;
long long f(){
    long long ans=n; long long i;
    for(i=2;i*i<=n;++i) if(n%i==0){
        int b=0;
        while(n%i==0) ++b,n/=i;
        ans/=i;
        ans*=b*i-b+i;
    } if(n>1) ans/=n, ans*=n+n-1; 
    return ans;
}
int main(){
    scanf("%lld",&n);
    printf("%lld",f());
    return 0;
}