hdu-Calculation 2(欧拉函数)

最新推荐文章于 2019-08-21 12:01:48 发布

双之城

最新推荐文章于 2019-08-21 12:01:48 发布

阅读量507

点赞数

分类专栏：数论

本文链接：https://blog.csdn.net/u011026037/article/details/23118885

版权

数论专栏收录该内容

46 篇文章 0 订阅

订阅专栏

众所周知若（n,x）= 1,则（n,n-x）=1;

所以只需求出n的欧拉函数值，则1~n-1中与n互素的和为m = Euler（n）* n / 2;

那么不互素的和为n*（n-1）/ 2 - m;

#include<stdio.h>
#include<math.h>
typedef __int64 LL;
const LL mod = 1000000007;
LL Euler(LL x)
{
    LL res = x;
    for(int i = 2;i <= sqrt(x);i++)
    {
        if(x%i==0)
        {
            while(x%i==0) x/=i;
            res = res/i*(i-1)%mod;
        }
    }
    if(x > 1) res = res/x*(x-1)%mod;
    return res;
}
int main()
{
    LL n;
    while(scanf("%I64d",&n),n)
    {
        LL m = Euler(n);
        printf("%I64d\n",(((n*(n-1)/2)%mod - m*n/2%mod)+mod)%mod);
    }
}