HDOJ 3501 Calculation 2

Calculation 2

这个题需要欧拉函数 

下面这个博客很好可以看看

欧拉函数最全总结_jiet07的博客-CSDN博客_欧拉函数

做此题需要使用欧拉函数的编程方法：

φ(N)=N∗∏(1−1/p)（P是数N的质因数）

φ(N)是小于N的数中与N互质的个数。这些N互质的数的和是φ(N)*N/2;

根据本题要求 就计算 1~N-1的和减去φ(N)*N/2再取余。

我写的时候没注意运算顺序什么的一些小问题。一直再最后计算出错

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll cacu(ll n)
{
    ll ans=n;
    for(ll i=2;i*i<=n;i++)
    {
        if(n%i==0)
        {
            ans=ans/i*(i-1);
            while (n%i==0)
                n/=i;
        }
    }
    if(n>1) ans=ans/n*(n-1);
    return ans;
}
int main()
{
    //freopen("C:\\Users\\Lenovo\\Desktop\\input.txt","r",stdin);
    ll n=0;
    while (cin>>n)
    {
        if(n==0)break;
        else
        {
            ll a=n*(n-1)/2;ll b=cacu(n)*n/2;
            ll output=(a-b)%1000000007 ;
            printf("%lld\n",output);
        }
    }
    return 0;
}