欧拉函数详解

 


欧拉函数

我们用ϕ(n)表示欧拉函数

定义:ϕ(n)表示对于整数n,小于等于n中与n互质的数的个数

性质

1.ϕ(n)为积性函数

证明:

此处证明需要用到下面计算方法1中的内容,建议先看后面再回过头来看这里

假设存在p,q,且pq=n

n,p,q进行质因数分解

n=a1p1a2p2...akpk

p=a1p1a2p2...ampm

q=am+1pm+1am+2m+2...akpk

那么

φ(n)=ni=1k(11pi)

φ(a)=ai=1m(11pi)

φ(b)=bi=m+1k(11pi)

因为n=ab

显然

φ(n)=φ(a)φ(b)

这种方法也是常见的证明一个函数是积性函数的方法

 

2.d|nϕ(d)=n

3.1n中与n互质的数的和为nϕ(n)2(n>1)

证明:若gcd(n,i)=1,那么gcd(n,ni)=1

因此与n互质的数都是成对出现的。且每一对的和都为n

这样最终答案为nϕ(n)2

4. aϕ(n)1(modn)

 

计算方法

(n)计算单值欧拉函数

假设我们需要计算ϕ(n)

分情况讨论

1.当n=1

很明显,答案为1

2.当n为质数时

根据素数的定义,答案为n1

(仅有nn不互质)

3.当n为合数时

我们已经知道了n为素数的情况

不妨对n进行质因数分解

n=a1p1a2p2...akpk

假设k=1

那么ϕ(pk)=pkpk1

证明:

考虑容斥,与一个数互素的数的个数就是这个数减去与它不互素的数的个数

因为p是素数,所以在pk中与其不互素的数为1p,2p....pk1p,有pk1

得证

 

k1

ϕ(n)

=φ(a1p1a2p2akPk)

=i=1kaPiaiPi1

=i=1kaiPi(11pi)

=ni=1k(11pi)

复制代码
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int MAXN = 1e7 + 10;
int p, ans = 1, N;
void GetPhi() {
    for(int i = 2; i * i <= p; i++) {
        if(p % i == 0) {
            int now = i - 1; p /= i;
            while(p % i == 0) now = now * i, p /= i;
            ans = ans * now;
        }
    }
    if(p != 1) ans *= (p - 1); 
}
int main() {
    cin >> p; N = p;
    GetPhi();
    cout << ans;
    return 0;
}
复制代码

 

线性筛

因为欧拉函数是积性函数

因此可以使用线性筛法

性质1

p为素数,则φ(p)=p1

证明:

1p中,只有(p,p)1

性质2

i mod p0,且p为素数

φ(ip)=φ(i)φ(p)

=φ(ip)=φ(i)(p1)

这一步同时利用了性质1和欧拉函数的积性

性质3

i mod p=0,且p为素数,

φ(ip)=φ(i)p

证明:

没怎么看懂,丢一个链接

http://blog.csdn.net/Lytning/article/details/24432651

复制代码
#include<bits/stdc++.h>
using namespace std;
const int MAXN = 3e5 + 10;
void GetPhi(int N) {
    static int phi[MAXN], vis[MAXN], prime[MAXN], tot = 0;
    for(int i = 2; i <= N; i++) {
        if(!vis[i]) prime[++tot] = i, phi[i] = i - 1;
        for(int j = 1; j <= tot && i * prime[j] <= N; j++) {
            vis[i * prime[j]] = 1;
            if(!(i % prime[j])) {phi[i * prime[j]] = phi[i] * prime[j]; break;}
            else phi[i * prime[j]] = phi[i] * (prime[j] - 1);
        }
    }
    while(cin >> N) cout << phi[N] << endl;
}
int main() {
    GetPhi(100);
    return 0;
}
复制代码

例题

放几道水题

http://poj.org/problem?id=2407

题解

http://poj.org/problem?id=2478

题解

https://www.luogu.org/problemnew/show/P2158

题解

参考资料

数论学习笔记 欧拉函数 (一些性质和运用)内置杜教筛

 

posted @   自为风月马前卒  阅读(2307)  评论(0编辑  收藏  举报
编辑推荐:
· go语言实现终端里的倒计时
· 如何编写易于单元测试的代码
· 10年+ .NET Coder 心语,封装的思维:从隐藏、稳定开始理解其本质意义
· .NET Core 中如何实现缓存的预热?
· 从 HTTP 原因短语缺失研究 HTTP/2 和 HTTP/3 的设计差异
阅读排行:
· 分享一个免费、快速、无限量使用的满血 DeepSeek R1 模型,支持深度思考和联网搜索!
· 基于 Docker 搭建 FRP 内网穿透开源项目(很简单哒)
· ollama系列1:轻松3步本地部署deepseek,普通电脑可用
· 按钮权限的设计及实现
· 25岁的心里话

Contact with me

点击右上角即可分享
微信分享提示