欧拉函数

欧拉函数

在数论,对正整数n,欧拉函数是小于n的正整数中与n互质的数的数目(φ(1)=1)。

其中p1, p2……pn为x的所有质因数,x是不为0的整数

 

分解n=p1q1 * p2q2 * p3q3 * ……* pkqk

φ(n)= n*(1 - 1/p1) *(1 - 1/p2) *(1 - 1/p3)*……. *(1 - 1/pk)

            = p1q1 * (1 - 1/p1) * p2q2 * (1 - 1/p2) * p3q3 * (1 - 1/p3) *…..* pkqk * (1 - 1/pk)

            = p1q1 * ( (p1 – 1)/p1 ) * p2q2 * ( (p2 – 1)/p2 ) * p3q3 * ( (p3- 1)/p3 ) *…..* pkqk * ( (pk – 1)/pk )

            = p1q1-1 * (p1 – 1) * p2q2-1 * (p2 – 1) * p3q3-1 * (p3- 1) *…..* pkqk-1 * (pk – 1)

 

那么对于一个素数,显然 φ(p)= p-1 

 

 

 

 

代码实现利用到了线性筛唉

g[n]表示n的因子个数

求一个数约数的个数,假设把它分解完了,那么g[i]=(1+q1)(1+q2)….(1+qn)

其实就是乘法原理啊,对于每个因子都可以选指数为0次,1次…q1次,一共(q1+1)种情况

 

phi[ ]是欧拉函数,欧拉函数是积性函数

证明: 假设n,m互质,那么他们一定没有相同的质因子

那么n=p1q1 * p2q2 * p3q3 * ……* pkqk

m=t1l1 * t2l2 * t3l3 * ……* tklk 

假设k=n*m

那么显然φ(k)=φ(n)*φ(m)

/*欧拉筛

筛数 i 
选取一个素数 p
把 p * i 筛掉
此时会检查 i 的最小素因子是否是p
所以有两种可能:
1、i 和 p 互素
    phi[i * p] = phi[i] * phi[p]
2、i 的最小素因子刚好是 p 
    phi[i * p] = phi[i] * p
*/


bool notprime[]
int tot, prime[], e[];
int phi[];// phi 是欧拉函数
int g[]; // g[n]表示n的因子个数
int f[]; // f[n]表示n的最小素因子的次数
         //    f[12] = 2;
         // f[81] = 4;
void sieve(int n) {
    for (int i = 2; i <= n; i++) {
        if (!notprime[i]) {
            prime[tot++] = i;
            e[i] = i;
            phi[i] = i - 1;
            g[i] = 2;
            f[i] = 1;
        }
        for (int j = 0; j < tot && (k = prime[j] * i) < n; j++) {
            e[k] = prime[j];
            notprime[k] = true;
            if (e[i] == prime[j]) {
                phi[k] = phi[i] * prime[j];  //i和prime[j]不互质 
                f[k] = f[i] + 1;
                g[k] = g[i] / (f[i] + 1) * (f[k] + 1);
                break;
            }
            else {
                phi[k] = phi[i] * phi[prime[j]]; //i和prime[j]互质 
                g[k] = g[i] * g[prime[j]];
                f[k] = 1;
            }
        }
    }
}

 

 

欧拉定理

 

posted @ 2019-04-09 10:27  晔子  阅读(2297)  评论(0编辑  收藏  举报