连续整数序列中计算三元互质子集的个数

/** 改题目是对于一道在线笔试题目的简化版本，问题描述有些不同，望读者以本文开头的问题描述为主，进行下文的观看**/

参考原文：https://blog.csdn.net/lov_vol/article/details/51174769 

Description

给出一个整数n，表示1，2，...，n。从这n个数中任意选择3个不同的数字x，y，z，问x，y，z的最大公约数等于m的方案有多少种？（注意：(1,2,3)，(1,3,2)，(2,1,3)，(2,3,1)，(3,1,2)，(3,2,1)属于同一种方案）

Input

第一行输入两个整数，n, m (1<= m<=n <=10^5) 

Output

输出一个整数表示答案

解题思路：

最大公约数(Greatest common divisor) ，问题求解情况 | gcd(x,y,z) ==m | 的个数，则容易得到x, y, z 均是m的倍数，故能够得到gcd==m 的数一定是在m的倍数中去寻找（1*m, 2*m, 3*m, .... d*m, | d = n/m)。故问题可以简化为: 求解 1~d 中任意三个数字，器最大公约数为1（这个转化是最关键的）。现在最初问题的 n == d，m == 1。根据排列组合，1~d中所以组合个数为C(d,3)，那么 | gcd(x,y,z | x,y,z  1~d ) == 1| 的个数为： C(d,3) - | gcd==(2,3,4,...,d) | 。

用dp[i] 表示 gcd == i 的个数情况，1~d中 i 的倍数有 b = d / i 个 ，故 dp[i] = C(b, 3) - dp[2*i, 3*i,....,b*i] .

public static void main(String[] args) {
		// TODO Auto-generated method stub
		Scanner sc=new Scanner(System.in);
		int[] num = new int[2];
		for (int i = 0; i < 2; i++) {
            num[i] = sc.nextInt();
        }
		int n = num[0];
		int m = num[1];
		n = n / m ;
		int ans = n*(n-1)*(n-2)/6;
		int[] dp = new int[10000];
		for(int i = n; i>=2; i--) {
			int a = n/i;
			dp[i] = a * (a-1) * (a-2) / 6;
			int b = 2*i;
			while(b<= n)
			{
				dp[i] -=dp[b];
				b += i;
			}
			
		}
		for(int i=2;i<=n;i++) {
			ans-=dp[i];
		}
		System.out.println("----------"+ans+"-------");
		sc.close();

		

	}