NYOJ 569 最大公约数之和

最大公约数之和

时间限制： 1000 ms  |  内存限制： 65535 KB

难度： 4

描述

题目很简单，求出：

输入

每行一个数n（n<2^31)，以文件结束。

输出

每个结果占一行。

样例输入

12

样例输出

40

考虑n=12时，最大公约数有1、2、3、4、12。
gcd（1，12）=1，gcd（5，12）=1，gcd（7，12）=1，gcd（11，12）=1，
gcd（2，12）=2，gcd（10，12）=2，
gcd（3，12）=3，gcd（9，12）=3，
gcd（4，12）=4，gcd（8，12）=4，
gcd（6，12）=6，
gcd（12，12）=12，
 
所以：
 
gcd（1，12）+gcd（2，12）+gcd（3，12）+gcd（4，12）+gcd（5，12）+gcd（6，12）+
 
gcd（7，12）+gcd（8，12）+gcd（9，12）+gcd（10，12）+gcd（11，12）+gcd（12，12）
 
=4*1+2*2+2*3+2*4+1*6+1*12=40
 
我们可以发现：上式等价与
 
                          φ（12）*1+φ（6）*2+φ（4）*3+φ（3）*4+φ（2）*6+φ（1）*12
 
                       =4*1+2*2+2*3+2*4+1*6+1*12
 
                        =40 
其中φ（x)表示不超过x且与x互 素的正整数的个数，称为x的欧拉函数值 。φ函数的值 通式：φ(x)=x(1-1/p1)(1-1/p2)(1-1/p3)(1-1/p4)…..(1-1/pn),其中p1, p2……pn为x的所有质因数，x是不为0的整数。具体实现见代码。

import java.util.Scanner;


public class Main{


	public static void main(String[] args) {
		Scanner in = new Scanner(System.in);
		int n, k;
		long sum;
		while (in.hasNext()) {
			n = in.nextInt();
			sum = 0;
			for (int i = 1; i * i <= n; i++) {
               //注意平方数的处理 只有一个因子i因为另一个因子也等于i
				if (i * i == n && n % i == 0) {
					sum += euler(i) * i;


					break;
				}
                //如果i是n的因子 那么K=/i也是n的因子 分别计算两个不同因子的欧拉函数值
				if (n % i == 0) {
					k = n / i;
					sum += euler(k) * i;


					sum += euler(i) * k;


				}


			}
			System.out.println(sum);


		}


	}


	//计算n的欧拉函数值 
	private static int euler(int n) {


		int ans = n;
		for (int i = 2; i * i <= n; i++) {
			if (n % i == 0) {


				while (n % i == 0) {
					n /= i;
				}
				//以上为找出n的质因子，然后根据公式求ans 注意先除防止两个大数相乘益处


				ans = ans / i * (i - 1);
			}


		}


		if (n > 1) {//此时n是质素 出去
			ans = ans / n * (n - 1);
		}


		return ans;


	}


}