/** 改题目是对于一道在线笔试题目的简化版本,问题描述有些不同,望读者以本文开头的问题描述为主,进行下文的观看**/
参考原文: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();
}