洛谷 P1072 Hankson的趣味题【GCD | LCM】

题目链接

题目描述：

给出四个正整数 $a_0,a_1,b_0,b_1$，设某正整数 x 满足：

1. $gcd(x,a_0)=a_1$
2. $lcm(x,b_0)=b_1$

输出满足条件的 x 的个数
其中 $1\le a_0,a_1,b_0,b_1\le 2\ast 10^9$ 且 n ≤ 2000。

---

题解：

朴素想法：

1. 由 $gcd(x,a_0)=a_1$ 我们可知，$a_1$ 是 x 的因子。
2. 由 $lcm(x,b_0)=b_1$ 我们可知，x 是 $b_1$ 的因子。
3. 所以我们想到用试除法求出 $b_1$ 所有的约数 x，逐一判断 x 是否满足这两个条件。
4. 但这样有一个点会超时，所以我们对 gcd 和 lcm 的过程进行如下优化：

gcd:

$gcd(x,a_0)=a_1$
—> $gcd(x/a_1,a_0/a_1)=1$

lcm:

$lcm(x,b_0)=b_1$
—> $x\ast b_0/gcd(x,b_0)=b_1$
—> $gcd(x,b_0)=x\ast b_0/b_1$
—> $gcd(b_1/b_0,b_1/x)=1$

这样就优化了gcd 和 lcm ，就成功通过最后一个点啦

---

AC code:

#include <iostream>
#include <vector>
#include <queue>
#include <stack>
#include <set>
#include <map>
#include <unordered_set>
#include <unordered_map>
#include <deque>
#include <list>
#include <iomanip>
#include <algorithm>
#include <fstream>
#include <cmath>
#include <cstdio>
#include <cstring>
#include <cstdlib>
//#pragma GCC optimize(2)
using namespace std;
typedef long long ll;
//cout << fixed << setprecision(2);
//cout << setw(2);
const int N = 2e5 + 6, M = 1e9 + 7;

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

int main() {
    //freopen("/Users/xumingfei/Desktop/ACM/test.txt", "r", stdin);
    ios::sync_with_stdio(false);
    cin.tie(0), cout.tie(0);
    int T;
    cin >> T;
    while (T--) {
        int a, b, c, d, ans = 0;
        cin >> a >> b >> c >> d;
        for (int i = 1; i * i <= d; i++) {
            if (d % i == 0) {
                if (i % b == 0 && gcd(a / b, i / b) == 1 && gcd(d / c, d / i) == 1) ans++;
                int j = d / i;
                if (j == i) continue;
                if (j % b == 0 && gcd(a / b, j / b) == 1 && gcd(d / c, d / j) == 1) ans++;
            }
        }
        cout << ans << '\n';
    }
    return 0;
}