算法训练 Hankson的趣味题
时间限制:1.0s 内存限制:64.0MB
提交此题 锦囊1 锦囊2
问题描述
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 的个数;
样例输入
2
41 1 96 288
95 1 37 1776
样例输出
6
2
样例说明
第一组输入数据,x 可以是9、18、36、72、144、288,共有6 个。
第二组输入数据,x 可以是48、1776,共有2 个。
数据规模和约定
对于 50%的数据,保证有1≤a0,a1,b0,b1≤10000 且n≤100。
对于 100%的数据,保证有1≤a0,a1,b0,b1≤2,000,000,000 且n≤2000。
// 挺好理解的题目,唯一分解定理
// 两个数的最小公约数 是 两个数的对应素数 位的 最小的指数
// 两个数的最大公倍数 是 两个数的对应素数 位的 最大的指数
// 所以做法为先打表找出素数,然后对素数表逐个遍历,判断求解
// 但是有细节没有处理好,结果不知道 错了多少回
#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<queue>
#include<algorithm>
#define ll long long
#define maxn 100007
using namespace std;
int vis[maxn];
int prime[maxn];
int cnt = 0;
int num[4];
int a0,a1,b0,b1;
void find_prime()
{
vis[1] = 1;
for(int i=2;i<maxn;i++)
{
if(vis[i] == 0)
{
prime[cnt++] = i;
for(int j=i+i;j<maxn;j+=i)
vis[j] = 1;
}
}
}
ll solve()
{
ll ans = 1;
for(int i=0;i<cnt;i++)
{
int cnt1=0,cnt2=0;
int cnt3=0,cnt4=0;
if(b1 % prime[i] == 0)
{
while(a1 % prime[i] == 0) { a1 /= prime[i]; cnt1 ++; }
while(a0 % prime[i] == 0) { a0 /= prime[i]; cnt2 ++; }
while(b1 % prime[i] == 0) { b1 /= prime[i]; cnt3 ++; }
while(b0 % prime[i] == 0) { b0 /= prime[i]; cnt4 ++; }
if( cnt1 > cnt2 || cnt3 < cnt4)
return 0;
// 分类讨论,我的讨论有合并
if(cnt1 == cnt2)
{
if(cnt1 > cnt3 ) return 0;
if(cnt3 == cnt4)
ans *= ( cnt3 - cnt1 + 1);
}
else if(cnt1 < cnt2 )
{
if(cnt3 == cnt4 && cnt1 > cnt3) return 0;
if(cnt3 > cnt4 && cnt1 != cnt3) return 0;
}
}
if(i==cnt-1 && b1 > 1) // 没处理完,放入继续处理,可能放进去后,改变结果(使结果为0)
prime[cnt++] = b1;
}
return ans;
}
int main()
{
memset(vis,0,sizeof(vis));
find_prime();
int n;
scanf("%d",&n);
for(int i=0;i<n;i++)
{
for(int j=0;j<4;j++)
scanf("%d",&num[j]);
a0 = num[0];
a1 = num[1];
b0 = num[2];
b1 = num[3];
ll ans = solve();
printf("%lld\n",ans);
}
return 0;
}