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.