题目链接:点击打开链接
#1284 : 机会渺茫
时间限制:
5000ms
单点时限:
1000ms
内存限制:
256MB
-
3 2
样例输出
-
4 1
描述
小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的几率)。
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<map>
#define LL long long
using namespace std;
map<LL,int> M;
LL n,m,cnt;
LL solve(LL x)
{
LL i,ans=0;
for(i=1;i*i<=x;i++)
{
if(x%i==0)
{
ans+=2;
if(M[i]!=1)
M[i]=1;
else
cnt++;
if(M[x/i]!=1)
M[x/i]=1;
else
cnt++;
}
}
if(i*i==x) ans--;
return ans;
}
LL GCD(LL x,LL y)
{
if(x%y==0) return y;
return GCD(y,x%y);
}
int main()
{
scanf("%lld%lld",&n,&m);
cnt=0;
LL sum=solve(n)*solve(m);
LL a=GCD(sum,cnt);
printf("%lld %lld\n",sum/a,cnt/a);
return 0;
}
/*********************** set实现 **********************/
<pre name="code" class="cpp">#include<cstdio>
#include<algorithm>
#include<cstring>
#include<set>
#define LL long long
using namespace std;
LL n,m,cnt;
set<LL> S;
LL solve(LL x)
{
LL i,ans=0;
for(i=1;i*i<=x;i++)
{
if(x%i==0)
{
ans+=2;
if(S.find(i)!=S.end())
cnt++;
else
S.insert(i);
if(S.find(x/i)!=S.end())
cnt++;
else
S.insert(x/i);
}
}
if(i*i==x) ans--;
return ans;
}
LL GCD(LL x,LL y)
{
if(x%y==0) return y;
return GCD(y,x%y);
}
int main()
{
scanf("%lld%lld",&n,&m);
cnt=0;
LL sum=solve(n)*solve(m);
LL a=GCD(sum,cnt);
printf("%lld %lld\n",sum/a,cnt/a);
return 0;
}