[offer收割]编程练习赛3-题目2：机会渺茫

题目2 : 机会渺茫

时间限制: 5000ms

单点时限: 1000ms

内存限制: 256MB

描述

小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的几率）。

样例输入

3 2

样例输出

4 1

思路：

第一步分别求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;
}