数论：最大公约数与最小公倍数的应用：最大公约数和最小公倍数问题

总结：

参考文献：https://sochiji.blog.luogu.org/solution-p1029

1.a * b == (a,b) * [a,b] (a,b的最大公约数 * a,b的最小公倍数)

2.如果a的素因数分解为 (x1 ^ p1) * (x2 ^ p2) * (x3 ^ p3)

b的素因数分解为 (x1 ^ p4) * (x5 ^ p5) * (x6 ^ p6)

可以增添一下素因数a : (x1 ^ p1) * (x2 ^ p2) * (x3 ^ p3) *(x5 ^ 0) * (x6 ^ 0)

b:(x1 ^ p4) * (x5 ^ p5) * (x6 ^ p6)* (x2 ^ 0) * (x3 ^ 0)

则a,b的最大公约数为:(x1 ^ min(p1,p4) ) * (x2 ^ min (p2 ,0) ) * (x3 ^ min(p3,0) ) * (x5 ^ min(0,p5)) * (x6 ^ min(0.p6) ).

a,b的最小公倍数为：(x1 ^ max(p1,p4) ) * (x2 ^ max (p2 ,0) ) * (x3 ^ max(p3,0) ) * (x5 ^ max(0,p5)) * (x6 ^ max(0.p6) ).

3.由2可以知道，给定一个最小公倍数 x , 和一个最大公约数 y. 如果 x / y 不能够整除的话，则不存在两个数a,b使得 (a,b) == y ,  [a,b] == x.

证明：a 能够整除 y , x 能够整除 a , 则 x必定能够整除y. 但是现在 x 不能整除 y 则矛盾。

题目链接：

https://www.luogu.com.cn/problem/P1029

题目：

题目描述

输入两个正整数 x_0, y_0x0​,y0​，求出满足下列条件的 P, QP,Q 的个数：

1. P,QP,Q 是正整数。

2. 要求 P, QP,Q 以 x_0x0​ 为最大公约数，以 y_0y0​ 为最小公倍数。

试求：满足条件的所有可能的 P, QP,Q 的个数。

输入格式

一行两个正整数 x_0, y_0x0​,y0​。

输出格式

一行一个数，表示求出满足条件的 P, QP,Q 的个数。

输入输出样例

输入 #1复制

3 60

输出 #1复制

4

说明/提示

P,QP,Q 有 44 种：

1. 3, 603,60。
2. 15, 1215,12。
3. 12, 1512,15。
4. 60, 360,3。

对于 100\%100% 的数据，2 \le x_0, y_0 \le {10}^52≤x0​,y0​≤105。

【题目来源】

NOIP 2001 普及组第二题

思路1：按照上面的总结2，3进行求解

 给了一个最大公因数x,最小公倍数y.

1.于是我们求出 y / x. 如果 y / x 不能整除，则由总结3可以知道，不存在这样一组a,b组合，得到最大公因数x,最小公倍数为y

2.如果能够整除， 求出其对应的质因数。

而这些质因数由总结2可以知道，可以随机分给a或者b。 但是质因数一次有且只能分给一个。即质因数x1给了a,则当下不能分给b.

代码实现：

# include <iostream>
using namespace std;

int main()
{
    int x,y;
    scanf("%d %d",&x,&y);
    int temp = y / x;
    
    if(x * temp != y)
    {
        printf("0\n");
        return 0;
    }
    
    int res = 1;
    for(int i = 2 ; i <= temp / i; i++)
    {
        if(temp % i == 0)
        {
            res *= 2;
            while(temp % i == 0)
            {
                temp /= i;
            }
        }
    }
    if(temp >1)
    {
        res *= 2;
    }
    printf("%d\n",res);
    return 0;
}

思路2：根据总结1 + 枚举求解：

我们知道，a,b最大公因数  * a,b最小公倍数 == a * b

现在已知最大公因数，最小公倍数。

可以求出temp = 最大公因数 * 最小公倍数 .

然后从1~temp值作为a值进行枚举。b值就是 temp / a.

则满足条件的情况为：

1.必须要temp / a 能够整除

2.必须gcd(a,b) == 对应输入的最大公因数。

但是数据范围是 10 ^ 5. 则temp 最大为1e10。

从1~temp枚举必定会超时。

这个时候就可以用到a * b 要 == temp,

则 a,b中最多只有一个值 大于 sqrt(temp)

所以我们只要枚举 1 ~ sqrt(temp)作为a即可。

而根据题目样例 3,60 和 60 , 3为两种不同的情况。 所以最后乘以2.

但是还要注意一点， 那就是 a,b 都为 sqrt(temp)的情况。

这个时候a,b和 b,a互换 实际上还是为一种，所以这种情况需要进行特判。

代码实现：

# include <iostream>
using namespace std;

long long gcd(long long a , long long b)
{
    return b ? gcd(b, a % b) : a;
}

int main()
{
    long long x,y;
    scanf("%lld %lld",&x,&y);
    long long temp = x * y;

    int res = 0;

    for(long long i = 1 ; i <= temp / i ; i++)
    {
        if(temp % i == 0 && gcd(temp / i , i) == x)
        {
            if(temp / i == i)
            {
                res += 1;
            }
            else
            {
                res += 2;
            }
        }
    }
    printf("%d\n",res);
    return 0;
}