基础数学1007 LightOJ 1236

题意:
给一个n,问有多少对lcm(i,j)=n|j>=i
思路:
因为要lcm==n,所以将n分解成质因子后
n=p1^k1 * p2^k2 * … pk^ek
(i,j)之间至少有一个数分解后也有pi^ei,剩下一个数可以是0~ei
所以对于每一个p,都有2*ei+1种情况
然后因为除了(n,n)之外都重复了一遍,所以Ans++,Ans/=2

#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
#include<math.h>
#include<queue>
#include<stack>
#include<string>
#include<vector>
#include<map>
#include<set>
using namespace std;
#define lowbit(x) (x&(-x))
typedef long long LL;
const int maxn = 100005;
const int inf=(1<<28)-1;
#define maxp 10000005
bool notprime[maxp];
int primes[800005];
void get_prime()
{
    notprime[1]=true;
    for(int i=2;i<maxp;++i)
    if(!notprime[i])
    {
        primes[++primes[0]]=i;
        for(LL j=(LL)i*i;j<maxp;j+=i)
        notprime[j]=true;
    }
}
int tot;
LL Fac[maxn];
int Num[maxn];
void GetFac(LL x)
{
    tot=0;
    memset(Num,0,sizeof(Num));
    LL k=sqrt(x);
    for(int i=1;i<=primes[0]&&primes[i]*primes[i]<=x;++i)
    if(x%primes[i]==0)
    {
        Fac[++tot]=primes[i];
        x/=primes[i];
        Num[tot]=1;
        while(x%primes[i]==0) 
        {
            x/=primes[i];
            Num[tot]++;
        }
    }
    if(x>1) 
    {
        Fac[++tot]=x;
        Num[tot]=1;
    }
}
int main()
{
    get_prime();
    //printf("%d\n",primes[0]);
    int Case=0,T;
    scanf("%d",&T);
    while(T--)
    {
        LL n;
        scanf("%lld",&n);
        LL Ans=1;
        GetFac(n);
        //for(int i=1;i<=tot;++i) printf("%lld(%d)  ",Fac[i],Num[i]);
        for(int i=1;i<=tot;++i)
        Ans*=Num[i]*2+1;
        Ans++;Ans/=2;
        printf("Case %d: %lld\n",++Case,Ans);
    }
    return 0;
}