poj 3604 Professor Ben

#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\n",sum);
}
return 0;
}