AcWing 200. Hankson的趣味题(最大公约数)

Acwing 200. Hankson的趣味题

扩展知识：
对于一个约数：

1. 一个整数 $N$ 的约数个数上限大约是 $\sqrt{n}$，但是 $1\sim N$ 中平均每个数的约数个数大约只有 $\log N$。一般达不到上限。根据实际计算 $10^9$ 之内的自然数中，约数的个数最多的自然数仅有1536个
2. 因为一个整数 $N$ 的约数可以用这个整数的质因子组成，那么就可以先进行质因数分解，然后用 $dfs$ 求出 $N$的约数。这样可以试除法的时间复杂度$O(\sqrt{N})$ 缩短到$O(\frac{N\sqrt{N}}{\log N})$。

思路：

1. 因为 $lcm(x,b_0)=b_1$,那么 $x$ 能取到的值都是 $b_1$ 的约数，那么就可以先进行求 $b_1$ 的约数，求约数的方法在扩展知识2中已经介绍过了。
2. 然后再进行对求得的约数进行遍历，看哪一个满足就 res++，那么得到的最终的值就是答案。

注意：
该题的代码不好实现，容易出错，建议读者自己按照思路敲一遍，不会的话再来看代码。

代码实现：

#include<bits/stdc++.h>
#define int long long
using namespace std;
typedef pair<int,int> PII;
const int N = 500005;

int primes[N],used[N],cnt;
vector<PII> sum;
unordered_set<int> st;

void get_primes(int x) {
    for(int i =2 ; i <= x; i++)
    {
        if(!used[i]) primes[++cnt] = i;
        for(int j = 1; i * primes[j] <= x; j++)
        {
            used[i * primes[j]] = true;
            if(i % primes[j] == 0) break;
        }
    }
}
void dfs(int u,int p)
{
    if(u >= sum.size())
    {
        st.insert(p);
        return ;
    }
    for(int i = 0; i <= sum[u].second; i++)
    {
        dfs(u+1,p);
        p *= sum[u].first;
    }
}
int gcd(int a,int b)
{
    return b ? gcd(b,a % b) : a;
}
signed main ()
{
    int T;
    cin >> T;
    get_primes(50000);
    while(T--)
    {
        int a1,a2,b1,b2;
        cin >> a1 >> a2 >> b1 >> b2;
        int t = b2;
        st.clear();
        sum.clear();
        for(int i = 1; primes[i] <= t/primes[i]; i++) {
            if (t % primes[i] == 0) {
                int len = 0;
                while (t % primes[i] == 0) {
                    len++;
                    t /= primes[i];
                }
                sum.emplace_back(primes[i], len);
            }
        }
        if(t > 1) sum.emplace_back(t,1);
        dfs(0,1);
        int res = 0;
        for(int i : st)
        {
            if(gcd(a1,i) == a2 && i * b1/gcd(i,b1) == b2)
                res ++;
        }
        cout << res << endl;
    }
    return 0;
}