hdu1299

///hdu1299   Diophantus of Alexandria

1/x+1/y=1/n，给定一个n，满足条件的x，y有多少对。。

x 、y、n都是正整数，并且 显然，x >= n ， y >= n ，现在假设 y = n +k (k为正整数) ，那么带入公式，可以得出 x = (n*(n+k))/k = n*n/k + n; 由于x 是正整数，现在的关键问题就是要求出 n*n/ k 有多少组正整数的可能，显然，所要求的就是 n*n 因子的个数// 问题已经非常接近答案了，但是最后还有一个问题，n<= 10^9 ， 那么n*n <= 10^18 ,对于一个这么大的数字怎样才能求出它因子的个数呢？

命题1： 一个正整数 n 可以用素因子唯一表示为 p1^r1 * p2^r2 * ... pk^rk (其中 pi 为素数) , 那么这个数的因子的个数就是，(r1+1)*(r2+1)*...*(rk+1).

如果一个数字 n = p1^r1 * p2^r2 * ... pk^rk ，那么 n*n = p1^r1 * p2^r2 * ... pk^rk   * p1^r1 * p2^r2 * ... pk^rk ，它的因子的个数就是 (2*r1+1)*(2*r2+1)*...*(2*rk+1).

这个问题就转化成了求 n <= 10^9 的素因子的问题，只需要先筛选出 sqrt(10^9) 内的素数，然后用n去试除 sqrt(n) 中的每一个即可，当然，n可能会有大于sqrt(n) 的因子呶//. 只试到 sqrt(n)，怎么找出所有的素因子呢，不妨想一想，如果n有大于sqrt(n) 的素因子，它会有几个，又是多少呢？ //（答案。此时除过剩下的n自己就是素数。。T-T

code：

#include<iostream>
#include<cmath>
using namespace std;
//const int MAX=3000000; //Çó[2,MAX]¼äµÄËØÊý
const int MAX=50000;
bool isprime[MAX+1];
int prime[MAX];
int pnum; //ËØÊý±íÔªËØ×ÜÊý
void getprime()
{  
   int i,j;
   memset(isprime,0,sizeof(isprime));
   pnum=0;
   for (i=2;i<=MAX;i++)
   {
       if(!isprime[i])    prime[pnum++]=i;
       for (j=0; j<pnum &&prime[j]*i<=MAX ; j++)
       {
             isprime[prime[j]*i] = 1;
             if (i%prime[j]==0) break;
       }
   }
}
int main()
{
    int T,n,ca=0,i,ans,up;
    getprime();
    cin>>T;
    while(T--)
    {
              cin>>n;
              printf("Scenario #%d:/n",++ca);
              ans=1;
              for(i=0;i<pnum;i++)
              {
                     up=sqrt(n*1.0)+1;
                     if(prime[i]>up)break;
                     int r=0;
                     while(n%prime[i]==0)
                     {
                              r++;
                              n/=prime[i];
                     }
                     ans*=(r*2+1);
              }
              if(n>1)ans*=3;
              printf("%d/n/n",(ans+1)/2);
    } return 0;
}