「NOIP 2009」Hankson 的趣味题

传送门

---

problem

给定 $a_0,a_1,b_0,b_1$，求满足下列条件的 $x$ 的数量：

• $\gcd (x,a_0)=a_1$。
• $\mathrm{l}\mathrm{c}\mathrm{m}(x,b_0)=b_1$。

数据范围：$1\le a_0,a_1,b_0,b_1\le 2\times 10^9$，$T\le 2000$，其中 $T$ 是数据组数。

---

solution

这道题不难，主要就是 $\gcd$ 的一个小性质：$\gcd (x,y)=k⟹\gcd (\frac{x}{k},\frac{y}{k})=1$。

那么第一个式子就是 $\gcd (\frac{x}{a_1},\frac{a_0}{a_1})=1$。

第二个式子，稍微化简一下可得：

$$\begin{array}{rl}\frac{x{b}_0}{\gcd (x,b_0)}& =b_1\\ \gcd (x,b_0)& =\frac{x{b}_0}{b_1}\\ \gcd (\frac{b_1}{b_0},\frac{b_1}{x})& =1\end{array}$$

由于 $x$ 是 $b_1$ 的因数，我们就可以 $O(\sqrt{b_1})$ 枚举因数来解决了。

时间复杂度 $O(T\sqrt{b_1})$。

---

code

#include<bits/stdc++.h>
using namespace std;
int T,a0,a1,b0,b1;
bool check(int x){
	if(x%a1)  return false;
	if(__gcd(b1/b0,b1/x)!=1)  return false;
	if(__gcd(x/a1,a0/a1)!=1)  return false;
	return true;
}
int main(){
	scanf("%d",&T);
	while(T--){
		int ans=0;
		scanf("%d%d%d%d",&a0,&a1,&b0,&b1);
		for(int i=1;i*i<=b1;++i){
			if(b1%i==0){
				ans+=check(i);
				if(i*i!=b1)  ans+=check(b1/i);
			}
		}
		printf("%d\n",ans);
	}
	return 0;
}