欧拉phi函数

#include <cstdio>
#include <cstring>
#include <cmath>
// 计算欧拉phi函数。phi(n)是不超过n且与n互素的正整数个数
int euler_phi(int n){
    int m = (int)sqrt(n + 0.5);
    int ans = n;
    // phi(n)=n(1-1/p1)(1-1/p2)...(1-1/pk)
    for(int i = 2 ; i <= m ; 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;
}
// 用类似筛法的方式计算phi(1)、phi(2)、...、phi(n)
#define maxn 1000
int phi[maxn];
void phi_table(int n){
    for(int i = 2 ; i <= n ; i++) phi[i] = 0;
    phi[1] = 1;
    for(int i = 2; i <= n ; i++) if(!phi[i])
        for(int j = i ; j <= n ; j += i){
            if(!phi[j]) phi[j] = j;
            phi[j] = phi[j] / i * (i - 1);
        }
}
int main(){
    phi_table(999);
    for(int i = 1 ; i <= 999 ; i++)
        if(phi[i] != euler_phi(i))
            printf("wrong at %d
",i);
    return 0;
}