HDU - 4497 - GCD and LCM （线性筛素数 + 计数）

题目链接：

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

题目大意：

给出一个GCD和一个LCM，求出所有满足条件的三个数的组合。

题解：

一个大神写的分析，贼牛逼

如果G%L != 0，说明一定无解。

把K = G / L质数分解，G / L = p1^t1 * p2^t2 * p3^t3 * ……；同时 x/= L, y/= L, z/=L，不影响结果。

假设三个数字的质数分解是：

x = p1^i1 * p2^i2 * p3^i3 * ……

y = p1^j1 * p2^j2 * p3^j3 * ……

z = p1^k1 * p2^k2 * p3^k3 * ……

因为lcm(x, y, z) = K, 那么对于p1来说，i1, j1, k1里面一定有一个是0，并且一定有一个是t1，所以有3种情况：

0 0 t1 有3种

t1 t1 0 有3种

t1 0 1~t1-1 有(t1-1)*6种

一共是6*t1种。
所以总数为 6 * t1 * 6 * t2 * 6 * t3

AC代码

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
using namespace std;

const int N = 100000;
bool a[N];
int primes[N];
int num;

void setPrime()
{
    memset(a, 0, sizeof(a));
    int  i, j;
    for(i = 2; i < N; i++)
    {
        if(!(a[i])) primes[num++] = i;
        for(j = 0; (j<num && i*primes[j]<N); j++)
        {
            a[i*primes[j]] = 1;
            if(!(i%primes[j])) break;
        }
    }
}
int main()
{
    int T,a,b,ans,t;
    setPrime();
    cin >> T;
    while(T--)
    {
        scanf("%d%d", &a, &b);
        if(b % a != 0)
            cout << 0 <<endl;
        else
        {
            int k = b/a;
            ans = 1, t = 0;
            cout << primes[num - 1] <<endl;
            for(int i = 0; i < num; i++)
            {
                if(k < primes[i])break;
                while(k % primes[i] == 0)
                {
                    k /= primes[i];
                    t++;
                    cout << k <<endl;
                }
                if(t != 0)
                    ans = ans * 6 * t;
                t = 0;
            }
            if(k > 1)//可能会出现比10万大的质数，但是比10万大的质数平方后会超int，不合题意
                ans = ans * 6;
            cout << ans <<endl;
        }
    }
    return 0;
}