#include <bits/stdc++.h>

using namespace std;
int prime[50010];
bool isprime[10000010];
int gcd(int a,int b)
{
	if(b==0)
		return a;
	return gcd(b,a%b);
}
int pf0[40][2],cnt0,pf1[40][2],cnt1;
int pf[40][2],cnt;
int a0,a1,b0,b1;
int x;
int ans=0;
void dfs(int p)
{
	for(int i=0;i<=pf[p][1];i++)
	{
		for(int j=1;j<=i;j++)
			x*=pf[p][0];
		if(p==cnt)
		{
			if(gcd(x,a0)==a1)
				ans++;			
		}
		else
			dfs(p+1);
		for(int j=1;j<=i;j++)
			x/=pf[p][0];
	}
}
int main()
{
	int t;
	scanf("%d",&t);
	int c=0;
	for(int i=2;i<=50000;i++)
	{
		if(!isprime[i])
			prime[++c]=i;
		for(int j=1;prime[j]*i<=50000;j++)
		{
			isprime[prime[j]*i]=1;
			//cout<<prime[j]<<" ";
			if(i%prime[j]==0)
				break;
		}
	}
	while(t--)
	{
		cnt0=0,cnt1=0;
		x=1,ans=0;
		cnt=0;
		memset(pf1,0,sizeof(pf1));
		memset(pf,0,sizeof(pf));
		memset(pf0,0,sizeof(pf0));
		scanf("%d%d%d%d",&a0,&a1,&b0,&b1);
		for(int i=1;i<=c&&prime[i]<=b0;i++)
		{
			if(b0%prime[i]==0)
				cnt0++,pf0[cnt0][0]=prime[i];
			while(b0%prime[i]==0)
				pf0[cnt0][1]++,b0/=prime[i];
		}
//		for(int i=1;i<=cnt0;i++)
//			cout<<pf0[i][0]<<" "<<pf0[i][1]<<endl;
		for(int i=1;i<=c&&prime[i]<=b1;i++)
		{
			if(b1%prime[i]==0)
				cnt1++,pf1[cnt1][0]=prime[i];
			while(b1%prime[i]==0)
				pf1[cnt1][1]++,b1/=prime[i];
		}
//		for(int i=1;i<=cnt1;i++)
//			cout<<pf1[i][0]<<" "<<pf1[i][1]<<endl;
		if(b0>1)
			pf0[++cnt0][0]=b0,pf0[cnt0][1]++;
		if(b1>1)
			pf1[++cnt1][0]=b1,pf1[cnt1][1]++;
		cnt=0;
		for(int i=1;i<=cnt1;i++)
		{
			for(int j=1;j<=cnt0;j++)
				if(pf0[j][0]==pf1[i][0])
				{
					if(pf0[j][1]==pf1[i][1])
						pf[++cnt][0]=pf0[j][0],pf[cnt][1]=pf0[j][1];
					else
						for(int k=1;k<=pf1[i][1];k++)
							x*=pf1[i][0];
				}
		}
//		cout<<x<<endl;
		dfs(1);
		printf("%d\n",ans);
	}
	return 0;
}