数论基础--欧拉函数，欧拉筛

带佬链接–>浅谈欧拉函数

欧拉函数的意义：

给你一个数n，求出1~n之间有多少个数和该数互质。

欧拉函数通式：

eula(n)=n*(1-1/p1)(1-1/p2)(1-1/p3)*(1-1/p4)……(1-1/pn)

模板：

求单个数的欧拉函数值：

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

进阶：

根据通式用埃氏筛原理优化打表1求出1~n之间每个数的欧拉函数（在处理多个数的欧拉函数时更快）：

int Eula[N];
void init()
{
    Eula[1]=1;
    for(int i=2;i<N;i++)
    {
        if(!Eula[i])
        {
            for(int j=i;j<N;j+=i)
            {
                if(!Eula[j])
                    Eula[j]=j;
                Eula[j]=Eula[j]/i*(i-1);
            }
        }
    }
}

用欧拉筛原理优化打表求出求出1~n之间每一个数的欧拉函数值（比埃氏筛更快，因为欧拉筛比埃氏筛快）：

int prime[N],vis[N],phi[N],cnt=0;
void init()
{
    phi[1]=1;
    for(int i=2;i<N;i++)
    {
        if(!vis[i])
        {
            cnt++;
            prime[cnt]=i;
            phi[i]=i-1;
        }
        for(int j=1;j<=cnt&&i*prime[j]<N;j++)
        {
            vis[i*prime[j]]=1;
            if(i%prime[j]==0)
            {
                phi[prime[j]*i]=phi[i]*prime[j];
                break;
            }
            else
            {
                phi[prime[j]*i]=phi[i]*(prime[j]-1);
            }
        }
    }
}

附带一个欧拉筛模板：

const int N=1e5+7;
int prime[N];
map<int,bool>flag;
void init()
{
    int cnt=0;
    for(int i=2;i<=N;i++)
    {
        if(!flag[i]) prime[++cnt]=i;
        for(int j=1;j<=cnt;j++)
        {
            if(i*prime[j]>N) break;
            flag[i*prime[j]]=true;
            if(i%prime[j]==0) break;
        }
    }
    flag[1]=true;
    flag[0]=true;
}