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;
}