Calculation 2

3
4
0

0
2

数论什么的最不懂了

代码：

#include <stdio.h>
#include <string.h>
#define LL __int64
#define INF 1000000007

LL oular(LL n){
    LL i, ans = n;
    for(i = 2; i*i <= n; i++){
        if(n%i == 0){
        n/=i;
        ans -= ans/i;
        while(n%i == 0){
            n/=i; 
        }
    }
    }    
    if(n > 1) ans -= ans/n;
    return ans;
}

int main(){
    LL n;
    while(scanf("%I64d", &n), n){
        LL ans = n*(n-1)/2;
        printf("%I64d\n", (ans-(n*oular(n)/2))%INF);
    }
    return 0;
}

hdoj 3501 Calculation 2 【欧拉函数】

原文：http://blog.csdn.net/shengweisong/article/details/41014079