机会渺茫 HihoCoder - 1284

 HihoCoder - 1284 

#1284 : 机会渺茫

Time Limit: 5000ms

Case Time Limit: 1000ms

Memory Limit: 256MB

Description

小Hi最近在追求一名学数学的女生小Z。小Z其实是想拒绝他的，但是找不到好的说辞，于是提出了这样的要求：对于给定的两个正整数N和M，小Hi随机选取一个N的约数N'，小Z随机选取一个M的约数M'，如果N'和M'相等，她就答应小Hi。

小Z让小Hi去编写这个随机程序，到时候她review过没有问题了就可以抽签了。但是小Hi写着写着，却越来越觉得机会渺茫。那么问题来了，小Hi能够追到小Z的几率是多少呢？

Input

每个输入文件仅包含单组测试数据。

每组测试数据的第一行为两个正整数N和M，意义如前文所述。

对于40%的数据，满足1<=N,M<=106

对于100%的数据，满足1<=N,M<=1012

Output

对于每组测试数据，输出两个互质的正整数A和B（以A分之B表示小Hi能够追到小Z的几率）。

Sample Input

3 2

Sample Output

4 1

   首先想到利用约数原理求出n,m的约数的个数，然后枚举求得n,m相同约数的个数。。。。。

这样显得比较麻烦，这里我们用到gcd...

令 d = gcd(a, b); d的约数就是n,m共有约数的个数。。。

代码：

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int maxn = 1e6+7;

int prim[maxn], cnt = 0;
bool vis[maxn];


ll gcd(ll a, ll b)
{
    if(!b) return a;
    else return gcd(b, a%b);
}

void getprim()
{
    memset(vis, 0, sizeof(vis));
    for(int i = 2; i*i < maxn; i++)
    {
        if(!vis[i]) {
            for(int j = i*i; j < maxn; j += i)
                vis[j] = true;
        }
    }
    for(int i = 2; i < maxn; i++)
        if(!vis[i]) prim[cnt++] = i;
}

ll gety(ll a){
    ll ans = 1;
    for(int i = 0; i < cnt&&a > 1; i++)
    {
        if(a%prim[i] == 0)
        {
            int c = 0;
            while(a%prim[i] == 0)
            {
                a /= prim[i];
                c++;
            }
            ans *= (c+1);
        }
    }
    if(a > 1) ans *= 2;
    return ans;
}

int main()
{
    ll n, m;
    getprim();
    scanf("%lld%lld", &n, &m);
    if(n < m) swap(n, m);
    ll d = gcd(n, m);
    ll a = gety(n)*gety(m), b = gety(d);
    ll k = gcd(a,b);
    printf("%lld %lld\n", a/k, b/k);
}