求助

#include<iostream>
using namespace std;

int n, ans, a0, a1, b0, b1;

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

int lcm(int a, int b) {
	return a/gcd(a, b)*b;
}

int main() {
	cin >> n;
	while (n--) {
		ans=0;
		int i;
		cin >> a0 >> a1 >> b0 >> b1;
		for (i=a1; i*i<=b1; i++) {
			if (b1%i==0) {
				if (gcd(i, a0)==a1 && lcm(i, b0)==b1) ans++;
				int tmp=b1/i;
				if (tmp==i) continue;
				if (gcd(tmp, a0)==a1 && lcm(tmp, b0)==b1) ans++;
				
			}
		}
		cout << ans << endl;
	}
	return 0;
}