UVA, 11064 Number Theory

题意：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;
}