Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 2548 Accepted Submission(s): 1064
3 4 0
0 2
#include<stdio.h>
#define mod 1000000007
__int64 eular(__int64 n)//欧拉函数模板
{
__int64 res=n,i;
for(i=2;i*i<=n;i++)
{
if(n%i==0)
res=res/i*(i-1);//先进行除法是为了防止中间数据的溢出
while(n%i==0)
n/=i;
}
if(n>1)
res=res/n*(n-1);//保证x一定是素数
return res;
}
int main()
{
__int64 n,ans;
while(scanf("%I64d",&n)&&n)
{
ans=n*(n-1)/2;//求出所有小于n的正整数和
printf("%I64d\n",(ans-n*eular(n)/2)%mod);
}// n*eular(n)/2为小于n且与n互质的数的和,因为若m与n互质,那么n-m也一定与n互质
return 0;
}HDOJ 3501 Calculation 2(欧拉函数拓展——求非互质数和)
原文:http://blog.csdn.net/zwj1452267376/article/details/46404765