求最大公因数最小公倍数的逆向问题

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;
}

测试输入：

问题得到解决