定义:1-N中与N互质的数的个数被称为欧拉函数,即为phi(n)
求欧拉函数主要由两类问题:
- 求单个数的欧拉函数O(sqrt(n))
- 求1-N的欧拉函数O(n)
对于第一类问题,我们可以使用分解质因数的方法来求,这个可以用容斥原理很容易证明。
int phi(int n) {
int ans = n;
for (int i = 2; i <= n / i; i++) {
if (n % i == 0) ans = ans / i * (i - 1);
while (n % i == 0) n /= i;
}
if (n > 1) ans = ans / n * (n - 1);
return ans;
}
对于第二类问题,我们可以通过对线性筛略加改进即可求。
void euler_linear_sieve {
for (int i = 2; i <= n; i++) {
if (!st[i]) prime[++cnt] = i, phi[i] = i - 1;
for (int j = 1; prime[j] <= n / i; j++) {
st[prime[j] * i] = true;
if (i % prime[j] == 0) {
phi[prime[j] * i] = prime[j] * phi[i];
break;
}
else phi[prime[j] * i] = (prime[j] - 1) * phi[i];
}
}
}
对于欧拉筛的正确性,在这里给出说明:
- 对于质数我们都知道phi[p] = p - 1,我们只需考虑合数,由线性筛知每个合数均会被筛掉且只会被其最小质因子筛掉,所以如果我们在筛掉合数时顺便把合数的欧拉函数求出,那么我们一定可以保证每个合数的欧拉函数都会被计算到且只会被计算一次。
- 筛合数时只会有两种情况,如果当前i % prime[j] == 0,那么phi[i * prime[j]] = prime[j] * phi[i],否则phi[i * prime[j]] = prime[j] * phi[i] * (1 - 1 / prime[j]) = (prime[j] - 1) * phi[i]。
欧拉函数的应用:欧拉定理
如果a与n互质,则有a ^ phi(n) = 1 (mod n)。
证明:
设与n互质的phi[n]个数为b1, b2, … , bphi[n],由于a与n互质所以ab1, ab2, … , abphi[n]也与n互质,并且者n个数在模n的意义下两两不相同, 由于与1-n中与n互质的数只有phi[n]个,所以b1b2*…b[phi[n]] = a^phi[n]b1b2…b[phi[n]] (mod n),又因为b1b2*…*b[phi[n]]与n互质,所以两边可以同时除掉,于是就得到了欧拉定理。