欧拉函数

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

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

n=p1a1p2a2...pkak4

那么可得

φ(n)=n(p11)/p1(p21)/p2...(pk1)/pk=nk1pi1k1pi φ ( 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 n 分解为素数幂乘积形式,即n=p1a1p2a2...pkak,设已经处理完前 k1 k − 1 个因子,设前 k1 k − 1 个因子所能形成的数的所有因子的欧拉函数和为 S S .加入素因子pk,则对答案的贡献为 S(pk1)ak1i=0pik S ∗ ( p k − 1 ) ∗ ∑ i = 0 a k − 1 p k i ,将其与原答案合并得: ANS=S((pk1)ak1i=0pik+1) A N S = S ∗ ( ( p k − 1 ) ∗ ∑ i = 0 a k − 1 p k i + 1 ) .化简即可得 ANS=Spakk A N S = S ∗ p k a k .推广一下,即 ANS=pa11pa22...pakk=n A N S = p 1 a 1 ∗ p 2 a 2 ∗ . . . ∗ p k a k = n

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

求欧拉函数
质因数分解法 O(n) 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) 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]沙拉公主的困惑
上帝与集合的正确用法

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值