POJ 3604 Professor Ben G++ 快筛素数 前N个自然数的立方和 背

#include <iostream>
#include <cstdio>
using namespace std;
//抄博友程序 看博友好分析 需复习  快筛素数   前N个自然数的立方和  背 
int da[5000008];
int prim[5000008];
int js;
int main()
{
	js=0;
	for(int i=2;i<=5000000;i++)
	{
		if(prim[i]==0)//抄博友程序  背 
		{
			da[js]=i;
			js++;
		} 
		for(int j=0;j<js && (da[j]*i)<=5000000;j++)//抄博友程序  背 
		{
			prim[da[j]*i]=1;
			if(i%da[j]==0)//抄博友程序 背 
			{
				break; 
			} 
		}	
	}
	int T;
	//cin>>T;
	scanf("%d",&T);
	for(int o=0;o<T;o++)
	{
		int n;
		//cin>>n;
		scanf("%d",&n);
		long long jg=1;
		for(int i=0;i<js && da[i]*da[i]<=n ;i++)//抄博友 && da[i]<=n超时 巧妙 
		{
			if(n%da[i]==0)
			{
				int cns=0;
				while(n%da[i]==0)
				{
					n=n/da[i];
					cns++;
				}
				jg=jg*(((cns+1)*(cns+2))/2)*(((cns+1)*(cns+2))/2);//前N个自然数的立方和
			}
		} 
		//cout<<jg<<endl;
		if(n>1)//n一般为1 n是素数时
		{
			jg=jg*9;//不能改为jg=9 
		}
		printf("%lld\n",jg);
	}
}