扩展知识:
对于一个约数:
- 一个整数 N N N 的约数个数上限大约是 n \sqrt{n} n,但是 1 ∼ N 1 \sim N 1∼N 中平均每个数的约数个数大约只有 log N \log{N} logN。一般达不到上限。根据实际计算 1 0 9 10^9 109 之内的自然数中,约数的个数最多的自然数仅有1536个
- 因为一个整数 N N N 的约数可以用这个整数的质因子组成,那么就可以先进行质因数分解,然后用 d f s dfs dfs 求出 N N N的约数。这样可以试除法的时间复杂度 O ( N ) O(\sqrt{N}) O(N) 缩短到 O ( N N log N ) O(\frac{N\sqrt{N}}{\log{N}}) O(logNNN)。
思路:
- 因为 l c m ( x , b 0 ) = b 1 lcm(x,b_0) = b_1 lcm(x,b0)=b1,那么 x x x 能取到的值都是 b 1 b_1 b1 的约数,那么就可以先进行求 b 1 b_1 b1 的约数,求约数的方法在扩展知识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;
}