HDU3501-Calculation 2-欧拉函数的性质的应用

HDU：​​http://acm.hdu.edu.cn/showproblem.php?pid=3501​​

性质①：N>1，不大于N且和N互素的所有正整数的和是 ( n*(n-1)/2 - n*eular(n)/2 ) % mod。

#include<iostream>
#include<cmath>
using namespace std;
//求单个数的欧拉函数
long long eular(long long n){
  long long ans=n;
  for(int i=2;i*i<=n;i++){
    if(n%i==0){
      ans-=ans*1.0/i;
      while(n%i==0){
        n/=i;
      }
    }
  }
  if(n>1) ans-=ans/n;
  return ans;
}
int main(){
  long long n;
  while(cin>>n){
    if(n==0) break;
    cout<<( n*(n-1)/2 - n*eular(n)/2 )%1000000007<<endl;
  }
}

注意：数据类型要用long long或者_int64.