欧拉函数

参考：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);
                        }
                }
        }
}