hihoCoder

来源:互联网 发布:无线传感器网络 编辑:程序博客网 时间:2024/06/16 15:38


时间限制: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 <cstdio>#include <cstring>#define LL long long// 求n和m的最大公约数LL GCD( LL n, LL m ) {    if( m == 0 )      return n;    return GCD( m, n % m );}// 求n和m公约数的个数LL num( LL s ) {    LL sum = 0;    for( LL i = 1; i * i <= s; i++ ) {        if( s % i == 0 )          sum++;        else continue;        if( s / i != i )          sum++;    }    return sum;}int main() {    LL n, m;    while( scanf( "%lld%lld", &n, &m ) != EOF ) {        LL a, b, c, p;        p = GCD( n, m );  // n和m的最大公约数        a = num( n );     // n的公约数个数        b = num( m );     // m的公约数个数        c = num( p );     // n和m公约数的个数,其实就是n和m最大公约数的约数个数,是分子        LL cnt = a * b;   // 分母        // 找到分子和分母的最大公约数,进行约分        LL x = GCD( cnt, c );        printf( "%lld %lld\n", cnt / x, c / x );    }    return 0;}


0 0
原创粉丝点击