欧拉函数模板

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;
}

 

posted @ 2013-07-14 11:21  calmound  阅读(267)  评论(0编辑  收藏  举报