hihocoder offer收割赛。。#1284

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

我想着蛋疼的做了质因数分解，然后再算的因子个数。。慢的一比，结果导致超时，还不如直接一个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;
}