UVa 10299 Relatives Solution

#include<iostream>
#include<list>
#include<string>
#include<cstring>
#include<sstream>
#include<cctype>
#include<string.h>
#include<algorithm>
#include<cmath>
#include<stack>
#include<fstream>
#include<cstdlib>
#include<vector>
#include<map>
#include<utility>
#include<iomanip>
#include<queue>
using namespace std;
#define clr(a) memset(a,0,sizeof(a))
#define fill(a,v) memset(a,v,sizeof(a))
#define PB push_back
#define pi acos(-1.0)
#define eps 1e-9

bool stats[32000+10];
vector<long>prime;
long N,n,primeSize,cnt;
map<long, long>factor;


void PrimeGenerate()
{
    long i,j,sqrtN;

    for(i=4;i<=32000;i+=2)
        stats[i]=true;

    sqrtN=sqrt(32000);

    for(i=3;i<=sqrtN;i+=2)
    {
        if(!stats[i])
            for(j=i*i;j<=32000;j+=i+i)
                stats[j]=true;
    }

    for(i=2;i<=32000;i++)
        if(!stats[i])
            prime.PB(i);

    primeSize=prime.size();
}


void PrimeFactorize()
{
    long i;
    cnt=0;

    for(i=0;i<primeSize;i++)
    {
        if(n%prime[i]==0)
        {
            factor[cnt]=prime[i];
            cnt++;

            while(n%prime[i]==0)
                n/=prime[i];
        }

        if(n==1)    break;
    }
}


long CoPrime()
{
    long i;

    for(i=0;i<cnt;i++)
    {
        N/=factor[i];
        N*=(factor[i]-1);
    }

    if(n!=1)
    {
        N/=n;
        N*=(n-1);
    }

    return N;

}


int main()
{
   long result;
   PrimeGenerate();
   while(cin>>n)
   {
       if(n==0) break;
       if(n==1)
       {
           cout<<"0"<<endl;
           continue;
       }
       N=n;
       factor.clear();
       PrimeFactorize();
       result=CoPrime();
       cout<<result<<endl;

   }
return 0;
}