POJ 2407 Relatives 【欧拉函数】

裸欧拉函数。

#include<stdio.h>
#include<string.h>
const int N=1e5+;
int p[N],pr[N],cnt;
void init(){
    for(int i=;i<N;i++){
        if(!p[i])pr[++cnt]=i;
        for(int j=;j<=cnt&&i*pr[j]<N;j++){
            p[i*pr[j]]=;
            if(i%pr[j]==)    break;
        }
    }
}
int er(int n){
    int ans=n;
    for(int i=;i<=cnt&&pr[i]*pr[i]<=n;i++){
        if(n%pr[i]==){
            ans=ans/pr[i]*(pr[i]-);//先除后乘，防止溢出
            while(n%pr[i]==)
                n/=pr[i];
        }
    }
    if(n>)    ans=ans/n*(n-);
    return ans;
}
int main(){
    int n;init();
    while(~scanf("%d",&n)&&n){
        if(n==){
            printf("0\n");
            continue ;
        }
        printf("%d\n",er(n));
    }
    return ;
}