POJ_2407 Relatives 【欧拉函数裸题】

#include <iostream>
#include <cmath>
#include <cstring>
using namespace std;
const int MAXN = 1e5+3;
int Prime[MAXN], nPrime;
bool isPrime[MAXN];

void getPrime()
{
    memset(isPrime, 1, sizeof(isPrime));
    isPrime[0] = isPrime[1] = 0;
    nPrime = 0;
    for(int i = 2; i < MAXN; i++)
    {
        if(isPrime[i])
            Prime[nPrime++] = i;
        for(int j = 0; j < nPrime && (long long)Prime[j]*i < MAXN; j++)
        {
            isPrime[Prime[j]*i] = 0;
            if(i%Prime[j] == 0)
                break;
        }
    }
}

int Euler(int n)
{
    int ans = n;
    for(int i = 0; Prime[i]*Prime[i] <= n; i++)
    {
        if(n%Prime[i] == 0)
        {
            ans = ans - ans/Prime[i];
            do
            {
                n/=Prime[i];
            }
            while(n%Prime[i] == 0);
        }
    }
    if(n > 1)
        ans = ans - ans/n;
    return ans;
}

int main()
{
    int N;
    getPrime();
    while(cin >> N && N)
    {
        cout << Euler(N) << endl;
    }
    return 0;
}