ACM第一课---数论之欧拉函数

原创 2013年11月05日 22:13:38

欧拉函数

定义

欧拉函数PHI(n)表示的是比n小,并且与n互质的正整数的个数(包括1)。
比如:PHI(1) = 1; PHI(2) = 1; PHI(3) = 2; PHI(4) = 2; ... PHI(9) = 6; ...

通式及其证明

要计算一个正整数n的欧拉函数的方法如下:
1. 将n表示成素数的乘积: n = p1 ^ k1 * p2 ^ k2 * ... * pn ^ kn(这里p1, p2, ..., pn是素数)

2. PHI(n) = (p1 ^ k1 - p1 ^ (k1 - 1)) * (p2 ^ k2 - p2 ^ (k2 - 1)) * ... * (pn ^ kn - pn ^ (kn - 1))
               =n*(p1-1)(p2-1)……(pi-1)/(p1*p2*……pi);
                  =n*(1-1/p1)*(1-1/p2)....(1-1/pn)



证明过程如下:
1. 容易想到:当n为素数时,PHI(n) = n - 1。因为每个比n小的正整数都和n互素。当n为素数p的k次方时,PHI(n) = p ^ k - p ^ (k - 1)。因为在1到n之间的正整数只有p的倍数和n不互素,这样的数有(p ^ k / p)个。
(第二步的证明没看太明白,还希望懂的人讲解一下)。。。。

2. 如果m和n互素,即GCD(m, n) = 1,那么PHI(m * n) = PHI(m) * PHI(n)。用中国剩余定理可以证明,证明的思路是建立这样一种一一对应的关系(a, b) <-> x,其中正整数a小于m并且gcd(a, m) = 1,正整数b小于n并且gcd(b, n) = 1,正整数x小于m*n并且gcd(m*n, x) = 1。证明过程如下:
    1)根据中国剩余定理,如果m和n互素,那么关于未知量x的方程组x % m = a, x % n = b(0 <= a < m, 0 <= b < n),当0 <= x < m * n时存在并且仅存在一个解。容易证明,如果两个这样的方程组有相同的m, n但是a, b不同,那么他们的解x一定不同。
    2)首先用反正法证明:gcd(m, a) = 1且gcd(n, b) = 1是gcd(m*n, x) = 1的必要条件:假设gcd(a, m) = k > 1,由此可得:a = a' * k; m = m' * k => x = k' * m + a = k' * k * m' + k * a' = k * (k' * m' + a'); 所以gcd(x, m) = k > 1。同理可证,如果gcd(b, n) > 1, 那么gcd(x, n) > 1。所以x和m * n互素的必要条件是a和m互诉且b和n互素。
    3)接下来我们证明充分性:由x % m = a 可以得到x = k * m + a;由欧几里德算法求最大公约数的过程(就不证明了,呵呵,还得想)可以知道gcd(x, m) = gcd(m, a) = 1;同理可得,如果gcd(n, b) = 1那么gcd(x, n) = 1。接下来很容易得到:gcd(m*n, x) = 1。从而证明了充分性。
    4)上面三步的结论表明,数对(a, b)是可以和x建立起一一对应的关系的,所以有多少个不同的(a, b),就有多少个不同的x。
3.将n分解成素数乘积后,显然对于任意的i, j(i != j)都满足 pi ^ ki和pj ^ kj是互素的,于是可以的到上面的公式。

跟据上面的公式,可以得到关于欧拉函数的递推关系:
假设数p能整除n,那么
如果p还能整除n / p, 即是若(N%p == 0 && (N/p) % p ==0 )PHI(n) = PHI(n / p) * p;
如果p不能整除n / p, 即是若 (N%p == 0 && (N/p) % p != 0)  PHI(n) = PHI(n / p) * (p - 1);

算法实现

先回顾筛法
(主要可参考上篇)
for(i=2;i<=1000000;i++)
    {
        if(!c[i])prime[len++]=i;
        for(j=0;j<len&&prime[j]*i<=1000000;j++)
        {
            c[prime[j]*i]=1;//不是质数
            if(i%prime[j]==0)break; //                                 
            }                  
        }    
    }

