GCD and LCM HDU - 4497

http://acm.hdu.edu.cn/showproblem.php?pid=4497

将gcd与lcm素因子分解 如果gcd某个素因子的幂次pi大于lcm对应素因子的幂次qi 那就是凑不出 因为gcd肯定是lcm因子

如果pi<=qi xyz三个数中必有一个幂次为pi 必有一个为qi 因为gcd就是取得三者最小 lcm取最大 只剩一个数的幂次ti可变 当ti=pi或ti=qi时 可以产生3!/2!种可能 意思是把这三个幂次全排列一下 当ti!=pi且ti!=qi时 则产生3!种可能 幂次内加法原理 幂次之间乘法原理即可

 

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=1e2+10;

int main()
{
    ll a1[maxn],a2[maxn],p1[maxn],p2[maxn];
    ll gcd,lcm,ans;
    int t,n1,n2,i,j,flag;
    scanf("%d",&t);
    while(t--){
        scanf("%lld%lld",&gcd,&lcm);

        n1=0;
        for(i=2;i*i<=gcd;i++){
            if(gcd%i==0){
                n1++;
                a1[n1]=i,p1[n1]=0;
                while(gcd%i==0){
                    gcd/=i;
                    p1[n1]++;
                }
            }
        }
        if(gcd!=1){
            n1++;
            a1[n1]=gcd,p1[n1]=1;
        }

        n2=0;
        for(i=2;i*i<=lcm;i++){
            if(lcm%i==0){
                n2++;
                a2[n2]=i,p2[n2]=0;
                while(lcm%i==0){
                    lcm/=i;
                    p2[n2]++;
                }
            }
        }
        if(lcm!=1){
            n2++;
            a2[n2]=lcm,p2[n2]=1;
        }
        /*
        printf("***%d***\n",n1);
        for(i=1;i<=n1;i++){
            printf("%lld %lld\n",a1[i],p1[i]);
        }
        printf("***%d***\n",n2);
        for(i=1;i<=n2;i++){
            printf("%lld %lld\n",a2[i],p2[i]);
        }
        */
        i=1,j=1,flag=1,ans=1;
        while(i<=n1&&j<=n2){
            if(a1[i]<a2[j]){
                flag=0;
                break;
            }
            else if(a1[i]==a2[j]){
                if(p1[i]>p2[j]){
                    flag=0;
                    break;
                }
                else if(p1[i]<p2[j]){
                    ans*=6*(p2[j]-p1[i]);
                }
                i++,j++;
            }
            else{
                ans*=6*p2[j];
                j++;
            }
        }
        if(i<=n1) flag=0;
        while(j<=n2){
            ans*=6*p2[j];
            j++;
        }
        if(flag) printf("%lld\n",ans);
        else printf("0\n");
    }
    return 0;
}