数论学习笔记之欧拉函数

　　　　　　　　　　　　数论学习笔记之欧拉函数

~1~.　　欧拉函数φ（n) 表示小于或等于 n 的正整数中与 n 互质的数的个数，即：φ(n) 为 集合{ a| 1 <= a <= n 且 gcd(a,n) = 1 }的元素个数，比如，φ（5） = 4，（1,2,3,4）。

~2~.　　欧拉函数的的推导（待续）

　　先看几个结论：

　　①，对于素数 p ，显然 $ 。(此式记为 Ⅰ）

　　②，对于素数 p，有 φ(p^k) = p^k -  p^k-1。　（此式记为 Ⅱ）

　　③，如果 gcd(m，n） = 1，则 φ(m*n) = φ(m)*φ(n) 。（此式记为 Ⅲ ）

　　④，算术基本定理：任意 m >= 2 都可以化成这种形式（且是唯一的）：

 m = p₁k₁*p₂k₂*p₃k₃*......*p_nk_n。（其中 p_i  都是素数，且互不相同，把此式记为 Ⅳ) 

　　⑤，要证明的式子： φ(m) = m*(1 - p₁) * (1 - p₂) * ......*(1 - p_n) 。(其中 p_i  都是素数，且互不相同,此式记为 Ⅴ）

　　

要证明 ⑤ 首先得证明前面的几个结论。下面依次证明：　　

①④不予证明。

②的证明（φ（p^k) = p^k - p^k-1)　　

　　这个可以根据定义来证明　　

自然数列：  1，2， 3， 4， 5， ，，p，，，2*p，，，3*p，，，，，，pⁿ    .

　　在 [1 , p] 中 满足1 <= a <= n 且 gcd(a,p) = 1 的数有 p-1 个，（除了 p之外都满足）

　　在 (p，2*p]中 满足1 <= a <= n 且 gcd(a,p) = 1 的数有 p-1 个，(除了 2*p)

　　。。。。。。

　　可以看出每个区间的长度都是 p ，一共有 p^k/p = p^k-1个区间 ，所以 φ(p^k) = p^k-1*(p - 1)。证明完成。、

③的证明 （如果 gcd(m，n） = 1，则 φ(m*n) = φ(m)*φ（n)）

这个是用中国剩余定理 证明的～待续。。。

 

~3~. 几个代码

 

　　1.单个数的欧拉函数值。

　　

LL eu(LL N)
{
    LL ans = N ;
    for (LL i = 2; i*i <= N; i ++) {
        if (N%i == 0) {
            ans -= ans/i ;
            while (N%i == 0) N /= i ;
        }
    }
    if (N > 1) ans -= ans/N ;
    return ans ;
}

 

　　2.筛法求欧拉函数

LL euler[maxn+1] ;
void get_euler()
{
    memset(euler,0,sizeof(euler)) ;
    euler[1] = 1 ;
    for (LL i = 2; i <= maxn; i ++) {
        if (!euler[i]) {
            for (LL j = i; j <= maxn; j += i) {
                if (!euler[j]) euler[j] = j ;
                euler[j] =euler[j]/i*(i-1) ;
            }
        }
    }
}

　　3.同时得到素数表和欧拉函数

bool ch_ch[maxn+1] ;
LL pr[maxn+1] ;
LL phi[maxn+1] ;
LL tot = 0 ;
void get_euler_prime()
{
    memset(ch_ch,false,sizeof(ch_ch)) ;
    phi[1] = 1;
    for (LL i = 2; i <= maxn; i ++) {
        if (!ch_ch[i]) {
            pr[tot++] = i ;
            phi[i] = i - 1 ;
        }
        for (LL j = 0; j < tot && pr[j] <= maxn/i; j ++) {
            ch_ch[i*pr[j]] = true ;
            if (i%pr[j] == 0) {
                phi[i*pr[j]] = phi[i]*pr[j] ;
                break ;
            }
            else {
                phi[i*pr[j]] = phi[i] * (pr[j] - 1) ;
            }
        }
    }
}

 　　感觉代码的第十七行十分巧妙，是这样的:

 i = t*pr[j]^k，i * pr[j] = t*pr[j]^(k+1)， phi(i*pr[j]) = phi(t*pr[j]^(k+1)) = phi(t) * phi(pr[j]^(k+1))

 = phi(t) * pr[j] * phi(pr[j]^k) = phi(t) * phi(pr[j]^k) * pr[j] = phi(i) * pr[j] ；

　　4.一个小结论。　　

　　　　设 b[i] 与 N 互质, 则 b[1] + b[2] + ... + b[φ(N)] = N*φ(N)/2 ；

倒序相加即可证明。若 gcd(a,N) = 1,那么 gcd (N-a,N) = 1 ;