欧拉函数（模版）

贴一下欧拉函数的模版，已备日后使用

int oula(int n) { //O(n)
    int rea = n;
    for (int i = 2; i <= n; i ++) {
        if(n % i == 0) {
            rea = rea - rea/i;
            do{
                n /= i;
            }while(n % i == 0);
        }
    }
    return rea;
}
int oula_youhua(int n){//O(√n)
    int rea = n;
    for (int i = 2; i * i <= n; i ++) {
        if(n % i == 0) {
            rea = rea - rea / i;
            do{
                n /= i;
            }while(n % i == 0);
        }
    }
    if(n > 1)
        rea = rea - rea / n;
    return rea;
}
bool Isprime[50000];
int p[maxn];
void prim() {//素数打表
    memset(Isprime, 0, sizeof(Isprime));
    Isprime[0] = Isprime[1] = 1;
    int k = 0;
    for (int i = 2; i < 50000; i ++) {
        if(!Isprime[i])
            p[k++] = i;
        for (int j = 0; j < k && i * p[j] < 50000; j ++) {
            Isprime[i *p[j]] = 1;
            if(!(i % p[j]))
                break;
        }
    }
}
//筛法打表
int phi(int n){
    int rea = n;
    for (int i = 0; p[i] * p[i] <= n; i ++) {
        if(n % p[i] == 0) {
            rea = rea - rea / n;
            do{
                n /= p[i];
            }while(n % p[i] == 0);
        }
    }
    if(n > 1)
        rea = rea - rea / n;
    return rea;
}

void ditui() {//递推打表
    int i,j;
    for(i = 1; i < maxn; i++)
        p[i] = i;
    for(i = 2;i < maxn; i++){
        if(p[i] == i){
            for(j = i; j < maxn; j += i){
                p[j] = (p[j] / i) * (i-1);
            }
        }
    }
}