GCD and LCM（hdu4497）

传送门

题意

给出$gcd(x,y,z)，lcm(x,y,z)$,求出所有满足条件的$(x,y,z)$组合的数量。注意：$(1,2,3)和(1,3,2)$是两个不同的组合。

思路

1. 明显地，若$lcm\%gcd\ne 0$,则无解。
2. 设$d=lcm/gcd,d=q_1^{c_1}\times q_2^{c_2}\times q_3^{c_3}\times \cdots \times q_m^{c_m}$
3. 求满足$gcd(x,y,z)，lcm(x,y,z)$的组合数等价于求满足$gcd(a,b,c)=1，lcm(a,b,c)=d$的组合数。此处a,b,c实质上是gcd的倍数($a=x/d,b=y/d,c=z/d$)。
4. $$a=q_1^{A_1}\times q_2^{A_2}\times \cdots \times q_m^{A_m}$$
   $$b=q_1^{B_1}\times q_2^{B_2}\times \cdots \times q_m^{B_m}$$
   $$c=q_1^{C_1}\times q_2^{C_2}\times \cdots \times q_m^{C_m}$$
   1. 因为$gcd(a,b,c)=1$,所以$min(A_i,B_i,C_i)=0,max(A_i,B_i,C_i)=c_i$。
   2. $(A_i,B_i,C_i)的组合数=3+3+(c_i-1)\times A_3^3=6c_i$
   3. $两个0一个c_i，两个c_i一个0，一个0一个c_i剩下一个任选[1,c_i-1]中的一个。$
5. 综上：$ans=6c_1\times 6c_2\times \cdots \times 6c_m$

#include<cstdio>
#include<queue>
#include<set>
#include<cstdlib>
#include<string.h>
#include<string>
#include<iostream>
#include<cmath>
#include<unordered_map>
#include<map>
#include<algorithm>
#define endl "\n"
#define IOS ios::sync_with_stdio(0), cin.tie(0),cout.tie(0)
#define ft first
#define sd second
#define pll pair<ll, ll>
#define pii pair<int, int>
#define ll long long int
#define ull unsigned long long int
#define mt(a,b) memset(a, b, sizeof a)
//#define int long long
const double PI = acos(-1.0);
const int inf = 0x3f3f3f3f;
const int INF = 0x7fffffff;
using namespace std;
const int N = 1e5 + 7, M = 1e6 + 10;
ll c[100], m;

void div(ll n)
{
	m = 0;
	int q = sqrt(n);
	for (int i = 2; i <= q; i++)
	{
		if (n % i == 0)
		{
			m++;
			c[m] = 0;
			while (n % i == 0) n /= i, c[m]++;
		}
	}
	if (n > 1) c[++m] = 1;
}

int main()
{
	IOS;
	int T; cin >> T;
	while (T--)
	{
		ll gcd, lcm; cin >> gcd >> lcm;

		if (lcm % gcd) cout << 0 << endl;
		else
		{
			ll ans = 1;
			ll d = lcm / gcd;
			div(d);
			
			for(int i = 1; i <= m; i++)
				ans *= 6 * c[i];
			cout << ans << endl;
		}
	}


	return 0;
}