2019-03-22 [蓝桥杯][算法提高VIP]欧拉函数

题目描述
给定一个大于1,不超过2000000的正整数n,输出欧拉函数,phi(n)的值。
如果你并不了解欧拉函数,那么请参阅提示。

提示
欧拉函数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)/(p1p2...*pk)

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

输入
在给定的输入文件中进行读入:
一行一个正整数n。
输出
将输出信息输出到指定的文件中:
一行一个整数表示phi(n)。
样例输入
17
样例输出
16
提示
C语言在线学习平台微信号dotcpp
来源
算法提高

使用素数筛

#include<iostream>
#include<algorithm>
#include<cstdio>
using namespace std;
const int N=2000010;
int prime[N],cnt=0;
int res[11],num[11],n,m,tmp;//m记录质数的个数 
bool p[N];
void FindPrime(int n)
{
    for(int i=2;i<=n;i++)
    {
        if(p[i]==false)
        {
            prime[++cnt]=i;
            for(int j=i+i;j<=n;j+=i) p[j]=true;
        }
    }
}
int main(void)
{
    cin>>n;
    tmp=n;
    FindPrime(n);
    //cout<<"cnt="<<cnt<<endl;
    //for(int i=1;i<=cnt;i++) cout<<prime[i]<<" ";
    cnt=1;
    while(n!=1)
    {
        if(n%prime[cnt]==0) 
        {
            res[++m]=prime[cnt];
            while(n%prime[cnt]==0)
            {
                num[m]++;
                n/=prime[cnt];  
            } 
        }
        cnt++;
    }
    for(int i=1;i<=m;i++)
    {
        tmp/=res[i];
        tmp*=res[i]-1;
    }
    cout<<tmp<<endl;
    return 0;
}
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容

  • 算法提高 欧拉函数 时间限制:1.0s 内存限制:512.0MB 提交此题 说明 2016.4.5 已更新试题...
    逍遥_9353阅读 221评论 0 1
  • 在C语言中,五种基本数据类型存储空间长度的排列顺序是: A)char B)char=int<=float C)ch...
    夏天再来阅读 3,417评论 0 2
  • 概述 RSA是"非对称加密算法",非对称加密算法需要两个密钥:公开密钥(publickey)和私有密钥(priva...
    倦飞知还阅读 742评论 0 1
  • 想了想自己好像做事都没什么计划,想到什么做什么。想想这样子还是不行的,所以还是补上今年的一个计划,自己监督一下自己...
    JenniferYe阅读 234评论 0 2
  • 为山方进篑, 遇水必观澜。 风起青蘋末, 翅展苍穹间。
    w5cp2518阅读 200评论 0 0