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

1.欧拉函数

代码:

#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 n;

void get_euler(int x,int res){
    for(int i=2;i<=x/i;++i){
    if(x%i==0){
       res=res/i*(i-1);
       while(x%i==0) x/=i;
    }
    }
    if((x>1))  res=res/x*(x-1);
    printf("%d\n",res);
}

int main(){
 ios::sync_with_stdio(false);
    cin>>n;
     while(n--){
    int a;
     cin>>a;
      int t=a;
       
      get_euler(a,t);
     }
  return 0;
}

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;
}