线性求欧拉函数

版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/ffgcc/article/details/81951954

仅记录模板

const int N=2e7;
const int maxn=N+10;
int primes,prime[maxn],phi[maxn];
bool vis[maxn];
void init()
{
    int i,j;
    phi[1]=1;
    for (i=2;i<=N;i++)
    {
        if (!vis[i])
        {
            prime[++primes]=i;
            phi[i]=i-1;
        }
        for (j=1;j<=primes&&i*prime[j]<=N;j++)
        {
            vis[i*prime[j]]=1;
            if (i%prime[j])
                phi[i*prime[j]]=phi[i]*(prime[j]-1);
            else
            {
                phi[i*prime[j]]=phi[i]*prime[j];
                break;
            }
        }
    }
}

猜你喜欢

转载自blog.csdn.net/ffgcc/article/details/81951954
今日推荐