唯一分解定理·UVA 10892·LCM的个数

题目大意：

给定整数n(n ≤2000000000),求有多少对(a, b) （a≤b）满足LCM(a, b) = n;

解题思路：

这种题不用唯一分解定理也不知道该从什么地方下手，我们应该已经知道根据唯一分解定理：

a = p1^k1*p2^k2*p3^k3.......；

b = q1^l1*q2^l2*q3^l3.......;

如果LCM(a, b )  = n的话，那么一定能够有：

n = m1^r1^m2^r2*m3^r3......；

其中 r1 = max(ki, li);

因此我们根据这个可以先行筛出素数并枚举素数，我们知道一种素数个数如果为s，那么a，b中这一素数的个数有如下三种情况：

1、a中个数为s，b中为0~s-1；

2、b中个数为s，a中为0~s-1；

3、a = s， b = s；

所以这一素数所组成的种类是2*s+1；

那么所有素数所组成的可能性，根据乘法原理，种类相乘即可；

但是有一个问题，因为由于时间的问题，我们只能枚举出上线为1e5的素数表，但是如果n本身有一个大于1e5的素数怎么办，也就是说枚举完素数表后，n不为1，也就是说还有三种情况（根据如上推断），那么将结果*3即可；

最后的结果因为要求a<b，所以除去赘余，ans/2+1即可；

AC代码：

#include <map>
#include <set>
#include <stack>
#include <cmath>
#include <queue>
#include <bitset>
#include <string>
#include <vector>
#include <cstdio>
#include <cctype>
#include <fstream>
#include <cstdlib>
#include <sstream>
#include <cstring>
#include <iostream>
#include <algorithm>
#pragma comment(linker, "/STACK:1024000000,1024000000")

using namespace std;
#define   maxn          1e5
#define   lson          l,m,rt<<1
#define   rson          m+1,r,rt<<1|1
#define   ms(x,y)      memset(x,y,sizeof(x))
#define   mc(x, y)     memcpy(x, y, sizeof(x))
#define   rep(i,n)      for(int i=0;i<(n);i++)
#define   repf(i,a,b)   for(int i=(a);i<=(b);i++)
#define   pii           pair<int,int>
#define   mp            make_pair
#define   FI            first
#define   SE            second
#define   IT            iterator
#define   PB            push_back
#define   Times         10
typedef   long long     ll;
typedef   unsigned long long ull;
typedef   long double   ld;

const double eps = 1e-10;
const double  pi = acos(-1.0);
const  ll    mod = 1e9+7;
const  int   inf = 0x3f3f3f3f;
const  ll    INF = (ll)1e18+300;
const int    maxd = 1e5;

int m, n;
int pri[maxd];
int vis[maxd];
int num;
void init() {
    ms(pri, 0);
    ms(vis, 0);
    num = 0;
    for (int i = 2; i < maxd; i++) {
        if(!vis[i]) {
            pri[num ++] = i;
            for (int j = 2*i; j <= maxd; j += i){
                pri[j] = 1;
            }
        }
    }
}

void solve (){

}

int main(){
    init();
    while (cin >> n && n) {
       long long ans = 1LL;
       m = n;
      for (int i = 0; n > 1 && i < num; ++ i) {
            int size = 0;
            while (n%pri[i] == 0) {
                n /= pri[i];
                size ++;
            }
            ans *= (1LL + (size<<1LL));
     }
      if(n > 1) {
          ans = 3LL * ans;
      }
     ans = ans/2LL + 1LL;
      printf("%d %lld\n",m,ans);
    }
    return 0;
}