题目描述
Hanks 博士是 BT(Bio-Tech，生物技术) 领域的知名专家，他的儿子名叫 Hankson。现在，刚刚放学回家的 Hankson 正在思考一个有趣的问题。
今天在课堂上，老师讲解了如何求两个正整数c_1c1 和 c_2c2 的最大公约数和最小公倍数。现在 Hankson 认为自己已经熟练地掌握了这些知识，他开始思考一个“求公约数”和“求公倍数”之类问题的“逆问题”，这个问题是这样的：已知正整数a_0,a_1,b_0,b_1a0,a1,b0,b1，设某未知正整数xx 满足：
1． x 和 a_0a0 的最大公约数是 a1；
2． x和 b0 的最小公倍数是b1。
Hankson 的“逆问题”就是求出满足条件的正整数x。但稍加思索之后，他发现这样的xx 并不唯一，甚至可能不存在。因此他转而开始考虑如何求解满足条件的 xx 的个数。请你帮助他编程求解这个问题。
输入输出格式
输入格式：
第一行为一个正整数 nn，表示有 n 组输入数据。接下来的n 行每行一组输入数据，为四个正整数 a0,a1,b0,b1，每两个整数之间用一个空格隔开。输入数据保证 a0 能被 a1 整除，b1 能被b_0b0整除。
输出格式：
共 nn行。每组输入数据的输出结果占一行，为一个整数。
对于每组数据：若不存在这样的 xx，请输出 00；
若存在这样的xx，请输出满足条件的xx 的个数；
#include<bits/stdc++.h> using namespace std; int a,b,c,d,n,cnt,ans; int prime[1000000],f[1000000]; void work(int p) { int na=0,nb=0,nc=0,nd=0; while(a%p==0){a/=p;na++;} while(b%p==0){b/=p;nb++;} while(c%p==0){c/=p;nc++;} while(d%p==0){d/=p;nd++;} if(na<nb||nc>nd||nb>nd) ans=0; if(na>nb&&nc<nd&&nd>nb) ans=0; if((na>nb||nc<nd)&&nd>=nb) ans*=1; if(nc==nd&&na==nb) ans*=(nd-nb+1); } int main() { cin>>n; for(int i=2;i<=45000;i++){if(!f[i])prime[++cnt]=i;for(int j=1;j*i<=45000;j++)f[j*i]=1;} for(int k=1;k<=n;k++) { ans=1; cin>>a>>b>>c>>d; for(int i=1;i<=cnt&&ans;i++)if(d%prime[i]==0)work(prime[i]); if(d>1)work(d); cout<<ans<<endl; } }