【欧拉函数】BZOJ2705(SDOI2012)[Longge的问题]题解

最新推荐文章于 2022-01-15 14:43:59 发布

ZigZagK

最新推荐文章于 2022-01-15 14:43:59 发布

阅读量321

点赞数

分类专栏： BZOJ题解欧拉函数及扩展定理

本文链接：https://blog.csdn.net/zzkksunboy/article/details/78764611

版权

BZOJ题解同时被 2 个专栏收录

117 篇文章 0 订阅

订阅专栏

欧拉函数及扩展定理

5 篇文章 0 订阅

订阅专栏

题目概述

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

解题报告

\sum i = 1 n (i, n) = \sum d | n d \sum i = 1 n d [(i, n d) = 1] = \sum d | n d φ (n d)

$\sum_{i=1}^{n}(i,n)=\sum_{d|n}d\sum_{i=1}^{{n\over d}}[(i,{n\over d})=1]=\sum_{d|n}d\varphi({n\over d})$

然后暴枚因子求欧拉函数就可以了，效率不知道怎么证明，反正可以秒过:P。

示例程序

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

LL n,ans;

inline LL phi(LL n)
{
    LL ans=n;
    for (int i=2,S=sqrt(n);n>1&&i<=S;i++)
        if (!(n%i)) {ans=ans/i*(i-1);while (!(n%i)) n/=i;}
    if (n>1) ans=ans/n*(n-1);return ans;
}
int main()
{
    freopen("program.in","r",stdin);
    freopen("program.out","w",stdout);
    scanf("%lld",&n);ans=0;
    for (int i=1,S=sqrt(n);i<=S;i++)
        if (!(n%i)) {ans+=phi(n/i)*i;if (n/i>S) ans+=phi(i)*n/i;}
    return printf("%lld\n",ans),0;
}