hdu 4497 GCD and LCM

思路：易知L不能整除G时为0；

将L/G质因数分解，对于其中的因子p,个数为cnt，则至少有一个包含p^cnt,至少有一个数不包含p；

只有一个数包含p^cnt时，有C(3,1);

有2个数包含p^cnt时，有C(3,1);

有2个数包含p因子，其中一个是p^cnt,另外一个有cnt-1种，总共有(cnt-1)A(3,2).

所以总共有6*cnt。

代码如下：

#include<iostream>
#include<stdio.h>
#include<algorithm>
#include<iomanip>
#include<cmath>
#include<cstring>
#include<vector>
#define ll __int64
#define pi acos(-1.0)
#define MAX 50000
using namespace std;
int prime[MAX],cnt,e[MAX],lena;
bool f[MAX];
void init()
{
    int i,j;
    cnt=0;
    for(i=2;i<MAX;i++){
        if(f[i]==0) prime[cnt++]=i;
        for(j=0;j<cnt&&i*prime[j]<MAX;j++){
            f[i*prime[j]]=1;
            if(i%prime[j]==0) break;
        }
    }
}
int solve(int n)
{
    int i,j=0,ans=1;
    for(i=0;i<cnt&&prime[i]*prime[i]<=n;i++){
        if(n%prime[i]==0){
            e[j]=1;
            n/=prime[i];
            while(n%prime[i]==0){
                e[j]++;
                n/=prime[i];
            }
            j++;
        }
    }
    if(n>1){
        e[j++]=1;
    }
    for(i=0;i<j;i++){
        ans*=6*e[i];
    }
    return ans;
}
int main(){
    int i,t,n,a,b,j,ans;
    init();
    scanf("%d",&t);
    while(t--){
        scanf("%d%d",&a,&b);
        if(b%a) printf("0
");
        else  printf("%d
",solve(b/a));
    }
    return 0;
}