題意：GCD（m，n）是指n和m的最大公約數，現在給你個數n，找所有比n小的m，m滿足GCD（m,n）！=1且GCD（m,n）！=m，就是讓你找所有既不是因子也不是互質的數。比如6的因子是1,2,3,6 而且5又和6互質，所以只有4滿足條件。
歐拉公式：假設p[t]爲所有n的素數因子，φ（n）=n/p[1]*(p[1]-1)/p[2]*(p[2]-1).../p[t]*(p[t]-1)。
素數因子的公式： 假如n=p[1]^(a1)*p[2]^(a2)*p[3]^(a3)*...p[t]*(at),那麼公約數就有：(a1+1)*(a2+1)*(a3+1)*...(at+1)個。
注意：因爲1既和n互質又是n的因子，所以多減的1要補回來。
思路：最大值是2的31次方，開平方根是4.6W多，保險點可以拿5W做爲素數的上限，求5W以內的所有素數以後，然後從第一個判斷是不是素數因子，然後除到沒有這個因子爲止，再判斷後面的素數是不是因子，直到大於5W還不是因子，說明最後這個數是大於5W的素數因子，而且這個素數因子只有1次，於是直接作爲p[t]處理。
上波代碼：
#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>
#include <vector>
using namespace std;
#define N 50000
vector <int> prime;
int IsNotPrime[N+10],prime_num=0;
void makeprime()
{
for(int i=2;i<N;i++)
{
if(!IsNotPrime[i])
{
prime.push_back(i);
prime_num++;
}
for(int j=0;j<prime_num&&prime[j]*i<=N;j++)
{
IsNotPrime[prime[j]*i]=1;
if(i%prime[j]==0)
break;
}
}
}
int NotRes(int n)
{
int res1=n,res2=1;
if(n==1)
return 1;
for(int i=0;i<prime_num&&n!=1;i++)
{
if(n%prime[i]==0)
{
res1=(res1/prime[i])*(prime[i]-1);n/=prime[i];
int c=1;
while(n%prime[i]==0)
{
n/=prime[i];
c++;
}
res2*=c+1;
}
}
if(n>1)
{
res1=(res1/n)*(n-1);
res2*=2;
}
return res1+res2-1;
}
int main()
{
int n;
makeprime();
cin.sync_with_stdio(false);
while(cin>>n&&n)
{
cout<<n-NotRes(n)<<"\n";
}
return 0;
}