集训第8天总结

• 今天学了一下欧拉函数，如果单独求一个数的欧拉函数我们可以用试触法而因为x的因子中能大于等于根号x的因子最多只能有一个所以我们可以用o(根号n)的复杂度取求出一个数的欧拉函数。
  要求n个数的欧拉函数的话我们可以用欧拉筛o(n)的实现，如果i为质数的话易知f(i)=i-1，如果i%p(j)==0的话那么由唯一分解定理可知i的质因数与ip(j)相同就推出f(ip(j))=p(j)f(i)，如果i%p(j)!=0时因为欧拉筛的每个数都是被他是最小因子所筛的从而可知ip(j)的质因数个数仅比i多一个p(j)，所以可以推出f(i*p(j))=f(i)p(j)(1-1/p(j))。（这里的/号是数学上的除号）

代码如下
求一个数的欧拉函数

#include<iostream>//(与n互质的数的个数)f(n)=n*(1-1/p1)(1-1/p2)...... 
using namespace std;
int main(){
    int x;
    scanf("%d",&x);
    int sum=x;
    for(int i=2;i<=x/i;i++){
        if(x%i == 0){
            sum=sum/i*(i-1);
            while(x%i==0){
                x/=i;
            }
        }
    }
    if(x>1){
        sum=sum/x*(x-1);
    }
    printf("%d\n",sum);
    return 0;
}

求前n个数的欧拉函数

#include<iostream>
using namespace std;
const int N=1000001;
int prime[N],vis[N],cnt=0,ola[N];
void Prime(int n){
    ola[1]=1;
    for(int i=2;i<=n;i++){
        if(!vis[i]){
            ola[i]=i-1;
            prime[++cnt]=i;
        }
        for(int j=1;i*prime[j]<=n;j++){
            vis[i*prime[j]]=1;
            if(i%prime[j]==0){
                ola[i*prime[j]]=ola[i]*prime[j];//如果prime[j]是i的因子 
                break;
            }
            ola[i*prime[j]]=ola[i]*(prime[j]-1);//如果prime[j]不是i的因子 
        }
    }
    long long ans=0;
    for(int i=1;i<=n;i++){
        ans+=ola[i];
    }
    printf("%lld\n",ans);
}
int main(){
    int n;
    scanf("%d",&n);
    Prime(n);
    return 0;
}