蓝桥杯C++——试题 算法提高 欧拉函数

蓝桥杯试题集试题总汇（C++）

问题描述

　　老师出了一道难题，小酱不会做，请你编个程序帮帮他，奖金一瓶酱油：
　　从1—n中有多少个数与n互质?
　　|||||╭══╮ ┌═════┐
　　╭╯让路║═║酱油专用车║
　　╰⊙═⊙╯ └══⊙═⊙═(坑爹的题面格式化，害得我用‘|’来代替空格，复制到记事本上看就变成正版的了)

输入格式

　　输入共一行，表示一个整数n。

输出格式

　　输出共一行，表示从1—n中与n互质的数的个数。

样例输入

30

样例输出

8

数据规模和约定

　　60%的数据≤10^6
　　100%的数据≤2*10^9

分析

欧拉函数的详细说明

欧拉函数phi(n)是数论中非常重要的一个函数，其表示1到n-1之间，与n互质的数的个数。显然的，我们可以通过定义直接计算phi(n)。
当然，phi(n)还有这么一种计算方法。
首先我们对n进行质因数分解，不妨设n=p1^a1 * p2^a2 * ... * pk^ak （这里a^b表示a的b次幂，p1到pk为k个互不相同的质数，a1到ak均为正整数），那么
phi(n)=n(1-(1/p1))(1-(1/p2))....(1-(1/pk))
稍稍化简一下就是
phi(n)=n(p1-1)(p2-1)...(pk-1)/(p1*p2*...*pk)

计算的时候小心中间计算结果超过int类型上界，可通过调整公式各项的计算顺序避免(比如先做除法)!

参考代码

#include <iostream>
using namespace std;

//欧拉函数直接筛选
long long direct_phi(long long x)
{
    long long ans = x;
    for (int i = 2; i * i <= x; i++)
    {
        if (x % i == 0)
        {
            ans = ans * (i-1) / i;	 这个对应的就是欧拉函数的通式
            while(x % i == 0){      // 完全筛去 i 因子，确保下一个得到的 i 是 n 的素因子
                x = x / i;
            }
        }
    }
    if (x > 1) {        // 为了保证防止未除掉的素因子出现
        ans = ans * (x-1) / x;
    }
    return ans;
}

int main(int argc, char** argv)
{
    long long n;
    cin>> n;
    cout<< direct_phi(n)<< endl;
    return 0;
}