51nod 1363 最小公倍数之和

51nod 1363

原题连接：
https://www.51nod.com/onlineJudge/questionCode.html#!problemId=1363

这个是一个不算很难的计算。

对于每次询问，给定$n$：

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

$$=\sum_{i=1}^n\frac{in}{gcd(i,n)}=n\sum_{i=1}^{n}\frac{i}{gcd(i,n)}$$

则有：

$$answer=n\sum_{i=1}^n\frac{i}{gcd(i,n)}\\=n\sum_{k|n}\big(\frac{1}{k}\sum_{i=1}^n[gcd(i,n)=k]*i\big)\\=n\sum_{k|n}\big(\sum_{i=1}^n[gcd(i,n)=k]*\frac{i}{k}\big)\\=n\sum_{k|n}\big(\sum_{i=1}^n[gcd(i,n)=k]*\frac{i}{k}\big)$$

因为这个时候 $\frac{i}{k}$ 与 $\frac{n}{k}$ 互素

则：

$$answer=n\sum_{k|n}\big(\sum_{i=1}^{\frac{n}{k}}[i\perp \frac{n}{k}]*i\big)\\=n\sum_{k|n}\big(\sum_{i=1}^{k}[i\perp k]*i\big)$$

记对于$1$ 到 $n$ ，与 $n$ 互素的 整数之和为$S(n)$：

$$S(n)=\sum_{i=1}^{n}[i\perp n]*i$$

当 有 $t\perp n$ 时， 则必然有： $(n-t)\perp n$

所以，当$n>1$ 时，每一个与$n$互素的整数 $t$ 都 必然存在一个整数 $u$ 有：

$$u\perp n\ \ 且\ \ t+u=n$$

所以$n>1$时：

$$S(n)=\frac{\phi(n)n}{2}$$

$n=1$时：

$$S(n)=1$$

所以：

$$answer=n\sum_{k|n}S(k)\\=n\big(1+\sum_{k|n\ ,k>1}\frac{\phi(k)k}{2}\big)\\=\frac{n}{2}\big(2+\sum_{k|n\ ,k>1}\phi(k)k\big)\\=\frac{n}{2}\big(1+\sum_{k|n}\phi(k)k\big)$$

积性函数的和函数依然是积性函数。

而积性函数的值可以由其 因式分解后的素数幂形式得来。

$f(n)=\phi(n)n$ 就是一个积性函数。

$f(p^k)=(p^k-p^{(k-1)})*p^k$ 其中，k>1.

那么其和函数素数幂形式是两个等比数列求和。

计算得：

$$F(p^k)=\sum_{k|p^k}f(k)\\=p\Big(\frac{p^{2k}-1}{p^2-1}p-\frac{p^{2k}-1}{p^2-1}\Big)+1$$

(注：貌似暴力分解不会超时。也算是暴力题了)

#include <stdio.h>
#include <algorithm>
#include <string.h>

using namespace std;
typedef long long LL;

const LL mod=1e9+7;

LL P(LL a,int b)
{
    LL tmp=1;
    while(b)
    {
        if(b&1)
            tmp=tmp*a%mod;
        a=a*a%mod;
        b>>=1;
    }
    return tmp;
}

LL F(LL p,int k)
{
    LL Iv=P((p*p%mod-1+mod)%mod,mod-2);
    LL S=(P(p,k<<1)-1+mod)%mod;
    S=S*Iv%mod;
    S=p*(S*p%mod-S+mod)%mod;
    return (S+1)%mod;
}
int main ()
{
    int T;
    LL Iv=P(2,mod-2);
    scanf("%d",&T);
    while(T--)
    {
        int n;
        scanf("%d",&n);
        LL N=(LL)n;
        LL ans=1;
        for(int i=2;i*i<=n;i++)
        {
            if(n%i) continue;
            int k=0;
            while(n%i==0)
            {
                k++;
                n/=i;
            }
            ans*=F((LL)i,k);
            ans%=mod;
        }

        if(n>1)
            ans=(ans*F((LL)n,1))%mod;
        printf("%lld\n",(((LL)1+ans)%mod)*N %mod * Iv % mod);
    }
    return 0;
}