求解欧拉函数

83 篇文章 10 订阅
40 篇文章 2 订阅

这里写图片描述
这里写图片描述
全文围绕以上两个公式展开

---------------------------------------------------------------------------------------------------------------------------

//直接对n进行素因子分解
#include <iostream>
#include <stdio.h>
#include <algorithm>

using namespace std;

int phi(int n){
    int ans=n;
    for(int i=2;i<=n;i++){
        if(n%i==0)
            ans=ans/i*(i-1);    //相当于 ans=ans-ans/i;
        while(n%i==0) n/=i;
    }
    return ans;
}

int main()
{
    int n;
    while(cin>>n){
        cout<<phi(n)<<endl;
    }
    return 0;
}

//对上面进行优化
#include <iostream>
#include <stdio.h>
#include <algorithm>

using namespace std;

int phi(int n){
    int ans=n;
    for(int i=2;i*i<=n;i++){
        if(n%i==0)
            ans=ans/i*(i-1);    //相当于 ans=ans-ans/i;
        while(n%i==0) n/=i;
    }
    if(n>1) ans=ans/n*(n-1);
    return ans;
}

int main()
{
    int n;
    while(cin>>n){
        cout<<phi(n)<<endl;
    }
    return 0;
}

//素数表实现
#include <iostream>
#include <stdio.h>
#include <algorithm>
#include<cstring>

using namespace std;

bool book[50000];
int p[2000];

//线性素数筛法
void prime(){       
    memset(book,0,sizeof(book));
    book[0]=book[1]=1;
    int k=0;
    for(int i=2;i<50000;i++){
        if(!book[i]) p[k++]=i;
        for(int j=0;j<k&&i*p[j]<50000;j++){
            book[i*p[j]]=1;
            if(i%p[j]==0) break;
        }
    }
}

int phi(int n){
    int ans=n;
    for(int i=0;p[i]*p[i]<=n;i++){
        if(n%p[i]==0)
            ans=ans/p[i]*(p[i]-1);   
        while(n%p[i]==0) n/=p[i];
    }
    if(n>1) ans=ans/n*(n-1);
    return ans;
}

int main()
{
    prime();
    int n;
    while(cin>>n){
        cout<<phi(n)<<endl;
    }
    return 0;
}

// 筛法求函数值
#include <iostream>
#include <stdio.h>
#include <algorithm>

using namespace std;

int phi[50010];

void phi_table(int n){
    for(int i=2;i<=n;i++){
        phi[i]=0;
    }
    phi[1]=1;
    for(int i=2;i<=n;i++){
        if(!phi[i]){   //如果i目前不存在,则对i以及i的倍数进行构造。
            for(int j=i;j<=n;j+=i){
                if(!phi[j]) phi[j]=j;
                phi[j]=phi[j]/i*(i-1);  //前面i不存在说明筛选下来i为素数,故对j进行素因子处理操作。
            }
        }
    }
}

int main()
{
    phi_table(50000);
    int n;
    while(cin>>n){
        cout<<phi[n]<<endl;
    }
    return 0;
}

线性筛法

#include<bits/stdc++.h>

using namespace std;

const int maxn=1e7+10;

bool vis[maxn];

int fai[maxn];

int prime[maxn];

void get_oula(int n){
    memset(vis,0,sizeof(vis));
    fai[1]=1;
    int tot=0;
    for(int i=2;i<n;i++){
        if(!vis[i]){
            prime[tot++]=i;
            fai[i]=i-1;
        }
        for(int j=0;j<tot;j++){
            if(i*prime[j]>n) break;
            vis[i*prime[j]]=true;
            if(i%prime[j]==0){
                fai[i*prime[j]]=fai[i]*prime[j];
                break;
            }else{
                fai[i*prime[j]]=fai[i]*(prime[j]-1);
            }
        }
    }
}

int main(){
    get_oula(5000000);
    for(int i=1;i<100;i++){
        cout<<fai[i]<<endl;
    }
    return 0;
}

  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值