欧拉函数模板

void eular()

{

    memset(vis,0,sizeof(vis));

    vis[0]=vis[1]=1;

    for(i=2;i*i<=N;i++)

    {

        if(vis[i]==0)

        {

            for(j=i*i;j<=N;j+=i)

                vis[j]=1;

        }

    } //这段求出了N内的所有素数 

    for(i=1;i<=N;i++)

        phi[i]=i;

    for(i=2;i<=N;i++)

    {

        if(vis[i]==0)

        {

            for(j=i;j<=N;j+=i)//这里从i开始,必定能整除i,其倍数也同理

                phi[j]=phi[j]/i*(i-1); //此处注意先/i再*(i-1),否则范围较大时会溢出

        }

    }

}

递归求欧拉函数

for (i = 1; i <= maxn; i++) phi[i] = i;

for (i = 2; i <= maxn; i += 2) phi[i] /= 2;

for (i = 3; i <= maxn; i += 2) if(phi[i] == i) {

for (j = i; j <= maxn; j += i)

phi[j] = phi[j] / i * (i - 1);

单独求欧拉函数

unsigned euler(unsigned x)

{

    // 就是公式

    unsigned i, res=x;

    for (i = 2; i < (int)sqrt(x * 1.0) + 1; i++)

        if(x%i==0)

        {

            res = res / i * (i - 1);

            while (x % i == 0) x /= i; // 保证i一定是素数

        }

    if (x > 1) res = res / x * (x - 1);

    return res;

}

 

你可能感兴趣的:(函数)