题目2 : 机会渺茫
时间限制:
5000ms
单点时限:
1000ms
内存限制:
256MB
-
3 2
样例输出
-
4 1
描述
小Hi最近在追求一名学数学的女生小Z。小Z其实是想拒绝他的,但是找不到好的说辞,于是提出了这样的要求:对于给定的两个正整数N和M,小Hi随机选取一个N的约数N',小Z随机选取一个M的约数M',如果N'和M'相等,她就答应小Hi。
小Z让小Hi去编写这个随机程序,到时候她review过没有问题了就可以抽签了。但是小Hi写着写着,却越来越觉得机会渺茫。那么问题来了,小Hi能够追到小Z的几率是多少呢?
输入
每个输入文件仅包含单组测试数据。
每组测试数据的第一行为两个正整数N和M,意义如前文所述。
对于40%的数据,满足1<=N,M<=106
对于100%的数据,满足1<=N,M<=1012
输出
对于每组测试数据,输出两个互质的正整数A和B(以A分之B表示小Hi能够追到小Z的几率)。
思路:
第一步分别求N,M的约数,第二步统计N,M约数中相同的个数,然后进行输出处理。
代码:
#include <iostream>
#include <cmath>
#include <set>
using namespace std;
/*
gcd算法:辗转相除法,
定理:gcd(a,b)=gcd(b,a mod b)
证明: a可以表示成a = kb+r,则r=a mod b
假设d是a,b的一个公约数,则 d整除a,d整除b,
而r = a - kb,所以d整除r,因此d是b,a mod b的公约数。
*/
int gcd(int a,int b) //求a,b的最大公约数
{
if(b == 0)
return a;
else
return gcd(b,a%b);
}
set<int> getdivisor(int a)
{
set<int> s; //使用set集合的目的是set的查找非常高效。
int len = static_cast<int>(sqrt(a*1.0));//提高查找效率
for(int i=1;i<=len;i++)
{
if(a%i==0)
{
s.insert(i);
s.insert(a/i);
}
}
return s;
}
int main()
{
int a,b;
cin >> a >> b;
set<int> s1,s2;
s1 = getdivisor(a);
s2 = getdivisor(b);
int cunt=0;
for(auto it=s1.begin();it!=s1.end();it++)//统计相同约数。
{
if(s2.find(*it)!=s2.end())
cunt++;
}
int len = s1.size()*s2.size();
int temp = gcd(len,cunt);//处理输出形式,将百分数转化为分数。
cout << len/temp << " " << cunt/temp << endl;
return 0;
}