参考:http://blog.csdn.net/leolin_/article/details/6642096
欧拉函数:对于一个正整数n,小于n且和n互质的正整数(包括1)的个数,记作φ(n) 。
通式:φ(x)=x*(1-1/p1)*(1-1/p2)*(1-1/p3)*(1-1/p4)…..(1-1/pn),其中p1, p2……pn为x的所有质因数,x是不为0的整数。φ(1)=1(唯一和1互质的数就是1本身)。
对于质数p,φ(p) = p - 1。注意φ(1)=1.
欧拉定理:对于互质的正整数a和n,有aφ(n) ≡ 1 mod n。
欧拉函数是积性函数——若m,n互质,φ(mn)=φ(m)φ(n)。
若n是质数p的k次幂,φ(n)=p^k-p^(k-1)=(p-1)p^(k-1),因为除了p的倍数外,其他数都跟n互质。
特殊性质:当n为奇数时,φ(2n)=φ(n)
欧拉函数还有这样的性质:
设a为N的质因数,若(N % a == 0 && (N / a) % a == 0) 则有E(N)=E(N / a) * a;若(N % a == 0 && (N / a) % a != 0) 则有:E(N) = E(N / a) * (a - 1)。
模板:
long long Euler(long long n) {
long long ans = n, i;
for(i = 2; i*i <= n; ++i) {
if(n%i == 0){
ans = ans - ans/i;
while(n%i == 0)n /= i;
}
}
if(n > 1) {
ans = ans - ans / n;
}
return ans;
}
加上素数筛选的欧拉函数实现:
typedef long long ll;
const ll MAXN = 500000;
ll p[MAXN];
bool vis[MAXN];
ll prime(ll n) {
memset(vis, false, sizeof(vis));
memset(p, 0, sizeof(p));
ll pos = 0;
for(ll i = 2; i <= n; ++i) {
if(!vis[i])
p[pos++] = i;
for(ll j = 0; j < pos && i*p[j] <= n; ++j) {
vis[p[j]*i] = true;
if((i % p[j] == 0)) break;
}
}
return pos;
}
ll Euler(ll n) {
ll ans = n;
for(ll i = 0; p[i]*p[i] <= n; ++i) {
if(n % p[i] == 0) {
ans = ans - ans / p[i];
while(n % p[i] == 0) {
n /= p[i];
}
}
}
if(n > 1) {
ans = ans - ans/n;
}
return ans;
}
递推求欧拉函数:
当需要频繁地使用欧拉函数时,可以先打一个表,复杂度约为O(nlogn)
typedef long long ll;
const ll MAXN = 500000;
ll phi[MAXN];
void Euler(ll n) {
for(ll i = 1; i <= n; ++i) phi[i] = i;
for(ll i = 2; i <= n; i += 2) phi[i] /= 2;
for(ll i = 3; i <= n; i += 2) {
if(phi[i] == i) {
for(ll j = i; j <= n; j += i) {
phi[j] = phi[j] / i * (i-1);
}
}
}
}