HDU3501

题意为求所以n的因子（除1外）的和，先求出欧拉函数p（n），p（n）*n/2求得所以小于n且与n互质的数之和，再用所有小于n的和减去它就为结果。
//欧拉函数：φ(x)=x(1-1/p1)(1-1/p2)(1-1/p3)(1-1/p4)…..(1-1/pn),其中p1, p2……pn为x的所有质因数，x是不为0的整数。
欧拉函数：P(N) 表示的是小于N且与N互质的数的个数；
小于N与N互质的数和 N*P(N)/2;

代码:

#include<iostream>
using namespace std;
#include<cstdio>
#define mod 1000000007
int main()
{
   __int64 ans,n,m,i;
    while(scanf("%I64d",&n)!=EOF)
 {
  if(n==0) break;
  ans=m=n;
  for(i=2;i*i<=n;i++)
  {
   if(n%i==0)
   {
    ans-=ans/i;
    while(n%i==0) n/=i;
    if(n==1) break;
   }
  }
  if(n!=1) ans-=ans/n;   //此时的ans就为欧拉函数值
  ans=(m*(m-1)/2-m*ans/2)%mod;
  printf("%I64dn",ans);
 }
 return 0;
}

欧拉函数
百科名片
在数论，对正整数n，欧拉函数是少于或等于n的数中与n互质的数的数目。此函数以其首名研究者欧拉命名，它又称为Euler’s totient function、φ函数、欧拉商数等。 例如φ(8)=4，因为1,3,5,7均和8互质。从欧拉函数引伸出来在环论方面的事实和拉格朗日定理构成了欧拉定理的证明。

Euler函数
通式：φ(x)=x(1-1/p1)(1-1/p2)(1-1/p3)(1-1/p4)…..(1-1/pn),其中p1, p2……pn为x的所有质因数，x是不为0的整数。φ(1)=1（唯一和1互质的数就是1本身）。 (注意：每种质因数只一个。比如12=2*2*3 　　那么φ（12）=12*（1-1/2）*(1-1/3)=4) 　　若n是质数p的k次幂，φ(n)=p^k-p^(k-1)=(p-1)p^(k-1)，因为除了p的倍数外，其他数都跟n互质.
欧拉函数是积性函数——若m,n互质，φ(mn)=φ(m)φ(n)。 　　特殊性质：当n为奇数时，φ(2n)=φ(n), 证明于上述类似。

　　设A, B, C是跟m, n, mn互质的数的集，据中国剩余定理，A*B和C可建立一一对应的关系。因此φ(n)的值使用算术基本定理便知， 　　
例如φ(72)=φ(2^3×3^2)=(2-1)2^(3-1)×(3-1)3^(2-1)=24 　　
与欧拉定理、费马小定理的关系 　　
对任何两个互质的正整数a, m, m>=2有 　
a^φ(m)≡1(mod m) 　即欧拉定理
　　当m是质数p时，此式则为： 　　
　　a^(p-1)≡1(mod m) 　　即费马小定理。