题目地址: 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;
}
}
}