[唯一分解定理] LightOJ 1236

题意：

求有多少组 （ i，j ）
使 lcm(i, j) = n and (i ≤ j).
(1 ≤ n ≤ 10^14)

思路

假设已找到一对（a，b）的lcm=n

有a=p1^d1+p2^d2+…pm^dm

b=p1^e1+p2^e2+…pm^em

易知max(di,ei)=ci

先考虑有序数对（a，b），由唯一分解定理知，a的每一个素因数的幂的大小都决定一个独一无二的数。

所以（a，b）的种数就是(di,ei)的种数，即2*(ci+1)-1(因为有序对(c1,c1)重复了一次所以-1)

所以有序对(a,b)的种数ans=(2*c1+1)(2*c2+1)(2*c3+1)…(2*cm+1)

但是要求求无序对的种数，已知(a,b)(b,a)重复，除了（n，n）只算了一个之外，所以ans = ans/2 +1

代码

#include <algorithm>
#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>

using namespace std;

typedef long long ll;

const int N = 1e7 + 100;

ll prime[ N / 10 ];
bool vis[ N ];
int tot;

void init () {

    memset ( vis, false, sizeof ( vis ) );

    for ( ll i = 2; i < N; ++i )
        if ( !vis[ i ] ) {
            prime[ tot++ ] = i;

            for ( ll j = i * i; j < N; j += i )
                vis[ j ] = true;
        }
}

int main () {
    init ();
    int T;
    scanf ( "%d", &T );

    for ( int k = 1; k <= T; ++k ) {
        ll n;
        scanf ( "%lld", &n );

        ll t = n; //因子
        ll ans = 1;
        //分解
        for ( int i = 0; i < tot && prime[ i ] * prime[ i ] <= n; ++i ) {
            ll cnt = 0;
            while ( t % prime[ i ] == 0 )
                ++cnt, t /= prime[ i ];

            ans *= ( 2 * cnt + 1 );
        }

        if ( t > 1 )
            ans *= 3;
    }
}