poj 3604 Professor Ben

质因数分解；牛人推导公式(1^3+2^3+……+(1+a1)^3)*……*(1^3+2^3+……+(1+ai)^3)……

链接http://poj.org/problem?id=3604

#include<iostream>
#include<stdio.h>
#include<algorithm>
#include<iomanip>
#include<cmath>
#include<string>
using namespace std;
int prime[3001],m;
bool f[3001];
void init()
{
    int i,j;
    m=0;
    memset(f,0,sizeof(f));
    for(i=2;i<3000;i++)
    {
        if(f[i]==0)
        {
            prime[m++]=i;
            for(j=i*i;j<3000;j+=i)
                f[i]=1;
        }
    }
}
int main()
{
    init();
    __int64 sum,temp;
    int i,j,a,n,t;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d",&n);
        sum=1;
        for(i=0;i<m&&prime[i]*prime[i]<=n;i++)
        {
            if(n%prime[i]==0)
            {
                j=0;
                while(n%prime[i]==0)
                {
                    n/=prime[i];
                    j++;
                }
                temp=(j+1)*(j+2)/2;
                sum*=temp*temp;
            }
        }
        if(n>1)
            sum*=9;
        printf("%I64d
",sum);
    }
    return 0;
}