您的位置:首页 > 其它

2013吉林通化邀请赛 1005 GCD and LCM

2013-08-24 21:16 281 查看
题目地址: http://acm.hdu.edu.cn/showproblem.php?pid=4497

本题就是整数的标准分解,然后结合一点组合计数, 我们知道每个素数对应指数的最大值和最小值,现在问原始的三个数有多少种。

首先做一个优化,将G变为1 ,将L 变为L/G (不整除就不谈了啊~ 无解) 、

然后对每一个素数,分三种情况:【1】 指数==0 (不整除,忽略) 【2】 指数==1 那么c[3][2]+ c[3][2](哪两个是最大的或最小的) 【3】指数>=2 情形2的6种+ 指定大的*指定小的*剩下一个取中间的(n-2)个。

分解10^9这样的数时,我们打了10^ 6以内的素数表后,用每一个素数去尽量除干净,那么如果剩下还有数,必定为一个素数的1次幂 (不可以容忍多于两个素数或者一个素数的高次幂),这样就是情形2 ,最后再乘6就可以了、

代码:

#include<iostream>
#include<algorithm>
#include<cmath>
#include<vector>
using namespace std;

#define N 10000000

bool  p[N+1]={0};

vector<int>  v;

void pre()
{
for(int i=2;i<=sqrt(N);i++)
if(p[i]==0)
for(int j=i*i;j<=N;j+=i)
p[j]=1;

for(int i=2;i<=N;i++)
if(p[i]==0) v.push_back(i);

}
int main()
{
pre();

//  for(int i=0;i<100;i++)
//   cout<<v[i]<<endl;
int T;
cin>>T;
int L,G;

while(T--)
{
cin>>G>>L;
if(L%G!=0)
{
cout<<"0"<<endl;
continue;

}

else
{
long long ans=1;
L=L/G;
G=1;

for(int i=0;i<v.size();i++)
{
if(L%v[i]!=0)  continue;

int exp=0;

while(L%v[i]==0)
{
exp++;
L/=v[i];
}

if(exp==1) ans*=6;
else if(exp>=2) ans*=(6*exp);
}

if(L!=1)   ans*=6;

cout<<ans<<endl;
}
}
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: