poj 2480 Longge's problem（积性函数 & 欧拉函数）

最新推荐文章于 2020-08-26 09:47:16 发布

theArcticOcean

最新推荐文章于 2020-08-26 09:47:16 发布

阅读量960

点赞数

分类专栏： algorithm_数论文章标签： Euler

本文链接：https://blog.csdn.net/thearcticocean/article/details/50539548

版权

algorithm_数论专栏收录该内容

81 篇文章 0 订阅

订阅专栏

http://poj.org/problem?id=2480

大意：求解 ∑gcd(i, N) 1<=i <=N。

对于最大公约数，它有这样的性质，gcd(n,m1*m2)=gcd(nm1)*gcd(n,m2) 比如gcd(6,12)=gcd(6,3)*gcd(6,4).

我们设f(N)=∑gcd(i, N)，那么它也应该是一个积性函数。（积性函数的积是积性函数，积性函数的和是积性函数）

设 $\\ F(n)=\sum\limits_{p=2}^{p\leq n}\phi(\frac{n}{p})p$ 其中p是和n互素的数字。

又设 $\\ n=p_1^{a_1}p_2^{a_2}\cdots p_k^{a_k}$ ，由积性函数的性质：

$\\ F(n)=F(p_1^{a_1})F(p_2^{a_2})\cdots F(p_k^{a_k})$ （*）

同时我们知道：

$\\ F(p_i^{a_i})=\phi(p_i^{a_i})+p_i\phi(p_i^{a_i-1})\cdots p_i^{a_i-1}\phi(p_i)+p_i^{a_i}\phi(1)$

欧拉函数有这样的计算方式： $\\ \phi(p_i^{a_i})=p_i^{a_i-1}(p_i-1)$

所以得到：

$\\ F(p_i^{a_i})=p_i^{a_i-1}(p_i-1)+p_ip_i^{a_i-2}(p_i-1)+\cdots+p_i^{a_i}\\ = (p_i^{a_i}-p_i^{a_i-1})a_i+p_i^{a_i}\\ = p_i^{a_i}(a_i+1-\frac{a_i}{p_i})$

带入（*）中：

$\\ F(n)=p_1^{a_1}p_2^{a_2}\cdots p_k^{a_k}(1+a_1-\frac{a_1}{p_1})(1+a_2-\frac{a_2}{p_2})\cdots (1+a_k-\frac{a_k}{p_k})\\=n\frac{1}{p_1}(p_1+a_1p_1-a_1)\frac{1}{p_2}(p_2+a_2p_2-a_2)\cdots \frac{1}{p_k}(p_k+a_kp_k-a_k)$

写代码：

#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
typedef long long LL;

int main()
{
    LL n,ans;
    while(cin>>n){
        int len=sqrt(n*1.0);
        ans=n;
        LL a=0,p=1;
        for(int i=2;i<=len;i++){
            a=0;
            if(n%i==0){
               a++;
               n/=i;
               p=i;
               while(n%i==0) {
                   a++;
                   n/=i;
               }
               ans=ans/p*(p+a*p-a);
            }
        }
        if(n>1){
            ans=ans/n*(n+n-1);
        }
        printf("%lld\n",ans);
    }
    return 0;
}