POJ2407求单个数的欧拉函数值

题意：
求n的欧拉函数值。
题外话：
今天本来打算学DP的，发现第一道求m个子序列就把我打趴了一整天，所以先学完数论再去战斗DP。。ORZ，好，下面来说一下欧拉函数吧。
欧拉函数
在数论中，对正整数n，欧拉函数是少于或等于n的数中与n互质的数的数目。（不是互质的素数数目，当初我看错了，以为是互质的素数数目，蒙了半天）。
n的所有素因子和就是euler[n]*n/2;
欧拉函数的求法就是euler[n]=n*(1-1/p1)(1-1/p2)……*（1-1/pk）p1,p2,pk是n的素因子，就是对于每个n的素因子都有n*1/pi个小于等于n的因子，要去除它们，但是因子之间可能会有重复，所以要连乘，每次去除的都是在上一次的基础上剩余的因子。
欧拉函数这个算法主要是看这位传说中的数论大佬的博客的：
http://blog.csdn.net/acdreamers/article/details/8007991

不多说，上我看了几个博客的模板。

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
int phi(int n)
{
    int res=n;//假设现有n个数与n互质，开始筛选剔除
    for(int i=2;i*i<=n;i++)//i*i应该是相当于sqrt(n),类似我们在求一个数是否为素数时候用的优化手段。 
    {
        if(n%i==0)//若这个数是n的因子，减去n以下含有这个因子的数个数，假设n=8，小于等于8，2为公因子的有8/2=4个
        {
            res-=res/i;
            while(n%i==0)
            n=n/i;
        }
    }
    if(n>1)   //若n大于1，则此时的n也是一个除1以外的因子
    res-=res/n;
    return res;
}
int main()
{
    int n;
    while(~scanf("%d",&n),n)
    {
        printf("%d\n",phi(n));
    }
}