LightOJ-1341 Aladdin and the Flying Carpet 因子个数

先附上大佬博客Orz：http://www.cnblogs.com/qq2424260747/p/4924216.html

这道题要用到一个结论：

设 N=p1^a1*p2^a2*p3^a3*...*pn^an(其中p1,p2,...,pn为N的因子，a1,a2,... ,an分别为因子的指数)（唯一分解定理）；则 N的因子个数 M=(1 + a1)*(1 + a2)*(1 + a3)*...*(1 + an) 

然后M/=2，再减去不合条件的个数即可。

附上AC代码：

#include<cstdio>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
using namespace std;
#define ll long long

const ll MAX=1e6+5;
bool isp[MAX];
ll pri[MAX];
void init()
{
    memset(isp,true,sizeof(isp));
    memset(pri,0,sizeof(pri));
    isp[0]=isp[1]=false;
    for(ll i=2;i<MAX;i++)
        if(isp[i])
        {
            pri[++pri[0]]=i;
            for(ll j=2;i*j<MAX;j++)
                isp[i*j]=false;
        }
}
ll fac[MAX],tot;//质因子
ll num[MAX];//每个质因子的幂
void decom(ll a)
{
    tot=0;
    memset(num,0,sizeof(num));
    ll n=a;
    for(ll j=1;j<=pri[0];j++)
    {
        if(pri[j]*pri[j]>n)
            break;
        if(n%pri[j]==0)
        {
            fac[tot]=pri[j];
            while(n%pri[j]==0)
            {
                num[tot]++;
                n/=pri[j];
            }
            tot++;
        }
    }
    if(n>1)
    {
        num[tot]++;
        fac[tot++]=n;
    }
    /*for(int i=0;i<tot;i++)
        cout<<"i="<<i<<" fac="<<fac[i]<<" num="<<num[i]<<endl;*/
}
ll a,b;
int main()
{
    int t;
    scanf("%d",&t);
    init();
    for(int tt=1;tt<=t;tt++)
    {
        scanf("%lld%lld",&a,&b);
        if(a<b*b)
        {
            printf("Case %d: 0\n",tt);
            continue;
        }
        decom(a);
        ll sum=1;
        for(ll i=0;i<tot;i++)
            sum*=(1+num[i]);
        sum/=2;
        //cout<<"sum="<<sum<<endl;
        for(ll i=1;i<b;i++)//此处不能用fac判断 eg:36,4
            if(a%i==0)
                sum--;
        printf("Case %d: %lld\n",tt,sum);
    }
    return 0;
}