数论---欧拉函数，筛法求欧拉函数

欧拉函数

欧拉函数：

互质：一个数不能被另一个数整除，表示互质，1，-1能与任何数互质

证明：欧拉函数表示的是1~N中与N互质的数

N能够被拆分为多个质因子，,并且有相应的指数

与分解质因数类似

#include <iostream>
#include <algorithm>
using namespace std;

int phi(int x)
{
    int res=x;
    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)    //最后一个数不为1，表示最大质因子
        res=res/x*(x-1);
    return res;
}

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

筛法求欧拉函数

筛法求欧拉函数直接套用线性筛的模板就可以了

分为以下三种情况考虑

• i是质数                        i的欧拉函数值=i-1（从定义出发）
• pj是i的质因子              pj*i的欧拉函数值=phi[i]*j (从欧拉函数展开式出发)
• pj不是i的质因子          pj*i的欧拉函数值=phi[i]*pj*(1-1/pj)

第二种情况证明：

 

第三种情况证明：

 

#include <iostream>
#include <algorithm>
using namespace std;
const int N=1e6+10;
typedef long long ll;
int phi[N];
int primes[N],cnt;
bool st[N];

ll get_oula(int n)
{
    ll res=0;
    phi[1]=1;
    for(int i=2;i<=n;i++)
    {
        if(!st[i]) 
        {
            primes[cnt++]=i;
            phi[i]=i-1; //质数前面所有数都互质
        }
        for(int j=0;primes[j]<=n/i;j++)
        {
            int t=primes[j]*i;
            st[t]=true;
            if(i%primes[j]==0) 
            {
                phi[t]=primes[j]*phi[i];
                break;
            }
            phi[t]=phi[i]*(primes[j]-1);
        }
    }
    for(int i=1;i<=n;i++)
        res+=phi[i];
    return res;
}

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