【hihoCode】-1284-机会渺茫（STL-set）

#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

题解：找到 n （保证n<m）的所有约数存入set。再看m的约数在set里可不可以找得到，找得到就不是s.end(),否则就是s.end().

这样就可以计算出 n 的约数个数num1，m 约数个数num2，和约数相同的个数ant， 计算即可

#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
#include<cstring> 
#include<set>
#define LL long long
LL GCD(LL n,LL m)
{
	if(n%m==0)
		return m;
	return GCD(m,n%m);
}
int main()
{
	LL n,m;
	LL num1=0,num2=0,ant=0;
	set<LL> s;
	scanf("%lld%lld",&n,&m);
	//if(n==m)
	//	printf("1 1\n");
	//else
	//{
		if(n>m)
		swap(n,m);
		LL t=sqrt(n);
		for(int i=1;i<=t;i++)
		{
			if(n%i==0)
			{
				num1+=2;
				s.insert(i);
				s.insert(n/i);
			}
		}
		if(t*t==n)
			num1--;
		t=sqrt(m);
		for(int i=1;i<=t;i++)
		{
			if(m%i==0)
			{
				num2+=2;
				if(s.find(i)!=s.end())
					ant++;
				if(s.find(m/i)!=s.end())
					ant++;	
			}
		}
		if(t*t==m&&s.find(t)!=s.end())
			ant--,num2--;
		LL num=num1*num2;
		LL g=GCD(ant,num);
		printf("%lld %lld\n",num/g,ant/g);
	//}	
	return 0;
}