【阶段1】【Hankson的趣味题】求约数，不含辗转相除法的做法

【题意】
Hanks博士是BT（Bio-Tech，生物技术）领域的知名专家，他的儿子名叫Hankson。
现在，刚刚放学回家的Hankson正在思考一个有趣的问题。
今天在课堂上，老师讲解了如何求两个正整数c1和c2的最大公约数和最小公倍数。
现在Hankson认为自己已经熟练地掌握了这些知识，他开始思考一个“求公约数”和“求公倍数”之类问题的“逆问题”，这个问题是这样的：
已知正整数a0,a1,b0,b1，设某未知正整数x满足：
1、 x和a0的最大公约数是a1；
2、 x和b0的最小公倍数是b1。
Hankson的“逆问题”就是求出满足条件的正整数x。
但稍加思索之后，他发现这样的x并不唯一，甚至可能不存在。
因此他转而开始考虑如何求解满足条件的x的个数。
请你帮助他编程求解这个问题。

【输入格式】
输入第一行为一个正整数n，表示有n组输入数据。
接下来的n行每行一组输入数据，为四个正整数a0，a1，b0，b1，每两个整数之间用一个空格隔开。
输入数据保证a0能被a1整除，b1能被b0整除。

【输出格式】
输出共n行。
每组输入数据的输出结果占一行，为一个整数。
对于每组数据：若不存在这样的x，请输出0；
若存在这样的x，请输出满足条件的x的个数；

【数据范围】
1≤n≤2000,
1≤a0,a1,b0,b1≤2∗10^9

【输入样例】
2
41 1 96 288
95 1 37 1776

【输出样例】
6
2

我们先把a0、a1、b0、b1都分解出质因数（xi表示第i个素数）

就会得到a0=,a1=

       以及b0=,b1=

假如我们要求的 X=

因为a1是a0和x的最大公因数，所以 (其中i取1~n)

因为b1是b0和x的最小公倍数，所以 (其中i取1~n)

• 由此可知，
• 当时，可以取大于等于的任意自然数
• 当时，说明的min取了的值，被固定
• 当时，可以取小于等于的任意自然数
• 当时，说明的max取了的值，被固定
• 因此，当被固定且合法时，ans在这一位上没有任何选择，依然等于本身
• 当没有被固定且合法时，ans在这一位上有选择，根据乘法原理，ans=ans*这一位的选择数
• 什么时候不合法呢？就是在被固定成两个不同的数时，这是无法满足的，此时ans=0，break

注意

• 我们用线性筛来预处理出质数，但是因为它们<=2*10^9，所以我们没办法求出它们所有的质因数
• 但是巧妙的地方来了，由于最小公倍数是<=2*10^9次方的，所以我们要求的X也是<=2*10^9的
• 线性筛没办法求2*10^9之内的所有素数，但是求之内的绰绰有余
• 但是X是<=2*10^9，所以它如果有大于的质因数，也有且只有1个
• 因为一旦有两个大于的质因数，它就>2*10^9了
• 因此我们放心地用线性筛出来的质数去分解它们，如果没分干净（最后的值不等于1）说明，剩下来的那个数就是大于的质因数，我们不用刻意去求，它自己就浮现出来了（太女少了）
• 然后就很简单啦

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
inline ll read()
{
	char c=getchar();ll s=0,f=1;
	while(!isdigit(c)){if(c=='-')f=-1;c=getchar();}
	while(isdigit(c)){s=s*10+c-'0';c=getchar();}
	return f*s;
}
ll a1,a2,b1,b2;
bool bk[51000000];
int tot,pri[4100];
inline void pre()
{
	memset(bk,true,sizeof(bk));
	for(int i=2;i<=200000;i++)
	{
		if(bk[i])pri[++tot]=i;
		if(tot>4000)break;
		for(int j=1;j<=tot && pri[j]*i<=200000;j++)
		{
			bk[i*pri[j]]=false;
			if(i%pri[j]==0)break;
		}
	}
}
ll ta1[4100],ta2[4100],tb1[4100],tb2[4100];
ll xia[4100],shang[4100];
int main()
{
	//freopen("input10.in","r",stdin);
	pre();int T=read();
	while(T--)
	{
		a1=read();a2=read();
		b1=read();b2=read();
		memset(ta1,0,sizeof(ta1));
		memset(ta2,0,sizeof(ta2));
		memset(tb1,0,sizeof(tb1));
		memset(tb2,0,sizeof(tb2));
		for(int i=1;i<=tot;i++)
		{
			if(a1==1)break;
			while(a1!=1 && a1%pri[i]==0)a1/=pri[i],ta1[i]++;
		}
		for(int i=1;i<=tot;i++)
		{
			if(a2==1)break;
			while(a2!=1 && a2%pri[i]==0)a2/=pri[i],ta2[i]++;
		}
		for(int i=1;i<=tot;i++)
		{
			if(b1==1)break;
			while(b1!=1 && b1%pri[i]==0)b1/=pri[i],tb1[i]++;
		}
		for(int i=1;i<=tot;i++)
		{
			if(b2==1)break;
			while(b2!=1 && b2%pri[i]==0)b2/=pri[i],tb2[i]++;
		}
		ll ans=1;
		for(int i=1;i<=tot;i++)
		{
			if(ta1[i]!=ta2[i])xia[i]=-1;
			else xia[i]=ta1[i];
			if(tb1[i]!=tb2[i])shang[i]=-1;
			else shang[i]=tb1[i];
			if(xia[i]==-1 && shang[i]==-1 && min(ta1[i],ta2[i])!=max(tb1[i],tb2[i])){ans=0;break;}
			if(xia[i]!=-1 && shang[i]!=-1)ans*=(shang[i]-xia[i]+1);
		}
		if(a1!=1 || a2!=1 || b1!=1 || b2!=1)
		{
			if(b1==b2 && b1!=1)ans<<=1;
		}
		printf("%lld\n",ans);
	}
	return 0;
}