欧拉函数及筛法求欧拉函数

2020-04-15 12:00:11 阅读：283 来源： 互联网

1.欧拉函数

代码:

 1 #include <iostream>
 2 #include <cstdio>
 3 #include <cstring>
 4 #include <string.h>
 5 #include <math.h>
 6 #include <algorithm>
 7 #include <stack>
 8 #include <queue>
 9 #include <vector>
10 #include <map>
11 #include <set>
12 #define ll long long
13 const int N=1e6+10;
14 using namespace std;
15 typedef pair<int,int>PII;
16 
17 int n;
18 
19 void get_euler(int x,int res){
20     for(int i=2;i<=x/i;++i){
21     if(x%i==0){
22        res=res/i*(i-1);
23        while(x%i==0) x/=i;
24     }
25     }
26     if((x>1))  res=res/x*(x-1);
27     printf("%d\n",res);
28 }
29 
30 int main(){
31  ios::sync_with_stdio(false);
32     cin>>n;
33      while(n--){
34     int a;
35      cin>>a;
36       int t=a;
37        
38       get_euler(a,t);
39      }
40   return 0;
41 }

2.筛法求1～n的欧拉函数之和

代码:

#include <iostream>
#include <cstdio>
#include <cstring>
#include <string.h>
#include <math.h>
#include <algorithm>
#include <stack>
#include <queue>
#include <vector>
#include <map>
#include <set>
#define ll long long
const int N=1e6+10;
using namespace std;
typedef pair<int,int>PII;

int primes[N],cnt;
int phi[N];
bool st[N];

ll get_eulers(int n){
    phi[1]=1;
    for(int i=2;i<=n;++i){
     if(!st[i]){
      primes[cnt++]=i;
      phi[i]=i-1;
    }
    for(int j=0;j<cnt && primes[j]<=n/i;++j){
       st[primes[j]*i]=true;
       if(i%primes[j]==0){
          phi[primes[j]*i]=phi[i]*primes[j];
          break;
       }
     else  phi[primes[j]*i]=phi[i]*(primes[j]-1);
    }
    }
    ll res=0;
    for(int i=1;i<=n;++i) res+=phi[i];
    return res;
}

int main(){
 ios::sync_with_stdio(false);
 int n;
  cin>>n;
   printf("%lld\n",get_eulers(n));
  return 0;
}

标签：函数,筛法,get,int,res,long,include,欧拉
来源： https://www.cnblogs.com/lr599909928/p/12704284.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

ICode9

欧拉函数及筛法求欧拉函数