[总结] 欧拉函数

欧拉函数 φ ( n ) 表示小于等于 n 的数中与 n 互质的数的个数
引理:
1. 如果 n 为素数 p ,则 φ ( p ) = p 1
2. 如果 n 为某个素数 p 的幂次 p a ,则 φ ( p a ) = ( p 1 ) p a 1
3. 如果 n 为任意两个互质的数 a , b 的乘积,则 φ ( a b ) = φ ( a ) φ ( b )

继续推导,可得:
4. 将 n 分解为素数幂乘积形式,即

n = p 1 a 1 p 2 a 2 . . . p k a k 4

那么可得
φ ( n ) = n ( p 1 1 ) / p 1 ( p 2 1 ) / p 2 . . . ( p k 1 ) / p k = n 1 k p i 1 1 k p i

我还能推:
5. n 的所有因子的欧拉函数和为 n

推导过程: n 分解为素数幂乘积形式,即 n = p 1 a 1 p 2 a 2 . . . p k a k ,设已经处理完前 k 1 个因子,设前 k 1 个因子所能形成的数的所有因子的欧拉函数和为 S .加入素因子 p k ,则对答案的贡献为 S ( p k 1 ) i = 0 a k 1 p k i ,将其与原答案合并得: A N S = S ( ( p k 1 ) i = 0 a k 1 p k i + 1 ) .化简即可得 A N S = S p k a k .推广一下,即 A N S = p 1 a 1 p 2 a 2 . . . p k a k = n

欧拉定理: a p 互质,则 a φ ( p ) 1 ( mod p )
其实费马小定理就是欧拉定理的一种特殊情况( p 为质数,则 φ ( p ) = p 1

求欧拉函数
质因数分解法 O ( n ) 推出单个欧拉函数
利用引理1,4直接推即可

代码实现如下:

int get_phi_single(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;
}

线性算法 O ( n ) 推出1~n的欧拉函数
利用引理1,4直接推即可,同时将质数筛了出来

代码实现如下:

int get_phi() {
    for(int i=2;i<size;i++) {
        if(!check[i]) {
            prime[++top]=i;
            phi[i]=i-1;
        }
        for(int j=1;j<=top && i*prime[j]<size;j++) {
            check[i*prime[j]]=1;
            if(i%prime[j]==0) {
                phi[i*prime[j]]=prime[j]*phi[i];
                break;
            }
            else phi[i*prime[j]]=(prime[j]-1)*phi[i];
        }
    }
}

题目:
[SDOI2012]Longge的问题
奇数国
Visible Lattice Points
[SDOI2008]仪仗队
[Sdoi2008]沙拉公主的困惑
上帝与集合的正确用法

猜你喜欢

转载自blog.csdn.net/force_chl/article/details/80543103