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
这个问题中我才用的是直接暴力法求解,找出求最大公因数时满足条件的x,再找出求最小公倍数时满足的x,最后找出两次查找中相同的x有多少个。
#include<iostream>
#include<cstdlib>
#include<vector>
using namespace std;
int main()
{
vector<int> unfun1(int a, int b, int max); //求x和a0的最大公因数等于a1的x
vector<int> unfun2(int c, int d, int max); //求x的b0的最小公倍数等于b1的x
int judge(vector<int> m, vector<int> n); //判断有多少x符合条件
int Fun_1(const int a1, const int b1); //求最大公因数
int Fun_2(int a1, int b1); //求最小公倍数
void swap(int& a, int& b); //a>b时,交换a,b
int Max(int n1, int n2, int n3, int n4); //找出输入的一组数据中的最大值
vector<int> mark_1;
vector<int> mark_2;
vector<int> n_count; //用来存放结果
int n_num = 0,m_max = 0;
int m_1, m_2, m_3, m_4;
cin >> n_num;
for (int i = 0; i < n_num; i++)
{
cin >> m_1 >> m_2 >> m_3 >> m_4;
m_max = Max(m_1, m_2, m_3, m_4);
mark_1 = unfun1(m_1, m_2, m_max);
mark_2 = unfun2(m_3, m_4, m_max);
n_count.push_back(judge(mark_1, mark_2));
mark_1.clear(); //清空mark_1
mark_2.clear(); //清空mark_2
}
for (unsigned int i = 0; i < n_num; i++) //输出结果
{
cout << n_count[i] << endl;
}
system("pause");
return 0;
}
void swap(int& a, int& b) /*找出参数中较大的数和较小的数*/
{
int t = 0;
if (a < b)
{
t = a;
a = b;
b = t;
}
}
int Fun_1(const int a1, const int b1) /*根据辗转相除法求出最大公因数*/
{
int x, y;
int t = 0;
x = a1;
y = b1;
swap(x, y);
while (y != 0)
{
t = x % y;
x = y;
y = t;
}
return x;
}
int Fun_2(int a1, int b1) /*求最小公倍数*/
{
int product = 0;
product = a1 * b1;
return product / Fun_1(a1, b1); /*最小公倍数数等于两数之积除以最大公因数*/
}
vector<int> unfun1(int a, int b, int max) //判断有哪些x是满足最大公因数规则的
{
vector<int> ma;
for (int i = 1; i <= max; i++)
{
if (b == Fun_1(i, a))
{
ma.push_back(i);
}
}
return ma;
}
vector<int> unfun2(int c, int d, int max) //判断有哪些x是满足最小公倍数规则的
{
vector<int> mb;
for (int i = 1; i <= max; i++)
{
if (d == Fun_2(i, c))
{
mb.push_back(i);
}
}
return mb;
}
int judge(vector<int> m, vector<int> n) //判断两次循环中的x有多少个是相等的
{
vector<int>::iterator l_1;
vector<int>::iterator l_2;
int count_1 = 0;
for (l_1 = m.begin(); l_1 != m.end(); l_1++)
{
for (l_2 = n.begin(); l_2 != n.end(); l_2++)
{
if ((*l_1) == (*l_2))
{
count_1++;
}
}
}
//cout << "总数是:" << count_1 << endl;
return count_1;
}
int Max(int n1, int n2, int n3, int n4) //找出输入的四个数中的最大值
{
n1 = n1 > n2 ? n1 : n2;
n1 = n1 > n3 ? n1 : n3;
n1 = n1 > n4 ? n1 : n4;
return n1;
}
测试输入:
问题得到解决