HDU 4497 GCD and LCM （分解质因数）

标签：

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

如果G不是L的约数 就不可能找到三个数。

L的所有素因子一定包含G的所有素因子 并且次方数一定大于等于G的。只需要三个数 对于每一个素因子的次方数 三个的最小值是G的，最大值是L的。考虑三个对应的次方数都不一样，那么其中两个是确定的 一个是G的一个是L的 剩下的一个在G和L的之间。算上排列 总共有6种。或者其中两个是一样的，那么也有6种情况。最后可以合并计算。

//#pragma comment(linker, "/STACK:10240000,10240000")
#include <algorithm>
#include <iostream>
#include <sstream>
#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <vector>
#include <cmath>
#include <queue>
#include <stack>
#include <set>
#include <map>
#define mod 4294967296
#define MAX 0x3f3f3f3f
#define lson o<<1, l, m
#define rson o<<1|1, m+1, r
#define SZ(x) ((int)ans.size())
#define MAKE make_pair
#define INFL 0x3f3f3f3f3f3f3f3fLL
#define mem(a) memset(a, 0, sizeof(a))
const double pi = acos(-1.0);
const double eps = 1e-9;
const int N = 200005;
const int M = 20005;
typedef __int64 ll;
using namespace std;

ll a, b;
struct C {
    ll num, cnt;
} s[N], t[N];
int T;

int Ini(ll a, C* f) {
    int tmp = sqrt(1.0*a + 0.5), e = 0;
    for(int i = 2; i <= tmp; i++) {
        if(a % i == 0) {
            int k = 0;
            while(a % i == 0) {
                k++;
                a /= i;
            }
            f[e].num = i;
            f[e++].cnt = k;
        }
    }
    if(a != 1) {
        f[e].cnt = 1;
        f[e++].num = a;
    }
    return e;

}

int main()  {

    //freopen("in.txt","r",stdin);
    cin >> T;
    while(T--) {
        cin >> a >> b;
        if(b % a) {
            puts("0");
            continue;
        }
        mem(s);
        mem(t);
        int n = Ini(a, s);
        int m = Ini(b, t);

        ll ans = 1, x;
        int fr = 0;
        for(int i = 0; i < m; i++) {
            if(t[i].num == s[fr].num) {
                x = t[i].cnt - s[fr].cnt;
                fr++;
                if(x) ans *= x * 6;
            } else ans *= t[i].cnt * 6;

        }
        cout << ans << endl;
    }
    return 0;
}

HDU 4497 GCD and LCM （分解质因数）

标签：

原文地址：http://blog.csdn.net/u013923947/article/details/46343869