最大公约数之和
-
描述
-
题目很简单,求出:
-
输入
- 每行一个数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; } }