码迷,mamicode.com
首页 > 其他好文 > 详细

欧拉函数

时间:2018-09-07 20:07:23      阅读:189      评论:0      收藏:0      [点我收藏+]

标签:c11   ddl   src   return   span   int   防止   i++   for   

phi[i]=从1到i与i互素的数的个数

公式:技术分享图片  ,pi为x的质因子,n为x的质因子个数

例如:

12=2*2*3;

12=12*(1-1/2)*(1-1/3)

 求单个的:

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

求很多个的:

特性:

1.若a为质数,phi[a]=a-1;

2.若a为质数,b与a不互素,即b%a==0,gcd(a,b)=a,那么phi[b*a]=phi[b]*a;

3.若a为质数,b与a互素,即b%a!=0,gcd(a,b)=1,phi[b*a]=phi[b]*phi[a]=phi[b]*(a-1);

由以上特性,我们需要先把质数找出来并存起来

利用找出来的质数和i把它们的乘积直接求出,类似于埃筛

break是为了防止一个数被求多次,我们总是让一个数被它最小的质因子筛去

const int maxn=1e5+10;
int phi[maxn],prime[30000];
int cnt=0;
void Euler()
{
    phi[1]=1;
    for(int i=2;i<maxn;i++)
    {
        if(!phi[i])
        {
            phi[i]=i-1;
            prime[cnt++]=i;
        }
        for(int j=0;j<cnt&&i*prime[j]<maxn;j++)
        {
            if(i%prime[j])phi[i*prime[j]]=phi[i]*(prime[j]-1);
            else
            {
                phi[i*prime[j]]=phi[i]*prime[j];
                break;
            }
        }
    }
}

 

欧拉函数

标签:c11   ddl   src   return   span   int   防止   i++   for   

原文地址:https://www.cnblogs.com/raincle/p/9606667.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!