欧拉函数

欧拉函数

这里先介绍一下欧拉函数：

       在数论中，对正整数n，欧拉函数表示的是小于等于n的数中与n互质的数的数目。此函数以其首名研究者欧拉命名，它又称为Euler's totient function、φ函数、欧拉商数等。 例如φ(8)=4，因为1,3,5,7均和8互质。

φ函数的值通式：φ(x)=x(1-1/p1)(1-1/p2)(1-1/p3)(1-1/p4)…..(1-1/pn),其中p1, p2……pn为x的所有质因数，x是不为0的整数。φ(1)=1（唯一和1互质的数(小于等于1)就是1本身）。 (注意：每种质因数只一个。比如12=2*2*3那么φ（12）=12*（1-1/2）*(1-1/3)=4若n是质数p的k次幂，φ(n)=p^k-p^(k-1)=(p-1)p^(k-1)，因为除了p的倍数外，其他数都跟n互质，设n为正整数，以 φ(n)表示不超过n且与n互素的正整数的个数，称为n的欧拉函数值，这里函数φ：N→N，n→φ(n)称为欧拉函数。欧拉函数是积性函数——若m,n互质，φ(mn)=φ(m)φ(n)。特殊性质：当n为奇数时，φ(2n)=φ(n), 证明与上述类似。

关系欧拉函数公式有几个模型：

1.p^k的欧拉函数

对于给定的一个素数 p ， φ(p) = p -1。则对于正整数 n = p^k ，

φ(n)= p^k - p^{k -1}

 证明：小于 p^k的正整数个数为 p^k- 1个，其中 和 p^k不互质的正整数有{p * 1,p* 2,...,p * (p^{k - 1}-1)}共计p^{k - 1} - 1 个

             所以 φ(n) = p^k- 1 - (p^{k - 1}- 1) = p^k- p^{k- 1} 

2. p * q的欧拉函数

假设 p, q是两个互质的正整数，则 p * q 的欧拉函数为

φ(p * q) = φ(p) * φ(q) ， gcd(p, q) =1 。

证明： 令 n = p * q，gcd(p,q) = 1，根据剩余定理，有Zn和 Zp ×Zq 之间存在一一映射( a ∈ Zp ,b ∈ Zq ⇔ b * p + a * q ∈ Zn )
       所以 n的完全余数集合的元素个数等于集合 Zp × Zq的元素个数, 而后者的元素个数为 φ(p) * φ(q) 

     所以有φ(p * q) = φ(p)* φ(q) 

3.任意正整数的欧拉函数

任意一个整数 n 都可以表示为其素因子的乘积(算术基本定理)为：

  m
  n= ∏  p_i^ki          (m 为 n 的素因子的个数)
       i=1

根据前面两个结论，很容易得出它的欧拉函数为：

             m                             m
  Φ(n) = ∏ p_i^ki-1(p_i-1) =  n * ∏ (1 - 1/pi)    对于任意 n > 2，2 | Φ(n) ,因为必存在  p_i -1 是偶数。
             i=1                          i=1        

下面讨论对于任意整数我们怎么用程序实现欧拉函数公式：

因为要分解质因数，我们先用素数筛得到一张素数表：

int const MAX = 1e6 + 2;
int prime[MAX];
void get_prime()
{
    memset(prime,1,sizeof(prime));
    prime[1] = 0;
    for(int i = 2; i * i <= MAX; i++)  
        if(prime[i])  
            for(int j = i * i; j <= MAX; j += i)  
                prime[j] = 0;
}

    然后就可利用素数表"分解质因数"，再用上述公式通过筛选得到欧拉函数表：

void euler()
{
    for(int i = 1; i <= MAX; i++)
        euler[i] = i;  //初始化值为本身
    for(int i = 2; i <= MAX; i++)
        if(prime[i])   //"分解质因数"
            for(int j = i; j <= MAX; j += i)
		//这步操作类似素数筛，用来筛选与j互素的数的个数
		//先进行除法是为了防止中间数据溢出   
                euler[j] = euler[j] / i * (i - 1);
}

这里举euler[18]的例子解释下euler[j] = euler[j] / i * (i - 1)这步操作的用途

   一开始euler[18] = 18，质数从2开始枚举，因为2是18的因子
所以第一次我们可以得到euler[18] = euler[18]  / 2 * (2 - 1)  = 9
此时euler[18] = 9，9表示小于18的数当中不是2的倍数的数的个数，即1,3,5,7,9,11,13,15,17。
质数到3，因为3也是18的因子，所以euler[18]  = euler[18] / 3 * (3 - 1) = 6 
此时euler[18] = 6，6表示小于18的数当中不是2且不是3的倍数的数的个数之和，表示的数为1,5,7,11,13,17。
质数到5,因为5不是18的因子，故不会影响euler[18]的值
以此类推最后我们得到euler[18] = 6 (即18与1,5,7,11,13,17这些数互素)