然后便是求解欧拉函数
phi[1] = 1;
    for (i = 2; i < 10000; i++) {
        if (!mark[i]) {
            phi[i] = i - 1;
            continue;
        }
        for (j = 0; j < size && prime[j] * prime[j] <= i; j++) {
            if (i % prime[j] == 0) {
                if (i / prime[j] % prime[j] == 0)
                    phi[i] = prime[j] * phi[i / prime[j]];
                else
                    phi[i] = (prime[j] - 1) * phi[i / prime[j]];
                break;
            }
        }
    }


总结

表示真心不好理解.所以,暂时只能先记忆,做做题目。
可能就会理解吧

【未完待续】........

欧拉函数求法与应用

欧拉函数在OI中是个非常重要的东西,不知道的话会吃大亏的. 欧拉函数用希腊字母φ表示,φ(N)表示N的欧拉函数. 对φ(N)的值,我们可以通俗地理解为小于N且与N互质的数的个数(包含1)....
  • sentimental_dog
  • sentimental_dog
  • 2016-07-23 10:34:04
  • 12814

欧拉函数及其证明

请思考以下问题:   任意给定正整数n,请问在小于等于n的正整数之中,有多少个与n构成互质关系?(比如,在1到8之中,有多少个数与8构成互质关系?) 计算这个值的方法就叫做欧拉函数,以φ(n)表示...
  • paxhujing
  • paxhujing
  • 2016-05-09 15:51:46
  • 9487

欧拉函数详解

欧拉函数详解 对一个正整数N,欧拉函数是小于N且与N互质的数的个数.。 例如φ(24)=8,因为1, 5, 7, 11, 13, 17, 19, 23均和 24 互质。 φ(n) = n*(1-1/...
  • w144215160044
  • w144215160044
  • 2016-04-15 09:28:07
  • 1656

欧拉函数

单个数的欧拉函数以及线性筛
  • HandsomeHow
  • HandsomeHow
  • 2016-05-05 21:41:38
  • 1126

算法-欧拉函数

欧拉函数简介——摘自百度百科 注意以下三个特殊性质 编程实现 利用欧拉函数和它本身不同质因数的关系,用筛法计算出某个范围内所有数的欧拉函数值。...
  • qq_36409190
  • qq_36409190
  • 2016-11-15 17:22:24
  • 1034

欧拉函数的一个性质及其证明

欧拉函数φ(n)指小于或等于n的正整数中与n互质的数的数目
  • chen1352
  • chen1352
  • 2016-02-19 10:02:59
  • 1856

欧拉函数的计算

直接上代码
  • WT_cnyali
  • WT_cnyali
  • 2016-11-17 08:21:37
  • 1121

欧拉函数及其求法

欧拉函数:对正整数n,欧拉函数是少于或等于n的数中与n互质的数的数目。例如euler(8)=4,因为1,3,5,7均和8互质。 Euler函数表达通式:euler(x)=x(1-1/p1)...
  • qq_27138357
  • qq_27138357
  • 2015-08-10 15:47:57
  • 2285

欧拉函数及其部分性质

欧拉函数是指:对于一个正整数n,小于n且和n互质的正整数(包括1)的个数,记作φ(n) 。 通式:φ(x)=x*(1-1/p1)*(1-1/p2)*(1-1/p3)*(1-1/p4)…..(1-...
  • leolin_
  • leolin_
  • 2011-07-28 20:41:04
  • 18437

欧拉线性筛法求素数(顺便实现欧拉函数的求值)

我们先来看一下最经典的埃拉特斯特尼筛法。时间复杂度为O(n loglog n) int ans[MAXN]; void Prime(int n) { int cnt=0; memset(prime...
  • NK_test
  • NK_test
  • 2015-05-29 23:02:16
  • 8503
收藏助手
不良信息举报
您举报文章:ACM第一课---数论之欧拉函数
举报原因:
原因补充:

(最多只允许输入30个字)