BZOJ2705: [SDOI2012]Longge的问题

最新推荐文章于 2018-09-04 20:02:01 发布

forezxl

最新推荐文章于 2018-09-04 20:02:01 发布

阅读量194

点赞数 1

分类专栏： BZOJ 数论---欧拉定理/函数蒟蒻zxl的Blog专栏文章标签： BZOJ SDOI2012 欧拉函数

本文链接：https://blog.csdn.net/a1799342217/article/details/80298531

版权

蒟蒻zxl的Blog专栏同时被 3 个专栏收录

460 篇文章 6 订阅

订阅专栏

BZOJ

281 篇文章 1 订阅

订阅专栏

数论---欧拉定理/函数

7 篇文章 0 订阅

订阅专栏

这篇博客介绍了BZOJ2705题目的解决方案，主要涉及欧拉函数的概念。通过枚举整数n的质因子，利用欧拉函数φ(n/d)来计算gcd(n,m)=d时m的个数，从而得出答案。算法的时间复杂度为O(n^(1/2))。" 117260105,7852161,Studio3T：MongoDB集群连接及使用指南,"['MongoDB', '数据库管理', '连接工具', '脚本编写']

摘要由CSDN通过智能技术生成

欧拉函数

题目传送门

我们枚举 $n$ 的质因子。答案转化为 $\sum_{d|n}d*f(d)$ ，其中 $f(d)$ 表示 $gcd(n,m)=d$ 的 $m$ 的个数。

$gcd(n,m)=d$ 等价于 $gcd(n/d,m/d)=1$ ，所以 $f(d)=\varphi(n/d)$ ， $O(\sqrt n)$ 枚举 n <script type="math/tex" id="MathJax-Element-10">n</script>的质因子即可。

代码：

#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long LL;
LL n,m,ans;
inline LL phi(LL x){
    LL ret=x;
    for (LL i=2;i<=m;i++)
        if (!(x%i)){
            ret=ret/i*(i-1);
            while (!(x%i)) x/=i;
        }
    if (x!=1) ret=ret/x*(x-1);
    return ret;
}
int main(){
    scanf("%lld",&n),m=sqrt(n);
    for (LL i=1;i<=m;i++)
        if (!(n%i)){
            if (i*i<n) ans+=(n/i)*phi(i);
            ans+=i*phi(n/i);
        }
    return printf("%lld\n",ans),0;
}