P1072 [NOIP2009 提高组] Hankson 的趣味题

知识点：质数，约数

我说实话，这个题感觉李煜东有点误导人了，明明有更简单的方法，为啥书上写的一坨，这个就只是一个绿题而已，结果书上搞得有多难多难，这个题其实用到后面挨着的一个概念会非常简单，就是互质，前面讲最大公约数的时候已经证明了一个结论，那就是两个数都除以他们的最大公约数，得到的两个数是互质的，也就是最大公约数是1，同理我们对最小公倍数那一组进行处理，能得到类似的结论，b1/b0与b1/x互质，关键点就在于这个式子，显然，x是b1的约数，那么我们就求b1的约数，看看能不能满足a1是x的约数，x/a1与a0/a1互质，b1/b0与b1/x互质这三个条件，满足了就加一，这样这个题就完成了，时间复杂度应该是根号乘以log，外层循环是寻找约数是根号，里面是求最大公约数时间复杂度是log，这个题就这么简单，李煜东的书上写的真的复杂，关键我去看了acwing，闫雪灿的也是复杂的一笔

#include <bits/stdc++.h>

using namespace std;

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

int main() {
	int T;
	scanf("%d", &T);
	while (T--) {
		int a0, a1, b0, b1;
		scanf("%d%d%d%d", &a0, &a1, &b0, &b1);
		int ans = 0, t1 = a0 / a1, t2 = b1 / b0;
		for (int i = 1; i <= (int) sqrt(b1); i++) {
			if (b1 % i) continue;
			if (!(i % a1) && gcd(i / a1, t1) == 1 && gcd(b1 / i, t2) == 1) ans++;
			if (i * i != b1 && !(b1 / i % a1) && gcd(b1 / i / a1, t1) == 1 && gcd(b1 / (b1 / i), t2) == 1) ans++;
		}
		printf("%d\n", ans);
	}
	return 0;
}