蓝桥杯 算法训练 Hankson的趣味题 （唯一分解定理）

算法训练 Hankson的趣味题
时间限制：1.0s 内存限制：64.0MB
提交此题 锦囊1 锦囊2
问题描述
　　Hanks 博士是BT (Bio-Tech，生物技术) 领域的知名专家，他的儿子名叫Hankson。现 在，刚刚放学回家的Hankson 正在思考一个有趣的问题。 今天在课堂上，老师讲解了如何求两个正整数c1 和c2 的最大公约数和最小公倍数。现 在Hankson 认为自己已经熟练地掌握了这些知识，他开始思考一个“求公约数”和“求公 倍数”之类问题的“逆问题”，这个问题是这样的：已知正整数a0,a1,b0,b1，设某未知正整 数x 满足： 1． x 和a0 的最大公约数是a1； 2． x 和b0 的最小公倍数是b1。 Hankson 的“逆问题”就是求出满足条件的正整数x。但稍加思索之后，他发现这样的 x 并不唯一，甚至可能不存在。因此他转而开始考虑如何求解满足条件的x 的个数。请你帮 助他编程求解这个问题。
输入格式
　　输入第一行为一个正整数n，表示有n 组输入数据。

　　接下来的n 行每 行一组输入数据，为四个正整数a0，a1，b0，b1，每两个整数之间用一个空格隔开。输入 数据保证a0 能被a1 整除，b1 能被b0 整除。
输出格式
　　输出共n 行。每组输入数据的输出结果占一行，为一个整数。
　　对于每组数据：若不存在这样的 x，请输出0； 若存在这样的 x，请输出满足条件的x 的个数；
样例输入
2
41 1 96 288
95 1 37 1776
样例输出
6
2
样例说明
　　第一组输入数据，x 可以是9、18、36、72、144、288，共有6 个。
　　第二组输入数据，x 可以是48、1776，共有2 个。
数据规模和约定
　　对于 50%的数据，保证有1≤a0，a1，b0，b1≤10000 且n≤100。
　　对于 100%的数据，保证有1≤a0，a1，b0，b1≤2,000,000,000 且n≤2000。

// 挺好理解的题目，唯一分解定理 
// 两个数的最小公约数 是 两个数的对应素数 位的 最小的指数 
// 两个数的最大公倍数 是 两个数的对应素数 位的 最大的指数
// 所以做法为先打表找出素数，然后对素数表逐个遍历，判断求解

// 但是有细节没有处理好，结果不知道 错了多少回 
#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<queue>
#include<algorithm>

#define ll long long

#define maxn 100007 

using namespace std;

int vis[maxn];
int prime[maxn];

int cnt = 0;
int num[4];
int a0,a1,b0,b1;

void find_prime()
{
    vis[1] = 1;
    for(int i=2;i<maxn;i++)
    {
        if(vis[i] == 0)
        {
            prime[cnt++] = i;
            for(int j=i+i;j<maxn;j+=i)
                vis[j] = 1;
        }
    }
} 

ll solve()
{   
    ll ans = 1;

    for(int i=0;i<cnt;i++)
    {
        int cnt1=0,cnt2=0;
        int cnt3=0,cnt4=0;

        if(b1 % prime[i] == 0)
        {
            while(a1 % prime[i] == 0)   { a1 /= prime[i]; cnt1 ++; }

            while(a0 % prime[i] == 0)   { a0 /= prime[i]; cnt2 ++; }

            while(b1 % prime[i] == 0)   { b1 /= prime[i]; cnt3 ++; }

            while(b0 % prime[i] == 0)   { b0 /= prime[i]; cnt4 ++; }


            if( cnt1 > cnt2 || cnt3 < cnt4)
                return 0;   

            // 分类讨论，我的讨论有合并 
            if(cnt1 == cnt2)
            {
                if(cnt1 > cnt3 ) return 0;
                if(cnt3 == cnt4)
                    ans *= ( cnt3 - cnt1 + 1);  
            }
            else if(cnt1 < cnt2 )
            {
                if(cnt3 == cnt4 && cnt1 > cnt3) return 0;
                if(cnt3 > cnt4 && cnt1 != cnt3) return 0;
            }
        }

        if(i==cnt-1 && b1 > 1) // 没处理完，放入继续处理，可能放进去后，改变结果（使结果为0） 
            prime[cnt++] = b1; 
    }

    return ans;
}

int main()
{
    memset(vis,0,sizeof(vis));
    find_prime();

    int n;
    scanf("%d",&n);

    for(int i=0;i<n;i++)
    {
        for(int j=0;j<4;j++)
            scanf("%d",&num[j]);

        a0 = num[0];
        a1 = num[1];
        b0 = num[2];
        b1 = num[3];

        ll ans = solve();
        printf("%lld\n",ans);
    }
    return 0;
}