好久没刷题，水一水，反正排不上名次。。这道题记下

我想着蛋疼的做了质因数分解，然后再算的因子个数。。慢的一比，结果导致超时，还不如直接一个for循环搞定。。也是醉了

最后其实就是算两个数的gcd，然后gcd的因子个数除两个数的因子个数乘积就是最后结果，再约分一下就行。

#1284 : 机会渺茫

时间限制:
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

#include <iostream>
#include <vector>
#include <algorithm>
#include<map>
#define N 100000
using namespace std;
typedef long long ll;
int plist[100000], pcount = 0;
int prime(int n) {
    int i;
    if ((n != 2 && !(n % 2)) || (n != 3 && !(n % 3)) || (n != 5 && !(n % 5)) || (n != 7 && !(n % 7)))
        return 0;</span><span style="color:#333333;">
    for (i = 0; plist[i] * plist[i] <= n; i++)
        if (!(n%plist[i]))
            return 0;
    return n>1;
}
void initprime() {
    int i;
    for (plist[pcount++] = 2, i = 3; i<100000; i++)
        if (prime(i))
            plist[pcount++] = i;
}
void factor(int a, map<int,int> &s) {//质因数分解
    int i = 2;
    while(a>1) {
        if (prime(i) && a%i == 0) {
            s[i] = 0;
            while (a%i == 0) {
                ++s[i];
                a /= i;
            }
        }
        i++;
    }
}
ll cal(ll x) {
    ll i, ans = 0;
    for (i = 1; i * i <= x; i++) {
        if (x % i == 0) ans += 2;
    }
    if ((i - 1) * (i - 1) == x) ans -= 1;
    return ans;
}
int min(int a, int b) {
    return a < b ? a : b;
}
ll gcd(ll a, ll b)
{
    ll t = a%b;
    while (t) {
        a = b;
        b = t;
        t = a%b;
    }
    return b;
}
int main()
{
    ll a, b;
    cin >> a >> b;
    ll s1,s2,s3,tmp;
    tmp = gcd(a,b);
    s1 = cal(a);
    s2 = cal(b);
    s3 = cal(tmp);
    //几率为(s3)/(s1*s2)
    ll m1,m2,m3;
    m1 = s3;
    m2 = s2*s1;
    m3 = gcd(m1,m2);//约分
    m1 = m1/m3;
    m2 = m2/m3;

    cout << m2 <<" "<<m1<<endl;
    return 0;
}