UVa 10179 -- Irreducable Basic Fractions

题目大意：题目大意：给一个正整数n，求出在[1, n]区间内和n互质的正整数的个数。Euler's Totient（欧拉函数）的直接应用。

欧拉函数：计算 1 到 n 的正整数当中，跟 n 互质（最大公因数是一）的数，总共有几个。1.先将n做质因数分解：n=p1^a1*p2^a2*...*pm^am;(p1,p2,...,pm为n的质因数)；2.以质因数计算欧拉数公式为n*(1-1/p1)*(1-1/p2)*...*(1-1/pm); 3.整理为：n/p1*(p1-1)/p2*(p2-1)/.../pm*(pm-1);

思路分析：

1.先建立素数表，运用线性筛选法；

2.将输入的数对素数表进行一一试除，找出n的质因数，此时应该注意对于每个质因数都应除尽（注意欧拉公式）

  对于 if(n>1) ans=ans/n*(n-1); 这部分如果筛了全部素数，那个if(n>1)是不需要判断的。通常，素数不会全部筛出来的，如果最大值是n，筛到sqrt(n)就可以，这时候，      sqrt(n)到n之间的范围内，会有一个素数。

代码实现：

#include<stdio.h>
#include<string.h>
#include<math.h>
const int maxn=100000;
int prime[maxn/3];
int flag[maxn];
int cnt;
void Prime(){
    cnt=0;
    memset(flag,0,sizeof(flag));
    for(int i=2;i<maxn;i++){
        if(!flag[i])
            prime[cnt++]=i;
        for(int j=0;i*prime[j]<maxn;j++){
            flag[i*prime[j]]=1;
            if(i%prime[j]==0)
                break;
        }
    }
    //printf("%d\n",cnt);
    //for(int i=0;i<30;i++)
        //printf("%d ",prime[i]);
}
void euler(int n){
    int ans=n;
    for(int i=0;i<cnt&&n>1;i++){
        if(n%prime[i]==0){
            ans=ans/prime[i]*(prime[i]-1);
            while(n%prime[i]==0)
                n/=prime[i];
        }
    }
    if(n>1) ans=ans/n*(n-1);
    printf("%d\n",ans);
}
int main(){
    int n;
    Prime();
    while(scanf("%d",&n)==1){
        if(n==0)
            break;
        euler(n);
    }
}