P1072 Hankson 的趣味题

题目链接https://www.luogu.com.cn/problem/P1072

题目描述

Hanks 博士是 BT（Bio-Tech，生物技术) 领域的知名专家，他的儿子名叫 Hankson。现在，刚刚放学回家的 Hankson 正在思考一个有趣的问题。

今天在课堂上，老师讲解了如何求两个正整数 c1c_1c1​ 和 c2c_2c2​ 的最大公约数和最小公倍数。现在 Hankson 认为自己已经熟练地掌握了这些知识，他开始思考一个“求公约数”和“求公倍数”之类问题的“逆问题”，这个问题是这样的：已知正整数a0,a1,b0,b1 a_0,a_1,b_0,b_1a0​,a1​,b0​,b1​，设某未知正整数 xxx 满足：

1． xxx 和 a0a_0a0​ 的最大公约数是 a1a_1a1​；

2． xxx 和 b0b_0b0​ 的最小公倍数是 b1b_1b1​。

Hankson 的“逆问题”就是求出满足条件的正整数 xxx。但稍加思索之后，他发现这样的 xxx 并不唯一，甚至可能不存在。因此他转而开始考虑如何求解满足条件的 xxx 的个数。请你帮助他编程求解这个问题。

输入格式

第一行为一个正整数 nnn，表示有 nnn 组输入数据。接下来的n nn 行每行一组输入数据，为四个正整数 a0,a1,b0,b1a_0,a_1,b_0,b_1a0​,a1​,b0​,b1​，每两个整数之间用一个空格隔开。输入数据保证 a0a_0a0​ 能被 a1a_1a1​ 整除，b1b_1b1​ 能被 b0b_0b0​ 整除。

输出格式

共 nnn 行。每组输入数据的输出结果占一行，为一个整数。

对于每组数据：若不存在这样的 xxx，请输出 000，若存在这样的 xxx，请输出满足条件的 xxx 的个数；

输入输出样例

输入 #1

2 
41 1 96 288 
95 1 37 1776 

输出 #1

6 
2

题意：题意很简单，就是x 和a0的最大公约数是a1， x和b0的最小公倍数是b1，求这样的x有多少个

思路：

引用大佬的推导过程，第一次了解到gcd里面的数字可以向里面做乘法或除法，例如第二行和第四五行的转化

知识点：

1.要做到整除就必须写上a%b==0的限制条件，尤其是在枚举的时候，否则纯粹的计算机上的除法做到的知识下取整，做不到整除所以在必要的时候必须要加上这一个限制条件。

2.第一次了解到gcd里面的数字可以向里面做乘法或除法，例如第二行和第四五行的转化

3.在枚举的时候，若有一个倍数上限，常采用乘积的方式来枚举，如，在很小的范围内枚举，得到一个数，再用除法得到与他配对的另一个数。

代码如下：

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
const int maxn=2010;
int a0,a1,b0,b1;
int p,q;
int n;
int gcd(int x,int y){
	if(!y) return x;
	else{
		return gcd(y,x%y);
	}
}
int lcm(int x,int y){
	return x*y/gcd(x,y);
}
int main(){
	cin>>n;
	while(n--){
	cin>>a0>>a1>>b0>>b1;
	p=a0/a1;
	q=b1/b0;
	int ans=0;
	for(int x=1;x*x<=b1;x++){
		if(b1%x==0){
		if(x%a1==0&&gcd(x/a1,p)==1&&gcd(b1/x,q)==1) ans++;
		int y=b1/x;
		if(y==x) continue;
		if(y%a1==0gcd(y/a1,p)==1&&gcd(b1/y,q)==1) ans++;
	}
}
    cout<<ans<<endl;
}
	return 0;
}