Pairs Forming LCM LightOJ - 1236（求n以内的最小公因数为n的数对的对数）

然后关键是如何讨论，主要是形成数学思维。
假设a,b的最小公倍数为n,有：
那么ak肯定是ek和dk中的最大值。
即对于两个数a和b，其中必有一个数是包含了某个素因子的最高次幂，而另一个数包含该素因子的幂的范围是0~ak。然后对n的所有素因子依次进行如此的讨论，那么就一共有(2* a1+1) * (2 * a2+1 ) … * (2*ak+1)个。
根据题目要求，a<=b，所以要除以2，但由于当a和b都含有每个素因子的最高次幂时，只有一种可能，而除2就把这种可能消除了，所以要加1，把这种可能加回来。
但仅仅这样写会t,由于n比较大，那么我们在素因子分解中，就不能一个数一个数的枚举，可以预先处理数1e7以内的素数，然后进行分解。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e5+5;
const int M=1e7+5;
ll n;
vector<int>prime;
bool vis[M];
int e[N],cnt;
void euler()
{
    memset(vis,0,sizeof(vis));
    prime.clear();
    for(int i=2;i<M;i++)
    {
        if(!vis[i])
            prime.push_back(i);
        for(int j=0;j<prime.size()&&i*prime[j]<M;j++)
        {
            vis[i*prime[j]]=1;
            if(i%prime[j]==0)
                break;
        }
    }
}
void divide()
{
    ll tn=n;
    cnt=0;
    memset(e,0,sizeof(e));
    for(int i=0;prime[i]*prime[i]<=n&&i<prime.size();i++)
    {
        if(tn%prime[i]==0)
        {
            cnt++;
            while(tn%prime[i]==0)
            {
                tn/=prime[i];
                e[cnt]++;
            }
        }
    }
    if(tn>1)
    {
        cnt++;
        e[cnt]++;
    }
}
int main()
{
    int t,cas=0;
    euler();
    scanf("%d",&t);
    while(t--)
    {
        scanf("%lld",&n);
        divide();
        ll ans=1;
        for(int i=1;i<=cnt;i++)
            ans*=(2*e[i]+1);
        ans=ans/2+1;
        printf("Case %d: %lld\n",++cas,ans);
    }
    return 0;
}
/*
15
2
3
4
6
8
10
12
15
18
20
21
24
25
27
29
*/